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

Timestamp:

Nov 18, 2015, 12:25:20 PM (9 years ago)

Author:

jyg

Message:

Add various intializations of arrays in lmdz1d.F90
and in the convection scheme. Add output variables
for boundary layer splitting.

File:

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

LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -71 +71 @@
          o_tSTDlevs, epsfra, o_t_oce_sic, &
          o_ale_bl, o_alp_bl, o_ale_wk, o_alp_wk, &
+         o_dtvdf_x    , o_dtvdf_w    , o_dqvdf_x    , o_dqvdf_w    , &
+         o_sens_x     , o_sens_w     , o_flat_x     , o_flat_w     , &
+         o_delta_tsurf, &
+         o_cdragh_x   , o_cdragh_w   , o_cdragm_x   , o_cdragm_w   , &
+         o_kh         , o_kh_x       , o_kh_w       , &
          o_ale, o_alp, o_cin, o_WAPE, o_wake_h, &
          o_wake_s, o_wake_deltat, o_wake_deltaq, &
@@ -167 +172 @@
          sollwdownclr, lwdn0, ftsol, ustar, u10m, &
          v10m, pbl_tke, wake_delta_pbl_TKE, &
+         delta_tsurf, &
          wstar, cape, ema_pcb, ema_pct, &
          ema_cbmf, Ma, fm_therm, ale_bl, alp_bl, ale, &
@@ -197 +203 @@
          s_pblh, s_pblt, s_lcl, s_therm, uwriteSTD, &
          vwriteSTD, wwriteSTD, phiwriteSTD, qwriteSTD, &
-         twriteSTD, ale_wake, alp_wake, wake_h, &
+         twriteSTD, ale_wake, alp_wake, &
+         dtvdf_x    ,dtvdf_w    ,dqvdf_x    ,dqvdf_w    , &
+         sens_x     ,sens_w     ,zxfluxlat_x,zxfluxlat_w, &
+         cdragh_x   ,cdragh_w   ,cdragm_x   ,cdragm_w   , &
+         kh         ,kh_x       ,kh_w       , &
+         wake_h, &
          wake_omg, d_t_wake, d_q_wake, Vprecip, &
          wdtrainA, wdtrainM, n2, s2, proba_notrig, &
@@ -754 +765 @@
              CALL histwrite_phy(o_ale_wk, ale_wake)
              CALL histwrite_phy(o_alp_wk, alp_wake)
+             IF (iflag_pbl_split>=1) THEN
+               IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=dtvdf_x(1:klon,1:klev)/pdtphys
+               CALL histwrite_phy(o_dtvdf_x    ,zx_tmp_fi3d) 
+               IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=dtvdf_w(1:klon,1:klev)/pdtphys
+               CALL histwrite_phy(o_dtvdf_w    ,zx_tmp_fi3d) 
+               IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=dqvdf_x(1:klon,1:klev)/pdtphys
+               CALL histwrite_phy(o_dqvdf_x    ,zx_tmp_fi3d) 
+               IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=dqvdf_w(1:klon,1:klev)/pdtphys
+               CALL histwrite_phy(o_dqvdf_w    ,zx_tmp_fi3d) 
+               CALL histwrite_phy(o_sens_x     ,sens_x     ) 
+               CALL histwrite_phy(o_sens_w     ,sens_w     ) 
+               CALL histwrite_phy(o_flat_x     ,zxfluxlat_x) 
+               CALL histwrite_phy(o_flat_w     ,zxfluxlat_w) 
+               CALL histwrite_phy(o_delta_tsurf,delta_tsurf) 
+               CALL histwrite_phy(o_cdragh_x   ,cdragh_x   ) 
+               CALL histwrite_phy(o_cdragh_w   ,cdragh_w   ) 
+               CALL histwrite_phy(o_cdragm_x   ,cdragm_x   ) 
+               CALL histwrite_phy(o_cdragm_w   ,cdragm_w   ) 
+               CALL histwrite_phy(o_kh         ,kh         ) 
+               CALL histwrite_phy(o_kh_x       ,kh_x       ) 
+               CALL histwrite_phy(o_kh_w       ,kh_w       ) 
+             ENDIF   ! (iflag_pbl_split>=1)
              CALL histwrite_phy(o_ale, ale)
              CALL histwrite_phy(o_alp, alp)