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

Last change on this file since 1565 was 1520, checked in by Ehouarn Millour, 13 years ago

Implementation of a different vertical discretization (from/for planets, but
can in principle also be used for Earth).
Choice of vertical discretization is set by flag 'disvert_type';
'disvert_type=1' is Earth standard (default; ie set to 1 if
planet_type=="earth") case.
With 'disvert_type=2', approximate altitude of layers and reference atmospheric
scale height must be given using an input file ("z2sig.def", first line
should give scale height, in km, following lines must specify the altitude,
in km above surface, of mid-layers, one per line; see disvert_noterre.F).

Checked that these changes do not impact on 'bench' results, on Vargas.

EM.

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.3 KB
Line 
1!
2! $Id: comvert.h 1520 2011-05-23 11:37:09Z jghattas $
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
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 !-----------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.