diff --git a/src/PHYEX/aux/mode_gradient_m_phy.f90 b/src/PHYEX/aux/mode_gradient_m_phy.f90
index aeb7f5d1d9c6ac052b5cfc51f62f707724dabb81..7a9d731e7ad8e15db20967830eba3d0ae5000677 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 89eef4a6023da4a4be1ea9a11f9d65aff426d3e6..e07a11db245758b1110f595a5aefaf03ccfb0982 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 69f793542be2f229603988822277fb711da29aed..42d8c9c3b6afcbc6f9966bbb1ebb20684ac80896 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 f94e9bf86341b1d7b8de792727d31604e5f55bfd..f347b97e2c56c59ae7a168a717494cdb79b23f48 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 177c8844bf722acf4b469bcbd2eb9aeee705745c..6abe76527513b4bb4ea7ecf1441f08ed61fbf408 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 91ddb2cbea511be6af5e1cc5ebe70958cfc4e790..8ada5dc1950b0d3ac348ac7b44d3d269561148de 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 73e59decbc9163924fe6aca076815cb285aef41c..9c787e74fa1aa9cc36d421b224c7d18ff8140c19 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 814455eb95f94099750bf83dbd882d87df97df4a..fa5cdbb2d88be5ef80aa43c1719bf6e7933ae669 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 ebf25507e71c3fb2d20f512948222501e62f3267..64949022f12364e94971235ba68e882bd356d429 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 1c32f4d4d1d4e3ce7b92c5064a730a1d7d8b8267..1165b4d750ad9738517d935394bca416f4733324 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 39ee732ff6bf91f529afe10b94b2024d5b1fa03f..a794a4bddb60932d41ec59c63c9d5f57a1902382 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 d99b011a296f0b6bd44b70c381c1a66cfc863831..22d69b2b0d27cbe0c693f06fe9709edb17603160 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 718da2b7bb108aeea6fbe32d9390bf34795cd636..e6750f4f2e421dcca290f6bc38f5969b5d53b549 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 c7e53eb86a025033c25684fe2e2db4d65c485f23..d21666562378d6b4f727b723b66a46311f74ab9e 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 19cf4725f1b4b3596a950ba43bb2e70b7dfb367e..aeba527082b251694c9448c358449d5b00f581ef 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 56abd00baf8c72a0cbe8896ee4574191a2b0bf19..d9e22715a171bd1c0e30444269f23773fd2d3dac 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 548e8bda5982cc59c3920f2241d94443a820de74..c82bf823f495309621ee3aae7dc6fbed8300f7ed 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 330ba0bbecd501e68e4703a3178a2550fb497a6b..783117451f94417023678152e876281c38e8ce11 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 43d7c194f0aec823f993e9946451827722d5fc16..e69453e14f6c301f04d95e09281006d21fb701d3 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 b144572c51801be9f045437ebf8323b2a33b8994..f8a82bc413ca88d3b41feb9750df1b2a0e794399 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 0989dc52dc8c64679eaaa4a718b81f54614f8f85..6b95d1fa4e590dc7a264e09c200db0cbaa961fa4 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 1a7fa798b429365740460f18850a2215c1a69bb3..18f437cd296c59bbff19f67f0d4a4a9667be44d3 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 b9be5a7aca535ad2b53341b25f47060c37f61b4e..d3a76799fcd804fcff1242aad384d39b5cddef85 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 4971d2215a15ab383b56047f7b555abd1defbdfe..2d9f4046359d6cfdb93f3ee7e6424c08c3d7b68c 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 364a7007f2c98d705bcfcfbb3d328b084f75678b..017c946ef3442f7d650b86653dbc595743cc74d1 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 ca1c22825215509ad1cffa4097e0b5145f69ac23..3b0d3f88ea27896471b31cc28fade61689112bff 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 251521427ddb0bcf6772c2bf1b0741e99dc82c6d..d928830de8fe55be010d5b3b1aef4f880eb13c93 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 47b59eb8066528c40a966b9b7e176aab7649e133..703af76523ba318abf2b7bca7f8dec4e186d5ce3 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 a2fb7d4ed358b40b14e2f32f95a528cf1a346fb3..f2265404656bcf908715f8244ad6fd3d6c489c97 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 78d415912605c585c35ffe75c47985db47629aa2..36059d58d2bc8811ef820ee64d133bb169c75315 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 caf74ddc2eea96bde0625f5047ee105d252222de..0b723d7327f5bad12af239ae030d5eb7fae69a2f 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 bac9042e65e5d46f2e3109f54f59b833e1d59d7e..9855d95153bc0c0e61f5dbf28a21eb387b46448a 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 7d374b464bc55aa26354582183ed144ebcc6bcd2..c2bfb1daf1727478c7e178efc373340af3b5b8a4 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 fc012c8f10967b765e86f625c4254681f5ad578b..3708626bfa7e7d72398305b4fe303f49d9411d9c 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 fd8f2297149c9eab2daf277cac9f4bf9baecfefa..4d0749feab90357be3f896a76d96e1e2c8cc6f60 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 db6aaa9587502e2c80bc9f6eb9be30cad41dade3..83b3442e34980dd4a44eb6ae938ed34ffb30be42 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 128c292e931d918d3e805371a080dfad70fbfe3f..130dbbbb24a5c6af6f31cf06f89ebe24c5f027fa 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 cc1e68531da53b478c04ce72ffdfc50a9bbea151..a1a0401c5a4e0e24d2c698d7c7dc966370fe7a5d 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 20f080c031a0494c221ec46b43f882902c6d620e..19506c9c9e29bf568b4237aec50fe011368390b3 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 1a903452b9a90973eb754f80e9a830f719875686..eeb3f9084165e2895e4681a4981c5547fdaa8da8 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 cc3e3adc1e5639286c37689f9569905947dc8501..e5ddccd64e76b6b3bf6bdaca08d49e625171909b 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 cabbebbd7f02b319e10869d5efa4bcfdc44772f4..652465f37c3444ab286ee7fcb443feeb4f49e36d 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 639abf3c138f37f8cc00f4c8524defdd4fee9d93..2e14cf96269d4b540eb1614ebdda3eedaac01729 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 af2a62376f1ebd296ca6c3b46f9d88c76e7a74ce..1cf805a41ecd2459c6cb490710c97fe28a646c59 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 3fba2eaabb5078c0426457590e71941afb8e1e3b..3229e1077c4f1069e4f67ead4349a3bb5c217f8c 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 7b55d1a92f675dbf860a9e904565656a4ba94808..147c35eae9430ec828542bec44482ca29d16783c 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 f3f8d60c1c12bee4bdaac51849dc01d9fa834d22..7d50b3ff43760b2504200881d3ca5f98cd00d928 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 d6920ff167d96a8bd1ef1b10ccccf9b50f7ff21a..a3aa65c682deb382228549e83dc4d0e141327f7f 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 3a20104d354d7229ea38f467a81db7bb10297e8a..9a4a943d375846e89de494dbbaa1741a5e3d90d7 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 f3b9dc125e3d7a7282392b769be65ea4950f0c7d..d87073fcb77ccafb58f44178b935070283a974e4 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 32537b585e8663dafafcc2499113e27100a67722..169de683be55bf4443183030f7c87b9795e4a540 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 8e0119a4380f264f99c14f7c3e76426f049c175e..6863d1d0f8d37f6e1b37e8c76e55332678a03884 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 fca3fbf804de1e6b60bc4fd54a9073a4b04eb369..0f361b967211f4192ef6babdfcf66b4cdb7bf2d2 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 687e161ee34a2f9c138f4510f06f81b4e15035ba..82b22d5b40c57498ff971a776b2bbe5a740e277b 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 6a4557c89610cb9cafac0c37c4c639926e4b9a90..81cc17aeacbade14f0a240a180cc38ffbe82ce9a 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 0c16b69b4d41c53ba69667de98731a79df209c38..d43bc6ad25d7d60c73902d30123444e136e5cf7c 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 84f56e43dc2e839cde629f0c423da950d8b6064c..1a62ad431ee2a01c494e55342d0f3992a9712b36 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 14e970084e17abc54a598dd2720b5e1e84156e84..d801f1d7aa207b31242b08ba47680b3510705f67 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 7c0231d89ae203c89aaedf4e014bb86e85b232ce..f7e138bd25ed9f259cb09b01fbb49b2b8f803c71 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 50eae03d08804182d5da09f0eefcf8fc30bed701..ee6e1656e3a37b3aafd86f6f59615a74aa2a0ab8 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 a706525749d8ca10a8ea1ab46481cc58927d683c..f00c023fe55fd8206af27facb75b0dedadd839d9 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 989c3562cbcac6400404794c726c1c6f7237fe2d..95677ea24f284ebb4e045624af581414d78830c1 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 41b8b5762cf9f83d9fbe762a21061fcfff7133ae..f6126aee035ad21014cf449c0bbec8a22c4c78ba 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 e7ee26f7e5bccf8d571eae3b7639011b2fe16042..fc54eea6fad2914322901b0a42a34e76bd5d69dc 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 86e59a8d7778eba2ba80b4263a1d9edea2980a3f..c0b08fe8fe25727edd9685c055b209eef8763862 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 b3bee2145b48b566fd0673e947a804fe0f37a8d4..2c2c3d1aeb2d0e8ca0d0a1b9533740b72fc3f220 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 1a3bdb581aedd1a875c96665989f8f6c1747a220..46cbc4abc5232880100e2344141b031f33dbdb07 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 5a3932bf461543d85ba4226d0c1d02531cd242fe..bacf3bb5ce979fd2ff48a6ef82d9e04c6cf9ab28 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 f3ef464948915356ab5f1884b0a0af8e18e604ef..b2474ffb490f4164efeebd762871ff3a42aa4a95 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 aaba695968f56ce9cff37a79ed2e5689e8b4f2fa..87b423dfdb412402bfa4fdf6c5b606ba94bfea6c 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 ae504233f2bf27919460e236fe796abf2021b6de..1850c78b9c5f9d9e9b658f4727380cd75167f8b8 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 6284e4f1fbbc787d4305cf1811b216db85401634..cf61ff365d540f1e7172099eb21b1870f4e5855b 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 bdd4de76281fe67a4fc3285f3f8e7ea957e2adf8..4211fccb83f60aa3efe7cd629d8377b56c1e8542 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 ef3ad4cbe1da3b08d6e7aca87a58f73b0079519e..2477954570186115dab7af4bbe5a29b81415af41 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 608af23415bcba280ec196b3d60deda7e21367cb..a6f55f1db1cee1d1eaec7104b19ebd172013096e 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 02c1a67500ade26748dec2fd8dfc661962931e26..5b24ef7e3c578e5e4b12fa94280d92bca05429f2 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 a5ed7627c1583b15f95ae79d8d771bcb716b8680..94398cc0f0faf976f29826cd02c103c754c94382 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 b91bca2b68fe7b2debfaa13368868f1b37e6b1e1..516df39ea9f54f62129b6c78e3bb2c35669ce403 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 d93fbfb87a71e076307a436dd14912812c083167..5b8fd64749487aa7df05244c40f64b478016c081 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 e895283b0fe679202672555b72d8bd15b05a8245..0493471bd9e62528d328176365a0771e9df5a2a7 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 5f3509398aa633335d26156660f8107ccbf351e1..dc39aea63adf36f78d5683f566cd9186901ae59a 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
!! -------------