From 0ff14efa8f1c9e1193d5566da9a9ec3f969eeeb7 Mon Sep 17 00:00:00 2001
From: Philippe WAUTELET <philippe.wautelet@cnrs.fr>
Date: Thu, 19 Sep 2024 11:43:38 +0200
Subject: [PATCH] Philippe 19/09/2024: PHYEX: remove tabulation characters (not
conformant to Fortran standard)
---
src/PHYEX/aux/mode_gradient_m_phy.f90 | 2 +-
src/PHYEX/aux/mode_ini_cst.f90 | 20 ++--
src/PHYEX/aux/shuman_phy.f90 | 32 +++----
src/PHYEX/conv/convect_closure.f90 | 12 +--
src/PHYEX/conv/convect_closure_adjust.f90 | 4 +-
.../conv/convect_closure_adjust_shal.f90 | 4 +-
src/PHYEX/conv/convect_closure_shal.f90 | 12 +--
src/PHYEX/conv/convect_downdraft.f90 | 16 ++--
src/PHYEX/conv/convect_mixing_funct.f90 | 28 +++---
src/PHYEX/conv/convect_precip_adjust.f90 | 28 +++---
src/PHYEX/conv/convect_trigger_funct.f90 | 4 +-
src/PHYEX/conv/convect_trigger_shal.f90 | 4 +-
src/PHYEX/conv/convect_tstep_pref.f90 | 16 ++--
src/PHYEX/conv/convect_updraft_shal.f90 | 12 +--
src/PHYEX/conv/modd_convpar.f90 | 4 +-
src/PHYEX/conv/modd_convpar_shal.f90 | 4 +-
src/PHYEX/conv/shallow_convection.f90 | 8 +-
src/PHYEX/micro/hypser.f90 | 6 +-
src/PHYEX/micro/ini_ice_c1r3.f90 | 16 ++--
src/PHYEX/micro/ini_param_elec.f90 | 92 +++++++++----------
src/PHYEX/micro/ini_rain_c2r2.f90 | 22 ++---
src/PHYEX/micro/ini_rain_ice_elec.f90 | 16 ++--
src/PHYEX/micro/lima_cold.f90 | 8 +-
src/PHYEX/micro/lima_cold_hom_nucl.f90 | 14 +--
src/PHYEX/micro/lima_cold_slow_processes.f90 | 6 +-
src/PHYEX/micro/lima_meyers.f90 | 6 +-
src/PHYEX/micro/lima_mixed.f90 | 14 +--
src/PHYEX/micro/lima_mixed_fast_processes.f90 | 44 ++++-----
src/PHYEX/micro/lima_phillips.f90 | 18 ++--
src/PHYEX/micro/lima_precip_scavenging.f90 | 14 +--
src/PHYEX/micro/lima_warm.f90 | 10 +-
src/PHYEX/micro/lima_warm_coal.f90 | 14 +--
src/PHYEX/micro/lima_warm_evap.f90 | 12 +--
src/PHYEX/micro/lima_warm_nucl.f90 | 14 +--
src/PHYEX/micro/modd_cloudparn.f90 | 4 +-
src/PHYEX/micro/modd_elec_descr.f90 | 30 +++---
src/PHYEX/micro/modd_elec_param.f90 | 16 ++--
src/PHYEX/micro/modd_param_c1r3.f90 | 8 +-
src/PHYEX/micro/modd_param_c2r2.f90 | 4 +-
src/PHYEX/micro/modd_param_lima.f90 | 4 +-
src/PHYEX/micro/modd_param_lima_cold.f90 | 4 +-
src/PHYEX/micro/modd_rain_c2r2_descr.f90 | 8 +-
src/PHYEX/micro/modd_rain_c2r2_khko_param.f90 | 40 ++++----
src/PHYEX/micro/mode_elec_tendencies.f90 | 54 +++++------
src/PHYEX/micro/mode_ini_lima.f90 | 6 +-
src/PHYEX/micro/mode_ini_lima_cold_mixed.f90 | 22 ++---
src/PHYEX/micro/mode_ini_lima_warm.f90 | 24 ++---
.../micro/mode_init_aerosol_properties.f90 | 4 +-
src/PHYEX/micro/mode_lima_ccn_activation.f90 | 14 +--
.../micro/mode_lima_ccn_hom_freezing.f90 | 6 +-
.../mode_lima_conversion_melting_snow.f90 | 4 +-
.../micro/mode_lima_droplets_accretion.f90 | 4 +-
.../mode_lima_droplets_autoconversion.f90 | 6 +-
.../micro/mode_lima_droplets_hom_freezing.f90 | 4 +-
.../mode_lima_droplets_self_collection.f90 | 4 +-
.../micro/mode_lima_drops_self_collection.f90 | 4 +-
src/PHYEX/micro/mode_lima_graupel.f90 | 18 ++--
.../micro/mode_lima_graupel_deposition.f90 | 4 +-
src/PHYEX/micro/mode_lima_hail.f90 | 18 ++--
src/PHYEX/micro/mode_lima_hail_deposition.f90 | 4 +-
...mode_lima_init_ccn_activation_spectrum.f90 | 12 +--
.../micro/mode_lima_meyers_nucleation.f90 | 6 +-
.../mode_lima_phillips_ifn_nucleation.f90 | 18 ++--
src/PHYEX/micro/mode_lima_rain_accr_snow.f90 | 18 ++--
.../micro/mode_lima_rain_evaporation.f90 | 4 +-
src/PHYEX/micro/mode_lima_rain_freezing.f90 | 4 +-
.../micro/mode_lima_snow_self_collection.f90 | 4 +-
src/PHYEX/micro/mode_nzcolx.f90 | 28 +++---
src/PHYEX/micro/mode_rzcolx.f90 | 4 +-
src/PHYEX/micro/prognos_lima.f90 | 4 +-
src/PHYEX/micro/radar_rain_ice.f90 | 14 +--
src/PHYEX/micro/rain_c2r2_khko.f90 | 12 +--
src/PHYEX/micro/rain_ice_elec.f90 | 56 +++++------
src/PHYEX/turb/modd_turbn.f90 | 4 +-
src/PHYEX/turb/mode_ibm_mixinglength.f90 | 4 +-
src/PHYEX/turb/mode_turb_ver.f90 | 22 ++---
src/PHYEX/turb/mode_turb_ver_dyn_flux.f90 | 22 ++---
src/PHYEX/turb/mode_turb_ver_sv_flux.f90 | 22 ++---
src/PHYEX/turb/mode_turb_ver_thermo_corr.f90 | 22 ++---
src/PHYEX/turb/mode_turb_ver_thermo_flux.f90 | 22 ++---
src/PHYEX/turb/mode_update_lm.f90 | 4 +-
81 files changed, 580 insertions(+), 580 deletions(-)
diff --git a/src/PHYEX/aux/mode_gradient_m_phy.f90 b/src/PHYEX/aux/mode_gradient_m_phy.f90
index aeb7f5d1d..7a9d731e7 100644
--- a/src/PHYEX/aux/mode_gradient_m_phy.f90
+++ b/src/PHYEX/aux/mode_gradient_m_phy.f90
@@ -41,7 +41,7 @@ CONTAINS
!!
!! AUTHOR
!! ------
-!! P. Hereil and J. Stein * Meteo France *
+!! P. Hereil and J. Stein * Meteo France *
!!
!! MODIFICATIONS
!! -------------
diff --git a/src/PHYEX/aux/mode_ini_cst.f90 b/src/PHYEX/aux/mode_ini_cst.f90
index 89eef4a60..e07a11db2 100644
--- a/src/PHYEX/aux/mode_ini_cst.f90
+++ b/src/PHYEX/aux/mode_ini_cst.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -36,7 +36,7 @@ CONTAINS
!!
!! AUTHOR
!! ------
-!! V. Ducrocq * Meteo France *
+!! V. Ducrocq * Meteo France *
!!
!! MODIFICATIONS
!! -------------
@@ -72,7 +72,7 @@ IF (LHOOK) CALL DR_HOOK('INI_CST',0,ZHOOK_HANDLE)
CALL CST_ASSOCIATE()
!
!* 1. FUNDAMENTAL CONSTANTS
-! ---------------------
+! ---------------------
!
XPI = 2.*ASIN(1.)
XKARMAN = 0.4
@@ -84,7 +84,7 @@ XAVOGADRO = 6.0221367E+23
!-------------------------------------------------------------------------------
!
!* 2. ASTRONOMICAL CONSTANTS
-! ----------------------
+! ----------------------
!
XDAY = 86400.
XSIYEA = 365.25*XDAY*2.*XPI/ 6.283076
@@ -96,7 +96,7 @@ NDAYSEC = 24*3600 ! Number of seconds in a day
!
!
!* 3. TERRESTRIAL GEOIDE CONSTANTS
-! ----------------------------
+! ----------------------------
!
XRADIUS = 6371229.
XG = 9.80665
@@ -104,7 +104,7 @@ XG = 9.80665
!-------------------------------------------------------------------------------
!
!* 4. REFERENCE PRESSURE
-! -------------------
+! -------------------
!
! Ocean model cst same as in 1D/CMO SURFEX
! values used in ini_cst to overwrite XP00 and XTH00
@@ -118,7 +118,7 @@ XTH00 = 300.
!-------------------------------------------------------------------------------
!
!* 5. RADIATION CONSTANTS
-! -------------------
+! -------------------
!
! Original: XSTEFAN = 2.* XPI**5 * XBOLTZ**4 / (15.* XLIGHTSPEED**2 * XPLANCK**3)
! Juan: XSTEFAN = ( 2.* XPI**5 / 15. ) * ( (XBOLTZ / XPLANCK) * XBOLTZ ) * (XBOLTZ/(XLIGHTSPEED*XPLANCK))**2
@@ -130,7 +130,7 @@ XI0 = 1370.
!-------------------------------------------------------------------------------
!
!* 6. THERMODYNAMIC CONSTANTS
-! -----------------------
+! -----------------------
!
XMD = 28.9644E-3
XMV = 18.0153E-3
@@ -162,14 +162,14 @@ XALPHAOC = 1.9E-4
XBETAOC= 7.7475E-4
!
!* 7. PRECOMPUTED CONSTANTS
-! ---------------------
+! ---------------------
!
RDSRV = XRD/XRV
RDSCPD = XRD/XCPD
RINVXP00 = 1./XP00
!
!* 8. MACHINE PRECISION VALUE DEPENDING of REAL4/8 USE
-! ---------------------
+! ---------------------
!
XMNH_EPSILON = EPSILON (XMNH_EPSILON )
XMNH_HUGE = HUGE (XMNH_HUGE )
diff --git a/src/PHYEX/aux/shuman_phy.f90 b/src/PHYEX/aux/shuman_phy.f90
index 69f793542..42d8c9c3b 100644
--- a/src/PHYEX/aux/shuman_phy.f90
+++ b/src/PHYEX/aux/shuman_phy.f90
@@ -39,7 +39,7 @@ CONTAINS
!!
!! AUTHOR
!! ------
-!! V. Ducrocq * Meteo France *
+!! V. Ducrocq * Meteo France *
!!
!! MODIFICATIONS
!! -------------
@@ -143,7 +143,7 @@ END SUBROUTINE MXM_PHY
!!
!! AUTHOR
!! ------
-!! V. Ducrocq * Meteo France *
+!! V. Ducrocq * Meteo France *
!!
!! MODIFICATIONS
!! -------------
@@ -244,7 +244,7 @@ END SUBROUTINE MXM2D_PHY
!!
!! AUTHOR
!! ------
-!! V. Ducrocq * Meteo France *
+!! V. Ducrocq * Meteo France *
!!
!! MODIFICATIONS
!! -------------
@@ -342,7 +342,7 @@ END SUBROUTINE MZM_PHY
!!
!! AUTHOR
!! ------
-!! V. Ducrocq * Meteo France *
+!! V. Ducrocq * Meteo France *
!!
!! MODIFICATIONS
!! -------------
@@ -438,7 +438,7 @@ END SUBROUTINE MYM2D_PHY
!!
!! AUTHOR
!! ------
-!! V. Ducrocq * Meteo France *
+!! V. Ducrocq * Meteo France *
!!
!! MODIFICATIONS
!! -------------
@@ -532,7 +532,7 @@ END SUBROUTINE MYM_PHY
!!
!! AUTHOR
!! ------
-!! V. Ducrocq * Meteo France *
+!! V. Ducrocq * Meteo France *
!!
!! MODIFICATIONS
!! -------------
@@ -628,7 +628,7 @@ END SUBROUTINE DZM_PHY
!!
!! AUTHOR
!! ------
-!! V. Ducrocq * Meteo France *
+!! V. Ducrocq * Meteo France *
!!
!! MODIFICATIONS
!! -------------
@@ -727,7 +727,7 @@ END SUBROUTINE MZF_PHY
!!
!! AUTHOR
!! ------
-!! V. Ducrocq * Meteo France *
+!! V. Ducrocq * Meteo France *
!!
!! MODIFICATIONS
!! -------------
@@ -832,7 +832,7 @@ END SUBROUTINE MXF_PHY
!!
!! AUTHOR
!! ------
-!! V. Ducrocq * Meteo France *
+!! V. Ducrocq * Meteo France *
!!
!! MODIFICATIONS
!! -------------
@@ -935,7 +935,7 @@ END SUBROUTINE MXF2D_PHY
!!
!! AUTHOR
!! ------
-!! V. Ducrocq * Meteo France *
+!! V. Ducrocq * Meteo France *
!!
!! MODIFICATIONS
!! -------------
@@ -1036,7 +1036,7 @@ END SUBROUTINE MYF_PHY
!!
!! AUTHOR
!! ------
-!! V. Ducrocq * Meteo France *
+!! V. Ducrocq * Meteo France *
!!
!! MODIFICATIONS
!! -------------
@@ -1132,7 +1132,7 @@ END SUBROUTINE MYF2D_PHY
!!
!! AUTHOR
!! ------
-!! V. Ducrocq * Meteo France *
+!! V. Ducrocq * Meteo France *
!!
!! MODIFICATIONS
!! -------------
@@ -1229,7 +1229,7 @@ END SUBROUTINE DZF_PHY
!!
!! AUTHOR
!! ------
-!! V. Ducrocq * Meteo France *
+!! V. Ducrocq * Meteo France *
!!
!! MODIFICATIONS
!! -------------
@@ -1335,7 +1335,7 @@ END SUBROUTINE DXF_PHY
!!
!! AUTHOR
!! ------
-!! V. Ducrocq * Meteo France *
+!! V. Ducrocq * Meteo France *
!!
!! MODIFICATIONS
!! -------------
@@ -1439,7 +1439,7 @@ END SUBROUTINE DXM_PHY
!!
!! AUTHOR
!! ------
-!! V. Ducrocq * Meteo France *
+!! V. Ducrocq * Meteo France *
!!
!! MODIFICATIONS
!! -------------
@@ -1541,7 +1541,7 @@ END SUBROUTINE DYM_PHY
!!
!! AUTHOR
!! ------
-!! V. Ducrocq * Meteo France *
+!! V. Ducrocq * Meteo France *
!!
!! MODIFICATIONS
!! -------------
diff --git a/src/PHYEX/conv/convect_closure.f90 b/src/PHYEX/conv/convect_closure.f90
index f94e9bf86..f347b97e2 100644
--- a/src/PHYEX/conv/convect_closure.f90
+++ b/src/PHYEX/conv/convect_closure.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -197,7 +197,7 @@ REAL, DIMENSION(KLON), INTENT(IN) :: PDXDY ! grid area (m^2)
REAL, DIMENSION(KLON,KLEV),INTENT(IN) :: PTHL ! grid scale enthalpy (J/kg)
REAL, DIMENSION(KLON,KLEV),INTENT(IN) :: PTH ! grid scale theta
REAL, DIMENSION(KLON,KLEV),INTENT(IN) :: PRW ! grid scale total water
- ! mixing ratio
+ ! mixing ratio
REAL, DIMENSION(KLON,KLEV),INTENT(IN) :: PRC ! grid scale r_c
REAL, DIMENSION(KLON,KLEV),INTENT(IN) :: PRI ! grid scale r_i
LOGICAL, DIMENSION(KLON), INTENT(IN) :: OTRIG1 ! logical to keep trace of
@@ -518,9 +518,9 @@ DO JITER = 1, 6 ! Enter adjustment loop to assure that all CAPE is
KFTSTEPS = MAXVAL( ITSTEP(:) )
DO JSTEP = 1, KFTSTEPS ! Enter the fractional time step loop here
!
- ICOUNT(:) = ICOUNT(:) + 1
+ ICOUNT(:) = ICOUNT(:) + 1
!
- GWORK3(:) = ITSTEP(:) >= ICOUNT(:) .AND. GWORK1(:)
+ GWORK3(:) = ITSTEP(:) >= ICOUNT(:) .AND. GWORK1(:)
!
!
!* 7. Assign enthalpy and r_w values at the top and bottom of each
@@ -662,12 +662,12 @@ DO JITER = 1, 6 ! Enter adjustment loop to assure that all CAPE is
CALL CONVECT_SATMIXRATIO( KLON, ZWORK1, ZTELCL, ZWORK3, ZLV, ZLS, ZCPH )
ZWORK3(:) = MIN( .1, MAX( 0., ZWORK3(:) ) )
!
- ! compute theta_e updraft undilute
+ ! compute theta_e updraft undilute
ZTHEUL(:) = ZTLCL(:) * ZPI(:) ** ( 1. - 0.28 * ZRVLCL(:) ) &
* EXP( ( 3374.6525 / ZTLCL(:) - 2.5403 ) &
* ZRVLCL(:) * ( 1. + 0.81 * ZRVLCL(:) ) )
!
- ! compute theta_e saturated environment at LCL
+ ! compute theta_e saturated environment at LCL
ZTHES1(:) = ZTELCL(:) * ZPI(:) ** ( 1. - 0.28 * ZWORK3(:) ) &
* EXP( ( 3374.6525 / ZTELCL(:) - 2.5403 ) &
* ZWORK3(:) * ( 1. + 0.81 * ZWORK3(:) ) )
diff --git a/src/PHYEX/conv/convect_closure_adjust.f90 b/src/PHYEX/conv/convect_closure_adjust.f90
index 177c8844b..6abe76527 100644
--- a/src/PHYEX/conv/convect_closure_adjust.f90
+++ b/src/PHYEX/conv/convect_closure_adjust.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -170,7 +170,7 @@ IKE = KLEV - JCVEXT
PTPR(:) = PZTPR(:) * PADJ(:)
!
DO JK = IKB + 1, IKE
- PUMF(:,JK) = PZUMF(:,JK) * PADJ(:)
+ PUMF(:,JK) = PZUMF(:,JK) * PADJ(:)
PUER(:,JK) = PZUER(:,JK) * PADJ(:)
PUDR(:,JK) = PZUDR(:,JK) * PADJ(:)
PDMF(:,JK) = PZDMF(:,JK) * PADJ(:)
diff --git a/src/PHYEX/conv/convect_closure_adjust_shal.f90 b/src/PHYEX/conv/convect_closure_adjust_shal.f90
index 91ddb2cbe..8ada5dc19 100644
--- a/src/PHYEX/conv/convect_closure_adjust_shal.f90
+++ b/src/PHYEX/conv/convect_closure_adjust_shal.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -126,7 +126,7 @@ IKE = KLEV - JCVEXT
! ----------------------------------------------------
!
DO JK = IKB + 1, IKE
- PUMF(:,JK) = PZUMF(:,JK) * PADJ(:)
+ PUMF(:,JK) = PZUMF(:,JK) * PADJ(:)
PUER(:,JK) = PZUER(:,JK) * PADJ(:)
PUDR(:,JK) = PZUDR(:,JK) * PADJ(:)
END DO
diff --git a/src/PHYEX/conv/convect_closure_shal.f90 b/src/PHYEX/conv/convect_closure_shal.f90
index 73e59decb..9c787e74f 100644
--- a/src/PHYEX/conv/convect_closure_shal.f90
+++ b/src/PHYEX/conv/convect_closure_shal.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -168,7 +168,7 @@ REAL, DIMENSION(KLON), INTENT(IN) :: PDXDY ! grid area (m^2)
REAL, DIMENSION(KLON,KLEV),INTENT(IN) :: PTHL ! grid scale enthalpy (J/kg)
REAL, DIMENSION(KLON,KLEV),INTENT(IN) :: PTH ! grid scale theta
REAL, DIMENSION(KLON,KLEV),INTENT(IN) :: PRW ! grid scale total water
- ! mixing ratio
+ ! mixing ratio
REAL, DIMENSION(KLON,KLEV),INTENT(IN) :: PRC ! grid scale r_c
REAL, DIMENSION(KLON,KLEV),INTENT(IN) :: PRI ! grid scale r_i
LOGICAL, DIMENSION(KLON), INTENT(IN) :: OTRIG1 ! logical to keep trace of
@@ -401,9 +401,9 @@ DO JITER = 1, 4 ! Enter adjustment loop to assure that all CAPE is
KFTSTEPS = MAXVAL( ITSTEP(:) )
DO JSTEP = 1, KFTSTEPS ! Enter the fractional time step loop here
!
- ICOUNT(:) = ICOUNT(:) + 1
+ ICOUNT(:) = ICOUNT(:) + 1
!
- GWORK3(:) = ITSTEP(:) >= ICOUNT(:) .AND. GWORK1(:)
+ GWORK3(:) = ITSTEP(:) >= ICOUNT(:) .AND. GWORK1(:)
!
!
!* 7. Assign enthalpy and r_w values at the top and bottom of each
@@ -534,12 +534,12 @@ DO JITER = 1, 4 ! Enter adjustment loop to assure that all CAPE is
CALL CONVECT_SATMIXRATIO( KLON, ZWORK1, ZTELCL, ZWORK3, ZLV, ZLS, ZCPH )
ZWORK3(:) = MIN( .1, MAX( 0., ZWORK3(:) ) )
!
- ! compute theta_e updraft undilute
+ ! compute theta_e updraft undilute
ZTHEUL(:) = ZTLCL(:) * ZPI(:) ** ( 1. - 0.28 * ZRVLCL(:) ) &
* EXP( ( 3374.6525 / ZTLCL(:) - 2.5403 ) &
* ZRVLCL(:) * ( 1. + 0.81 * ZRVLCL(:) ) )
!
- ! compute theta_e saturated environment at LCL
+ ! compute theta_e saturated environment at LCL
ZTHES1(:) = ZTELCL(:) * ZPI(:) ** ( 1. - 0.28 * ZWORK3(:) ) &
* EXP( ( 3374.6525 / ZTELCL(:) - 2.5403 ) &
* ZWORK3(:) * ( 1. + 0.81 * ZWORK3(:) ) )
diff --git a/src/PHYEX/conv/convect_downdraft.f90 b/src/PHYEX/conv/convect_downdraft.f90
index 814455eb9..fa5cdbb2d 100644
--- a/src/PHYEX/conv/convect_downdraft.f90
+++ b/src/PHYEX/conv/convect_downdraft.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -35,12 +35,12 @@ REAL, DIMENSION(KLON,KLEV), INTENT(IN) :: PRC ! grid scale r_c (cloud water)
REAL, DIMENSION(KLON,KLEV), INTENT(IN) :: PRI ! grid scale r_i (cloud ice)
REAL, DIMENSION(KLON,KLEV), INTENT(IN) :: PPRES ! pressure (Pa)
REAL, DIMENSION(KLON,KLEV), INTENT(IN) :: PDPRES! pressure difference between
- ! bottom and top of layer (Pa)
+ ! bottom and top of layer (Pa)
REAL, DIMENSION(KLON,KLEV), INTENT(IN) :: PZ ! level height (m)
INTEGER, DIMENSION(KLON), INTENT(IN) :: KLCL ! contains vert. index of LCL
INTEGER, DIMENSION(KLON), INTENT(IN) :: KCTL ! contains vert. index of CTL
INTEGER, DIMENSION(KLON), INTENT(IN) :: KETL ! contains vert. index of
- ! equilibrium (zero buoyancy) level
+ ! equilibrium (zero buoyancy) level
INTEGER, DIMENSION(KLON), INTENT(IN) :: KML ! " vert. index of melting level
REAL, DIMENSION(KLON,KLEV), INTENT(IN) :: PUTHL ! updraft enthalpy (J/kg)
REAL, DIMENSION(KLON,KLEV), INTENT(IN) :: PURW ! updraft total water (kg/kg)
@@ -164,12 +164,12 @@ REAL, DIMENSION(KLON,KLEV), INTENT(IN) :: PRC ! grid scale r_c (cloud water)
REAL, DIMENSION(KLON,KLEV), INTENT(IN) :: PRI ! grid scale r_i (cloud ice)
REAL, DIMENSION(KLON,KLEV), INTENT(IN) :: PPRES ! pressure (Pa)
REAL, DIMENSION(KLON,KLEV), INTENT(IN) :: PDPRES! pressure difference between
- ! bottom and top of layer (Pa)
+ ! bottom and top of layer (Pa)
REAL, DIMENSION(KLON,KLEV), INTENT(IN) :: PZ ! level height (m)
INTEGER, DIMENSION(KLON), INTENT(IN) :: KLCL ! contains vert. index of LCL
INTEGER, DIMENSION(KLON), INTENT(IN) :: KCTL ! contains vert. index of CTL
INTEGER, DIMENSION(KLON), INTENT(IN) :: KETL ! contains vert. index of
- ! equilibrium (zero buoyancy) level
+ ! equilibrium (zero buoyancy) level
INTEGER, DIMENSION(KLON), INTENT(IN) :: KML ! " vert. index of melting level
REAL, DIMENSION(KLON,KLEV), INTENT(IN) :: PUTHL ! updraft enthalpy (J/kg)
REAL, DIMENSION(KLON,KLEV), INTENT(IN) :: PURW ! updraft total water (kg/kg)
@@ -308,8 +308,8 @@ END DO
! temperature perturbation due to melting at LFS
ZWORK3(:) = 0.
WHERE( KML(:) > IKB + 2 )
- ZWORK3(:) = ZWORK1(:) * ( ZLS(:) - ZLV(:) ) / ZCPH(:)
- ZDT(:) = ZDT(:) - ZWORK3(:) * REAL(KICE)
+ ZWORK3(:) = ZWORK1(:) * ( ZLS(:) - ZLV(:) ) / ZCPH(:)
+ ZDT(:) = ZDT(:) - ZWORK3(:) * REAL(KICE)
END WHERE
!
!
@@ -353,7 +353,7 @@ JKM = MAXVAL( KLFS(:) )
DO JK = JKM - 1, IKB + 1, -1
DO JI = 1, IIE
IF ( JK < KLFS(JI) .AND. ZWORK2(JI) > PTHES(JI,JK) .AND. GWORK1(JI) ) THEN
- KDBL(JI) = JK
+ KDBL(JI) = JK
GWORK1(JI) = .FALSE.
END IF
END DO
diff --git a/src/PHYEX/conv/convect_mixing_funct.f90 b/src/PHYEX/conv/convect_mixing_funct.f90
index ebf25507e..64949022f 100644
--- a/src/PHYEX/conv/convect_mixing_funct.f90
+++ b/src/PHYEX/conv/convect_mixing_funct.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -117,25 +117,25 @@ IF( KMF == 1 ) THEN
ZW1(:) = 1. / ( 1.+ ZP * ABS ( ZX(:) ) )
ZY(:) = EXP( -0.5 * ZX(:) * ZX(:) )
ZW2(:) = ZA1 * ZW1(:) + ZA2 * ZW1(:) * ZW1(:) + &
- ZA3 * ZW1(:) * ZW1(:) * ZW1(:)
+ ZA3 * ZW1(:) * ZW1(:) * ZW1(:)
ZW11 = ZA1 * ZT1 + ZA2 * ZT1 * ZT1 + ZA3 * ZT1 * ZT1 * ZT1
ENDIF
!
WHERE ( KMF == 1 .AND. ZX(:) >= 0. )
- PER(:) = ZSIGMA * ( 0.5 * ( ZSQRTP - ZE45 * ZW11 &
- - ZY(:) * ZW2(:) ) + ZSIGMA * ( ZE45 - ZY(:) ) ) &
- - 0.5 * ZE45 * PMIXC(:) * PMIXC(:)
- PDR(:) = ZSIGMA*( 0.5 * ( ZY(:) * ZW2(:) - ZE45 * ZW11 ) &
- + ZSIGMA * ( ZE45 - ZY(:) ) ) &
- - ZE45 * ( 0.5 + 0.5 * PMIXC(:) * PMIXC(:) - PMIXC(:) )
+ PER(:) = ZSIGMA * ( 0.5 * ( ZSQRTP - ZE45 * ZW11 &
+ - ZY(:) * ZW2(:) ) + ZSIGMA * ( ZE45 - ZY(:) ) ) &
+ - 0.5 * ZE45 * PMIXC(:) * PMIXC(:)
+ PDR(:) = ZSIGMA*( 0.5 * ( ZY(:) * ZW2(:) - ZE45 * ZW11 ) &
+ + ZSIGMA * ( ZE45 - ZY(:) ) ) &
+ - ZE45 * ( 0.5 + 0.5 * PMIXC(:) * PMIXC(:) - PMIXC(:) )
END WHERE
WHERE ( KMF == 1 .AND. ZX(:) < 0. )
- PER(:) = ZSIGMA*( 0.5 * ( ZY(:) * ZW2(:) - ZE45 * ZW11 ) &
- + ZSIGMA * ( ZE45 - ZY(:) ) ) &
- - 0.5 * ZE45 * PMIXC(:) * PMIXC(:)
- PDR(:) = ZSIGMA * ( 0.5 * ( ZSQRTP - ZE45 * ZW11 - ZY(:) &
- * ZW2(:) ) + ZSIGMA * ( ZE45 - ZY(:) ) ) &
- - ZE45 * ( 0.5 + 0.5 * PMIXC(:) * PMIXC(:) - PMIXC(:) )
+ PER(:) = ZSIGMA*( 0.5 * ( ZY(:) * ZW2(:) - ZE45 * ZW11 ) &
+ + ZSIGMA * ( ZE45 - ZY(:) ) ) &
+ - 0.5 * ZE45 * PMIXC(:) * PMIXC(:)
+ PDR(:) = ZSIGMA * ( 0.5 * ( ZSQRTP - ZE45 * ZW11 - ZY(:) &
+ * ZW2(:) ) + ZSIGMA * ( ZE45 - ZY(:) ) ) &
+ - ZE45 * ( 0.5 + 0.5 * PMIXC(:) * PMIXC(:) - PMIXC(:) )
END WHERE
!
PER(:) = PER(:) * ZFE
diff --git a/src/PHYEX/conv/convect_precip_adjust.f90 b/src/PHYEX/conv/convect_precip_adjust.f90
index 1c32f4d4d..1165b4d75 100644
--- a/src/PHYEX/conv/convect_precip_adjust.f90
+++ b/src/PHYEX/conv/convect_precip_adjust.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -33,7 +33,7 @@ REAL, DIMENSION(KLON), INTENT(IN) :: PMIXF ! critical mixed fraction at LCL
INTEGER, DIMENSION(KLON), INTENT(IN) :: KLCL ! contains vert. index of LCL
INTEGER, DIMENSION(KLON), INTENT(IN) :: KCTL ! contains vert. index of CTL
INTEGER, DIMENSION(KLON), INTENT(IN) :: KETL ! contains vert. index of equilibrium
- ! (zero buoyancy) level
+ ! (zero buoyancy) level
INTEGER, DIMENSION(KLON), INTENT(INOUT) :: KLFS ! contains vert. index of LFS
INTEGER, DIMENSION(KLON), INTENT(INOUT) :: KDBL ! contains vert. index of DBL
!
@@ -136,7 +136,7 @@ REAL, DIMENSION(KLON), INTENT(IN) :: PMIXF ! critical mixed fraction at LCL
INTEGER, DIMENSION(KLON), INTENT(IN) :: KLCL ! contains vert. index of LCL
INTEGER, DIMENSION(KLON), INTENT(IN) :: KCTL ! contains vert. index of CTL
INTEGER, DIMENSION(KLON), INTENT(IN) :: KETL ! contains vert. index of equilibrium
- ! (zero buoyancy) level
+ ! (zero buoyancy) level
INTEGER, DIMENSION(KLON), INTENT(INOUT) :: KLFS ! contains vert. index of LFS
INTEGER, DIMENSION(KLON), INTENT(INOUT) :: KDBL ! contains vert. index of DBL
!
@@ -164,7 +164,7 @@ INTEGER :: JI ! horizontal loop index
!
INTEGER, DIMENSION(KLON) :: IPRL
REAL, DIMENSION(KLON) :: ZWORK1, ZWORK2, ZWORK3, &
- ZWORK4, ZWORK5, ZWORK6 ! work arrays
+ ZWORK4, ZWORK5, ZWORK6 ! work arrays
!
!
!-------------------------------------------------------------------------------
@@ -189,7 +189,7 @@ PTPR(:) = 0.
!
WHERE ( PDTEVR(:) == 0. )
PTPR(:) = PUTPR(:) ! no downdraft evaporation => no downdraft, all
- ! precipitation occurs in updraft
+ ! precipitation occurs in updraft
END WHERE
IF ( COUNT( PDTEVR(:) > 0. ) == 0 ) THEN ! exit routine if no downdraft exists
RETURN
@@ -228,8 +228,8 @@ END DO
PTPR(:) = 0.
DO JK = JKT3, JKT2
WHERE ( JK >= KLCL(:) .AND. PPRES(:,JK) >= ZWORK5(:) - XUSRDPTH )
- PTPR(:) = PTPR(:) + PUPR(:,JK)
- IPRL(:) = JK
+ PTPR(:) = PTPR(:) + PUPR(:,JK)
+ IPRL(:) = JK
END WHERE
END DO
IPRL(:) = MIN( KETL(:), IPRL(:) )
@@ -251,7 +251,7 @@ ZWORK4(:) = PUTPR(:) - PTPR(:)
ZWORK5(:) = 0.
DO JK = JKT3, JKT1
WHERE ( JK >= KLCL(:) .AND. JK <= KLFS(:) )
- ZWORK5(:) = ZWORK5(:) + PUPR(:,JK)
+ ZWORK5(:) = ZWORK5(:) + PUPR(:,JK)
END WHERE
END DO
!
@@ -304,8 +304,8 @@ ZWORK4(:) = ZWORK2(:)
DO JI = 1, IIE
JK = KLFS(JI)
ZWORK2(JI) = ( 1. - ZWORK6(JI) ) + ZWORK6(JI) * &
- ( PUMF(JI,JK) - ( 1. - PMIXF(JI) ) * ZWORK1(JI) ) / &
- MAX( 1.E-1, PUMF(JI,JK) )
+ ( PUMF(JI,JK) - ( 1. - PMIXF(JI) ) * ZWORK1(JI) ) / &
+ MAX( 1.E-1, PUMF(JI,JK) )
END DO
!
!
@@ -313,10 +313,10 @@ JKT1 = MAXVAL( KLFS(:) ) ! value of KLFS might have been reset to IKB above
DO JK = IKB, JKT1
DO JI = 1, IIE
IF ( JK <= KLFS(JI) ) THEN
- PUMF(JI,JK) = PUMF(JI,JK) * ZWORK2(JI)
- PUER(JI,JK) = PUER(JI,JK) * ZWORK2(JI)
- PUDR(JI,JK) = PUDR(JI,JK) * ZWORK2(JI)
- PUPR(JI,JK) = PUPR(JI,JK) * ZWORK2(JI)
+ PUMF(JI,JK) = PUMF(JI,JK) * ZWORK2(JI)
+ PUER(JI,JK) = PUER(JI,JK) * ZWORK2(JI)
+ PUDR(JI,JK) = PUDR(JI,JK) * ZWORK2(JI)
+ PUPR(JI,JK) = PUPR(JI,JK) * ZWORK2(JI)
END IF
END DO
END DO
diff --git a/src/PHYEX/conv/convect_trigger_funct.f90 b/src/PHYEX/conv/convect_trigger_funct.f90
index 39ee732ff..a794a4bdd 100644
--- a/src/PHYEX/conv/convect_trigger_funct.f90
+++ b/src/PHYEX/conv/convect_trigger_funct.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -272,7 +272,7 @@ DO JKK = JKP, JKT
ZTHLCL(:) = ZTHLCL(:) / ZDPTHMIX(:) +.3 ! add small perturbation
ZRVLCL(:) = ZRVLCL(:) / ZDPTHMIX(:) +1.e-4
ZTHVLCL(:) = ZTHLCL(:) * ( 1. + ZEPSA * ZRVLCL(:) ) &
- / ( 1. + ZRVLCL(:) )
+ / ( 1. + ZRVLCL(:) )
!
!* 4.1 Use an empirical direct solution ( Bolton formula )
! to determine temperature and pressure at LCL.
diff --git a/src/PHYEX/conv/convect_trigger_shal.f90 b/src/PHYEX/conv/convect_trigger_shal.f90
index d99b011a2..22d69b2b0 100644
--- a/src/PHYEX/conv/convect_trigger_shal.f90
+++ b/src/PHYEX/conv/convect_trigger_shal.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -276,7 +276,7 @@ DO JKK = JKP, JKT
& (XATPERT * MAX(3.,PTKECLS(:))/XCPD +XBTPERT) * XDTPERT ! add small Temp Perturb.
ZRVLCL(:) = ZRVLCL(:) / ZDPTHMIX(:)
ZTHVLCL(:) = ZTHLCL(:) * ( 1. + ZEPSA * ZRVLCL(:) ) &
- / ( 1. + ZRVLCL(:) )
+ / ( 1. + ZRVLCL(:) )
!
!* 4.1 Use an empirical direct solution ( Bolton formula )
! to determine temperature and pressure at LCL.
diff --git a/src/PHYEX/conv/convect_tstep_pref.f90 b/src/PHYEX/conv/convect_tstep_pref.f90
index 718da2b7b..e6750f4f2 100644
--- a/src/PHYEX/conv/convect_tstep_pref.f90
+++ b/src/PHYEX/conv/convect_tstep_pref.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -144,18 +144,18 @@ END DO
!* 2. Compute convective time step
! ----------------------------
!
- ! compute wind speed at LCL, 500 hPa, CTL
+ ! compute wind speed at LCL, 500 hPa, CTL
DO JI = 1, IIE
JKLC = KLCL(JI)
JKP5 = IP500(JI)
JKCT = KCTL(JI)
ZWORK1(JI) = SQRT( PU(JI,JKLC) * PU(JI,JKLC) + &
- PV(JI,JKLC) * PV(JI,JKLC) )
+ PV(JI,JKLC) * PV(JI,JKLC) )
ZWORK2(JI) = SQRT( PU(JI,JKP5) * PU(JI,JKP5) + &
- PV(JI,JKP5) * PV(JI,JKP5) )
+ PV(JI,JKP5) * PV(JI,JKP5) )
ZWORK3(JI) = SQRT( PU(JI,JKCT) * PU(JI,JKCT) + &
- PV(JI,JKCT) * PV(JI,JKCT) )
+ PV(JI,JKCT) * PV(JI,JKCT) )
END DO
!
ZWORK2(:) = MAX( 0.1, 0.5 * ( ZWORK1(:) + ZWORK2(:) ) )
@@ -178,11 +178,11 @@ DO JI = 1, IIE
( PV(JI,JKCT) - PV(JI,JKLC) ) * &
( PV(JI,JKCT) - PV(JI,JKLC) )
ZWORK1(JI) = 1.E3 * ZWORK2(JI) * SQRT( ZWORK1(JI) ) / &
- MAX( 1.E-2, PZ(JI,JKCT) - PZ(JI,JKLC) )
+ MAX( 1.E-2, PZ(JI,JKCT) - PZ(JI,JKLC) )
END DO
!
PPREF(:) = 1.591 + ZWORK1(:) * ( -.639 + ZWORK1(:) * ( &
- 9.53E-2 - ZWORK1(:) * 4.96E-3 ) )
+ 9.53E-2 - ZWORK1(:) * 4.96E-3 ) )
PPREF(:) = MAX( .4, MIN( PPREF(:), .9 ) )
!
!* 3.2 Precipitation efficiency as a function of cloud base height
@@ -193,7 +193,7 @@ DO JI = 1, IIE
ZCBH(JI) = MAX( 3., ( PZ(JI,JKLC) - PZ(JI,IKB) ) * 3.281E-3 )
END DO
ZWORK1(:) = .9673 + ZCBH(:) * ( -.7003 + ZCBH(:) * ( .1622 + &
- ZCBH(:) * ( -1.2570E-2 + ZCBH(:) * ( 4.2772E-4 - &
+ ZCBH(:) * ( -1.2570E-2 + ZCBH(:) * ( 4.2772E-4 - &
ZCBH(:) * 5.44E-6 ) ) ) )
ZWORK1(:) = MAX( .4, MIN( .9, 1./ ( 1. + ZWORK1(:) ) ) )
!
diff --git a/src/PHYEX/conv/convect_updraft_shal.f90 b/src/PHYEX/conv/convect_updraft_shal.f90
index c7e53eb86..d21666562 100644
--- a/src/PHYEX/conv/convect_updraft_shal.f90
+++ b/src/PHYEX/conv/convect_updraft_shal.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -50,7 +50,7 @@ INTEGER, DIMENSION(KLON), INTENT(IN) :: KPBL ! " vert. index of source layer
!
INTEGER, DIMENSION(KLON), INTENT(OUT):: KCTL ! contains vert. index of CTL
INTEGER, DIMENSION(KLON), INTENT(OUT):: KETL ! contains vert. index of &
- !equilibrium (zero buoyancy) level
+ !equilibrium (zero buoyancy) level
REAL, DIMENSION(KLON,KLEV), INTENT(OUT):: PUMF ! updraft mass flux (kg/s)
REAL, DIMENSION(KLON,KLEV), INTENT(OUT):: PUER ! updraft entrainment (kg/s)
REAL, DIMENSION(KLON,KLEV), INTENT(OUT):: PUDR ! updraft detrainment (kg/s)
@@ -182,7 +182,7 @@ INTEGER, DIMENSION(KLON), INTENT(IN) :: KPBL ! " vert. index of source layer
!
INTEGER, DIMENSION(KLON), INTENT(OUT):: KCTL ! contains vert. index of CTL
INTEGER, DIMENSION(KLON), INTENT(OUT):: KETL ! contains vert. index of &
- !equilibrium (zero buoyancy) level
+ !equilibrium (zero buoyancy) level
REAL, DIMENSION(KLON,KLEV), INTENT(OUT):: PUMF ! updraft mass flux (kg/s)
REAL, DIMENSION(KLON,KLEV), INTENT(OUT):: PUER ! updraft entrainment (kg/s)
REAL, DIMENSION(KLON,KLEV), INTENT(OUT):: PUDR ! updraft detrainment (kg/s)
@@ -283,9 +283,9 @@ JKM = MINVAL( KDPL(:) )
DO JK = JKM, JKP
DO JI = 1, IIE
IF ( JK >= KDPL(JI) .AND. JK < KLCL(JI) ) THEN
- PUMF(JI,JK) = PMFLCL(JI)
- PUTHL(JI,JK) = ZWORK1(JI)
- PUTHV(JI,JK) = PTHLCL(JI) * ( 1. + ZEPSA * PRVLCL(JI) ) / &
+ PUMF(JI,JK) = PMFLCL(JI)
+ PUTHL(JI,JK) = ZWORK1(JI)
+ PUTHV(JI,JK) = PTHLCL(JI) * ( 1. + ZEPSA * PRVLCL(JI) ) / &
( 1. + PRVLCL(JI) )
PURW(JI,JK) = PRVLCL(JI)
END IF
diff --git a/src/PHYEX/conv/modd_convpar.f90 b/src/PHYEX/conv/modd_convpar.f90
index 19cf4725f..aeba52708 100644
--- a/src/PHYEX/conv/modd_convpar.f90
+++ b/src/PHYEX/conv/modd_convpar.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -56,7 +56,7 @@ REAL, SAVE :: XWTRIG ! constant in vertical velocity trigger
!
!
REAL, SAVE :: XNHGAM ! accounts for non-hydrost. pressure
- ! in buoyancy term of w equation
+ ! in buoyancy term of w equation
! = 2 / (1+gamma)
REAL, SAVE :: XTFRZ1 ! begin of freezing interval
REAL, SAVE :: XTFRZ2 ! end of freezing interval
diff --git a/src/PHYEX/conv/modd_convpar_shal.f90 b/src/PHYEX/conv/modd_convpar_shal.f90
index 56abd00ba..d9e22715a 100644
--- a/src/PHYEX/conv/modd_convpar_shal.f90
+++ b/src/PHYEX/conv/modd_convpar_shal.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -63,7 +63,7 @@ REAL, SAVE :: XWTRIG ! constant in vertical velocity trigger
!
!
REAL, SAVE :: XNHGAM ! accounts for non-hydrost. pressure
- ! in buoyancy term of w equation
+ ! in buoyancy term of w equation
! = 2 / (1+gamma)
REAL, SAVE :: XTFRZ1 ! begin of freezing interval
REAL, SAVE :: XTFRZ2 ! end of freezing interval
diff --git a/src/PHYEX/conv/shallow_convection.f90 b/src/PHYEX/conv/shallow_convection.f90
index 548e8bda5..c82bf823f 100644
--- a/src/PHYEX/conv/shallow_convection.f90
+++ b/src/PHYEX/conv/shallow_convection.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1996-2019 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1996-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -486,7 +486,7 @@ ALLOCATE( ILCL(ICONV) )
ALLOCATE( ICTL(ICONV) )
ALLOCATE( IETL(ICONV) )
!
- ! grid scale variables
+ ! grid scale variables
!
ALLOCATE( ZZ(ICONV,IKS) ) ; ZZ = 0.0
ALLOCATE( ZPRES(ICONV,IKS) ); ZPRES = 0.0
@@ -788,8 +788,8 @@ ENDIF
END DO
END DO
!
- ! execute a "scatter"= pack command to store the tendencies in
- ! the final 2D tables
+ ! execute a "scatter"= pack command to store the tendencies in
+ ! the final 2D tables
!
DO JK = IKB, IKE
DO JI = 1, ICONV
diff --git a/src/PHYEX/micro/hypser.f90 b/src/PHYEX/micro/hypser.f90
index 330ba0bbe..783117451 100644
--- a/src/PHYEX/micro/hypser.f90
+++ b/src/PHYEX/micro/hypser.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1996-2019 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1996-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -105,8 +105,8 @@ SERIE: DO JN = 1,5000
ZFAC = ZFAC * ZXH / REAL(JN)
PHYP = ZTEMP + ZFAC
IF (ABS(PHYP-ZTEMP).LE.ZPREC) THEN
- JFLAG = 1
- EXIT SERIE
+ JFLAG = 1
+ EXIT SERIE
END IF
ZTEMP = PHYP
ZZA = ZZA + 1.
diff --git a/src/PHYEX/micro/ini_ice_c1r3.f90 b/src/PHYEX/micro/ini_ice_c1r3.f90
index 43d7c194f..e69453e14 100644
--- a/src/PHYEX/micro/ini_ice_c1r3.f90
+++ b/src/PHYEX/micro/ini_ice_c1r3.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2000-2019 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2000-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -188,7 +188,7 @@ REAL :: ZFAC_ZRNIC ! Zrnic factor used to decrease Long Kernels
!
!
!* 0. FUNCTION STATEMENTS
-! -------------------
+! -------------------
!
!
!* 0.1 G(p) for p_moment of the Generalized GAMMA function
@@ -224,7 +224,7 @@ IF (ALLOCATED(XRTMIN)) RETURN ! In case of nesting microphysics constants of
!-------------------------------------------------------------------------------
!
!* 2. CHARACTERISTICS OF THE SPECIES
-! ------------------------------
+! ------------------------------
!
!
!* 2.1 Raindrops characteristics
@@ -307,7 +307,7 @@ XC1G = 1./2.
!-------------------------------------------------------------------------------
!
!* 3. DIMENSIONAL DISTRIBUTIONS OF THE SPECIES
-! ----------------------------------------
+! ----------------------------------------
!
!
!* 3.2 Ice crystal distribution
@@ -366,7 +366,7 @@ XCTMIN(6) = 1.0E-3
!-------------------------------------------------------------------------------
!
!* 4. CONSTANTS FOR THE SEDIMENTATION
-! -------------------------------
+! -------------------------------
!
!
!* 4.1 Exponent of the fall-speed air density correction
@@ -407,7 +407,7 @@ XFSEDG = XCG*XAG*XCCG*MOMG(XALPHAG,XNUG,XBG+XDG)* &
!-------------------------------------------------------------------------------
!
!* 5. CONSTANTS FOR THE SLOW COLD PROCESSES
-! -------------------------------------
+! -------------------------------------
!
!
!* 5.1 Constants for ice nucleation
@@ -609,7 +609,7 @@ END IF
!-------------------------------------------------------------------------------
!
!* 7. CONSTANTS FOR THE FAST COLD PROCESSES FOR THE AGGREGATES
-! --------------------------------------------------------
+! --------------------------------------------------------
!
!
!* 7.1 Constants for the riming of the aggregates
@@ -821,7 +821,7 @@ END IF
!-------------------------------------------------------------------------------
!
!* 8. CONSTANTS FOR THE FAST COLD PROCESSES FOR THE GRAUPELN
-! ------------------------------------------------------
+! ------------------------------------------------------
!
!
!* 8.1 Constants for the rain contact freezing
diff --git a/src/PHYEX/micro/ini_param_elec.f90 b/src/PHYEX/micro/ini_param_elec.f90
index b144572c5..f8a82bc41 100644
--- a/src/PHYEX/micro/ini_param_elec.f90
+++ b/src/PHYEX/micro/ini_param_elec.f90
@@ -1,11 +1,11 @@
-!MNH_LIC Copyright 1994-2019 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
!-----------------------------------------------------------------
-! ##########################
+! ##########################
MODULE MODI_INI_PARAM_ELEC
-! ##########################
+! ##########################
!
IMPLICIT NONE
INTERFACE
@@ -41,7 +41,7 @@ END MODULE MODI_INI_PARAM_ELEC
!! PURPOSE
!! -------
!! The purpose of this routine is to initialize the constants used to
-!! resolve the electrical scheme.
+!! resolve the electrical scheme.
!!
!!** METHOD
!! ------
@@ -56,19 +56,19 @@ END MODULE MODI_INI_PARAM_ELEC
!!
!! REFERENCE
!! ---------
-!! Helsdon and Farley, 1987: A numerical study of a Montana thunderstorm:
-!! 2. Model results versus observations involving electrical aspects.
-!! J. Geophys. Res., 92, 5661-5675.
+!! Helsdon and Farley, 1987: A numerical study of a Montana thunderstorm:
+!! 2. Model results versus observations involving electrical aspects.
+!! J. Geophys. Res., 92, 5661-5675.
!!
-!! Takahashi, 1978: Riming electrification as a charge generation
-!! mechanism in thunderstorms. J. Atmos. Sci., 35, 1536-1548.
+!! Takahashi, 1978: Riming electrification as a charge generation
+!! mechanism in thunderstorms. J. Atmos. Sci., 35, 1536-1548.
!!
-!! Gardiner et al., 1985: Measurements of initial potential gradient and
-!! particles charges in a Montana supercell thunderstorm.
-!! J. Geophys. Res., 90, 6079-6086.
+!! Gardiner et al., 1985: Measurements of initial potential gradient and
+!! particles charges in a Montana supercell thunderstorm.
+!! J. Geophys. Res., 90, 6079-6086.
!!
-!! Saunders et al., 1991: The effect of liquid water on thunderstorm
-!! charging. J. Geophys. Res., 96, 11007-11017.
+!! Saunders et al., 1991: The effect of liquid water on thunderstorm
+!! charging. J. Geophys. Res., 96, 11007-11017.
!!
!! AUTHOR
!! ------
@@ -97,8 +97,8 @@ END MODULE MODI_INI_PARAM_ELEC
!
!-------------------------------------------------------------------------------
!
-!* 0. DECLARATIONS
-! ------------
+!* 0. DECLARATIONS
+! ------------
!
USE MODD_CST
USE MODD_ELEC_n
@@ -158,7 +158,7 @@ USE MODI_VQZCOLX
!
IMPLICIT NONE
!
-!* 0.1 Declaration of dummy arguments
+!* 0.1 Declaration of dummy arguments
!
TYPE(TFILEDATA), INTENT(IN) :: TPINIFILE ! Initial file
CHARACTER (LEN=*), DIMENSION(:),INTENT(IN) :: HGETSVT
@@ -170,7 +170,7 @@ INTEGER, INTENT(IN) :: IIU ! Upper dimension in x direction (local)
INTEGER, INTENT(IN) :: IJU ! Upper dimension in y direction (local)
INTEGER, INTENT(IN) :: IKU ! Upper dimension in z direction
!
-!* 0.2 Declaration of local variables
+!* 0.2 Declaration of local variables
!
INTEGER :: IND ! Number of intervals to integrate kernels
INTEGER :: J1, JLWC, JTEMP
@@ -205,8 +205,8 @@ REAL, DIMENSION(:), ALLOCATABLE :: ZT, ZLWCC, ZEW
!
!-------------------------------------------------------------------------------
!
-!* 1. PRELIMINARIES
-! -------------
+!* 1. PRELIMINARIES
+! -------------
!
!* 1.1 Constants for electricity
!
@@ -361,8 +361,8 @@ END IF
!
!-------------------------------------------------------------------------------
!
-!* 2. SHAPE PARAMETERS
-! ----------------
+!* 2. SHAPE PARAMETERS
+! ----------------
!
XCXR = -1.0 ! Raindrop characteristic : XCXR (not initialized in ini_rain_ice.f90)
!
@@ -421,8 +421,8 @@ XJCURR_FW = -2.7E-12
!
!-------------------------------------------------------------------------------
!
-!* 3. COEFFICIENTS FOR CHARGE TRANSFERS
-! ---------------------------------
+!* 3. COEFFICIENTS FOR CHARGE TRANSFERS
+! ---------------------------------
!
! proportionality coefficient between mass transfer and charge transfer rates
! the mixing ratio is proportional to the volume of the particle
@@ -439,8 +439,8 @@ XCOEF_RQ_H = XFH / ZBH
!
!-------------------------------------------------------------------------------
!
-!* 3. HOMOGENEOUS NUCLEATION
-! ----------------------
+!* 3. HOMOGENEOUS NUCLEATION
+! ----------------------
!
XALPHACQ = 3. !>
XNUCQ = 1. ! >--- generic values
@@ -454,8 +454,8 @@ XQHON = XQHON / (XLBDACQ**XFC)
!
!-------------------------------------------------------------------------------
!
-!* 4. SEDIMENTATION
-! -------------
+!* 4. SEDIMENTATION
+! -------------
!
IF (ALLOCATED(XQTMIN)) DEALLOCATE(XQTMIN)
IF (ALLOCATED(XRTMIN_ELEC)) DEALLOCATE(XRTMIN_ELEC)
@@ -555,8 +555,8 @@ END IF
!
!-------------------------------------------------------------------------------
!
-!* 5. EVAPORATION OF RAINDROPS
-! ------------------------
+!* 5. EVAPORATION OF RAINDROPS
+! ------------------------
!
!XQREVAV1 = (2. / XPI) * MOMG(ZALPHAR,ZNUR,XFR) / MOMG(ZALPHAR,ZNUR,2.)
!XQREVAV2 = (XPI / ZAR) * (MOMG(ZALPHAR,ZNUR,2.) / MOMG(ZALPHAR,ZNUR,ZBR)) * &
@@ -565,8 +565,8 @@ END IF
!
!-------------------------------------------------------------------------------
!
-!* 6. RIMING OF CLOUD DROPLETS ON SNOW
-! --------------------------------
+!* 6. RIMING OF CLOUD DROPLETS ON SNOW
+! --------------------------------
!
IF (HELEC == 'ELE4') THEN
XEXQSRIMCG = -XFS
@@ -588,8 +588,8 @@ END DO
!
!-------------------------------------------------------------------------------
!
-!* 7. CONTACT FREEZING BETWEEN RAINDROPS AND PRISTINE ICE
-! ---------------------------------------------------
+!* 7. CONTACT FREEZING BETWEEN RAINDROPS AND PRISTINE ICE
+! ---------------------------------------------------
!
IF (HELEC == 'ELE4') THEN
XEXQRCFRIG = - ZDR - XFR - 2.0
@@ -604,8 +604,8 @@ END IF
!
!-------------------------------------------------------------------------------
!
-!* 8. INITIALIZATIONS FOR THE NON INDUCTIVE PROCESSES
-! -----------------------------------------------
+!* 8. INITIALIZATIONS FOR THE NON INDUCTIVE PROCESSES
+! -----------------------------------------------
!
! arrays allocation for NI charging rate
!
@@ -1062,8 +1062,8 @@ END IF
!
!-------------------------------------------------------------------------------
!
-!* 9. NON_INDUCTIVE PROCESS: AGGREGATION OF ICE ON SNOW
-! -------------------------------------------------
+!* 9. NON_INDUCTIVE PROCESS: AGGREGATION OF ICE ON SNOW
+! -------------------------------------------------
!
!* 9.1 Helsdon and Farley (1987) parameterization
!
@@ -1109,8 +1109,8 @@ END IF
!
!-------------------------------------------------------------------------------
!
-!* 10. ACCRETION OF RAINDROPS ON SNOW
-! ------------------------------
+!* 10. ACCRETION OF RAINDROPS ON SNOW
+! ------------------------------
!
IF( .NOT.ALLOCATED(XKER_Q_RACCSS)) ALLOCATE( XKER_Q_RACCSS(IACCLBDAS,IACCLBDAR) )
IF( .NOT.ALLOCATED(XKER_Q_RACCS)) ALLOCATE( XKER_Q_RACCS (IACCLBDAS,IACCLBDAR) )
@@ -1152,8 +1152,8 @@ CALL RSCOLRG (IND, ZALPHAS, ZNUS, ZALPHAR, ZNUR, &
!
!-------------------------------------------------------------------------------
!
-!* 11. DRY GROWTH OF GRAUPELN BY CAPTURE OF SNOW OR ICE
-! ------------------------------------------------
+!* 11. DRY GROWTH OF GRAUPELN BY CAPTURE OF SNOW OR ICE
+! ------------------------------------------------
!
!* 11.1 charge transfer associated to mass transfer
!
@@ -1325,8 +1325,8 @@ ENDIF
!
!-------------------------------------------------------------------------------
!
-!* 12. DRY GROWTH OF GRAUPELN BY CAPTURE OF RAINDROP
-! ---------------------------------------------
+!* 12. DRY GROWTH OF GRAUPELN BY CAPTURE OF RAINDROP
+! ---------------------------------------------
!
IF( .NOT.ALLOCATED(XKER_Q_RDRYG)) ALLOCATE( XKER_Q_RDRYG(IDRYLBDAG,IDRYLBDAR) )
!
@@ -1350,8 +1350,8 @@ CALL RZCOLX (IND, ZALPHAG, ZNUG, ZALPHAR, ZNUR, &
!
!-------------------------------------------------------------------------------
!
-!* 13. UPDATE THE Q=f(D) RELATION
-! --------------------------
+!* 13. UPDATE THE Q=f(D) RELATION
+! --------------------------
!
IF (HCLOUD(1:3) == 'ICE') THEN
XFQUPDC = 400.E6 * MOMG(XALPHACQ,XNUCQ,XFC) / XLBDACQ**XFC ! Nc~400E6 m-3 as
diff --git a/src/PHYEX/micro/ini_rain_c2r2.f90 b/src/PHYEX/micro/ini_rain_c2r2.f90
index 0989dc52d..6b95d1fa4 100644
--- a/src/PHYEX/micro/ini_rain_c2r2.f90
+++ b/src/PHYEX/micro/ini_rain_c2r2.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1996-2019 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1996-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -154,7 +154,7 @@ LOGICAL :: GFLAG ! Logical flag for printing the constatnts on the output
!
!
!* 0. FUNCTION STATEMENTS
-! -------------------
+! -------------------
!
!
!* 0.1 G(p) for p_moment of the Generalized GAMMA function
@@ -187,7 +187,7 @@ IF (ALLOCATED(XRTMIN)) RETURN ! In case of nesting microphysics constants of
!-------------------------------------------------------------------------------
!
!* 2. CHARACTERISTICS OF THE SPECIES
-! ------------------------------
+! ------------------------------
!
!
!* 2.1 Cloud droplet characteristics
@@ -280,7 +280,7 @@ XLBEXR = 1.0/XBR
!-------------------------------------------------------------------------------
!
!* 4. CONSTANTS FOR THE SEDIMENTATION
-! -------------------------------
+! -------------------------------
!
!* 4.1 Exponent of the fall-speed air density correction
!
@@ -304,10 +304,10 @@ XFSEDCC = XCC*GAMMA(XNUC+XDC/XALPHAC)/GAMMA(XNUC)* &
!-------------------------------------------------------------------------------
!
!* 5. CONSTANTS FOR THE NUCLEATION PROCESS
-! -------------------------------------
+! -------------------------------------
!
!
-! Compute CCN spectra parameters from CCN characteristics
+! Compute CCN spectra parameters from CCN characteristics
!
IF (HPARAM_CCN == 'CPB' .AND. HINI_CCN == 'AER') THEN
SELECT CASE (HTYPE_CCN)
@@ -345,7 +345,7 @@ XTHCO = 24.3E-3 ! Air thermal conductivity
ZSURF_TEN = 76.1E-3 ! Surface tension of a water drop at T=0 C
XCSTDCRIT = (XPI/6.)*XRHOLW*( (8.0*ZSURF_TEN )/( 3.0*XRV*XRHOLW ) )**3
!
-! Tabulation of the hypergeometric functions
+! Tabulation of the hypergeometric functions
!
! F(mu,k/2, k/2+1 ,-Beta s**2) and
! F(mu,k/2,(k+3)/2,-Beta s**2) as a function of s
@@ -378,13 +378,13 @@ ELSE ! other cases (but not used)
! defined by CPB98, are unbounded
END IF
!
-! Compute the tabulation of function of T :
+! Compute the tabulation of function of T :
!
! (Psi1)**(3/2)
! XAHENG = -----------------------
! G**(3/2)
!
-! XAHENY = a2 C**p k**q as given by Feingold
+! XAHENY = a2 C**p k**q as given by Feingold
!
NAHEN = 81 ! Tabulation for each Kelvin degree in the range XTT-40 to XTT+40
XAHENINTP1 = 1.0
@@ -538,7 +538,7 @@ XEX2EVAR = -0.5*XCEXVT
!
X0EVAR = (2.0*XPI)*XF0R*GAMMA(XNUR+1./XALPHAR)/GAMMA(XNUR)
X1EVAR = (2.0*XPI)*XF1R*((ZRHO00)**(XCEXVT)*(XCR/0.15E-4))**0.5* &
- GAMMA(XNUR+(XDR+3.0)/(2.0*XALPHAR))/GAMMA(XNUR)
+ GAMMA(XNUR+(XDR+3.0)/(2.0*XALPHAR))/GAMMA(XNUR)
!
XEX0EVAR = 2.0
XEX1EVAR = 2.0 - (XDR+1.0)*0.5
@@ -546,7 +546,7 @@ XEX2EVAR = -0.5*XCEXVT
!
X0EVAR = (12.0)*XF0R*GAMMA(XNUR+1./XALPHAR)/GAMMA(XNUR+3./XALPHAR)
X1EVAR = (12.0)*XF1R*((ZRHO00)**(XCEXVT)*(XCR/0.15E-4))**0.5* &
- GAMMA(XNUR+(XDR+3.0)/(2.0*XALPHAR))/GAMMA(XNUR+3./XALPHAR)
+ GAMMA(XNUR+(XDR+3.0)/(2.0*XALPHAR))/GAMMA(XNUR+3./XALPHAR)
!
!-------------------------------------------------------------------------------
!
diff --git a/src/PHYEX/micro/ini_rain_ice_elec.f90 b/src/PHYEX/micro/ini_rain_ice_elec.f90
index 1a7fa798b..18f437cd2 100644
--- a/src/PHYEX/micro/ini_rain_ice_elec.f90
+++ b/src/PHYEX/micro/ini_rain_ice_elec.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2002-2019 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2002-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -178,7 +178,7 @@ REAL :: PWETLBDAH_MAX, PWETLBDAH_MIN
!-------------------------------------------------------------------------------
!
!* 0. FUNCTION STATEMENTS
-! -------------------
+! -------------------
!
!* 0.1 p_moment of the Generalized GAMMA function
!
@@ -219,7 +219,7 @@ IF (HCLOUD == 'ICE4') XRTMIN(7) = 1.0E-15
!-------------------------------------------------------------------------------
!
!* 2. CHARACTERISTICS OF THE SPECIES
-! ------------------------------
+! ------------------------------
!
!* 2.1 Cloud droplet and Raindrop characteristics
!
@@ -324,7 +324,7 @@ XC1H = 1./2.
!-------------------------------------------------------------------------------
!
!* 3. DIMENSIONAL DISTRIBUTIONS OF THE SPECIES
-! ----------------------------------------
+! ----------------------------------------
!
! 3.1 Cloud droplet distribution
!
@@ -395,7 +395,7 @@ XCONC_URBAN = 5.E8 ! 500/cm3
!-------------------------------------------------------------------------------
!
!* 4. CONSTANTS FOR THE SEDIMENTATION
-! -------------------------------
+! -------------------------------
!
!* 4.1 Exponent of the fall-speed air density correction
!
@@ -448,7 +448,7 @@ XFSEDH = XCH * XAH * XCCH * MOMG(XALPHAH,XNUH,XBH+XDH) * &
!-------------------------------------------------------------------------------
!
!* 5. CONSTANTS FOR THE SLOW COLD PROCESSES
-! -------------------------------------
+! -------------------------------------
!
!* 5.1 Constants for ice nucleation
!
@@ -558,7 +558,7 @@ END IF
!-------------------------------------------------------------------------------
!
!* 6. CONSTANTS FOR THE SLOW WARM PROCESSES
-! -------------------------------------
+! -------------------------------------
!
!* 6.1 Constants for the cloud droplets autoconversion
!
@@ -588,7 +588,7 @@ XEX1EVAR = -1.0 - 0.5 * (XDR + 3.0)
!-------------------------------------------------------------------------------
!
!* 7. CONSTANTS FOR THE FAST COLD PROCESSES FOR THE AGGREGATES
-! --------------------------------------------------------
+! --------------------------------------------------------
!
!
!* 7.1 Constants for the riming of the aggregates
diff --git a/src/PHYEX/micro/lima_cold.f90 b/src/PHYEX/micro/lima_cold.f90
index b9be5a7ac..d3a76799f 100644
--- a/src/PHYEX/micro/lima_cold.f90
+++ b/src/PHYEX/micro/lima_cold.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2013-2020 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2013-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -222,7 +222,7 @@ REAL, DIMENSION(:,:,:), ALLOCATABLE :: PNHS ! Haze homogeneous activation
!
!
!* 0. 3D MICROPHYSCAL VARIABLES
-! -------------------------
+! -------------------------
!
!
! Prepare 3D water mixing ratios
@@ -327,7 +327,7 @@ END IF
!
!
!* 1. COMPUTE THE SEDIMENTATION (RS) SOURCE
-! -------------------------------------
+! -------------------------------------
!
if ( lbu_enable ) then
if ( lbudget_ri .and. osedi ) call Budget_store_init( tbudgets(NBUDGET_RI), 'SEDI', pris(:, :, :) * prhodj(:, :, :) )
@@ -362,7 +362,7 @@ end if
!
!
! COMPUTE THE NUCLEATION PROCESS SOURCES
-! --------------------------------------
+! --------------------------------------
!
IF (LNUCL) THEN
!
diff --git a/src/PHYEX/micro/lima_cold_hom_nucl.f90 b/src/PHYEX/micro/lima_cold_hom_nucl.f90
index 4971d2215..2d9f40463 100644
--- a/src/PHYEX/micro/lima_cold_hom_nucl.f90
+++ b/src/PHYEX/micro/lima_cold_hom_nucl.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2013-2020 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2013-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -221,7 +221,7 @@ INTEGER :: JL, JMOD_CCN, JMOD_IMM ! Loop index
INTEGER :: INEGT ! Case number of hom. nucleation
integer :: idx
LOGICAL, DIMENSION(SIZE(PRHODREF,1),SIZE(PRHODREF,2),SIZE(PRHODREF,3)) &
- :: GNEGT ! Test where to compute the hom. nucleation
+ :: GNEGT ! Test where to compute the hom. nucleation
INTEGER , DIMENSION(SIZE(GNEGT)) :: I1,I2,I3 ! Used to replace the COUNT
!
REAL :: ZEPS ! molar mass ratio
@@ -230,7 +230,7 @@ REAL :: ZEPS ! molar mass ratio
!
!
!* 1. PRELIMINARY COMPUTATIONS
-! ------------------------
+! ------------------------
!
!
! Physical domain
@@ -347,7 +347,7 @@ IF (INEGT.GT.0) THEN
!
!
!* 2. Haze homogeneous freezing
-! ------------------------
+! ------------------------
!
if ( nbumod == kmi .and. lbu_enable .and. ohhoni .and. nmod_ccn > 0 ) then
if ( lbudget_th ) call Budget_store_init( tbudgets(NBUDGET_TH), 'HONH', pths(:, :, :) * prhodj(:, :, :) )
@@ -481,7 +481,7 @@ IF (INEGT.GT.0) THEN
!
!
!* 3. Cloud droplets homogeneous freezing
-! -----------------------------------
+! -----------------------------------
!
!
! Compute the droplet homogeneous nucleation source: RCHONI
@@ -544,7 +544,7 @@ END IF
!
!
!* 4. Rain drops homogeneous freezing
-! -------------------------------
+! -------------------------------
!
!
! Compute the drop homogeneous nucleation source: RRHONG
@@ -590,7 +590,7 @@ END IF
!
!
!* 4. Unpack variables, clean
-! -----------------------
+! -----------------------
!
!
! End of homogeneous nucleation processes
diff --git a/src/PHYEX/micro/lima_cold_slow_processes.f90 b/src/PHYEX/micro/lima_cold_slow_processes.f90
index 364a7007f..017c946ef 100644
--- a/src/PHYEX/micro/lima_cold_slow_processes.f90
+++ b/src/PHYEX/micro/lima_cold_slow_processes.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2013-2020 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2013-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -150,7 +150,7 @@ REAL, DIMENSION(:,:,:),OPTIONAL, INTENT(INOUT) :: PCSS ! Snow/aggregates C.
!* 0.2 Declarations of local variables :
!
LOGICAL, DIMENSION(SIZE(PRHODREF,1),SIZE(PRHODREF,2),SIZE(PRHODREF,3)) &
- :: GMICRO ! Computations only where necessary
+ :: GMICRO ! Computations only where necessary
INTEGER :: IMICRO
INTEGER , DIMENSION(SIZE(GMICRO)) :: I1,I2,I3 ! Used to replace PACK
INTEGER :: JL ! and PACK intrinsics
@@ -344,7 +344,7 @@ IF( IMICRO >= 1 ) THEN
ZCST(:) = (XNS*ZRST(:)*ZLBDAS(:)**XBS)
ZCSS(:) = (XNS*ZRSS(:)*ZLBDAS(:)**XBS)
ELSE IF (NMOM_S.GE.2) THEN
- WHERE (ZRST(:)>XRTMIN(5) .AND. ZCST(:)>XCTMIN(5))
+ WHERE (ZRST(:)>XRTMIN(5) .AND. ZCST(:)>XCTMIN(5))
ZLBDAS(:) = ( XLBS*ZCST(:) / ZRST(:) )**XLBEXS
END WHERE
ELSE
diff --git a/src/PHYEX/micro/lima_meyers.f90 b/src/PHYEX/micro/lima_meyers.f90
index ca1c22825..3b0d3f88e 100644
--- a/src/PHYEX/micro/lima_meyers.f90
+++ b/src/PHYEX/micro/lima_meyers.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2013-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2013-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -174,7 +174,7 @@ INTEGER :: JL ! Loop index
INTEGER :: INEGT ! Case number of nucleation
!
LOGICAL, DIMENSION(SIZE(PRHODREF,1),SIZE(PRHODREF,2),SIZE(PRHODREF,3)) &
- :: GNEGT ! Test where to compute the nucleation
+ :: GNEGT ! Test where to compute the nucleation
!
INTEGER, DIMENSION(SIZE(PRHODREF)) :: I1,I2,I3 ! Indexes for PACK replacement
!
@@ -220,7 +220,7 @@ REAL, DIMENSION(:), ALLOCATABLE :: ZTCELSIUS
!
!
!* 1. PRELIMINARY COMPUTATIONS
-! ------------------------
+! ------------------------
!
!
! Physical domain
diff --git a/src/PHYEX/micro/lima_mixed.f90 b/src/PHYEX/micro/lima_mixed.f90
index 251521427..d928830de 100644
--- a/src/PHYEX/micro/lima_mixed.f90
+++ b/src/PHYEX/micro/lima_mixed.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2013-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2013-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -271,7 +271,7 @@ INTEGER :: JMOD_IFN ! Loop index
!
!
!* 0. 3D MICROPHYSCAL VARIABLES
-! -------------------------
+! -------------------------
!
!
! Prepare 3D water mixing ratios
@@ -376,7 +376,7 @@ END IF
!
!
!* 1. Pack variables, computations only where necessary
-! -------------------------------------------------
+! -------------------------------------------------
!
! Physical domain
!
@@ -548,7 +548,7 @@ IF( IMICRO >= 0 ) THEN
ZCST(:) = (XNS*ZRST(:)*ZLBDAS(:)**XBS)
ZCSS(:) = (XNS*ZRSS(:)*ZLBDAS(:)**XBS)
ELSE IF (NMOM_S.GE.2) THEN
- WHERE (ZRST(:)>XRTMIN(5) .AND. ZCST(:)>XCTMIN(5))
+ WHERE (ZRST(:)>XRTMIN(5) .AND. ZCST(:)>XCTMIN(5))
ZLBDAS(:) = ( XLBS*ZCST(:) / ZRST(:) )**XLBEXS
END WHERE
ELSE
@@ -592,7 +592,7 @@ IF( IMICRO >= 0 ) THEN
!
!
!* 2. Compute the slow processes involving cloud water and graupel
-! ------------------------------------------------------------
+! ------------------------------------------------------------
!
CALL LIMA_MIXED_SLOW_PROCESSES(ZRHODREF, ZZT, ZSSI, PTSTEP, &
ZLSFACT, ZLVFACT, ZAI, ZCJ, &
@@ -608,7 +608,7 @@ IF( IMICRO >= 0 ) THEN
!
!
! 3. Compute the fast RS and RG processes
-! ------------------------------------
+! ------------------------------------
!
IF (NMOM_S.GE.1) THEN
CALL LIMA_MIXED_FAST_PROCESSES(ZRHODREF, ZZT, ZPRES, PTSTEP, &
@@ -628,7 +628,7 @@ END IF
!
!
! 4. Unpack variables
-! ----------------
+! ----------------
!
!
ZW(:,:,:) = PRVS(:,:,:)
diff --git a/src/PHYEX/micro/lima_mixed_fast_processes.f90 b/src/PHYEX/micro/lima_mixed_fast_processes.f90
index 47b59eb80..703af7652 100644
--- a/src/PHYEX/micro/lima_mixed_fast_processes.f90
+++ b/src/PHYEX/micro/lima_mixed_fast_processes.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2013-2020 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2013-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -378,10 +378,10 @@ IF( IGRIM>0 ) THEN
!
WHERE ( GRIM(:) )
ZZW1(:,1) = MIN( PRCS1D(:), &
- XCRIMSS * ZZW(:) * PRCT1D(:) * PCST1D(:) &
- * PLBDAS(:)**XEXCRIMSS &
+ XCRIMSS * ZZW(:) * PRCT1D(:) * PCST1D(:) &
+ * PLBDAS(:)**XEXCRIMSS &
* (1+(XFVELOS/PLBDAS(:))**XALPHAS)**(-XNUS+XEXCRIMSS/XALPHAS) &
- * PRHODREF(:)**(-XCEXVT) )
+ * PRHODREF(:)**(-XCEXVT) )
PRCS1D(:) = PRCS1D(:) - ZZW1(:,1)
PRSS1D(:) = PRSS1D(:) + ZZW1(:,1)
PTHS1D(:) = PTHS1D(:) + ZZW1(:,1)*(PLSFACT(:)-PLVFACT(:)) ! f(L_f*(RCRIMSS))
@@ -400,13 +400,13 @@ IF( IGRIM>0 ) THEN
!
WHERE ( GRIM(:) .AND. (PRSS1D(:)>XRTMIN(5)/PTSTEP) .AND. (PCSS1D(:)>XCTMIN(5)/PTSTEP))
ZZW1(:,2) = MIN( PRCS1D(:), &
- XCRIMSG * PRCT1D(:) * PCST1D(:) & ! RCRIMSG
- * PLBDAS(:)**XEXCRIMSG*(1+(XFVELOS/PLBDAS(:))**XALPHAS)**(-XNUS+XEXCRIMSG/XALPHAS) &
- * PRHODREF(:)**(-XCEXVT) &
- - ZZW1(:,1) )
+ XCRIMSG * PRCT1D(:) * PCST1D(:) & ! RCRIMSG
+ * PLBDAS(:)**XEXCRIMSG*(1+(XFVELOS/PLBDAS(:))**XALPHAS)**(-XNUS+XEXCRIMSG/XALPHAS) &
+ * PRHODREF(:)**(-XCEXVT) &
+ - ZZW1(:,1) )
ZZW1(:,3) = MIN( PRSS1D(:), PCST1D(:) * &
XSRIMCG * PLBDAS(:)**XEXSRIMCG & ! RSRIMCG
- * (1.0 - ZZW(:) )/PTSTEP)
+ * (1.0 - ZZW(:) )/PTSTEP)
PRCS1D(:) = PRCS1D(:) - ZZW1(:,2)
PRSS1D(:) = PRSS1D(:) - ZZW1(:,3)
PRGS1D(:) = PRGS1D(:) + ZZW1(:,2) + ZZW1(:,3)
@@ -906,10 +906,10 @@ IF( IGACC>0 .AND. NMOM_R.GE.2) THEN
DO JJ = 1,IGACC
ZVEC3(JJ) = ( XKER_N_RACCSS(IVEC1(JJ)+1,IVEC2(JJ)+1)* ZVEC2(JJ) &
- XKER_N_RACCSS(IVEC1(JJ)+1,IVEC2(JJ) )*(ZVEC2(JJ) - 1.0) ) &
- * ZVEC1(JJ) &
+ * ZVEC1(JJ) &
- ( XKER_N_RACCSS(IVEC1(JJ) ,IVEC2(JJ)+1)* ZVEC2(JJ) &
- XKER_N_RACCSS(IVEC1(JJ) ,IVEC2(JJ) )*(ZVEC2(JJ) - 1.0) ) &
- * (ZVEC1(JJ) - 1.0)
+ * (ZVEC1(JJ) - 1.0)
END DO
ZZNW(:) = UNPACK( VECTOR=ZVEC3(:),MASK=GACC,FIELD=0.0 )
!
@@ -981,10 +981,10 @@ IF( IGACC>0 .AND. NMOM_R.GE.2) THEN
DO JJ = 1,IGACC
ZVEC3(JJ) = ( XKER_N_SACCRG(IVEC2(JJ)+1,IVEC1(JJ)+1)* ZVEC1(JJ) &
- XKER_N_SACCRG(IVEC2(JJ)+1,IVEC1(JJ) )*(ZVEC1(JJ) - 1.0) ) &
- * ZVEC2(JJ) &
+ * ZVEC2(JJ) &
- ( XKER_N_SACCRG(IVEC2(JJ) ,IVEC1(JJ)+1)* ZVEC1(JJ) &
- XKER_N_SACCRG(IVEC2(JJ) ,IVEC1(JJ) )*(ZVEC1(JJ) - 1.0) ) &
- * (ZVEC2(JJ) - 1.0)
+ * (ZVEC2(JJ) - 1.0)
END DO
ZZNW(:) = UNPACK( VECTOR=ZVEC3(:),MASK=GACC,FIELD=0.0 )
!
@@ -1355,7 +1355,7 @@ IF( IGDRY>0 ) THEN
DO JJ = 1,IGDRY
ZVEC3(JJ) = ( XKER_N_SDRYG(IVEC1(JJ)+1,IVEC2(JJ)+1)* ZVEC2(JJ) &
- XKER_N_SDRYG(IVEC1(JJ)+1,IVEC2(JJ) )*(ZVEC2(JJ) - 1.0) ) &
- * ZVEC1(JJ) &
+ * ZVEC1(JJ) &
- ( XKER_N_SDRYG(IVEC1(JJ) ,IVEC2(JJ)+1)* ZVEC2(JJ) &
- XKER_N_SDRYG(IVEC1(JJ) ,IVEC2(JJ) )*(ZVEC2(JJ) - 1.0) ) &
* (ZVEC1(JJ) - 1.0)
@@ -1442,10 +1442,10 @@ IF( IGDRY>0 ) THEN
DO JJ = 1,IGDRY
ZVEC3(JJ) = ( XKER_N_RDRYG(IVEC1(JJ)+1,IVEC2(JJ)+1)* ZVEC2(JJ) &
- XKER_N_RDRYG(IVEC1(JJ)+1,IVEC2(JJ) )*(ZVEC2(JJ) - 1.0) ) &
- * ZVEC1(JJ) &
+ * ZVEC1(JJ) &
- ( XKER_N_RDRYG(IVEC1(JJ) ,IVEC2(JJ)+1)* ZVEC2(JJ) &
- XKER_N_RDRYG(IVEC1(JJ) ,IVEC2(JJ) )*(ZVEC2(JJ) - 1.0) ) &
- * (ZVEC1(JJ) - 1.0)
+ * (ZVEC1(JJ) - 1.0)
END DO
ZZNW(:) = UNPACK( VECTOR=ZVEC3(:),MASK=GDRY,FIELD=0.0 )
!
@@ -1619,7 +1619,7 @@ WHERE( PRGT1D(:)>XRTMIN(6) .AND. PZT(:)<XTT &
PRRS1D(:) = PRRS1D(:) - ZZW1(:,4)
PRGS1D(:) = PRGS1D(:) + ZRDRYG(:)
PTHS1D(:) = PTHS1D(:) + (ZZW1(:,1)+ZZW1(:,4))*(PLSFACT(:)-PLVFACT(:)) !
- ! f(L_f*(RCDRYG+RRDRYG))
+ ! f(L_f*(RCDRYG+RRDRYG))
!
PCCS1D(:) = MAX( PCCS1D(:)-ZZNW1(:,1),0.0 )
PCIS1D(:) = MAX( PCIS1D(:)-ZZNW1(:,2),0.0 )
@@ -1883,17 +1883,17 @@ IF( IHAIL>0 ) THEN
DO JJ = 1,IGWET
ZVEC3(JJ) = ( XKER_N_SWETH(IVEC1(JJ)+1,IVEC2(JJ)+1)* ZVEC2(JJ) &
- XKER_N_SWETH(IVEC1(JJ)+1,IVEC2(JJ) )*(ZVEC2(JJ) - 1.0) ) &
- * ZVEC1(JJ) &
+ * ZVEC1(JJ) &
- ( XKER_N_SWETH(IVEC1(JJ) ,IVEC2(JJ)+1)* ZVEC2(JJ) &
- XKER_N_SWETH(IVEC1(JJ) ,IVEC2(JJ) )*(ZVEC2(JJ) - 1.0) ) &
- * (ZVEC1(JJ) - 1.0)
+ * (ZVEC1(JJ) - 1.0)
END DO
ZZNW(:) = UNPACK( VECTOR=ZVEC3(:),MASK=GWET,FIELD=0.0 )
!
WHERE( GWET(:) )
ZZNW1(:,3) = MIN( PCSS1D(:),XFNSWETH*ZZNW(:) * PCST1D(:) & ! NSWETH
*( PCHT1D(:) ) &
- *( PRHODREF(:)**(-XCEXVT+1.) ) &
+ *( PRHODREF(:)**(-XCEXVT+1.) ) &
*( XLBNSWETH1/( PLBDAH(:)**2 ) + &
XLBNSWETH2/( PLBDAH(:) * PLBDAS(:) ) + &
XLBNSWETH3/( PLBDAS(:)**2) ) )
@@ -1966,10 +1966,10 @@ IF( IHAIL>0 ) THEN
DO JJ = 1,IGWET
ZVEC3(JJ) = ( XKER_N_GWETH(IVEC1(JJ)+1,IVEC2(JJ)+1)* ZVEC2(JJ) &
- XKER_N_GWETH(IVEC1(JJ)+1,IVEC2(JJ) )*(ZVEC2(JJ) - 1.0) ) &
- * ZVEC1(JJ) &
+ * ZVEC1(JJ) &
- ( XKER_N_GWETH(IVEC1(JJ) ,IVEC2(JJ)+1)* ZVEC2(JJ) &
- XKER_N_GWETH(IVEC1(JJ) ,IVEC2(JJ) )*(ZVEC2(JJ) - 1.0) ) &
- * (ZVEC1(JJ) - 1.0)
+ * (ZVEC1(JJ) - 1.0)
END DO
ZZNW(:) = UNPACK( VECTOR=ZVEC3(:),MASK=GWET,FIELD=0.0 )
!
diff --git a/src/PHYEX/micro/lima_phillips.f90 b/src/PHYEX/micro/lima_phillips.f90
index a2fb7d4ed..f22654046 100644
--- a/src/PHYEX/micro/lima_phillips.f90
+++ b/src/PHYEX/micro/lima_phillips.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2013-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2013-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -201,7 +201,7 @@ INTEGER :: INEGT ! Case number of sedimentation, nucleation,
integer :: idx
!
LOGICAL, DIMENSION(SIZE(PRHODREF,1),SIZE(PRHODREF,2),SIZE(PRHODREF,3)) &
- :: GNEGT ! Test where to compute the nucleation
+ :: GNEGT ! Test where to compute the nucleation
!
INTEGER, DIMENSION(SIZE(PRHODREF)) :: I1,I2,I3 ! Indexes for PACK replacement
!
@@ -257,7 +257,7 @@ REAL, DIMENSION(:), ALLOCATABLE :: ZTCELSIUS, ZZT_SI0_BC
!
!
!* 1. PRELIMINARY COMPUTATIONS
-! ------------------------
+! ------------------------
!
!
! Physical domain
@@ -290,7 +290,7 @@ ZW(:,:,:) = PRVT(:,:,:)*( PPABST(:,:,:)-ZW(:,:,:) ) / ( (CST%XMV/CST%XMD) * ZW(:
!
!
!* 2. COMPUTATIONS ONLY WHERE NECESSARY : PACK
-! ----------------------------------------
+! ----------------------------------------
!
!
GNEGT(:,:,:) = .FALSE.
@@ -386,7 +386,7 @@ ALLOCATE( ZZY (INEGT) ) ; ZZY(:) = 0.0
!
!
!* 3. COMPUTE THE SATURATION OVER WATER AND ICE
-! -----------------------------------------
+! -----------------------------------------
!
!
ZTCELSIUS(:) = ZZT(:)-CST%XTT ! T [°C]
@@ -423,7 +423,7 @@ END IF
!
!
!* 4. COMPUTE THE ACTIVABLE FRACTION OF EACH IFN SPECIE
-! -------------------------------------------------
+! -------------------------------------------------
!
!
! Computation of the reference activity spectrum ( ZZY = N_{IN,1,*} )
@@ -440,7 +440,7 @@ CALL LIMA_PHILLIPS_INTEG(CST, ZZT, ZSI, ZSI0, ZSW, ZZY, Z_FRAC_ACT)
!
!
!* 5. COMPUTE THE HETEROGENEOUS NUCLEATION OF INSOLUBLE IFN
-! -----------------------------------------------------
+! -----------------------------------------------------
!
if ( nbumod == kmi .and. lbu_enable ) then
if ( lbudget_th ) call Budget_store_init( tbudgets(NBUDGET_TH), 'HIND', pths(:, :, :) * prhodj(:, :, :) )
@@ -512,7 +512,7 @@ end if
!
!
!* 6. COMPUTE THE HETEROGENEOUS NUCLEATION OF COATED IFN
-! --------------------------------------------------
+! --------------------------------------------------
!
if ( nbumod == kmi .and. lbu_enable ) then
if ( lbudget_th ) call Budget_store_init( tbudgets(NBUDGET_TH), 'HINC', &
@@ -602,7 +602,7 @@ end if
!
!
!* 7. UNPACK VARIABLES AND CLEAN
-! --------------------------
+! --------------------------
!
!
! End of the heterogeneous nucleation following Phillips 08
diff --git a/src/PHYEX/micro/lima_precip_scavenging.f90 b/src/PHYEX/micro/lima_precip_scavenging.f90
index 78d415912..36059d58d 100644
--- a/src/PHYEX/micro/lima_precip_scavenging.f90
+++ b/src/PHYEX/micro/lima_precip_scavenging.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2013-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2013-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -247,7 +247,7 @@ if ( BUCONF%lbudget_sv ) then
end if
!
!* 1. PRELIMINARY COMPUTATIONS
-! ------------------------
+! ------------------------
!
!
IIB=1+JPHEXT
@@ -278,7 +278,7 @@ CALL GAULAG(ZABSCISSR,ZWEIGHTR,NDIAMR,ZNUM)
!
!
!* 2. NUMERICAL OPTIMIZATION
-! ----------------------
+! ----------------------
!
!
! Optimization : compute in advance parameters depending on rain particles and
@@ -328,7 +328,7 @@ END DO
!
!
!* 3. AEROSOL SCAVENGING
-! ------------------
+! ------------------
!
!
! Iteration over the aerosol type and mode
@@ -568,7 +568,7 @@ end if
!
!
!* 3. SUBROUTINE AND FUNCTION
-! -----------------------
+! -----------------------
!
!
CONTAINS
@@ -676,7 +676,7 @@ INTEGER, SAVE :: ISPLITR
!-------------------------------------------------------------------------------
!
!* 1. COMPUTE THE LOOP BOUNDS
-! -----------------------
+! -----------------------
!
IIB=1+JPHEXT
IIE=SIZE(PZZ,1) - JPHEXT
@@ -688,7 +688,7 @@ IKE=SIZE(PZZ,3) - JPVEXT
!-------------------------------------------------------------------------------
!
!* 2. COMPUTE THE SEDIMENTATION (RS) SOURCE
-! -------------------------------------
+! -------------------------------------
!
!* 2.1 splitting factor for high Courant number C=v_fall*(del_Z/del_T)
!
diff --git a/src/PHYEX/micro/lima_warm.f90 b/src/PHYEX/micro/lima_warm.f90
index caf74ddc2..0b723d732 100644
--- a/src/PHYEX/micro/lima_warm.f90
+++ b/src/PHYEX/micro/lima_warm.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2013-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2013-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -252,7 +252,7 @@ LOGICAL, DIMENSION(SIZE(PRHODREF,1),SIZE(PRHODREF,2)) :: GDEP
!
!
!* 0. 3D MICROPHYSCAL VARIABLES
-! -------------------------
+! -------------------------
!
!
! Prepare 3D water mixing ratios
@@ -297,7 +297,7 @@ END IF
!
!
!* 1. COMPUTE THE SLOPE PARAMETERS ZLBDC,ZLBDR
-! ----------------------------------------
+! ----------------------------------------
!
!
ZWLBDC3(:,:,:) = XMNH_HUGE
@@ -320,7 +320,7 @@ ZT(:,:,:) = PTHT(:,:,:) * (PPABST(:,:,:)/XP00)**(XRD/XCPD)
!
!
!* 2. COMPUTE THE SEDIMENTATION (RS) SOURCE
-! -------------------------------------
+! -------------------------------------
!
!
if ( lbudget_rc .and. osedc ) call Budget_store_init( tbudgets(NBUDGET_RC), 'SEDI', prcs(:, :, :) * prhodj(:, :, :) )
@@ -368,7 +368,7 @@ END IF
!-------------------------------------------------------------------------------
!
!* 2. COMPUTES THE NUCLEATION PROCESS SOURCES
-! --------------------------------------
+! --------------------------------------
!
!
IF ( LACTI .AND. NMOD_CCN > 0 .AND. .NOT. LSPRO ) THEN
diff --git a/src/PHYEX/micro/lima_warm_coal.f90 b/src/PHYEX/micro/lima_warm_coal.f90
index bac9042e6..9855d9515 100644
--- a/src/PHYEX/micro/lima_warm_coal.f90
+++ b/src/PHYEX/micro/lima_warm_coal.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2013-2020 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2013-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -186,7 +186,7 @@ INTEGER :: IIB, IIE, IJB, IJE, IKB, IKE ! Physical domain
!
!
!* 1. PREPARE COMPUTATIONS - PACK
-! ---------------------------
+! ---------------------------
!
!
IIB=1+JPHEXT
@@ -249,7 +249,7 @@ IF( IMICRO >= 0 ) THEN
IF (NMOM_R.GE.2) THEN
!
!* 2. Self-collection of cloud droplets
-! ------------------------------------
+! ------------------------------------
!
!
if ( lbudget_sv ) call Budget_store_init( tbudgets(NBUDGET_SV1 - 1 + nsv_lima_nc), 'SELF', pccs(:, :, :) * prhodj(:, :, :) )
@@ -270,7 +270,7 @@ IF (NMOM_R.GE.2) THEN
!
!
!* 3. Autoconversion of cloud droplets (Berry-Reinhardt parameterization)
-! ----------------------------------------------------------------------
+! ----------------------------------------------------------------------
!
!
!
@@ -333,7 +333,7 @@ IF (NMOM_R.GE.2) THEN
!
!
!* 4. Accretion sources
-! --------------------
+! --------------------
!
!
GACCR(:) = ZRRT(:)>XRTMIN(3) .AND. ZCRT(:)>XCTMIN(3)
@@ -405,7 +405,7 @@ IF (NMOM_R.GE.2) THEN
!
!
!* 5. Self collection - Coalescence/Break-up
-! -----------------------------------------
+! -----------------------------------------
!
!
IF( IACCR >= 0 ) THEN
@@ -454,7 +454,7 @@ END IF
!
!
!* 6. Unpack and clean
-! -------------------
+! -------------------
!
!
ZW(:,:,:) = PRCS(:,:,:)
diff --git a/src/PHYEX/micro/lima_warm_evap.f90 b/src/PHYEX/micro/lima_warm_evap.f90
index 7d374b464..c2bfb1daf 100644
--- a/src/PHYEX/micro/lima_warm_evap.f90
+++ b/src/PHYEX/micro/lima_warm_evap.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2013-2019 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2013-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -165,7 +165,7 @@ INTEGER :: IIB, IIE, IJB, IJE, IKB, IKE ! Physical domain
!
!
!* 1. PREPARE COMPUTATIONS - PACK
-! ---------------------------
+! ---------------------------
!
!
IIB=1+JPHEXT
@@ -239,7 +239,7 @@ IF( IEVAP >= 0 ) THEN
!
!
!* 2. compute the evaporation of rain drops
-! ----------------------------------------
+! ----------------------------------------
!
!
ZZW3(:) = MAX((1.0 - ZRVT(:)/ZZW1(:)),0.0) ! Subsaturation
@@ -274,7 +274,7 @@ IF( IEVAP >= 0 ) THEN
!
!
!* 3. Unpack and clean
-! -------------------
+! -------------------
!
!
ZW(:,:,:) = PRVS(:,:,:)
@@ -315,12 +315,12 @@ IF( IEVAP >= 0 ) THEN
!
!
!* 4. Update Nr if: 80 microns < Dr < D_h
-! ---------------------------------------
+! ---------------------------------------
!
!
IF (LKHKO) THEN
!* correct negative values for rain
-! --------------------------------
+! --------------------------------
!
WHERE (PRRS(:,:,:)<0.)
PRCS(:,:,:) = PRCS(:,:,:)+PRRS(:,:,:)
diff --git a/src/PHYEX/micro/lima_warm_nucl.f90 b/src/PHYEX/micro/lima_warm_nucl.f90
index fc012c8f1..3708626bf 100644
--- a/src/PHYEX/micro/lima_warm_nucl.f90
+++ b/src/PHYEX/micro/lima_warm_nucl.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2013-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2013-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -226,7 +226,7 @@ TYPE(TFIELDMETADATA) :: TZFIELD
ILUOUT = TLUOUT%NLU
!
!* 1. PREPARE COMPUTATIONS - PACK
-! ---------------------------
+! ---------------------------
!
!
IIB=1+JPHEXT
@@ -368,7 +368,7 @@ ELSE
!
!
!* 2. compute the constant term (ZZW3) relative to smax
-! ----------------------------------------------------
+! ----------------------------------------------------
!
! Remark : in LIMA's nucleation parameterization, Smax=0.01 for a supersaturation of 1% !
!
@@ -442,7 +442,7 @@ ELSE
!
!
!* 3. Compute the maximum of supersaturation
-! -----------------------------------------
+! -----------------------------------------
!
!
! estimate S_max for the CPB98 parameterization with SEVERAL aerosols mode
@@ -462,7 +462,7 @@ ELSE
!
!
!* 4. Compute the nucleus source
-! -----------------------------
+! -----------------------------
!
!
! Again : Smax=0.01 for a 1% supersaturation
@@ -545,7 +545,7 @@ END IF ! CH_AER_ACTIVATION
!
!
!* 5. Cleaning
-! -----------
+! -----------
!
!
DEALLOCATE(IVEC1)
@@ -615,7 +615,7 @@ END IF
!
!
!* 6. Functions used to compute the maximum of supersaturation
-! -----------------------------------------------------------
+! -----------------------------------------------------------
!
!
CONTAINS
diff --git a/src/PHYEX/micro/modd_cloudparn.f90 b/src/PHYEX/micro/modd_cloudparn.f90
index fd8f22971..4d0749fea 100644
--- a/src/PHYEX/micro/modd_cloudparn.f90
+++ b/src/PHYEX/micro/modd_cloudparn.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -31,7 +31,7 @@
!!
!! AUTHOR
!! ------
-!! E. Richard *Laboratoire d'Aerologie*
+!! E. Richard *Laboratoire d'Aerologie*
!!
!! MODIFICATIONS
!! -------------
diff --git a/src/PHYEX/micro/modd_elec_descr.f90 b/src/PHYEX/micro/modd_elec_descr.f90
index db6aaa958..83b3442e3 100644
--- a/src/PHYEX/micro/modd_elec_descr.f90
+++ b/src/PHYEX/micro/modd_elec_descr.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2002-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2002-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -9,23 +9,23 @@
!
!!**** *MODD_ELEC_DESCR* - declaration of the electrical descriptive constants
!!
-!! PURPOSE
-!! -------
+!! PURPOSE
+!! -------
!
-!!** IMPLICIT ARGUMENTS
-!! ------------------
-!! None
+!!** IMPLICIT ARGUMENTS
+!! ------------------
+!! None
!!
-!! REFERENCE
-!! ---------
+!! REFERENCE
+!! ---------
!!
-!! AUTHOR
-!! ------
+!! AUTHOR
+!! ------
!! Gilles Molinie * Laboratoire d'Aerologie *
!!
-!! MODIFICATIONS
-!! -------------
-!! Original 14/11/02
+!! MODIFICATIONS
+!! -------------
+!! Original 14/11/02
!! M. Chong 26/01/10 Small ions parameters
!! +Option for Fair weather field from
!! Helsdon-Farley (JGR, 1987, 5661-5675)
@@ -36,8 +36,8 @@
!!
!-------------------------------------------------------------------------------
!
-!* 0. DECLARATIONS
-! ------------
+!* 0. DECLARATIONS
+! ------------
!
USE MODD_PARAMETERS, ONLY: JPSVNAMELGTMAX
diff --git a/src/PHYEX/micro/modd_elec_param.f90 b/src/PHYEX/micro/modd_elec_param.f90
index 128c292e9..130dbbbb2 100644
--- a/src/PHYEX/micro/modd_elec_param.f90
+++ b/src/PHYEX/micro/modd_elec_param.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2002-2023 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2002-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -8,16 +8,16 @@
! #######################
!
!!**** *MODD_ELEC_PARAM* - declaration of some electrical factors
-!! extensively used in the electrical scheme.
+!! extensively used in the electrical scheme.
!!
!! PURPOSE
!! -------
-! The purpose of this declarative module is to declare some precomputed
-! electrical parameters directly used in routines related to cloud electricity
+! The purpose of this declarative module is to declare some precomputed
+! electrical parameters directly used in routines related to cloud electricity
!!
!!** IMPLICIT ARGUMENTS
!! ------------------
-!! None
+!! None
!!
!! REFERENCE
!! ---------
@@ -37,8 +37,8 @@
!!
!-------------------------------------------------------------------------------
!
-!* 0. DECLARATIONS
-! ------------
+!* 0. DECLARATIONS
+! ------------
!
USE MODD_PARAMETERS, ONLY: JPMODELMAX
!
@@ -51,7 +51,7 @@ REAL :: XCOEF_RQ_V, XCOEF_RQ_C, & ! Constants for proportionality
XCOEF_RQ_S, XCOEF_RQ_G, & ! charge transfer
XCOEF_RQ_H
!
-REAL :: XQHON ! Constant for spontaneous freezing of droplets if T<-35°
+REAL :: XQHON ! Constant for spontaneous freezing of droplets if T<-35�
!
REAL, DIMENSION(:), ALLOCATABLE :: XFQSED ! Constant for sedimentation of
! electric charge in LIMA
diff --git a/src/PHYEX/micro/modd_param_c1r3.f90 b/src/PHYEX/micro/modd_param_c1r3.f90
index cc1e68531..a1a0401c5 100644
--- a/src/PHYEX/micro/modd_param_c1r3.f90
+++ b/src/PHYEX/micro/modd_param_c1r3.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -33,7 +33,7 @@
!!
!! AUTHOR
!! ------
-!! J.-P. Pinty *Laboratoire d'Aerologie*
+!! J.-P. Pinty *Laboratoire d'Aerologie*
!!
!! MODIFICATIONS
!! -------------
@@ -47,8 +47,8 @@
!
IMPLICIT NONE
REAL,SAVE :: XALPHAI,XNUI, & ! Pristine ice distribution parameters
- XALPHAS,XNUS, & ! Snow/aggregate distribution parameters
- XALPHAG,XNUG ! Graupel distribution parameters
+ XALPHAS,XNUS, & ! Snow/aggregate distribution parameters
+ XALPHAG,XNUG ! Graupel distribution parameters
REAL,SAVE :: XFACTNUC_DEP,XFACTNUC_CON ! Amplification factor for IN conc.
! DEP refers to DEPosition mode
! CON refers to CONtact mode
diff --git a/src/PHYEX/micro/modd_param_c2r2.f90 b/src/PHYEX/micro/modd_param_c2r2.f90
index 20f080c03..19506c9c9 100644
--- a/src/PHYEX/micro/modd_param_c2r2.f90
+++ b/src/PHYEX/micro/modd_param_c2r2.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -39,7 +39,7 @@
!!
!! AUTHOR
!! ------
-!! J.-P. Pinty *Laboratoire d'Aerologie*
+!! J.-P. Pinty *Laboratoire d'Aerologie*
!!
!! MODIFICATIONS
!! -------------
diff --git a/src/PHYEX/micro/modd_param_lima.f90 b/src/PHYEX/micro/modd_param_lima.f90
index 1a903452b..eeb3f9084 100644
--- a/src/PHYEX/micro/modd_param_lima.f90
+++ b/src/PHYEX/micro/modd_param_lima.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2013-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2013-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -24,7 +24,7 @@
!!
!! AUTHOR
!! ------
-!! J.-P. Pinty *Laboratoire d'Aerologie*
+!! J.-P. Pinty *Laboratoire d'Aerologie*
!! S. Berthet * Laboratoire d'Aerologie*
!! B. Vié * Laboratoire d'Aerologie*
!!
diff --git a/src/PHYEX/micro/modd_param_lima_cold.f90 b/src/PHYEX/micro/modd_param_lima_cold.f90
index cc3e3adc1..e5ddccd64 100644
--- a/src/PHYEX/micro/modd_param_lima_cold.f90
+++ b/src/PHYEX/micro/modd_param_lima_cold.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2013-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2013-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -12,7 +12,7 @@
!! the LIMA cold scheme.
!! AUTHOR
!! ------
-!! J.-P. Pinty *Laboratoire d'Aerologie*
+!! J.-P. Pinty *Laboratoire d'Aerologie*
!! S. Berthet * Laboratoire d'Aerologie*
!! B. Vié * Laboratoire d'Aerologie*
!!
diff --git a/src/PHYEX/micro/modd_rain_c2r2_descr.f90 b/src/PHYEX/micro/modd_rain_c2r2_descr.f90
index cabbebbd7..652465f37 100644
--- a/src/PHYEX/micro/modd_rain_c2r2_descr.f90
+++ b/src/PHYEX/micro/modd_rain_c2r2_descr.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2000-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2000-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -37,7 +37,7 @@
!!
!! AUTHOR
!! ------
-!! J.-P. Pinty *Laboratoire d'Aerologie*
+!! J.-P. Pinty *Laboratoire d'Aerologie*
!!
!! MODIFICATIONS
!! -------------
@@ -54,14 +54,14 @@ IMPLICIT NONE
REAL,SAVE :: XCEXVT ! air density fall speed correction
!
REAL,SAVE :: XAR,XBR,XCR,XDR,XF0R,XF1R, & ! Raindrop charact.
- XAC,XBC,XCC,XDC,XF0C,XF2C,XC1C ! Cloud droplet charact.
+ XAC,XBC,XCC,XDC,XF0C,XF2C,XC1C ! Cloud droplet charact.
!
REAL,DIMENSION(:),SAVE,ALLOCATABLE :: XRTMIN
! Min values of the mixing ratios
REAL,DIMENSION(:),SAVE,ALLOCATABLE :: XCTMIN
! Min values of the drop concentrations
REAL,SAVE :: XLBC, XLBEXC, & ! shape parameters of the cloud droplets
- XLBR, XLBEXR ! shape parameters of the raindrops
+ XLBR, XLBEXR ! shape parameters of the raindrops
!
CHARACTER(LEN=JPSVNAMELGTMAX),DIMENSION(4),PARAMETER &
:: C2R2NAMES=(/'CCCN ','CCLOUD','CRAIN ','SUPSAT'/)
diff --git a/src/PHYEX/micro/modd_rain_c2r2_khko_param.f90 b/src/PHYEX/micro/modd_rain_c2r2_khko_param.f90
index 639abf3c1..2e14cf962 100644
--- a/src/PHYEX/micro/modd_rain_c2r2_khko_param.f90
+++ b/src/PHYEX/micro/modd_rain_c2r2_khko_param.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -30,8 +30,8 @@
!!
!! AUTHOR
!! ------
-!! J.-P. Pinty *Laboratoire d'Aerologie*
-!! O.Geoffroy (GMEI)
+!! J.-P. Pinty *Laboratoire d'Aerologie*
+!! O.Geoffroy (GMEI)
!!
!! MODIFICATIONS
!! -------------
@@ -51,43 +51,43 @@ REAL,SAVE :: XFSEDRR,XFSEDCR, & ! Constants for sedimentation
!
!
REAL,SAVE :: XDIVA, & ! Diffusivity of water vapor
- XTHCO ! Thermal conductivity
+ XTHCO ! Thermal conductivity
REAL,SAVE :: XWMIN ! Min value of updraft velocity
- ! to enable nucleation process
+ ! to enable nucleation process
REAL,SAVE :: XTMIN ! Min value of
! temperature evolution
- ! to enable nucleation process
+ ! to enable nucleation process
REAL,SAVE :: XCSTHEN,XCSTDCRIT ! Cst for HEN precalculations
INTEGER, SAVE :: NHYP ! Number of value of the HYP
- ! functions
+ ! functions
REAL,SAVE :: XHYPINTP1, XHYPINTP2 ! Factors defining the
- ! supersaturation log scale
+ ! supersaturation log scale
REAL, DIMENSION(:), SAVE, ALLOCATABLE & ! Tabulated HYPgeometric
- :: XHYPF12, XHYPF32 ! functions used in HEN
+ :: XHYPF12, XHYPF32 ! functions used in HEN
INTEGER, SAVE :: NAHEN ! Number of value of the AHEN
- ! functions
+ ! functions
REAL,SAVE :: XAHENINTP1, XAHENINTP2 ! Factors defining the
- ! temperatures in lin scale
+ ! temperatures in lin scale
REAL, DIMENSION(:), SAVE, ALLOCATABLE & !
:: XAHENG,XPSI1, XPSI3, & ! Twomey-CPB98 and
- XAHENF,XAHENY ! Feingold-Heymsfield
- ! parameterization to compute Smax
+ XAHENF,XAHENY ! Feingold-Heymsfield
+ ! parameterization to compute Smax
REAL,SAVE :: XWCOEF_F1, XWCOEF_F2, XWCOEF_F3, & ! COEF_F of the polynomial temp.
XWCOEF_Y1, XWCOEF_Y2, XWCOEF_Y3 ! COEF_Y of the polynomial temp.
- ! function powering W
+ ! function powering W
!
!
REAL,SAVE :: XKERA1, XKERA2 ! Constants to define the lin
- ! and parabolic kernel param.
+ ! and parabolic kernel param.
REAL,SAVE :: XSELFC ! Constants for cloud droplet
! selfcollection : SELF
!
REAL,SAVE :: XAUTO1, XAUTO2, XCAUTR, & ! Constants for cloud droplet
- XLAUTR, XLAUTR_THRESHOLD, & ! autoconversion : AUT
- XITAUTR, XITAUTR_THRESHOLD
+ XLAUTR, XLAUTR_THRESHOLD, & ! autoconversion : AUT
+ XITAUTR, XITAUTR_THRESHOLD
!
REAL,SAVE :: XACCR1, XACCR2, XACCR3, & ! Constants for the accretion
- XACCR4, XACCR5, XACCR6, & ! process
+ XACCR4, XACCR5, XACCR6, & ! process
XACCR_CLARGE1, XACCR_CLARGE2, XACCR_RLARGE1, XACCR_RLARGE2, &
XACCR_CSMALL1, XACCR_CSMALL2, XACCR_RSMALL1, XACCR_RSMALL2
!
@@ -98,7 +98,7 @@ REAL,SAVE :: XSPONBUD1,XSPONBUD2,XSPONBUD3, & ! Spontaneous Break-up
XSPONCOEF2 ! (drop size limiter)
!
REAL,SAVE :: X0EVAR, X1EVAR, & ! Constants for raindrop
- XEX0EVAR, XEX1EVAR, XEX2EVAR ! evaporation: EVA
+ XEX0EVAR, XEX1EVAR, XEX2EVAR ! evaporation: EVA
!
REAL,SAVE :: XCONCC_INI, XCONCR_PARAM_INI ! Used to initialize the
! concentrations from mixing ratios
@@ -112,7 +112,7 @@ REAL,SAVE :: XCREC, XCRER
! Factors to compute reff when cloud and rain are present
!
REAL,SAVE :: XR0 ! new drizzle drops radius
- ! autoconversion
+ ! autoconversion
!
REAL,SAVE :: XCEVAP ! Constants for raindrop
! evaporation
diff --git a/src/PHYEX/micro/mode_elec_tendencies.f90 b/src/PHYEX/micro/mode_elec_tendencies.f90
index af2a62376..1cf805a41 100644
--- a/src/PHYEX/micro/mode_elec_tendencies.f90
+++ b/src/PHYEX/micro/mode_elec_tendencies.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2022-2023 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2022-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -2265,8 +2265,8 @@ CONTAINS
! following Gardiner et al. (1985)
!
!
-!* 0. DECLARATIONS
-! ------------
+!* 0. DECLARATIONS
+! ------------
!
IMPLICIT NONE
!
@@ -2316,8 +2316,8 @@ END SUBROUTINE ELEC_INIT_NOIND_GARDI
! following Saunders et al. (1991)
!
!
-!* 0. DECLARATIONS
-! ------------
+!* 0. DECLARATIONS
+! ------------
!
IMPLICIT NONE
!
@@ -2415,8 +2415,8 @@ END SUBROUTINE ELEC_INIT_NOIND_EWC
! following Saunders and Peck (1998)
!
!
-!* 0. DECLARATIONS
-! ------------
+!* 0. DECLARATIONS
+! ------------
!
IMPLICIT NONE
!
@@ -2636,8 +2636,8 @@ END SUBROUTINE ELEC_INIT_NOIND_RAR
! following Takahashi (1978)
!
!
-!* 0. DECLARATIONS
-! ------------
+!* 0. DECLARATIONS
+! ------------
!
IMPLICIT NONE
!
@@ -2681,8 +2681,8 @@ END SUBROUTINE ELEC_INIT_NOIND_TAKAH
! following Tsenova and Mitzeva (2009)
!
!
-!* 0. DECLARATIONS
-! ------------
+!* 0. DECLARATIONS
+! ------------
!
IMPLICIT NONE
!
@@ -2761,8 +2761,8 @@ END SUBROUTINE ELEC_INIT_NOIND_TEEWC
! Purpose : compute dQ(RAR,T) in the parameterization of Saunders and Peck (1998)
!
!
-!* 0. DECLARATIONS
-! ------------
+!* 0. DECLARATIONS
+! ------------
!
IMPLICIT NONE
!
@@ -2774,7 +2774,7 @@ REAL, DIMENSION(KMICRO), INTENT(INOUT) :: PDQ ! interpolated dQ
!
!
!* 1. COMPUTE dQ(RAR, T)
-! ------------------
+! ------------------
!
PDQ(:) = 0.
!
@@ -2803,8 +2803,8 @@ END SUBROUTINE ELEC_INI_NI_SAP98
! Purpose : interpolate dQ from a lookup table at each gridpoint
!
!
-!* 0. DECLARATIONS
-! ------------
+!* 0. DECLARATIONS
+! ------------
!
IMPLICIT NONE
!
@@ -2823,8 +2823,8 @@ REAL, DIMENSION(:), ALLOCATABLE :: ZVECT1, ZVECT2
INTEGER, DIMENSION(:), ALLOCATABLE :: IVECT1, IVECT2
!
!
-!* 1. FIND THE INDEXES FOR RAR/EW AND T
-! ---------------------------------
+!* 1. FIND THE INDEXES FOR RAR/EW AND T
+! ---------------------------------
!
PDQ(:) = 0.
!
@@ -2881,8 +2881,8 @@ IF (IGAUX > 0) THEN
ENDWHERE
!
!
-!* 2. INTERPOLATE dQ(RAR or EW,T)
-! ---------------------------
+!* 2. INTERPOLATE dQ(RAR or EW,T)
+! ---------------------------
!
ZDQ_INTERP(:) = BI_LIN_INTP_V( PTABLE, IVECT2, IVECT1, ZVECT2, ZVECT1, &
IGAUX )
@@ -3197,7 +3197,7 @@ END SUBROUTINE ELEC_IDRYG_B
! between snow and graupel
!
!
-!* 0. DECLARATIONS
+!* 0. DECLARATIONS
! ------------
!
IMPLICIT NONE
@@ -3424,12 +3424,12 @@ END SUBROUTINE COMPUTE_CHARGE_TRANSFER
! | PDX |x1 |
! | |
!
-!* 0. DECLARATIONS
-! ------------
+!* 0. DECLARATIONS
+! ------------
!
IMPLICIT NONE
!
-!* 0.1 Declaration of dummy arguments
+!* 0.1 Declaration of dummy arguments
!
INTEGER, INTENT(IN) :: KN ! Size of the result vector
INTEGER, INTENT(IN), DIMENSION(KN) :: KI ! Tabulated coordinate
@@ -3439,13 +3439,13 @@ REAL, INTENT(IN), DIMENSION(KN) :: PDX, PDY !
!
REAL, DIMENSION(KN) :: PY ! Interpolated value
!
-!* 0.2 Declaration of local variables
+!* 0.2 Declaration of local variables
!
INTEGER :: JJ ! Loop index
!
!
-!* 1. INTERPOLATION
-! -------------
+!* 1. INTERPOLATION
+! -------------
!
DO JJ = 1, KN
PY(JJ) = (1.0 - PDX(JJ)) * (1.0 - PDY(JJ)) * ZT(KI(JJ), KJ(JJ)) + &
diff --git a/src/PHYEX/micro/mode_ini_lima.f90 b/src/PHYEX/micro/mode_ini_lima.f90
index 3fba2eaab..3229e1077 100644
--- a/src/PHYEX/micro/mode_ini_lima.f90
+++ b/src/PHYEX/micro/mode_ini_lima.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2013-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2013-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -68,7 +68,7 @@ INTEGER :: JI
!
!
!* 1. INIT OUTPUT LISTING, COMPUTE KSPLITR AND KSPLITG
-! ------------------------------------------------
+! ------------------------------------------------
!
!
! Init output listing
@@ -145,7 +145,7 @@ XCEXVT = 0.4
!
!
!* 2. DEFINE SPECIES CHARACTERISTICS AND PROCESSES CONSTANTS
-! ------------------------------------------------------
+! ------------------------------------------------------
!
!
CALL INI_LIMA_WARM(PTSTEP, PDZMIN)
diff --git a/src/PHYEX/micro/mode_ini_lima_cold_mixed.f90 b/src/PHYEX/micro/mode_ini_lima_cold_mixed.f90
index 7b55d1a92..147c35eae 100644
--- a/src/PHYEX/micro/mode_ini_lima_cold_mixed.f90
+++ b/src/PHYEX/micro/mode_ini_lima_cold_mixed.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2013-2019 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2013-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -141,7 +141,7 @@ CALL PARAM_LIMA_MIXED_ASSOCIATE()
!
!
!* 1. CHARACTERISTICS OF THE SPECIES
-! ------------------------------
+! ------------------------------
!
!CALL RAIN_ICE_PARAM_ASSOCIATE()
!
@@ -302,7 +302,7 @@ XC1H = 1./2.
!
!
!* 2. DIMENSIONAL DISTRIBUTIONS OF THE SPECIES
-! ----------------------------------------
+! ----------------------------------------
!
!
!* 2.1 Ice, snow, graupel and hail distribution
@@ -401,7 +401,7 @@ XREFFI = (3*XAI/(2*ZRHOIW*10**3*XGAMMAI)*MOMG(XALPHAI,XNUI,XBI)/MOMG(XALPHAI,XNU
!
!
!* 3. CONSTANTS FOR THE SEDIMENTATION
-! -------------------------------
+! -------------------------------
!
!
!* 3.1 Exponent of the fall-speed air density correction
@@ -512,7 +512,7 @@ XFSEDC(7) = XCH*GAMMA_X0D(XNUH+XDH/XALPHAH)/GAMMA_X0D(XNUH)*
!
!
!* 4. CONSTANTS FOR HETEROGENEOUS NUCLEATION
-! --------------------------------------
+! --------------------------------------
!
!
! ***************
@@ -664,7 +664,7 @@ XBETA2 = 0.639
!
!
!* 5. CONSTANTS FOR THE SLOW COLD PROCESSES
-! -------------------------------------
+! -------------------------------------
!
!
!* 5.1.2 Constants for homogeneous nucleation from haze particules
@@ -764,7 +764,7 @@ XEX1DEPH = -0.5*(XDH+3.0)
!
!
!* 6. CONSTANTS FOR THE COALESCENCE PROCESSES
-! ---------------------------------------
+! ---------------------------------------
!
!
!* 6.0 Precalculation of the gamma function momentum
@@ -835,7 +835,7 @@ XAGGS_RLARGE2 = XKER_ZRNIC_A2*ZGAMI(7)*ZGAMS(2)*XAI
!
!
!* 7. CONSTANTS FOR THE FAST COLD PROCESSES FOR THE AGGREGATES
-! --------------------------------------------------------
+! --------------------------------------------------------
!
!
!* 7.1 Constants for the riming of the aggregates
@@ -1278,7 +1278,7 @@ XMOMGR_RDSF = MOMG(XALPHAR,XNUR,6.0+XDR)
!
!
!* 8. CONSTANTS FOR THE FAST COLD PROCESSES FOR THE GRAUPELN
-! --------------------------------------------------------
+! --------------------------------------------------------
!
!
!* 8.1 Constants for the rain contact freezing
@@ -1782,7 +1782,7 @@ END IF
!-------------------------------------------------------------------------------
!
!* 10. SET-UP RADIATIVE PARAMETERS
-! ---------------------------
+! ---------------------------
!
!
! R_eff_i = XFREFFI * (rho*r_i/N_i)**(1/3)
@@ -1793,7 +1793,7 @@ XFREFFI = 0.5 * ZGAMI(8) * (1.0/XLBI)**XLBEXI
!
!
!* 11. SOME PRINTS FOR CONTROL
-! -----------------------
+! -----------------------
!
!
!!$GFLAG = .TRUE.
diff --git a/src/PHYEX/micro/mode_ini_lima_warm.f90 b/src/PHYEX/micro/mode_ini_lima_warm.f90
index f3f8d60c1..7d50b3ff4 100644
--- a/src/PHYEX/micro/mode_ini_lima_warm.f90
+++ b/src/PHYEX/micro/mode_ini_lima_warm.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2013-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2013-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -84,7 +84,7 @@ REAL :: ZSMIN, ZSMAX ! Minimal and maximal supersaturation used to
CALL PARAM_LIMA_WARM_ASSOCIATE()
!
!* 1. CHARACTERISTICS OF THE SPECIES
-! ------------------------------
+! ------------------------------
!
!
!* 1.1 Cloud droplet characteristics
@@ -117,7 +117,7 @@ XF1R = 0.308
!
!
!* 2. DIMENSIONAL DISTRIBUTIONS OF THE SPECIES
-! ----------------------------------------
+! ----------------------------------------
!
!
!* 2.1 Cloud droplet distribution
@@ -169,7 +169,7 @@ END IF
!
!
!* 3. CONSTANTS FOR THE SEDIMENTATION
-! -------------------------------
+! -------------------------------
!
!
!* 4.1 Exponent of the fall-speed air density correction
@@ -207,7 +207,7 @@ XFSEDC(3) = XFSEDCR
!
!
!* 4. CONSTANTS FOR THE NUCLEATION PROCESS
-! ------------------------------------
+! ------------------------------------
!
!
XWMIN = 0.01 ! Minimal positive vertical velocity required
@@ -226,7 +226,7 @@ XCSTDCRIT = (XPI/6.)*XRHOLW*( (8.0*ZSURF_TEN )/( 3.0*XRV*XRHOLW ) )**3
!
!
!
-! 4.1 Tabulation of the hypergeometric functions in 'no units'
+! 4.1 Tabulation of the hypergeometric functions in 'no units'
! --------------------------------------------------------
!
! In LIMA's nucleation parameterization,
@@ -261,7 +261,7 @@ NAHEN = 81 ! Tabulation for each Kelvin degree in the range XTT-40 to XTT+40
XAHENINTP1 = 1.0
XAHENINTP2 = 0.5*REAL(NAHEN-1) - XTT
!
-! Compute the tabulation of function of T :
+! Compute the tabulation of function of T :
!
! 1
! XAHENG = -----------------------
@@ -312,7 +312,7 @@ XCONCR_PARAM_INI = (1.E7)**3/(XPI*XRHOLW) ! MP law with N_O=1.E7 m-1 is assumed
!
!
!* 5. CONSTANTS FOR THE COALESCENCE PROCESSES
-! ---------------------------------------
+! ---------------------------------------
!
!
!* 6.1 Csts for the coalescence processes
@@ -370,7 +370,7 @@ XSCBUEXP1 = -2500.0
!
!
!* 6. CONSTANTS FOR THE "SONTANEOUS" BREAK-UP
-! ---------------------------------------
+! ---------------------------------------
!
!
XSPONBUD1 = 3.0E-3
@@ -383,7 +383,7 @@ XSPONCOEF2 = ((XSPONBUD3/XSPONBUD2)**3 - 1.0)/(XSPONBUD3-XSPONBUD1)**2
!
!
!* 7. CONSTANTS FOR EVAPORATION PROCESS
-! ---------------------------------------
+! ---------------------------------------
!
!
X0CNDC = (4.0*XPI)*XC1C*XF0C*MOMG(XALPHAC,XNUC,1.)
@@ -416,7 +416,7 @@ XCEVAP = 0.86
!
!
!* 8. SET-UP RADIATIVE PARAMETERS
-! ---------------------------
+! ---------------------------
!
!
! R_eff_c = XFREFFC * (rho*r_c/N_c)**(1/3)
@@ -435,7 +435,7 @@ XCRER = 1.0/ (ZGAMR(6) * XAR**(2.0/3.0))
!
!
!* 9. SOME PRINTS FOR CONTROL
-! -----------------------
+! -----------------------
!
!
!!$GFLAG = .TRUE.
diff --git a/src/PHYEX/micro/mode_init_aerosol_properties.f90 b/src/PHYEX/micro/mode_init_aerosol_properties.f90
index d6920ff16..a3aa65c68 100644
--- a/src/PHYEX/micro/mode_init_aerosol_properties.f90
+++ b/src/PHYEX/micro/mode_init_aerosol_properties.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2013-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2013-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -216,7 +216,7 @@ END IF
END IF
!
!
-! Compute CCN spectra parameters from CCN characteristics
+! Compute CCN spectra parameters from CCN characteristics
!
!* INPUT : XBETAHEN_TEST is in 'percent' and XBETAHEN_MULTI in 'no units',
! XK... and XMU... are invariant
diff --git a/src/PHYEX/micro/mode_lima_ccn_activation.f90 b/src/PHYEX/micro/mode_lima_ccn_activation.f90
index 3a20104d3..9a4a943d3 100644
--- a/src/PHYEX/micro/mode_lima_ccn_activation.f90
+++ b/src/PHYEX/micro/mode_lima_ccn_activation.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2013-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2013-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -173,7 +173,7 @@ INTEGER :: IIB, IIE, IJB, IJE, IKB, IKE ! Physical domain
!ILUOUT = TLUOUT%NLU
!
!* 1. PREPARE COMPUTATIONS - PACK
-! ---------------------------
+! ---------------------------
!
IIB=1+JPHEXT
IIE=SIZE(PRHODREF,1) - JPHEXT
@@ -333,7 +333,7 @@ ELSE ! CCN activation from Cohard-Pinty
!
!
!* 2. compute the constant term (ZZW3) relative to smax
-! ----------------------------------------------------
+! ----------------------------------------------------
!
! Remark : in LIMA's nucleation parameterization, Smax=0.01 for a supersaturation of 1% !
!
@@ -406,7 +406,7 @@ ELSE ! CCN activation from Cohard-Pinty
!
!
!* 3. Compute the maximum of supersaturation
-! -----------------------------------------
+! -----------------------------------------
!
!
! estimate S_max for the CPB98 parameterization with SEVERAL aerosols mode
@@ -430,7 +430,7 @@ ELSE ! CCN activation from Cohard-Pinty
!
!
!* 4. Compute the nucleus source
-! -----------------------------
+! -----------------------------
!
!
! Again : Smax=0.01 for a 1% supersaturation
@@ -518,7 +518,7 @@ END IF ! AER_ACTIVATION
!
!
!* 5. Cleaning
-! -----------
+! -----------
!
!
DEALLOCATE(IVEC1)
@@ -583,7 +583,7 @@ END IF ! INUCT
!
!
!* 6. Functions used to compute the maximum of supersaturation
-! -----------------------------------------------------------
+! -----------------------------------------------------------
!
!
CONTAINS
diff --git a/src/PHYEX/micro/mode_lima_ccn_hom_freezing.f90 b/src/PHYEX/micro/mode_lima_ccn_hom_freezing.f90
index f3b9dc125..d87073fcb 100644
--- a/src/PHYEX/micro/mode_lima_ccn_hom_freezing.f90
+++ b/src/PHYEX/micro/mode_lima_ccn_hom_freezing.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2013-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2013-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -137,7 +137,7 @@ REAL :: ZEPS ! molar mass ratio
!
!
!* 1. Preliminary computations and packing
-! ------------------------------------
+! ------------------------------------
!
!
! Physical domain
@@ -234,7 +234,7 @@ IF (INEGT.GT.0) THEN
!
!
!* 2. Haze homogeneous freezing
-! ------------------------
+! ------------------------
!
!
! Compute the haze homogeneous nucleation source: RHHONI
diff --git a/src/PHYEX/micro/mode_lima_conversion_melting_snow.f90 b/src/PHYEX/micro/mode_lima_conversion_melting_snow.f90
index 32537b585..169de683b 100644
--- a/src/PHYEX/micro/mode_lima_conversion_melting_snow.f90
+++ b/src/PHYEX/micro/mode_lima_conversion_melting_snow.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2018-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2018-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -69,7 +69,7 @@ REAL, DIMENSION(SIZE(PRST)) :: ZW ! work arrays
!
!
!* 1. Conversion-melting of snow
-! --------------------------
+! --------------------------
!
!
P_RS_CMEL(:)=0.
diff --git a/src/PHYEX/micro/mode_lima_droplets_accretion.f90 b/src/PHYEX/micro/mode_lima_droplets_accretion.f90
index 8e0119a43..6863d1d0f 100644
--- a/src/PHYEX/micro/mode_lima_droplets_accretion.f90
+++ b/src/PHYEX/micro/mode_lima_droplets_accretion.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2018-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2018-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -73,7 +73,7 @@ LOGICAL, DIMENSION(SIZE(PRCT)) :: GACCR
!
!
!* 1. Accretion of cloud droplets on rain drops
-! --------------------------------------------
+! --------------------------------------------
!
P_RC_ACCR(:) = 0.0
P_CC_ACCR(:) = 0.0
diff --git a/src/PHYEX/micro/mode_lima_droplets_autoconversion.f90 b/src/PHYEX/micro/mode_lima_droplets_autoconversion.f90
index fca3fbf80..0f361b967 100644
--- a/src/PHYEX/micro/mode_lima_droplets_autoconversion.f90
+++ b/src/PHYEX/micro/mode_lima_droplets_autoconversion.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2018-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2018-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -77,7 +77,7 @@ IF (NMOM_C.EQ.1 .AND. LKESSLERAC) THEN
ELSE IF (LKHKO) THEN
!
! 1. Autoconversion of cloud droplets
-! -----------------------------------
+! -----------------------------------
!
WHERE ( PRCT(:)>XRTMIN(2) .AND. PCCT(:)>XCTMIN(2) .AND. LDCOMPUTE(:) )
!
@@ -96,7 +96,7 @@ ELSE IF (LKHKO) THEN
ELSE
!
! 2. Autoconversion of cloud droplets (Berry-Reinhardt parameterization)
-! ----------------------------------------------------------------------
+! ----------------------------------------------------------------------
!
WHERE( PRCT(:)>XRTMIN(2) .AND. PCCT(:)>XCTMIN(2) .AND. PLBDC(:)>0. .AND. LDCOMPUTE(:) )
ZW2(:) = MAX( 0.0, &
diff --git a/src/PHYEX/micro/mode_lima_droplets_hom_freezing.f90 b/src/PHYEX/micro/mode_lima_droplets_hom_freezing.f90
index 687e161ee..82b22d5b4 100644
--- a/src/PHYEX/micro/mode_lima_droplets_hom_freezing.f90
+++ b/src/PHYEX/micro/mode_lima_droplets_hom_freezing.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2013-2018 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2013-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -65,7 +65,7 @@ REAL, DIMENSION(SIZE(PT)) :: ZZW, ZZX, ZZY, ZTCELSIUS
!
!
!* 1. Cloud droplets homogeneous freezing
-! -----------------------------------
+! -----------------------------------
!
!
P_TH_HONC(:) = 0.
diff --git a/src/PHYEX/micro/mode_lima_droplets_self_collection.f90 b/src/PHYEX/micro/mode_lima_droplets_self_collection.f90
index 6a4557c89..81cc17aea 100644
--- a/src/PHYEX/micro/mode_lima_droplets_self_collection.f90
+++ b/src/PHYEX/micro/mode_lima_droplets_self_collection.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2018-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2018-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -58,7 +58,7 @@ REAL, DIMENSION(SIZE(PCCT)) :: ZW ! work arrays
!
!
!* 1. Cloud droplets self collection
-! ------------------------------
+! ------------------------------
!
!
P_CC_SELF(:)=0.
diff --git a/src/PHYEX/micro/mode_lima_drops_self_collection.f90 b/src/PHYEX/micro/mode_lima_drops_self_collection.f90
index 0c16b69b4..d43bc6ad2 100644
--- a/src/PHYEX/micro/mode_lima_drops_self_collection.f90
+++ b/src/PHYEX/micro/mode_lima_drops_self_collection.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2018-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2018-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -65,7 +65,7 @@ REAL, DIMENSION(SIZE(PCRT)) :: &
!
!
!* 1. Rain drops self-collection and break-up
-! ---------------------------------------
+! ---------------------------------------
!
!
P_CR_SCBU(:)=0.
diff --git a/src/PHYEX/micro/mode_lima_graupel.f90 b/src/PHYEX/micro/mode_lima_graupel.f90
index 84f56e43d..1a62ad431 100644
--- a/src/PHYEX/micro/mode_lima_graupel.f90
+++ b/src/PHYEX/micro/mode_lima_graupel.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2018-2019 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2018-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -254,10 +254,10 @@ WHERE( GDRY )
Z4(:) = GET_XKER_SDRYG(IVEC1(:) ,IVEC2(:) )
ZVEC3(:) = ( Z1(:)* ZVEC2(:) &
- Z2(:)*(ZVEC2(:) - 1.0) ) &
- * ZVEC1(:) &
+ * ZVEC1(:) &
- ( Z3(:)* ZVEC2(:) &
- Z4(:)*(ZVEC2(:) - 1.0) ) &
- * (ZVEC1(:) - 1.0)
+ * (ZVEC1(:) - 1.0)
ZZW(:) = ZVEC3(:)
ZZW3(:) = XFSDRYG * ZZW(:) * EXP( XCOLEXSG*(PT(:)-XTT) ) & ! RSDRYG - rs collected by graupel in dry mode
* PRST(:) * PCGT(:) &
@@ -272,10 +272,10 @@ WHERE( GDRY )
Z4(:) = GET_XKER_N_SDRYG(IVEC1(:) ,IVEC2(:) )
ZVEC3(:) = ( Z1(:)* ZVEC2(:) &
- Z2(:)*(ZVEC2(:) - 1.0) ) &
- * ZVEC1(:) &
+ * ZVEC1(:) &
- ( Z3(:)* ZVEC2(:) &
- Z4(:)*(ZVEC2(:) - 1.0) ) &
- * (ZVEC1(:) - 1.0)
+ * (ZVEC1(:) - 1.0)
ZZW(:) = ZVEC3(:)
ZZW3N(:) = XFNSDRYG * ZZW(:) * EXP( XCOLEXSG*(PT(:)-XTT) ) & ! NSDRYG - Ns collected by graupel in dry mode
* PCST(:) * PCGT(:) &
@@ -320,10 +320,10 @@ WHERE( GDRY )
Z4(:) = GET_XKER_RDRYG(IVEC1(:) ,IVEC2(:) )
ZVEC3(:) = ( Z1(:)* ZVEC2(:) &
- Z2(:)*(ZVEC2(:) - 1.0) ) &
- * ZVEC1(:) &
+ * ZVEC1(:) &
- ( Z3(:)* ZVEC2(:) &
- Z4(:)*(ZVEC2(:) - 1.0) ) &
- * (ZVEC1(:) - 1.0)
+ * (ZVEC1(:) - 1.0)
ZZW(:) = ZVEC3(:)
ZZW4(:) = XFRDRYG * ZZW(:) & ! RRDRYG
* PRRT(:) * PCGT(:) &
@@ -338,10 +338,10 @@ WHERE( GDRY )
Z4(:) = GET_XKER_N_RDRYG(IVEC1(:) ,IVEC2(:) )
ZVEC3(:) = ( Z1(:)* ZVEC2(:) &
- Z2(:)*(ZVEC2(:) - 1.0) ) &
- * ZVEC1(:) &
+ * ZVEC1(:) &
- ( Z3(:)* ZVEC2(:) &
- Z4(:)*(ZVEC2(:) - 1.0) ) &
- * (ZVEC1(:) - 1.0)
+ * (ZVEC1(:) - 1.0)
ZZW(:) = ZVEC3(:)
ZZW4N(:) = XFNRDRYG * ZZW(:) & ! NRDRYG
* PCRT(:) * PCGT(:) &
diff --git a/src/PHYEX/micro/mode_lima_graupel_deposition.f90 b/src/PHYEX/micro/mode_lima_graupel_deposition.f90
index 14e970084..d801f1d7a 100644
--- a/src/PHYEX/micro/mode_lima_graupel_deposition.f90
+++ b/src/PHYEX/micro/mode_lima_graupel_deposition.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2018-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2018-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -60,7 +60,7 @@ REAL, DIMENSION(:), INTENT(OUT) :: P_RG_DEPG
!
!
!* 1. Deposition of vapour on graupel
-! -------------------------------
+! -------------------------------
!
P_TH_DEPG(:) = 0.0
P_RG_DEPG(:) = 0.0
diff --git a/src/PHYEX/micro/mode_lima_hail.f90 b/src/PHYEX/micro/mode_lima_hail.f90
index 7c0231d89..f7e138bd2 100644
--- a/src/PHYEX/micro/mode_lima_hail.f90
+++ b/src/PHYEX/micro/mode_lima_hail.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2018-2019 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2018-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -228,10 +228,10 @@ WHERE( GWET )
Z4(:) = GET_XKER_SWETH(IVEC1(:) ,IVEC2(:) )
ZVEC3(:) = ( Z1(:)* ZVEC2(:) &
- Z2(:)*(ZVEC2(:) - 1.0) ) &
- * ZVEC1(:) &
+ * ZVEC1(:) &
- ( Z3(:)* ZVEC2(:) &
- Z4(:)*(ZVEC2(:) - 1.0) ) &
- * (ZVEC1(:) - 1.0)
+ * (ZVEC1(:) - 1.0)
ZZW(:) = ZVEC3(:)
ZZW3(:) = XFSWETH * ZZW(:) & ! RSWETH
* PRST(:) * PCHT(:) &
@@ -246,10 +246,10 @@ WHERE( GWET )
Z4(:) = GET_XKER_N_SWETH(IVEC1(:) ,IVEC2(:) )
ZVEC3(:) = ( Z1(:)* ZVEC2(:) &
- Z2(:)*(ZVEC2(:) - 1.0) ) &
- * ZVEC1(:) &
+ * ZVEC1(:) &
- ( Z3(:)* ZVEC2(:) &
- Z4(:)*(ZVEC2(:) - 1.0) ) &
- * (ZVEC1(:) - 1.0)
+ * (ZVEC1(:) - 1.0)
ZZW(:) = ZVEC3(:)
ZZW3N(:) = XFNSWETH * ZZW(:) & ! NSWETH
* PCST(:) * PCHT(:) &
@@ -295,10 +295,10 @@ WHERE( GWET )
Z4(:) = GET_XKER_GWETH(IVEC1(:) ,IVEC2(:) )
ZVEC3(:) = ( Z1(:)* ZVEC2(:) &
- Z2(:)*(ZVEC2(:) - 1.0) ) &
- * ZVEC1(:) &
+ * ZVEC1(:) &
- ( Z3(:)* ZVEC2(:) &
- Z4(:)*(ZVEC2(:) - 1.0) ) &
- * (ZVEC1(:) - 1.0)
+ * (ZVEC1(:) - 1.0)
ZZW(:) = ZVEC3(:)
ZZW4(:) = XFGWETH * ZZW(:) & ! RGWETH
* PRGT(:) * PCHT(:) &
@@ -313,10 +313,10 @@ WHERE( GWET )
Z4(:) = GET_XKER_N_GWETH(IVEC1(:) ,IVEC2(:) )
ZVEC3(:) = ( Z1(:)* ZVEC2(:) &
- Z2(:)*(ZVEC2(:) - 1.0) ) &
- * ZVEC1(:) &
+ * ZVEC1(:) &
- ( Z3(:)* ZVEC2(:) &
- Z4(:)*(ZVEC2(:) - 1.0) ) &
- * (ZVEC1(:) - 1.0)
+ * (ZVEC1(:) - 1.0)
ZZW(:) = ZVEC3(:)
ZZW4N(:) = XFNGWETH * ZZW(:) & ! NGWETH
* PCGT(:) * PCHT(:) &
diff --git a/src/PHYEX/micro/mode_lima_hail_deposition.f90 b/src/PHYEX/micro/mode_lima_hail_deposition.f90
index 50eae03d0..ee6e1656e 100644
--- a/src/PHYEX/micro/mode_lima_hail_deposition.f90
+++ b/src/PHYEX/micro/mode_lima_hail_deposition.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2018-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2018-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -60,7 +60,7 @@ REAL, DIMENSION(:), INTENT(INOUT) :: P_RH_DEPH
!
!
!* 1. Deposition of vapour on graupel
-! -------------------------------
+! -------------------------------
!
P_TH_DEPH(:) = 0.0
P_RH_DEPH(:) = 0.0
diff --git a/src/PHYEX/micro/mode_lima_init_ccn_activation_spectrum.f90 b/src/PHYEX/micro/mode_lima_init_ccn_activation_spectrum.f90
index a70652574..f00c023fe 100644
--- a/src/PHYEX/micro/mode_lima_init_ccn_activation_spectrum.f90
+++ b/src/PHYEX/micro/mode_lima_init_ccn_activation_spectrum.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2007-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2007-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -82,7 +82,7 @@ REAL, DIMENSION(1) :: XT ! temperature
!-------------------------------------------------------------------------------
!
!* 1. Select kappa value based on CTYPE_CCN
-! ---------------------------------
+! ---------------------------------
!
! Kappa values are from Petters and Kreidenweis (2007), table 1.
!
@@ -122,7 +122,7 @@ END DO
DO IJ=1, SIZE(XT)
!
!* 2. Compute Nccn(s) for several supersaturation values
-! --------------------------------------------------
+! --------------------------------------------------
!
! Get the value of Scrit at Ddry=0.1 micron
!
@@ -150,7 +150,7 @@ DO IJ=1, SIZE(XT)
!-------------------------------------------------------------------------------
!
!* 3. Compute C, k, mu, beta, using the Levenberg-Marquardt algorithm
-! ---------------------------------------------------------------
+! ---------------------------------------------------------------
!
PARAMS(1:3) = (/ 1., 1., 1000. /)
IFLAG = 1
@@ -166,7 +166,7 @@ END DO ! loop on temperatures
!-------------------------------------------------------------------------------
!
!* 6. Functions used to compute Scrit at Ddry=0.1 micron
-! --------------------------------------------------
+! --------------------------------------------------
!
CONTAINS
!
@@ -358,7 +358,7 @@ END FUNCTION DSDD
!-------------------------------------------------------------------------------
!
!* 7. Functions used to fit the CCN activation spectra with C s**k F()
-! ----------------------------------------------------------------
+! ----------------------------------------------------------------
!
SUBROUTINE DISTANCE(M,N,X,FVEC,IFLAG)
!!
diff --git a/src/PHYEX/micro/mode_lima_meyers_nucleation.f90 b/src/PHYEX/micro/mode_lima_meyers_nucleation.f90
index 989c3562c..95677ea24 100644
--- a/src/PHYEX/micro/mode_lima_meyers_nucleation.f90
+++ b/src/PHYEX/micro/mode_lima_meyers_nucleation.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2018-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2018-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -87,7 +87,7 @@ INTEGER :: JL ! Loop index
INTEGER :: INEGT ! Case number of nucleation
!
LOGICAL, DIMENSION(SIZE(PRHODREF,1),SIZE(PRHODREF,2),SIZE(PRHODREF,3)) &
- :: GNEGT ! Test where to compute the nucleation
+ :: GNEGT ! Test where to compute the nucleation
!
INTEGER, DIMENSION(SIZE(PRHODREF)) :: I1,I2,I3 ! Indexes for PACK replacement
!
@@ -127,7 +127,7 @@ REAL, DIMENSION(:), ALLOCATABLE :: ZTCELSIUS
!
!
!* 1. PRELIMINARY COMPUTATIONS
-! ------------------------
+! ------------------------
!
P_TH_HIND(:,:,:) = 0.
P_RI_HIND(:,:,:) = 0.
diff --git a/src/PHYEX/micro/mode_lima_phillips_ifn_nucleation.f90 b/src/PHYEX/micro/mode_lima_phillips_ifn_nucleation.f90
index 41b8b5762..f6126aee0 100644
--- a/src/PHYEX/micro/mode_lima_phillips_ifn_nucleation.f90
+++ b/src/PHYEX/micro/mode_lima_phillips_ifn_nucleation.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2018-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2018-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -127,7 +127,7 @@ INTEGER :: JL, JMOD_CCN, JMOD_IFN, JSPECIE, JMOD_IMM ! Loop index
INTEGER :: INEGT ! Case number of sedimentation, nucleation,
!
LOGICAL, DIMENSION(SIZE(PRHODREF,1),SIZE(PRHODREF,2),SIZE(PRHODREF,3)) &
- :: GNEGT ! Test where to compute the nucleation
+ :: GNEGT ! Test where to compute the nucleation
!
INTEGER, DIMENSION(SIZE(PRHODREF)) :: I1,I2,I3 ! Indexes for PACK replacement
!
@@ -173,7 +173,7 @@ REAL, DIMENSION(:), ALLOCATABLE :: ZTCELSIUS, ZZT_SI0_BC
!
!
!* 1. PRELIMINARY COMPUTATIONS
-! ------------------------
+! ------------------------
!
P_TH_HIND(:,:,:) = 0.
P_RI_HIND(:,:,:) = 0.
@@ -205,7 +205,7 @@ ZW(:,:,:) = PRVT(:,:,:)*( PPABST(:,:,:)-ZW(:,:,:) ) / ( (CST%XMV/CST%XMD) * ZW(:
!
!
!* 2. COMPUTATIONS ONLY WHERE NECESSARY : PACK
-! ----------------------------------------
+! ----------------------------------------
!
!
GNEGT(:,:,:) = .FALSE.
@@ -284,7 +284,7 @@ IF (INEGT > 0) THEN
!
!
!* 3. COMPUTE THE SATURATION OVER WATER AND ICE
-! -----------------------------------------
+! -----------------------------------------
!
!
ZTCELSIUS(:) = ZZT(:)-CST%XTT ! T [°C]
@@ -321,7 +321,7 @@ IF (INEGT > 0) THEN
!
!
!* 4. COMPUTE THE ACTIVABLE FRACTION OF EACH IFN SPECIE
-! -------------------------------------------------
+! -------------------------------------------------
!
!
! Computation of the reference activity spectrum ( ZZY = N_{IN,1,*} )
@@ -338,7 +338,7 @@ IF (INEGT > 0) THEN
!
!
!* 5. COMPUTE THE HETEROGENEOUS NUCLEATION OF INSOLUBLE IFN
-! -----------------------------------------------------
+! -----------------------------------------------------
!
!
!
@@ -379,7 +379,7 @@ IF (INEGT > 0) THEN
!
!
!* 6. COMPUTE THE HETEROGENEOUS NUCLEATION OF COATED IFN
-! --------------------------------------------------
+! --------------------------------------------------
!
!
! Heterogeneous nucleation by immersion of the activated CCN
@@ -426,7 +426,7 @@ IF (INEGT > 0) THEN
!
!
!* 7. CLEAN
-! -----
+! -----
!
!
DEALLOCATE(ZRVT)
diff --git a/src/PHYEX/micro/mode_lima_rain_accr_snow.f90 b/src/PHYEX/micro/mode_lima_rain_accr_snow.f90
index e7ee26f7e..fc54eea6f 100644
--- a/src/PHYEX/micro/mode_lima_rain_accr_snow.f90
+++ b/src/PHYEX/micro/mode_lima_rain_accr_snow.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2018-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2018-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -170,10 +170,10 @@ WHERE( GACC )
Z4(:) = GET_XKER_RACCSS(IVEC1(:) ,IVEC2(:) )
ZVEC3(:) = ( Z1(:)* ZVEC2(:) &
- Z2(:)*(ZVEC2(:) - 1.0) ) &
- * ZVEC1(:) &
+ * ZVEC1(:) &
- ( Z3(:)* ZVEC2(:) &
- Z4(:)*(ZVEC2(:) - 1.0) ) &
- * (ZVEC1(:) - 1.0)
+ * (ZVEC1(:) - 1.0)
ZZW1(:) = ZVEC3(:)
!
! 1.3.3b perform the bilinear interpolation of the normalized
@@ -185,10 +185,10 @@ WHERE( GACC )
Z4(:) = GET_XKER_N_RACCSS(IVEC1(:) ,IVEC2(:) )
ZVEC3(:) = ( Z1(:)* ZVEC2(:) &
- Z2(:)*(ZVEC2(:) - 1.0) ) &
- * ZVEC1(:) &
+ * ZVEC1(:) &
- ( Z3(:)* ZVEC2(:) &
- Z4(:)*(ZVEC2(:) - 1.0) ) &
- * (ZVEC1(:) - 1.0)
+ * (ZVEC1(:) - 1.0)
ZZWC1(:) = ZVEC3(:)
!
! 1.3.4 perform the bilinear interpolation of the normalized
@@ -236,10 +236,10 @@ WHERE( GACC )
Z4(:) = GET_XKER_SACCRG(IVEC2(:) ,IVEC1(:) )
ZVEC3(:) = ( Z1(:)* ZVEC1(:) &
- Z2(:)*(ZVEC1(:) - 1.0) ) &
- * ZVEC2(:) &
+ * ZVEC2(:) &
- ( Z3(:)* ZVEC1(:) &
- Z4(:)*(ZVEC1(:) - 1.0) ) &
- * (ZVEC2(:) - 1.0)
+ * (ZVEC2(:) - 1.0)
ZZW3(:) = ZVEC3(:)
!
! 1.3.5b perform the bilinear interpolation of the normalized
@@ -251,10 +251,10 @@ WHERE( GACC )
Z4(:) = GET_XKER_N_SACCRG(IVEC2(:) ,IVEC1(:) )
ZVEC3(:) = ( Z1(:)* ZVEC1(:) &
- Z2(:)*(ZVEC1(:) - 1.0) ) &
- * ZVEC2(:) &
+ * ZVEC2(:) &
- ( Z3(:)* ZVEC1(:) &
- Z4(:)*(ZVEC1(:) - 1.0) ) &
- * (ZVEC2(:) - 1.0)
+ * (ZVEC2(:) - 1.0)
ZZWC3(:) = ZVEC3(:)
!
! 1.3.4 raindrop accretion on the small sized aggregates
diff --git a/src/PHYEX/micro/mode_lima_rain_evaporation.f90 b/src/PHYEX/micro/mode_lima_rain_evaporation.f90
index 86e59a8d7..c0b08fe8f 100644
--- a/src/PHYEX/micro/mode_lima_rain_evaporation.f90
+++ b/src/PHYEX/micro/mode_lima_rain_evaporation.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2018-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2018-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -78,7 +78,7 @@ REAL, DIMENSION(SIZE(PRHODREF)) :: ZZW1, ZZW2
!
!
!* 1. PREPARE COMPUTATIONS - PACK
-! ---------------------------
+! ---------------------------
!
P_TH_EVAP(:) = 0.
P_RR_EVAP(:) = 0.
diff --git a/src/PHYEX/micro/mode_lima_rain_freezing.f90 b/src/PHYEX/micro/mode_lima_rain_freezing.f90
index b3bee2145..2c2c3d1ae 100644
--- a/src/PHYEX/micro/mode_lima_rain_freezing.f90
+++ b/src/PHYEX/micro/mode_lima_rain_freezing.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2018-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2018-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -69,7 +69,7 @@ REAL, DIMENSION(SIZE(PRRT)) :: ZW1, ZW2 ! work arrays
!
!
!* 1. PRELIMINARY COMPUTATIONS
-! ------------------------
+! ------------------------
!
!
P_TH_CFRZ(:)=0.
diff --git a/src/PHYEX/micro/mode_lima_snow_self_collection.f90 b/src/PHYEX/micro/mode_lima_snow_self_collection.f90
index 1a3bdb581..46cbc4abc 100644
--- a/src/PHYEX/micro/mode_lima_snow_self_collection.f90
+++ b/src/PHYEX/micro/mode_lima_snow_self_collection.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2018-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2018-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -68,7 +68,7 @@ REAL, DIMENSION(:), ALLOCATABLE :: ZVEC1, ZVEC3 ! Work vectors
!
!
!* 1. Snow self-collection and break-up
-! ---------------------------------
+! ---------------------------------
!
!
P_CS_SSC(:)=0.
diff --git a/src/PHYEX/micro/mode_nzcolx.f90 b/src/PHYEX/micro/mode_nzcolx.f90
index 5a3932bf4..bacf3bb5c 100644
--- a/src/PHYEX/micro/mode_nzcolx.f90
+++ b/src/PHYEX/micro/mode_nzcolx.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -101,13 +101,13 @@ INTEGER, INTENT(IN) :: KND ! Number of discrete size intervals in DX and DZ
!
!
REAL, INTENT(IN) :: PALPHAX ! First shape parameter of the specy X
- ! size distribution (generalized gamma law)
+ ! size distribution (generalized gamma law)
REAL, INTENT(IN) :: PNUX ! Second shape parameter of the specy X
- ! size distribution (generalized gamma law)
+ ! size distribution (generalized gamma law)
REAL, INTENT(IN) :: PALPHAZ ! First shape parameter of the specy Z
- ! size distribution (generalized gamma law)
+ ! size distribution (generalized gamma law)
REAL, INTENT(IN) :: PNUZ ! Second shape parameter of the specy Z
- ! size distribution (generalized gamma law)
+ ! size distribution (generalized gamma law)
REAL, INTENT(IN) :: PEXZ ! Efficiency of specy X collecting specy Z
REAL, INTENT(IN) :: PFALLX ! Fall speed constant of specy X
REAL, INTENT(IN) :: PEXFALLX ! Fall speed exponent of specy X
@@ -120,11 +120,11 @@ REAL, INTENT(IN) :: PLBDAZMAX ! Maximun slope of size distribution of specy Z
REAL, INTENT(IN) :: PLBDAXMIN ! Minimun slope of size distribution of specy X
REAL, INTENT(IN) :: PLBDAZMIN ! Minimun slope of size distribution of specy Z
REAL, INTENT(IN) :: PDINFTY ! Factor to define the largest diameter up to
- ! which the diameter integration is performed
+ ! which the diameter integration is performed
!
REAL, DIMENSION(:,:), INTENT(INOUT) :: PNZCOLX ! Scaled fall speed difference in
- ! the mass collection kernel as a
- ! function of LAMBDAX and LAMBDAZ
+ ! the mass collection kernel as a
+ ! function of LAMBDAX and LAMBDAZ
!
!
!* 0.2 Declarations of local variables
@@ -146,13 +146,13 @@ REAL :: ZDDZ ! Integration step of the diameter of specy Z
REAL :: ZDX ! Current diameter of the particle specy X
REAL :: ZDZ ! Current diameter of the particle specy Z
REAL :: ZCOLLZ ! Single integral of the mass weighted fall speed difference
- ! over the spectrum of specy Z
+ ! over the spectrum of specy Z
REAL :: ZCOLLXZ ! Double integral of the mass weighted fall speed difference
- ! over the spectra of specy X and specy Z
+ ! over the spectra of specy X and specy Z
REAL :: ZSCALZ ! Single integral of the scaling factor over
- ! the spectrum of specy Z
+ ! the spectrum of specy Z
REAL :: ZSCALXZ ! Double integral of the scaling factor over
- ! the spectra of specy X and specy Z
+ ! the spectra of specy X and specy Z
REAL :: ZFUNC ! Ancillary function
!
!
@@ -200,8 +200,8 @@ DO JLBDAX = 1,SIZE(PNZCOLX(:,:),1)
!* 1.6 Compute the normalization factor by integration over the
! dimensional spectrum of specy Z
!
- ZFUNC = (ZDX+ZDZ)**2 * GENERAL_GAMMA(PALPHAZ,PNUZ,ZLBDAZ,ZDZ)
- ZSCALZ = ZSCALZ + ZFUNC
+ ZFUNC = (ZDX+ZDZ)**2 * GENERAL_GAMMA(PALPHAZ,PNUZ,ZLBDAZ,ZDZ)
+ ZSCALZ = ZSCALZ + ZFUNC
!
!* 1.7 Compute the scaled fall speed difference by integration over
! the dimensional spectrum of specy Z
diff --git a/src/PHYEX/micro/mode_rzcolx.f90 b/src/PHYEX/micro/mode_rzcolx.f90
index f3ef46494..b2474ffb4 100644
--- a/src/PHYEX/micro/mode_rzcolx.f90
+++ b/src/PHYEX/micro/mode_rzcolx.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1995-2019 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1995-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -12,7 +12,7 @@ CONTAINS
SUBROUTINE RZCOLX( KND, PALPHAX, PNUX, PALPHAZ, PNUZ, &
PEXZ, PEXMASSZ, PFALLX, PEXFALLX, PFALLEXPX, &
PFALLZ, PEXFALLZ, PFALLEXPZ, &
- PLBDAXMAX, PLBDAZMAX, PLBDAXMIN, PLBDAZMIN, &
+ PLBDAXMAX, PLBDAZMAX, PLBDAXMIN, PLBDAZMIN, &
PDINFTY, PRZCOLX )
USE YOMHOOK , ONLY : LHOOK, DR_HOOK, JPHOOK
! ########################################################################
diff --git a/src/PHYEX/micro/prognos_lima.f90 b/src/PHYEX/micro/prognos_lima.f90
index aaba69596..87b423dfd 100644
--- a/src/PHYEX/micro/prognos_lima.f90
+++ b/src/PHYEX/micro/prognos_lima.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2012-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2012-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -140,7 +140,7 @@ INUCT = SIZE(PTT,1)
! 'ZCHEN_MULTI(1,1)=',ZCHEN_MULTI(1,1)
!
!* . Compute the nucleus source
-! -----------------------------
+! -----------------------------
!
!
! Modified values for Beta and C (see in init_aerosol_properties) account for that
diff --git a/src/PHYEX/micro/radar_rain_ice.f90 b/src/PHYEX/micro/radar_rain_ice.f90
index ae504233f..1850c78b9 100644
--- a/src/PHYEX/micro/radar_rain_ice.f90
+++ b/src/PHYEX/micro/radar_rain_ice.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1996-2019 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1996-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -161,18 +161,18 @@ REAL, DIMENSION(:,:,:), INTENT(IN),OPTIONAL :: PCHT ! hail concentration at
INTEGER :: IKB ! Coordinates of the first physical points along z
INTEGER :: IND ! Number of interval to integrate the kernels
REAL :: ZALPHA, ZNU, ZP ! Parameters to compute the value of the p_moment
- ! of the generalized Gamma function
+ ! of the generalized Gamma function
REAL :: ZDINFTY ! Factor used to define the "infinite" diameter
!
REAL :: ZCXR=-1.0 ! for rain N ~ 1/N_0
! (in Kessler parameterization)
REAL :: ZSLOPE, ZINTERCEPT, ZEXPONENT ! parameters defining the mean axis ratio
- ! functionnal
+ ! functionnal
REAL :: ZDMELT_FACT ! factor used to compute the equivalent
- ! melted diameter
+ ! melted diameter
REAL :: ZEQICE ! factor used to convert the ice crystals
- ! reflectivity into an equivalent liquid
- ! water reflectivity (from Smith, JCAM 84)
+ ! reflectivity into an equivalent liquid
+ ! water reflectivity (from Smith, JCAM 84)
REAL :: ZEXP ! anciliary parameter
REAL :: ZRHO00 ! Surface reference air density
!
@@ -205,7 +205,7 @@ REAL :: ZLB,ZLBEX , ZCC,ZCX,ZC,ZD
!
!
!* 1. FUNCTION STATEMENTS
-! -------------------
+! -------------------
!
!
!* 1.1 p_moment of the Generalized GAMMA function
diff --git a/src/PHYEX/micro/rain_c2r2_khko.f90 b/src/PHYEX/micro/rain_c2r2_khko.f90
index 6284e4f1f..cf61ff365 100644
--- a/src/PHYEX/micro/rain_c2r2_khko.f90
+++ b/src/PHYEX/micro/rain_c2r2_khko.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1996-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1996-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -424,7 +424,7 @@ TYPE(TFIELDMETADATA) :: TZFIELD
!-------------------------------------------------------------------------------
!
!* 1. COMPUTE THE SLOPE PARAMETERS ZLBDC,ZLBDR
-! ----------------------------------------
+! ----------------------------------------
!
CALL GET_INDICE_ll (IIB,IJB,IIE,IJE)
IKB=1+JPVEXT
@@ -458,7 +458,7 @@ ZT(:,:,:) = PTHT(:,:,:) * (PPABST(:,:,:)/XP00)**(XRD/XCPD)
!
!* 2. COMPUTES THE NUCLEATION PROCESS SOURCES
-! --------------------------------------
+! --------------------------------------
!
IF ((HACTCCN == 'ABRK').AND.((LORILAM).OR.(LDUST).OR.(LSALT))) THEN
CALL AER_NUCLEATION
@@ -518,7 +518,7 @@ IF (ORAIN) THEN
ENDIF
!-------------------------------------------------------------------------------
!* 6. COMPUTE THE SEDIMENTATION (RS) SOURCE
-! -------------------------------------
+! -------------------------------------
!
!* 6.1 Calculation of the mean volumic radius (ZRAY) and
! the terminal vertical velocity ZCC for precipitating clouds
@@ -527,7 +527,7 @@ ZTSPLITR = PTSTEP / REAL(KSPLITR) ! Small time step
!
!
!* 6.2 compute the sedimentation velocities for rain
-! --------------------------------------------
+! --------------------------------------------
!
ZMVRR(:,:,:) = 0.
ZVRR(:,:,:) = 0.
@@ -1689,7 +1689,7 @@ IF (HCLOUD == 'C2R2'.OR. HCLOUD=='C3R5') THEN
ELSE ! KHKO
!* correct negative values for rain
-! --------------------------------
+! --------------------------------
!
WHERE (PRRS(:,:,:)<0.)
PRCS(:,:,:) = PRCS(:,:,:)+PRRS(:,:,:)
diff --git a/src/PHYEX/micro/rain_ice_elec.f90 b/src/PHYEX/micro/rain_ice_elec.f90
index bdd4de762..4211fccb8 100644
--- a/src/PHYEX/micro/rain_ice_elec.f90
+++ b/src/PHYEX/micro/rain_ice_elec.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2002-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2002-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -541,7 +541,7 @@ REAL :: ZRHO00, ZCOR00 ! Surface reference air density
!-------------------------------------------------------------------------------
!
!* 1. COMPUTE THE LOOP BOUNDS
-! -----------------------
+! -----------------------
!
CALL GET_INDICE_ll (IIB,IJB,IIE,IJE)
IKB = 1 + JPVEXT
@@ -552,7 +552,7 @@ ZCOR00 = ZRHO00**XCEXVT
!
!
!* 2. COMPUTES THE SLOW COLD PROCESS SOURCES
-! --------------------------------------
+! --------------------------------------
!
!* 2.1 compute the ice nucleation
!
@@ -862,19 +862,19 @@ IF (IMICRO > 0) THEN
END IF
!
!
-!* 2.4 Initialization for the non-inductive charging process
+!* 2.4 Initialization for the non-inductive charging process
!
CALL ELEC_INI_NI_PROCESS
!
!
-!* 2.5 Compute the slow cold process sources
+!* 2.5 Compute the slow cold process sources
!
CALL RAIN_ICE_ELEC_SLOW
!
!-------------------------------------------------------------------------------
!
!* 3. COMPUTES THE SLOW WARM PROCESS SOURCES
-! --------------------------------------
+! --------------------------------------
!
IF( OWARM ) THEN ! Check if the formation of the raindrops by the slow
! warm processes is allowed
@@ -916,7 +916,7 @@ IF (IMICRO > 0) THEN
!-------------------------------------------------------------------------------
!
!* 8. UPDATE MIXING 3D RATIOS AND VOLUMETRIC CHARGE CONCENTRATIONS
-! ------------------------------------------------------------
+! ------------------------------------------------------------
!
!* 8.1 Update the mixing ratio
!
@@ -937,7 +937,7 @@ IF (IMICRO > 0) THEN
END IF
!
!
-!* 8.2 Compute the volumetric charge concentration
+!* 8.2 Compute the volumetric charge concentration
!
DO JL=1,IMICRO
PQPIS(I1(JL),I2(JL),I3(JL)) = ZQPIS(JL)
@@ -1076,7 +1076,7 @@ END IF
!-------------------------------------------------------------------------------
!
!* 8. COMPUTE THE SEDIMENTATION (RS) SOURCE
-! -------------------------------------
+! -------------------------------------
!
!* 8.1 time splitting loop initialization
!
@@ -2929,8 +2929,8 @@ IMPLICIT NONE
ZZW(:) = UNPACK( VECTOR=ZVEC1(:),MASK=GRIM,FIELD=0.0 )
!
!
-!* 5.1.6 perform the linear interpolation of the normalized
-!* "XFS"-moment of the incomplete gamma function
+!* 5.1.6 perform the linear interpolation of the normalized
+!* "XFS"-moment of the incomplete gamma function
!
ZVEC1(1:IGRIM) = XGAMINC_RIM3( IVEC2(1:IGRIM)+1 ) * ZVEC2(1:IGRIM) &
- XGAMINC_RIM3( IVEC2(1:IGRIM) ) * (ZVEC2(1:IGRIM) - 1.0)
@@ -3070,7 +3070,7 @@ IMPLICIT NONE
ZWQ1(:,5) = UNPACK( VECTOR=ZVECQ5(:), MASK=GACC, FIELD=0.0 )
!
! 5.2.4 raindrop accretion on the small sized aggregates:
-! RRACCSS & QRACCSS
+! RRACCSS & QRACCSS
!
WHERE ( GACC(:) )
ZZW1(:,2) = & !! coef of RRACCS
@@ -3469,7 +3469,7 @@ IMPLICIT NONE
ZWQ4(:) = UNPACK( VECTOR=ZVECQ6(:), MASK=GDRY, FIELD=0.0 ) ! Dvqsgmn if charge<0
ENDIF
!
-! 6.2.3.5 compute RSDRYG and QSDRYG = QSDRYG_coal + QSDRYG_boun
+! 6.2.3.5 compute RSDRYG and QSDRYG = QSDRYG_coal + QSDRYG_boun
!
WHERE( GDRY(:) )
ZZW1(:,3) = MIN( ZRSS(:),XFSDRYG*ZZW(:) & ! RSDRYG
@@ -4313,8 +4313,8 @@ IMPLICIT NONE
!
SUBROUTINE COMPUTE_LBDA(ZRR, ZRS, ZRG, ZRH)
!
-!* 0. DECLARATIONS
-! ------------
+!* 0. DECLARATIONS
+! ------------
!
IMPLICIT NONE
!
@@ -4357,12 +4357,12 @@ SUBROUTINE ELEC_UPDATE_QD(ZDUM, ZER, ZEI, ZES, ZEG, ZQR, ZQI, ZQS, ZQG, &
ZRR, ZRI, ZRS, ZRG, &
ZEH, ZQH, ZRH, ZEC, ZQC, ZRC)
!
-! Purpose : update the parameter e_x in the relation q_x = e_x d**f_x
+! Purpose : update the parameter e_x in the relation q_x = e_x d**f_x
! e_x = q_x/(N_x * M(f_x))
!
!
-!* 0. DECLARATIONS
-! ------------
+!* 0. DECLARATIONS
+! ------------
!
IMPLICIT NONE
!
@@ -4446,7 +4446,7 @@ END SUBROUTINE ELEC_UPDATE_QD
!
SUBROUTINE ELEC_INI_NI_PROCESS
!
-! Purpose : initialization for the non-inductive charging process
+! Purpose : initialization for the non-inductive charging process
!
! GELEC(:,1) : logical variable for Ice-Snow process --> ELEC_IAGGS_B
! from RAIN_ICE_ELEC_SLOW routine
@@ -4456,8 +4456,8 @@ END SUBROUTINE ELEC_UPDATE_QD
! from RAIN_ICE_ELEC_FAST_RG
!
!
-!* 0. DECLARATIONS
-! ------------
+!* 0. DECLARATIONS
+! ------------
!
IMPLICIT NONE
!
@@ -5381,7 +5381,7 @@ END SUBROUTINE ELEC_IAGGS_B
! Purpose : compute charge separation process during the dry collision
! between ice and graupeln
!
-!* 0. DECLARATIONS
+!* 0. DECLARATIONS
! ------------
!
IMPLICIT NONE
@@ -5557,7 +5557,7 @@ END SUBROUTINE ELEC_IDRYG_B
! Purpose : compute the charge separation during the dry collision
! between snow and graupeln
!
-!* 0. DECLARATIONS
+!* 0. DECLARATIONS
! ------------
!
IMPLICIT NONE
@@ -5819,12 +5819,12 @@ END SUBROUTINE INDUCTIVE_PROCESS
! | PDX |x1 |
! | |
!
-!* 0. DECLARATIONS
-! ------------
+!* 0. DECLARATIONS
+! ------------
!
IMPLICIT NONE
!
-!* 0.2 Declaration of local variables
+!* 0.2 Declaration of local variables
!
INTEGER, INTENT(IN) :: KN ! Size of the result vector
INTEGER, INTENT(IN), DIMENSION(KN) :: KI ! Tabulated coordinate
@@ -5835,8 +5835,8 @@ REAL, DIMENSION(KN) :: Y ! Interpolated value
!
INTEGER :: JJ ! Loop index
!
-!* 1. INTERPOLATION
-! -------------
+!* 1. INTERPOLATION
+! -------------
!
DO JJ = 1, KN
Y(JJ) = (1.0 - PDX(JJ)) * (1.0 - PDY(JJ)) * ZT(KI(JJ), KJ(JJ)) + &
diff --git a/src/PHYEX/turb/modd_turbn.f90 b/src/PHYEX/turb/modd_turbn.f90
index ef3ad4cbe..247795457 100644
--- a/src/PHYEX/turb/modd_turbn.f90
+++ b/src/PHYEX/turb/modd_turbn.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1995-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1995-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -25,7 +25,7 @@
!!
!! AUTHOR
!! ------
-!! J. Cuxart and J. Stein * I.N.M. and Meteo France*
+!! J. Cuxart and J. Stein * I.N.M. and Meteo France*
!!
!! MODIFICATIONS
!! -------------
diff --git a/src/PHYEX/turb/mode_ibm_mixinglength.f90 b/src/PHYEX/turb/mode_ibm_mixinglength.f90
index 608af2341..a6f55f1db 100644
--- a/src/PHYEX/turb/mode_ibm_mixinglength.f90
+++ b/src/PHYEX/turb/mode_ibm_mixinglength.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2019-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2019-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -32,7 +32,7 @@ SUBROUTINE IBM_MIXINGLENGTH(D,PLM,PLEPS,PMU,PHI,PTKE)
!
! AUTHOR
! ------
- !
+ !
! Franck Auguste * CERFACS(AE) *
!
! MODIFICATIONS
diff --git a/src/PHYEX/turb/mode_turb_ver.f90 b/src/PHYEX/turb/mode_turb_ver.f90
index 02c1a6750..5b24ef7e3 100644
--- a/src/PHYEX/turb/mode_turb_ver.f90
+++ b/src/PHYEX/turb/mode_turb_ver.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2022 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -37,12 +37,12 @@ SUBROUTINE TURB_VER(D,CST,CSTURB,TURBN,NEBN,TLES, &
!! -------
! The purpose of this routine is to compute the vertical turbulent
! fluxes of the evolutive variables and give back the source
-! terms to the main program. In the case of large horizontal meshes,
+! terms to the main program. In the case of large horizontal meshes,
! the divergence of these vertical turbulent fluxes represent the whole
! effect of the turbulence but when the three-dimensionnal version of
! the turbulence scheme is activated (CTURBDIM="3DIM"), these divergences
! are completed in the next routine TURB_HOR.
-! An arbitrary degree of implicitness has been implemented for the
+! An arbitrary degree of implicitness has been implemented for the
! temporal treatment of these diffusion terms.
! The vertical boundary conditions are as follows:
! * at the bottom, the surface fluxes are prescribed at the same
@@ -60,12 +60,12 @@ SUBROUTINE TURB_VER(D,CST,CSTURB,TURBN,NEBN,TLES, &
!! implicit scheme (a Crank-Nicholson type with coefficients different
!! than 0.5), which allows to vary the degree of implicitness of the
!! formulation.
-!! The different prognostic variables are treated one by one.
+!! The different prognostic variables are treated one by one.
!! The contributions of each turbulent fluxes are cumulated into the
!! tendency PRvarS, and into the dynamic and thermal production of
!! TKE if necessary.
!!
-!! In section 2 and 3, the thermodynamical fields are considered.
+!! In section 2 and 3, the thermodynamical fields are considered.
!! Only the turbulent fluxes of the conservative variables
!! (Thetal and Rnp stored in PRx(:,:,:,1)) are computed.
!! Note that the turbulent fluxes at the vertical
@@ -76,14 +76,14 @@ SUBROUTINE TURB_VER(D,CST,CSTURB,TURBN,NEBN,TLES, &
!! a function ETHETA or EMOIST, which preform the transformation from the
!! conservative variables to the virtual potential temperature.
!!
-!! In section 4, the variance of the statistical variable
+!! In section 4, the variance of the statistical variable
!! s indicating presence or not of condensation, is determined in function
!! of the turbulent moments of the conservative variables and its
!! squarred root is stored in PSIGS. This information will be completed in
!! the horizontal turbulence if the turbulence dimensionality is not
!! equal to "1DIM".
!!
-!! In section 5, the x component of the stress tensor is computed.
+!! In section 5, the x component of the stress tensor is computed.
!! The surface flux <u'w'> is computed from the value of the surface
!! fluxes computed in axes linked to the orography ( i", j" , k"):
!! i" is parallel to the surface and in the direction of the maximum
@@ -102,7 +102,7 @@ SUBROUTINE TURB_VER(D,CST,CSTURB,TURBN,NEBN,TLES, &
!! in the surface layer.
!!
!! In section 6, the same steps are repeated but for the y direction
-!! and in section 7, a diagnostic computation of the W variance is
+!! and in section 7, a diagnostic computation of the W variance is
!! performed.
!!
!! In section 8, the turbulent fluxes for the scalar variables are
@@ -112,11 +112,11 @@ SUBROUTINE TURB_VER(D,CST,CSTURB,TURBN,NEBN,TLES, &
!! EXTERNAL
!! --------
!! GX_U_M, GY_V_M, GZ_W_M : cartesian gradient operators
-!! GX_U_UW,GY_V_VW (X,Y,Z) represent the direction of the gradient
+!! GX_U_UW,GY_V_VW (X,Y,Z) represent the direction of the gradient
!! _(M,U,...)_ represent the localization of the
!! field to be derivated
!! _(M,UW,...) represent the localization of the
-!! field derivated
+!! field derivated
!!
!! SUBROUTINE TRIDIAG : to compute the split implicit evolution
!! of a variable located at a mass point
@@ -176,7 +176,7 @@ SUBROUTINE TURB_VER(D,CST,CSTURB,TURBN,NEBN,TLES, &
!! Modifications: October 10, 1995 (J.Cuxart and J. Stein)
!! Psi for scal var and LES tools
!! Modifications: November 10, 1995 (J. Stein)
-!! change the surface relations
+!! change the surface relations
!! Modifications: February 20, 1995 (J. Stein) optimization
!! Modifications: May 21, 1996 (J. Stein)
!! bug in the vertical flux of the V wind
diff --git a/src/PHYEX/turb/mode_turb_ver_dyn_flux.f90 b/src/PHYEX/turb/mode_turb_ver_dyn_flux.f90
index a5ed7627c..94398cc0f 100644
--- a/src/PHYEX/turb/mode_turb_ver_dyn_flux.f90
+++ b/src/PHYEX/turb/mode_turb_ver_dyn_flux.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -27,12 +27,12 @@ SUBROUTINE TURB_VER_DYN_FLUX(D,CST,CSTURB,TURBN,TLES,KSV,O2D,OFLAT, &
!! -------
! The purpose of this routine is to compute the vertical turbulent
! fluxes of the evolutive variables and give back the source
-! terms to the main program. In the case of large horizontal meshes,
+! terms to the main program. In the case of large horizontal meshes,
! the divergence of these vertical turbulent fluxes represent the whole
! effect of the turbulence but when the three-dimensionnal version of
! the turbulence scheme is activated (CTURBDIM="3DIM"), these divergences
! are completed in the next routine TURB_HOR.
-! An arbitrary degree of implicitness has been implemented for the
+! An arbitrary degree of implicitness has been implemented for the
! temporal treatment of these diffusion terms.
! The vertical boundary conditions are as follows:
! * at the bottom, the surface fluxes are prescribed at the same
@@ -50,12 +50,12 @@ SUBROUTINE TURB_VER_DYN_FLUX(D,CST,CSTURB,TURBN,TLES,KSV,O2D,OFLAT, &
!! implicit scheme (a Crank-Nicholson type with coefficients different
!! than 0.5), which allows to vary the degree of implicitness of the
!! formulation.
-!! The different prognostic variables are treated one by one.
+!! The different prognostic variables are treated one by one.
!! The contributions of each turbulent fluxes are cumulated into the
!! tendency PRvarS, and into the dynamic and thermal production of
!! TKE if necessary.
!!
-!! In section 2 and 3, the thermodynamical fields are considered.
+!! In section 2 and 3, the thermodynamical fields are considered.
!! Only the turbulent fluxes of the conservative variables
!! (Thetal and Rnp stored in PRx(:,:,:,1)) are computed.
!! Note that the turbulent fluxes at the vertical
@@ -66,14 +66,14 @@ SUBROUTINE TURB_VER_DYN_FLUX(D,CST,CSTURB,TURBN,TLES,KSV,O2D,OFLAT, &
!! a function ETHETA or EMOIST, which preform the transformation from the
!! conservative variables to the virtual potential temperature.
!!
-!! In section 4, the variance of the statistical variable
+!! In section 4, the variance of the statistical variable
!! s indicating presence or not of condensation, is determined in function
!! of the turbulent moments of the conservative variables and its
!! squarred root is stored in PSIGS. This information will be completed in
!! the horizontal turbulence if the turbulence dimensionality is not
!! equal to "1DIM".
!!
-!! In section 5, the x component of the stress tensor is computed.
+!! In section 5, the x component of the stress tensor is computed.
!! The surface flux <u'w'> is computed from the value of the surface
!! fluxes computed in axes linked to the orography ( i", j" , k"):
!! i" is parallel to the surface and in the direction of the maximum
@@ -92,7 +92,7 @@ SUBROUTINE TURB_VER_DYN_FLUX(D,CST,CSTURB,TURBN,TLES,KSV,O2D,OFLAT, &
!! in the surface layer.
!!
!! In section 6, the same steps are repeated but for the y direction
-!! and in section 7, a diagnostic computation of the W variance is
+!! and in section 7, a diagnostic computation of the W variance is
!! performed.
!!
!! In section 8, the turbulent fluxes for the scalar variables are
@@ -102,11 +102,11 @@ SUBROUTINE TURB_VER_DYN_FLUX(D,CST,CSTURB,TURBN,TLES,KSV,O2D,OFLAT, &
!! EXTERNAL
!! --------
!! GX_U_M, GY_V_M, GZ_W_M : cartesian gradient operators
-!! GX_U_UW,GY_V_VW (X,Y,Z) represent the direction of the gradient
+!! GX_U_UW,GY_V_VW (X,Y,Z) represent the direction of the gradient
!! _(M,U,...)_ represent the localization of the
!! field to be derivated
!! _(M,UW,...) represent the localization of the
-!! field derivated
+!! field derivated
!!
!!
!! MXM,MXF,MYM,MYF,MZM,MZF
@@ -170,7 +170,7 @@ SUBROUTINE TURB_VER_DYN_FLUX(D,CST,CSTURB,TURBN,TLES,KSV,O2D,OFLAT, &
!! Modifications: October 10, 1995 (J.Cuxart and J. Stein)
!! Psi for scal var and LES tools
!! Modifications: November 10, 1995 (J. Stein)
-!! change the surface relations
+!! change the surface relations
!! Modifications: February 20, 1995 (J. Stein) optimization
!! Modifications: May 21, 1996 (J. Stein)
!! bug in the vertical flux of the V wind
diff --git a/src/PHYEX/turb/mode_turb_ver_sv_flux.f90 b/src/PHYEX/turb/mode_turb_ver_sv_flux.f90
index b91bca2b6..516df39ea 100644
--- a/src/PHYEX/turb/mode_turb_ver_sv_flux.f90
+++ b/src/PHYEX/turb/mode_turb_ver_sv_flux.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2023 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -26,12 +26,12 @@ SUBROUTINE TURB_VER_SV_FLUX(D,CST,CSTURB,TURBN,TLES,ONOMIXLG, &
!! -------
! The purpose of this routine is to compute the vertical turbulent
! fluxes of the evolutive variables and give back the source
-! terms to the main program. In the case of large horizontal meshes,
+! terms to the main program. In the case of large horizontal meshes,
! the divergence of these vertical turbulent fluxes represent the whole
! effect of the turbulence but when the three-dimensionnal version of
! the turbulence scheme is activated (CTURBDIM="3DIM"), these divergences
! are completed in the next routine TURB_HOR.
-! An arbitrary degree of implicitness has been implemented for the
+! An arbitrary degree of implicitness has been implemented for the
! temporal treatment of these diffusion terms.
! The vertical boundary conditions are as follows:
! * at the bottom, the surface fluxes are prescribed at the same
@@ -49,12 +49,12 @@ SUBROUTINE TURB_VER_SV_FLUX(D,CST,CSTURB,TURBN,TLES,ONOMIXLG, &
!! implicit scheme (a Crank-Nicholson type with coefficients different
!! than 0.5), which allows to vary the degree of implicitness of the
!! formulation.
-!! The different prognostic variables are treated one by one.
+!! The different prognostic variables are treated one by one.
!! The contributions of each turbulent fluxes are cumulated into the
!! tendency PRvarS, and into the dynamic and thermal production of
!! TKE if necessary.
!!
-!! In section 2 and 3, the thermodynamical fields are considered.
+!! In section 2 and 3, the thermodynamical fields are considered.
!! Only the turbulent fluxes of the conservative variables
!! (Thetal and Rnp stored in PRx(:,:,:,1)) are computed.
!! Note that the turbulent fluxes at the vertical
@@ -65,14 +65,14 @@ SUBROUTINE TURB_VER_SV_FLUX(D,CST,CSTURB,TURBN,TLES,ONOMIXLG, &
!! a function ETHETA or EMOIST, which preform the transformation from the
!! conservative variables to the virtual potential temperature.
!!
-!! In section 4, the variance of the statistical variable
+!! In section 4, the variance of the statistical variable
!! s indicating presence or not of condensation, is determined in function
!! of the turbulent moments of the conservative variables and its
!! squarred root is stored in PSIGS. This information will be completed in
!! the horizontal turbulence if the turbulence dimensionality is not
!! equal to "1DIM".
!!
-!! In section 5, the x component of the stress tensor is computed.
+!! In section 5, the x component of the stress tensor is computed.
!! The surface flux <u'w'> is computed from the value of the surface
!! fluxes computed in axes linked to the orography ( i", j" , k"):
!! i" is parallel to the surface and in the direction of the maximum
@@ -91,7 +91,7 @@ SUBROUTINE TURB_VER_SV_FLUX(D,CST,CSTURB,TURBN,TLES,ONOMIXLG, &
!! in the surface layer.
!!
!! In section 6, the same steps are repeated but for the y direction
-!! and in section 7, a diagnostic computation of the W variance is
+!! and in section 7, a diagnostic computation of the W variance is
!! performed.
!!
!! In section 8, the turbulent fluxes for the scalar variables are
@@ -101,11 +101,11 @@ SUBROUTINE TURB_VER_SV_FLUX(D,CST,CSTURB,TURBN,TLES,ONOMIXLG, &
!! EXTERNAL
!! --------
!! GX_U_M, GY_V_M, GZ_W_M : cartesian gradient operators
-!! GX_U_UW,GY_V_VW (X,Y,Z) represent the direction of the gradient
+!! GX_U_UW,GY_V_VW (X,Y,Z) represent the direction of the gradient
!! _(M,U,...)_ represent the localization of the
!! field to be derivated
!! _(M,UW,...) represent the localization of the
-!! field derivated
+!! field derivated
!!
!!
!! MXM,MXF,MYM,MYF,MZM,MZF
@@ -172,7 +172,7 @@ SUBROUTINE TURB_VER_SV_FLUX(D,CST,CSTURB,TURBN,TLES,ONOMIXLG, &
!! Modifications: October 10, 1995 (J.Cuxart and J. Stein)
!! Psi for scal var and LES tools
!! Modifications: November 10, 1995 (J. Stein)
-!! change the surface relations
+!! change the surface relations
!! Modifications: February 20, 1995 (J. Stein) optimization
!! Modifications: May 21, 1996 (J. Stein)
!! bug in the vertical flux of the V wind
diff --git a/src/PHYEX/turb/mode_turb_ver_thermo_corr.f90 b/src/PHYEX/turb/mode_turb_ver_thermo_corr.f90
index d93fbfb87..5b8fd6474 100644
--- a/src/PHYEX/turb/mode_turb_ver_thermo_corr.f90
+++ b/src/PHYEX/turb/mode_turb_ver_thermo_corr.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2021 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -31,12 +31,12 @@ SUBROUTINE TURB_VER_THERMO_CORR(D,CST,CSTURB,TURBN,NEBN,TLES, &
!! -------
! The purpose of this routine is to compute the vertical turbulent
! fluxes of the evolutive variables and give back the source
-! terms to the main program. In the case of large horizontal meshes,
+! terms to the main program. In the case of large horizontal meshes,
! the divergence of these vertical turbulent fluxes represent the whole
! effect of the turbulence but when the three-dimensionnal version of
! the turbulence scheme is activated (CTURBDIM="3DIM"), these divergences
! are completed in the next routine TURB_HOR.
-! An arbitrary degree of implicitness has been implemented for the
+! An arbitrary degree of implicitness has been implemented for the
! temporal treatment of these diffusion terms.
! The vertical boundary conditions are as follows:
! * at the bottom, the surface fluxes are prescribed at the same
@@ -54,12 +54,12 @@ SUBROUTINE TURB_VER_THERMO_CORR(D,CST,CSTURB,TURBN,NEBN,TLES, &
!! implicit scheme (a Crank-Nicholson type with coefficients different
!! than 0.5), which allows to vary the degree of implicitness of the
!! formulation.
-!! The different prognostic variables are treated one by one.
+!! The different prognostic variables are treated one by one.
!! The contributions of each turbulent fluxes are cumulated into the
!! tendency PRvarS, and into the dynamic and thermal production of
!! TKE if necessary.
!!
-!! In section 2 and 3, the thermodynamical fields are considered.
+!! In section 2 and 3, the thermodynamical fields are considered.
!! Only the turbulent fluxes of the conservative variables
!! (Thetal and Rnp stored in PRx(:,:,:,1)) are computed.
!! Note that the turbulent fluxes at the vertical
@@ -70,14 +70,14 @@ SUBROUTINE TURB_VER_THERMO_CORR(D,CST,CSTURB,TURBN,NEBN,TLES, &
!! a function ETHETA or EMOIST, which preform the transformation from the
!! conservative variables to the virtual potential temperature.
!!
-!! In section 4, the variance of the statistical variable
+!! In section 4, the variance of the statistical variable
!! s indicating presence or not of condensation, is determined in function
!! of the turbulent moments of the conservative variables and its
!! squarred root is stored in PSIGS. This information will be completed in
!! the horizontal turbulence if the turbulence dimensionality is not
!! equal to "1DIM".
!!
-!! In section 5, the x component of the stress tensor is computed.
+!! In section 5, the x component of the stress tensor is computed.
!! The surface flux <u'w'> is computed from the value of the surface
!! fluxes computed in axes linked to the orography ( i", j" , k"):
!! i" is parallel to the surface and in the direction of the maximum
@@ -96,7 +96,7 @@ SUBROUTINE TURB_VER_THERMO_CORR(D,CST,CSTURB,TURBN,NEBN,TLES, &
!! in the surface layer.
!!
!! In section 6, the same steps are repeated but for the y direction
-!! and in section 7, a diagnostic computation of the W variance is
+!! and in section 7, a diagnostic computation of the W variance is
!! performed.
!!
!! In section 8, the turbulent fluxes for the scalar variables are
@@ -106,11 +106,11 @@ SUBROUTINE TURB_VER_THERMO_CORR(D,CST,CSTURB,TURBN,NEBN,TLES, &
!! EXTERNAL
!! --------
!! GX_U_M, GY_V_M, GZ_W_M : cartesian gradient operators
-!! GX_U_UW,GY_V_VW (X,Y,Z) represent the direction of the gradient
+!! GX_U_UW,GY_V_VW (X,Y,Z) represent the direction of the gradient
!! _(M,U,...)_ represent the localization of the
!! field to be derivated
!! _(M,UW,...) represent the localization of the
-!! field derivated
+!! field derivated
!!
!!
!! MXM,MXF,MYM,MYF,MZM,MZF
@@ -171,7 +171,7 @@ SUBROUTINE TURB_VER_THERMO_CORR(D,CST,CSTURB,TURBN,NEBN,TLES, &
!! Modifications: October 10, 1995 (J.Cuxart and J. Stein)
!! Psi for scal var and LES tools
!! Modifications: November 10, 1995 (J. Stein)
-!! change the surface relations
+!! change the surface relations
!! Modifications: February 20, 1995 (J. Stein) optimization
!! Modifications: May 21, 1996 (J. Stein)
!! bug in the vertical flux of the V wind
diff --git a/src/PHYEX/turb/mode_turb_ver_thermo_flux.f90 b/src/PHYEX/turb/mode_turb_ver_thermo_flux.f90
index e895283b0..0493471bd 100644
--- a/src/PHYEX/turb/mode_turb_ver_thermo_flux.f90
+++ b/src/PHYEX/turb/mode_turb_ver_thermo_flux.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2022 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -34,12 +34,12 @@ SUBROUTINE TURB_VER_THERMO_FLUX(D,CST,CSTURB,TURBN,TLES, &
!! -------
! The purpose of this routine is to compute the vertical turbulent
! fluxes of the evolutive variables and give back the source
-! terms to the main program. In the case of large horizontal meshes,
+! terms to the main program. In the case of large horizontal meshes,
! the divergence of these vertical turbulent fluxes represent the whole
! effect of the turbulence but when the three-dimensionnal version of
! the turbulence scheme is activated (CTURBDIM="3DIM"), these divergences
! are completed in the next routine TURB_HOR.
-! An arbitrary degree of implicitness has been implemented for the
+! An arbitrary degree of implicitness has been implemented for the
! temporal treatment of these diffusion terms.
! The vertical boundary conditions are as follows:
! * at the bottom, the surface fluxes are prescribed at the same
@@ -57,12 +57,12 @@ SUBROUTINE TURB_VER_THERMO_FLUX(D,CST,CSTURB,TURBN,TLES, &
!! implicit scheme (a Crank-Nicholson type with coefficients different
!! than 0.5), which allows to vary the degree of implicitness of the
!! formulation.
-!! The different prognostic variables are treated one by one.
+!! The different prognostic variables are treated one by one.
!! The contributions of each turbulent fluxes are cumulated into the
!! tendency PRvarS, and into the dynamic and thermal production of
!! TKE if necessary.
!!
-!! In section 2 and 3, the thermodynamical fields are considered.
+!! In section 2 and 3, the thermodynamical fields are considered.
!! Only the turbulent fluxes of the conservative variables
!! (Thetal and Rnp stored in PRx(:,:,:,1)) are computed.
!! Note that the turbulent fluxes at the vertical
@@ -73,14 +73,14 @@ SUBROUTINE TURB_VER_THERMO_FLUX(D,CST,CSTURB,TURBN,TLES, &
!! a function ETHETA or EMOIST, which preform the transformation from the
!! conservative variables to the virtual potential temperature.
!!
-!! In section 4, the variance of the statistical variable
+!! In section 4, the variance of the statistical variable
!! s indicating presence or not of condensation, is determined in function
!! of the turbulent moments of the conservative variables and its
!! squarred root is stored in PSIGS. This information will be completed in
!! the horizontal turbulence if the turbulence dimensionality is not
!! equal to "1DIM".
!!
-!! In section 5, the x component of the stress tensor is computed.
+!! In section 5, the x component of the stress tensor is computed.
!! The surface flux <u'w'> is computed from the value of the surface
!! fluxes computed in axes linked to the orography ( i", j" , k"):
!! i" is parallel to the surface and in the direction of the maximum
@@ -99,7 +99,7 @@ SUBROUTINE TURB_VER_THERMO_FLUX(D,CST,CSTURB,TURBN,TLES, &
!! in the surface layer.
!!
!! In section 6, the same steps are repeated but for the y direction
-!! and in section 7, a diagnostic computation of the W variance is
+!! and in section 7, a diagnostic computation of the W variance is
!! performed.
!!
!! In section 8, the turbulent fluxes for the scalar variables are
@@ -109,11 +109,11 @@ SUBROUTINE TURB_VER_THERMO_FLUX(D,CST,CSTURB,TURBN,TLES, &
!! EXTERNAL
!! --------
!! GX_U_M, GY_V_M, GZ_W_M : cartesian gradient operators
-!! GX_U_UW,GY_V_VW (X,Y,Z) represent the direction of the gradient
+!! GX_U_UW,GY_V_VW (X,Y,Z) represent the direction of the gradient
!! _(M,U,...)_ represent the localization of the
!! field to be derivated
!! _(M,UW,...) represent the localization of the
-!! field derivated
+!! field derivated
!!
!!
!! MXM,MXF,MYM,MYF,MZM,MZF
@@ -180,7 +180,7 @@ SUBROUTINE TURB_VER_THERMO_FLUX(D,CST,CSTURB,TURBN,TLES, &
!! Modifications: October 10, 1995 (J.Cuxart and J. Stein)
!! Psi for scal var and LES tools
!! Modifications: November 10, 1995 (J. Stein)
-!! change the surface relations
+!! change the surface relations
!! Modifications: February 20, 1995 (J. Stein) optimization
!! Modifications: May 21, 1996 (J. Stein)
!! bug in the vertical flux of the V wind
diff --git a/src/PHYEX/turb/mode_update_lm.f90 b/src/PHYEX/turb/mode_update_lm.f90
index 5f3509398..dc39aea63 100644
--- a/src/PHYEX/turb/mode_update_lm.f90
+++ b/src/PHYEX/turb/mode_update_lm.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2006-2019 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2006-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
@@ -29,7 +29,7 @@ SUBROUTINE UPDATE_LM(D,HLBCX,HLBCY,PLM,PLEPS)
!!
!! AUTHOR
!! ------
-!! V. Masson * Meteo France *
+!! V. Masson * Meteo France *
!!
!! MODIFICATIONS
!! -------------
--
GitLab