Changeset 2992 for LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

Timestamp:

Sep 21, 2017, 11:51:25 AM (7 years ago)

Author:

oboucher

Message:

Adding p, z, and t of tropopause in output files.
Stratosphere_mask is streamlined a bit.
This includes a treatment of missing values.

File:

• LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90 (modified) (6 diffs)

LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -184 +184 @@
          o_map_emis_Anv, o_map_pcld_Anv, o_map_tcld_Anv, &
          o_map_ntot, o_map_hc,o_map_hist,o_map_Cb,o_map_ThCi,o_map_Anv, &
-         o_alt_tropo
+         o_alt_tropo, & 
+! Tropopause
+         o_p_tropopause, o_z_tropopause, o_t_tropopause
 
 #ifdef CPP_StratAer
@@ -196 +198 @@
          o_budg_ocs_to_so2, o_budg_so2_to_h2so4, o_budg_h2so4_to_part, &
          o_surf_PM25_sulf, o_ext_strat_550, o_tau_strat_550, &
-         o_p_tropopause, o_vsed_aer, o_tau_strat_1020, o_ext_strat_1020, o_f_r_wet
+         o_vsed_aer, o_tau_strat_1020, o_ext_strat_1020, o_f_r_wet
 #endif
 
@@ -287 +289 @@
          east_gwstress, west_gwstress, &
          d_q_ch4, pmfd, pmfu, ref_liq, ref_ice, rhwriteSTD, &
-         ep, epmax_diag ! epmax_cape
+         ep, epmax_diag, &  ! epmax_cape
+         p_tropopause, t_tropopause, z_tropopause
 
 #ifdef CPP_StratAer 
@@ -298 +301 @@
          budg_emi_ocs, budg_emi_so2, budg_emi_h2so4, budg_emi_part, & 
          budg_ocs_to_so2, budg_so2_to_h2so4, budg_h2so4_to_part, & 
-         surf_PM25_sulf, tau_strat_550, p_tropopause, tausum_strat, &
+         surf_PM25_sulf, tau_strat_550, tausum_strat, &
          vsed_aer, tau_strat_1020, f_r_wet
 #endif
@@ -1220 +1223 @@
           ENDIF
        ENDIF
+
+       CALL histwrite_phy(o_p_tropopause, p_tropopause)
+       CALL histwrite_phy(o_t_tropopause, t_tropopause)
+       CALL histwrite_phy(o_z_tropopause, z_tropopause)
+
 #ifdef CPP_StratAer
        IF (type_trac=='coag') THEN
@@ -1244 +1252 @@
           CALL histwrite_phy(o_budg_sed_part,        budg_sed_part)
           CALL histwrite_phy(o_surf_PM25_sulf, surf_PM25_sulf)
-          CALL histwrite_phy(o_p_tropopause, p_tropopause)
           CALL histwrite_phy(o_vsed_aer, vsed_aer)
           CALL histwrite_phy(o_f_r_wet, f_r_wet)