Changeset 3853 for LMDZ6/branches/IPSL-CM6A-MR/libf/phylmd

Timestamp:

Feb 23, 2021, 4:01:34 PM (4 years ago)

Author:

lmdz-users

Message:

Modifications to code to introduce a new vertical axis on llm+1 levels so that fluxes can
be output on their proper levels as per the CMIP6 DataRequest?
LF

Location:

LMDZ6/branches/IPSL-CM6A-MR/libf/phylmd

Files:

• iophy.F90 (modified) (1 diff)
• phys_output_mod.F90 (modified) (2 diffs)

LMDZ6/branches/IPSL-CM6A-MR/libf/phylmd/iophy.F90

@@ -1227 +1227 @@
 
     nlev=SIZE(field,2)
-    IF (nlev.EQ.klev+1) THEN
-        nlevx=klev
-    ELSE
-        nlevx=nlev
-    ENDIF
+    nlevx=nlev
+!    IF (nlev.EQ.klev+1) THEN
+!        nlevx=klev
+!    ELSE
+!        nlevx=nlev
+!    ENDIF
 
     IF (SIZE(field,1) == klon) then

LMDZ6/branches/IPSL-CM6A-MR/libf/phylmd/phys_output_mod.F90

@@ -44 +44 @@
     USE mod_grid_phy_lmdz, only: klon_glo,nbp_lon,nbp_lat
     USE print_control_mod, ONLY: prt_level,lunout
-    USE vertical_layers_mod, ONLY: ap,bp,preff,presnivs, aps, bps, pseudoalt
+    USE vertical_layers_mod, ONLY: ap,bp,preff,presnivs, aps, bps, pseudoalt, presinter
     USE time_phylmdz_mod, ONLY: day_ini, itau_phy, start_time, annee_ref, day_ref
 #ifdef CPP_XIOS
@@ -348 +348 @@
     CALL wxios_add_vaxis("presnivs", &
             levmax(iff) - levmin(iff) + 1, presnivs(levmin(iff):levmax(iff)))
+    CALL wxios_add_vaxis("presinter", &
+            klev + 1, presinter(1:klev+1))
     CALL wxios_add_vaxis("Ahyb", &
             levmax(iff) - levmin(iff) + 1, aps(levmin(iff):levmax(iff)), positif='down', &