Changeset 2786 for LMDZ5/trunk/libf/dyn3d_common

Timestamp:

Jan 30, 2017, 11:24:47 AM (7 years ago)

Author:

Ehouarn Millour

Message:

Further convergence with planetary models:

• move vertical_layers_mod.F90 to phy_common and call init_vertical_layers in inigeomphy_mod
• add aps and bps (mid-layer coordinates) to vertical_layers_mod.F90
• same adaptations for the 1D case

EM

File:

• LMDZ5/trunk/libf/dyn3d_common/disvert.F90 (modified) (2 diffs)

LMDZ5/trunk/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))