source: LMDZ5/trunk/libf/dyn3d/comvert.h @ 1625

Last change on this file since 1625 was 1625, checked in by lguez, 12 years ago

Created logical variable "pressure_exner". "pressure_exner" replaces
"disvert_type" to choose between calls to "exner_hyb" and
"exner_milieu". If "pressure_exner" is true, pressure inside layers is
computed from Exner function ("exner_hyb"), else it is the mean of
pressure values at interfaces ("exner_milieu"). "disvert_type" now
only chooses between "disvert" and "disvert_noterre". Default value
for "pressure_exner" is true if "disvert_type" equals 1.
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.5 KB
Line 
1!
2! $Id: comvert.h 1625 2012-05-09 13:14:48Z lguez $
3!
4!-----------------------------------------------------------------------
5!   INCLUDE 'comvert.h'
6
7      COMMON/comvertr/ap(llm+1),bp(llm+1),presnivs(llm),dpres(llm),     &
8     &               pa,preff,nivsigs(llm),nivsig(llm+1),               &
9     &               aps(llm),bps(llm),scaleheight
10
11      common/comverti/disvert_type, pressure_exner
12
13      real ap     ! hybrid pressure contribution at interlayers
14      real bp     ! hybrid sigma contribution at interlayer
15      real presnivs ! (reference) pressure at mid-layers
16      real dpres
17      real pa     ! reference pressure (Pa) at which hybrid coordinates
18                  ! become purely pressure
19      real preff  ! reference surface pressure (Pa)
20      real nivsigs
21      real nivsig
22      real aps    ! hybrid pressure contribution at mid-layers
23      real bps    ! hybrid sigma contribution at mid-layers
24      real scaleheight ! atmospheric (reference) scale height (km)
25
26      integer disvert_type ! type of vertical discretization:
27                           ! 1: Earth (default for planet_type==earth),
28                           !     automatic generation
29                           ! 2: Planets (default for planet_type!=earth),
30                           !     using 'z2sig.def' (or 'esasig.def) file
31
32      logical pressure_exner
33!     compute pressure inside layers using Exner function, else use mean
34!     of pressure values at interfaces
35
36 !-----------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.