Changeset 3003 for LMDZ5/trunk/libf/phylmd/phys_output_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_mod.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/phys_output_mod.F90

@@ -126 +126 @@
                                                                       -90.,  -90.,  -90.,  -90.,  -90. /)
     REAL, DIMENSION(nfiles), SAVE     :: phys_out_latmax        = (/   90.,   90.,   90.,   90.,   90., &
-                                                                       90.,   90.,   90.,   90.,   90. /)                        
+                                                                       90.,   90.,   90.,   90.,   90. /)
+    REAL, DIMENSION(klev,2) :: Ahyb_bounds, Bhyb_bounds
+    REAL, DIMENSION(klev)   :: lev_index
+                
 #ifdef CPP_XIOS
     ! ug Variables utilis\'ees pour r\'ecup\'erer le calendrier pour xios
@@ -137 +140 @@
     WRITE(lunout,*) 'Debut phys_output_mod.F90'
     ! Initialisations (Valeurs par defaut
+
+    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)
+      lev_index(ilev) = REAL(ilev)
+    END DO
 
     IF (.NOT. ALLOCATED(o_trac)) ALLOCATE(o_trac(nqtot))
@@ -327 +338 @@
       print*,'phys_output_open: Declare vertical axes for each file'
     ENDIF
+
    IF (iff.LE.6.OR.iff.EQ.10) THEN
     CALL wxios_add_vaxis("presnivs", &
             levmax(iff) - levmin(iff) + 1, presnivs(levmin(iff):levmax(iff)))
     CALL wxios_add_vaxis("Ahyb", &
-            levmax(iff) - levmin(iff) + 1, aps)
+            levmax(iff) - levmin(iff) + 1, aps(levmin(iff):levmax(iff)), positif='down', &
+            bnds=Ahyb_bounds(levmin(iff):levmax(iff),:))
     CALL wxios_add_vaxis("Bhyb", &
-            levmax(iff) - levmin(iff) + 1, bps)
-    CALL wxios_add_vaxis("Alt", &
+            levmax(iff) - levmin(iff) + 1, bps(levmin(iff):levmax(iff)), positif='down', &
+            bnds=Bhyb_bounds(levmin(iff):levmax(iff),:))
+    CALL wxios_add_vaxis("klev", levmax(iff) - levmin(iff) + 1, &
+                          lev_index(levmin(iff):levmax(iff)))
+    CALL wxios_add_vaxis("bnds", 2, (/1.,2./))
+
+     CALL wxios_add_vaxis("Alt", &
             levmax(iff) - levmin(iff) + 1, pseudoalt)
    ELSE