Changeset 3629 for LMDZ5/branches/IPSLCM5A2.1/libf/phylmd/vertical_layers_mod.F90

Timestamp:

Feb 10, 2020, 9:54:26 AM (4 years ago)

Author:

acozic

Message:

Add new grid, new axis and new variables for cmip protocole and dr2xml

• field_group id="coord_hyb"
• grid_ref="klevp1_bnds"
• grid_ref="klev_bnds"
• domain id="greordered"
• axis id="axis_lat"
• axis id="bnds"
• axis id="klevp1"
• axis id="klev"

File:

• LMDZ5/branches/IPSLCM5A2.1/libf/phylmd/vertical_layers_mod.F90 (modified) (4 diffs)

LMDZ5/branches/IPSLCM5A2.1/libf/phylmd/vertical_layers_mod.F90

@@ -9 +9 @@
    REAL,SAVE,ALLOCATABLE :: bp(:) ! hybrid (sigma contribution) coordinate 
                                   ! at layer interfaces (Pa)
+   REAL,SAVE,ALLOCATABLE :: aps(:) ! hybrid (pressure contribution) coordinate 
+                                   ! at mid-layer (Pa)
+   REAL,SAVE,ALLOCATABLE :: bps(:) ! hybrid (sigma contribution) coordinate 
+                                   ! at mid-layer
    REAL,SAVE,ALLOCATABLE :: presnivs(:) ! reference pressure at mid-layer (Pa),
                                         ! based on preff, ap and bp
@@ -14 +18 @@
                                          ! based on preff and scaleheight
    
-!$OMP THREADPRIVATE(preff,scaleheight,ap,bp,presnivs,pseudoalt)
-   REAL, SAVE, ALLOCATABLE :: Ahyb(:), Bhyb(:)
-!$OMP THREADPRIVATE(Ahyb, Bhyb) 
+!$OMP THREADPRIVATE(preff,scaleheight,ap,bp,aps,bps,presnivs,pseudoalt)
+
 
 CONTAINS
 
   SUBROUTINE init_vertical_layers(nlayer,preff_,scaleheight_,ap_,bp_,&
-                                 presnivs_, pseudoalt_)
+                                 aps_,bps_,presnivs_, pseudoalt_)
     IMPLICIT NONE
     INTEGER,INTENT(IN) :: nlayer ! number of atmospheric layers
@@ -28 +31 @@
     REAL,INTENT(IN)    :: ap_(nlayer+1) ! hybrid coordinate at interfaces
     REAL,INTENT(IN)    :: bp_(nlayer+1) ! hybrid coordinate at interfaces
+    REAL,INTENT(IN)    :: aps_(nlayer) ! hybrid coordinate at mid-layer
+    REAL,INTENT(IN)    :: bps_(nlayer) ! hybrid coordinate at mid-layer
     REAL,INTENT(IN)    :: presnivs_(nlayer) ! Appproximative pressure of atm. layers (Pa)
     REAL,INTENT(IN)    :: pseudoalt_(nlayer) ! pseudo-altitude of atm. layers (km)
@@ -33 +38 @@
     ALLOCATE(ap(nlayer+1))
     ALLOCATE(bp(nlayer+1))
+    ALLOCATE(aps(nlayer))
+    ALLOCATE(bps(nlayer))
     ALLOCATE(presnivs(nlayer))
     ALLOCATE(pseudoalt(nlayer))
-    ALLOCATE(Ahyb(nlayer)) 
-    ALLOCATE(Bhyb(nlayer)) 
-
-
+  
     preff = preff_
     scaleheight=scaleheight_
     ap(:) = ap_(:)
     bp(:) = bp_(:)
+    aps(:) = aps_(:)
+    bps(:) = bps_(:)
     presnivs(:) = presnivs_(:)
     pseudoalt(:) = pseudoalt_(:)