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

Timestamp:

Jan 21, 2022, 3:50:54 PM (2 years ago)

Author:

oboucher

Message:

Audran Borella's parametrisation for ice supersaturation
activated with flag_ice_sursat (FALSE by default)

File:

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

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -209 +209 @@
          o_p_tropopause, o_z_tropopause, o_t_tropopause,  &
          o_col_O3_strato, o_col_O3_tropo,                 & 
+!--aviation & supersaturation
+         o_oclr, o_ocld, o_oss, o_ovc, o_rnebss, o_rnebclr, o_rnebseri, o_gammass, &
+         o_N1_ss, o_N2_ss, o_qsatl, o_qsats, &
+         o_drnebsub, o_drnebcon, o_drnebtur, o_drnebavi, o_flight_m, o_flight_h2o, &
+         o_Tcontr, o_qcontr, o_qcontr2, o_fcontrN, o_fcontrP, &
 !--interactive CO2
          o_flx_co2_ocean, o_flx_co2_ocean_cor, &
@@ -230 +235 @@
          o_vsed_aer, o_tau_strat_1020, o_ext_strat_1020, o_f_r_wet
 #endif
+
+    USE ice_sursat_mod, ONLY: flight_m, flight_h2o
     
     USE phys_output_ctrlout_mod, ONLY: o_heat_volc, o_cool_volc !NL
@@ -280 +287 @@
          cldq, flwp, fiwp, ue, ve, uq, vq, &
          uwat, vwat, &
+         rneb_seri, d_rneb_dyn, &
          plcl, plfc, wbeff, convoccur, upwd, dnwd, dnwd0, prw, prlw, prsw, &
          s_pblh, s_pblt, s_lcl, s_therm, uwriteSTD, &
@@ -293 +301 @@
          wdtrainA, wdtrainS, wdtrainM, n2, s2, proba_notrig, &
          random_notrig, &
+         qclr, qcld, qss, qvc, rnebclr, rnebss, gamma_ss, &
+         N1_ss, N2_ss, zqsatl, zqsats, &
+         Tcontr, qcontr, qcontr2, fcontrN, fcontrP, &
+         drneb_sub, drneb_con, drneb_tur, drneb_avi, &
          alp_bl_det, alp_bl_fluct_m, alp_bl_conv, &
          alp_bl_stat, alp_bl_fluct_tke, slab_wfbils, &
@@ -1813 +1825 @@
        ENDIF
 
+!--aviation & supersaturation
+       CALL histwrite_phy(o_oclr, qclr)
+       CALL histwrite_phy(o_ocld, qcld)
+       CALL histwrite_phy(o_oss, qss)
+       CALL histwrite_phy(o_ovc, qvc)
+       CALL histwrite_phy(o_rnebclr, rnebclr)
+       CALL histwrite_phy(o_rnebss, rnebss)
+       CALL histwrite_phy(o_rnebseri, rneb_seri)
+       CALL histwrite_phy(o_gammass, gamma_ss)
+       CALL histwrite_phy(o_N1_ss, N1_ss)
+       CALL histwrite_phy(o_N2_ss, N2_ss)
+       CALL histwrite_phy(o_drnebsub, drneb_sub)
+       CALL histwrite_phy(o_drnebcon, drneb_con)
+       CALL histwrite_phy(o_drnebtur, drneb_tur)
+       CALL histwrite_phy(o_drnebavi, drneb_avi)
+       CALL histwrite_phy(o_qsatl, zqsatl)
+       CALL histwrite_phy(o_qsats, zqsats)
+       CALL histwrite_phy(o_flight_m, flight_m)
+       CALL histwrite_phy(o_flight_h2o, flight_h2o)
+       CALL histwrite_phy(o_Tcontr, Tcontr)
+       CALL histwrite_phy(o_qcontr, qcontr)
+       CALL histwrite_phy(o_qcontr2, qcontr2)
+       CALL histwrite_phy(o_fcontrN, fcontrN)
+       CALL histwrite_phy(o_fcontrP, fcontrP)
        
        IF (vars_defined) zx_tmp_fi3d = wo(:, :, 1) * dobson_u * 1e3 / zmasse / rmo3 * rmd