source: LMDZ5/branches/LF-private/libf/dyn3dmem/comvert.h @ 3893

Last change on this file since 3893 was 1793, checked in by Ehouarn Millour, 11 years ago

Improved sponge layer:

  • Sponge tendencies are now computed analytically, instead than using a Forward Euler approximation.
  • Sponge tendencies are added within top_bound, and the sponge is applied after physics tendencies have been taken into account.

These changes imply that GCM results (when using sponge layer) will be differentwrt bench test cases using previous revisions.
EM

File size: 1.7 KB
Line 
1!
2! $Id: comvert.h 1654 2012-09-24 15:07:18Z aslmd $
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,pseudoalt(llm)
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      real pseudoalt ! pseudo-altitude of model levels (km), based on presnivs(),
26                     ! preff and scaleheight
27
28      integer disvert_type ! type of vertical discretization:
29                           ! 1: Earth (default for planet_type==earth),
30                           !     automatic generation
31                           ! 2: Planets (default for planet_type!=earth),
32                           !     using 'z2sig.def' (or 'esasig.def) file
33
34      logical pressure_exner
35!     compute pressure inside layers using Exner function, else use mean
36!     of pressure values at interfaces
37
38 !-----------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.