Changeset 3629 for LMDZ5/branches/IPSLCM5A2.1/libf/dyn3d_common/disvert.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/dyn3d_common/disvert.F90 (modified) (2 diffs)

LMDZ5/branches/IPSLCM5A2.1/libf/dyn3d_common/disvert.F90

@@ -10 +10 @@
   use new_unit_m, only: new_unit
   use assert_m, only: assert
-  USE comvert_mod, ONLY: ap, bp, nivsigs, nivsig, dpres, presnivs, &
-                         pa, preff, scaleheight
+  USE comvert_mod, ONLY: ap, bp, aps, bps, nivsigs, nivsig, dpres, presnivs, &
+                         pseudoalt, pa, preff, scaleheight
   USE logic_mod, ONLY: ok_strato
 
@@ -346 +346 @@
   DO l = 1, llm
      dpres(l) = bp(l) - bp(l+1)
+     aps(l) =  0.5 *( ap(l) +ap(l+1))
+     bps(l) =  0.5 *( bp(l) +bp(l+1))
      presnivs(l) = 0.5 *( ap(l)+bp(l)*preff + ap(l+1)+bp(l+1)*preff )
+     pseudoalt(l) = log(preff/presnivs(l))*scaleheight
      write(lunout, *)'PRESNIVS(', l, ')=', presnivs(l), ' Z ~ ', &
-          log(preff/presnivs(l))*scaleheight &
+          pseudoalt(l) &
           , ' DZ ~ ', scaleheight*log((ap(l)+bp(l)*preff)/ &
           max(ap(l+1)+bp(l+1)*preff, 1.e-10))