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

Timestamp:

Jun 16, 2014, 4:33:38 PM (10 years ago)

Author:

Ehouarn Millour

Message:

Preparatory stuff to fix and improve the slab ocean model.
FC

File:

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

LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -24 +24 @@
     ! defined and initialised in phys_output_mod.F90
 
-    USE dimphy, only: klon, klev, klevp1
+    USE dimphy, only: klon, klev, klevp1, nslay
     USE control_mod, only: day_step, iphysiq
     USE phys_output_ctrlout_mod, only: o_phis, o_aire, is_ter, is_lic, is_oce, &
@@ -77 +77 @@
          o_ale_bl_trig, o_alp_bl_det, &
          o_alp_bl_fluct_m, o_alp_bl_fluct_tke, &
-         o_alp_bl_conv, o_alp_bl_stat, o_slab_bils, &
+         o_alp_bl_conv, o_alp_bl_stat, &
+         o_slab_qflux, o_tslab, o_slab_bils, &
          o_weakinv, o_dthmin, o_cldtau, &
          o_cldemi, o_pr_con_l, o_pr_con_i, &
@@ -219 +220 @@
          bils_ec,bils_ech, bils_tke, bils_kinetic, bils_latent, bils_enthalp, &
          itau_con, nfiles, clef_files, nid_files, zvstr_gwd_rando
+    USE ocean_slab_mod, only: tslab, slab_bils
     USE indice_sol_mod, only: nbsrf
     USE infotrac, only: nqtot
@@ -726 +728 @@
        ENDIF  !(iflag_clos_bl>=1)
 !!! fin nrlmd le 10/04/2012
+       ! Output of slab ocean variables
        IF (type_ocean=='slab ') THEN
-          CALL histwrite_phy(o_slab_bils, slab_wfbils)
+          CALL histwrite_phy(o_slab_qflux, slab_wfbils)
+          CALL histwrite_phy(o_slab_bils, slab_bils)
+          IF (nslay.EQ.1) THEN
+              zx_tmp_fi2d(:)=tslab(:,1)
+              CALL histwrite_phy(o_tslab, zx_tmp_fi2d)
+          ELSE
+              CALL histwrite_phy(o_tslab, tslab)
+          END IF
        ENDIF !type_ocean == force/slab
        CALL histwrite_phy(o_weakinv, weak_inversion)