Changeset 5813 for LMDZ6/trunk/libf/phylmdiso

Timestamp:

Sep 17, 2025, 10:23:08 AM (2 months ago)

Author:

musat

Message:

Ajout iflag_tropo pour choisir le calcul de la tropopause
pour l'ozone:
iflag_tropo=0 : tropopause "dynamique"=> valeur par defaut pour

retrocompatibilite

iflag_tropo=1 : tropopause "lapse-rate" comme pour les aerosols

stratospheriques ==> valeur recommandee

Ionela Musat

File:

• LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

@@ -192 +192 @@
        ! [Variables internes non sauvegardees de la physique]
        ! Variables locales pour effectuer les appels en serie
+       p_tropopause, dyntropo, &
        t_seri,q_seri,ql_seri,qs_seri,qbs_seri, &
        u_seri,v_seri,cf_seri,rvc_seri,tr_seri, &
@@ -2992 +2993 @@
           IF(SIZE(time_climoz)==360.AND..NOT.ok_daily_climoz) ro3i=ro3i*360./year_len
           IF(adjust_tropopause) THEN
+             CALL stratosphere_mask(missing_val, pphis, t_seri, pplay, latitude_deg)
+             IF(iflag_tropo.EQ.0) THEN
+               dyntropo=dyn_tropopause(t_seri, ztsol, paprs, pplay, rot)
+             ELSEIF(iflag_tropo.EQ.1) THEN
+               dyntropo=p_tropopause
+             ENDIF
              CALL regr_pr_time_av(ncid_climoz, vars_climoz(1:read_climoz),   &
                       ro3i, 'C', press_cen_climoz, pplay, wo, paprs(:,1),    &
                       time_climoz ,  longitude_deg,   latitude_deg,          &
-                      dyn_tropopause(t_seri, ztsol, paprs, pplay, rot))
+                      dyntropo)
           ELSE
              CALL regr_pr_time_av(ncid_climoz,  vars_climoz(1:read_climoz),  &