Changeset 5614 for LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

Timestamp:

Apr 14, 2025, 9:21:07 PM (8 weeks ago)

Author:

evignon

Message:

Commission liée à un update majeur de la routine de condensation grande echelle suite au travail
de Lea, Audran et Etienne
Elle inclue une restructuration des routines pour clarifier le role "moniteur" de la routine lscp_main,
une mise à jour de la parametrisation de partitionnement de phase de Lea pour inclure les nuages de couche limite,
ainsi que des corrections des routines de precipitations "poprecip".
Convergence numerique verifiee en prod et debug pour les physiques NPv6.3 et 7.0.1c

File:

• LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -145 +145 @@
          o_zfull, o_zhalf, o_rneb, o_rnebjn, o_rnebcon, &
          o_rnebls, o_rneblsvol, o_rhum, o_rhl, o_rhi, o_ozone, o_ozone_light, &
+         o_distcltop, o_tempcltop,  &
          o_pfraclr, o_pfracld, o_cldfraliq, o_sigma2_icefracturb, o_mean_icefracturb,  &
+         o_cldfraliqth, o_sigma2_icefracturbth, o_mean_icefracturbth,  & 
          o_qrainlsc, o_qsnowlsc, o_dqreva, o_dqrauto, o_dqrcol, o_dqrmelt, o_dqrfreez, &
          o_dqssub, o_dqsauto, o_dqsagg, o_dqsrim, o_dqsmelt, o_dqsfreez, &
@@ -269 +271 @@
          o_SAD_sulfate, o_reff_sulfate, o_sulfmmr, o_nd_mode, o_sulfmmr_mode,o_SO2_chlm
 
-    USE lmdz_lscp_ini, ONLY: ok_poprecip
+    USE lmdz_lscp_ini, ONLY: ok_poprecip, iflag_icefrac
 
     USE phys_output_ctrlout_mod, ONLY: o_heat_volc, o_cool_volc !NL
@@ -377 +379 @@
          zphi, u_seri, v_seri, omega, cldfra, &
          rneb, rnebjn, rneblsvol,  &
-         zx_rh, zx_rhl, zx_rhi, &
+         zx_rh, zx_rhl, zx_rhi, distcltop, temp_cltop, &
          pfraclr, pfracld, cldfraliq, sigma2_icefracturb, mean_icefracturb, &
+         cldfraliqth, sigma2_icefracturbth, mean_icefracturbth, &
          qraindiag, qsnowdiag, dqreva, dqssub, &
          dqrauto,dqrcol,dqrmelt,dqrfreez, &
@@ -2087 +2090 @@
            CALL histwrite_phy(o_pfraclr, pfraclr)
            CALL histwrite_phy(o_pfracld, pfracld)
+           IF (iflag_icefrac .GT. 0) THEN
            CALL histwrite_phy(o_cldfraliq, cldfraliq)
            CALL histwrite_phy(o_sigma2_icefracturb, sigma2_icefracturb)
            CALL histwrite_phy(o_mean_icefracturb, mean_icefracturb)
+           CALL histwrite_phy(o_cldfraliqth, cldfraliqth)
+           CALL histwrite_phy(o_sigma2_icefracturbth, sigma2_icefracturbth)
+           CALL histwrite_phy(o_mean_icefracturbth, mean_icefracturbth)
+           ELSE
+           CALL histwrite_phy(o_distcltop, distcltop)
+           CALL histwrite_phy(o_tempcltop, temp_cltop)
+           ENDIF
            IF (ok_poprecip) THEN
            CALL histwrite_phy(o_qrainlsc, qraindiag)