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

Timestamp:

Oct 3, 2017, 4:45:26 PM (7 years ago)

Author:

Laurent Fairhead

Message:

Modifications to the code and xml files to output Ap and B, the coefficients
of the hybrid coordinates as requested by the CMIP6 DataRequest?
LF (with guidance from A. Caubel and S. Senesi)

File:

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

LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -29 +29 @@
     USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
     USE time_phylmdz_mod, ONLY: day_step_phy, start_time, itau_phy
+    USE vertical_layers_mod, ONLY : ap, bp, aps, bps
     USE phys_output_ctrlout_mod, ONLY: o_phis, o_aire, is_ter, is_lic, is_oce, &
+         o_Ahyb, o_Bhyb,o_Ahyb_inter, o_Bhyb_inter, & 
          is_ave, is_sic, o_contfracATM, o_contfracOR, &
          o_aireTER, o_flat, o_slp, o_ptstar, o_pt0, o_tsol, &
@@ -392 +394 @@
     REAL, PARAMETER :: dobson_u = 2.1415e-05 ! Dobson unit, in kg m-2
 !   REAL, PARAMETER :: missing_val=nf90_fill_real
+    REAL, DIMENSION(klev,2) :: Ahyb_bounds, Bhyb_bounds
+    INTEGER :: ilev
 #ifndef CPP_XIOS
     REAL :: missing_val
@@ -418 +422 @@
        iinitend = 1
     ENDIF
+
+    DO ilev=1,klev
+      Ahyb_bounds(ilev,1) = ap(ilev)
+      Ahyb_bounds(ilev,2) = ap(ilev+1)
+      Bhyb_bounds(ilev,1) = bp(ilev)
+      Bhyb_bounds(ilev,2) = bp(ilev+1)
+    END DO
 
 #ifdef CPP_XIOS
@@ -531 +542 @@
        CALL histwrite_phy("R_incl",R_incl)
        CALL histwrite_phy("solaire",solaire)
+       CALL histwrite_phy(o_Ahyb, aps)
+       CALL histwrite_phy(o_Bhyb, bps)
+       CALL histwrite_phy(o_Ahyb_inter, Ahyb_bounds)
+       CALL histwrite_phy(o_Bhyb_inter, Bhyb_bounds)
 !
 #ifdef CPP_RRTM