Changeset 5825 for LMDZ6

Timestamp:

Sep 22, 2025, 3:06:26 PM (2 months ago)

Author:

musat

Message:

Changement du nom du flag iflag_tropo qui devient iflag_tropopause_height;
et de la valeur défaut à 1 :
iflag_tropopause_height=0 : PV tropopause
iflag_tropopause_height=1 ! Lapse-rate tropopause (valeur par défaut)
Ionela Musat

Location:

LMDZ6/trunk/libf/phylmd

Files:

• clesphys_mod_h.f90 (modified) (3 diffs)
• physiq_mod.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/clesphys_mod_h.f90

@@ -51 +51 @@
           , iflag_thermals, nsplit_thermals              &
           , iflag_physiq, ok_3Deffect, ok_water_mass_fixer             &
-          , ok_mass_dtcon, ok_mass_dqcon, ok_mass_duvcon, iflag_tropo
+          , ok_mass_dtcon, ok_mass_dqcon, ok_mass_duvcon, iflag_tropopause_height
 
 
@@ -158 +158 @@
 
   ! flag pour le choix du calcul de tropopause pour l'ozone
-  INTEGER :: iflag_tropo
+  INTEGER :: iflag_tropopause_height
 
 
@@ -208 +208 @@
   !$OMP      , iflag_thermals, nsplit_thermals              &
   !$OMP      , iflag_physiq, ok_3Deffect, ok_water_mass_fixer             &
-  !$OMP      , ok_mass_dtcon, ok_mass_dqcon, ok_mass_duvcon, iflag_tropo )
+  !$OMP      , ok_mass_dtcon, ok_mass_dqcon, ok_mass_duvcon, iflag_tropopause_height )
 
 END MODULE clesphys_mod_h

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -1521 +1521 @@
        ok_water_mass_fixer=.FALSE.  ! OB: by default we do not apply the mass fixer
        CALL getin_p('ok_water_mass_fixer',ok_water_mass_fixer)
-       iflag_tropo = 0 ! IM: iflag_tropo = 0 => PV tropopause; iflag_tropo=1 => LR tropopause
-       CALL getin_p('iflag_tropo', iflag_tropo)
+       iflag_tropopause_height = 1 ! IM: iflag_tropopause_height = 0 => PV tropopause; iflag_tropopause_height=1 => LR tropopause
+       CALL getin_p('iflag_tropopause_height', iflag_tropopause_height)
 IF (CPPKEY_DUST) THEN
        IF (iflag_phytrac.EQ.0) THEN
@@ -2688 +2688 @@
           IF(adjust_tropopause) THEN
              CALL stratosphere_mask(missing_val, pphis, t_seri, pplay, latitude_deg)
-             IF(iflag_tropo.EQ.0) THEN
+             IF(iflag_tropopause_height.EQ.0) THEN
                dyntropo=dyn_tropopause(t_seri, ztsol, paprs, pplay, rot)
-             ELSEIF(iflag_tropo.EQ.1) THEN
+             ELSEIF(iflag_tropopause_height.EQ.1) THEN
                dyntropo=p_tropopause
              ENDIF