source: LMDZ6/branches/DYNAMICO-conv-GC/libf/dyn3d_common/conf_planete.F90

Last change on this file was 2600, checked in by Ehouarn Millour, 8 years ago

Cleanup in the dynamics: turn comvert.h into module comvert_mod.F90
EM

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
File size: 1.5 KB
Line 
1!
2! $Id$
3!
4SUBROUTINE conf_planete
5!
6#ifdef CPP_IOIPSL
7USE IOIPSL
8#else
9! if not using IOIPSL, we still need to use (a local version of) getin
10USE ioipsl_getincom
11#endif
12USE comconst_mod, ONLY: pi, g, molmass, kappa, cpp, omeg, rad, &
13                        year_day, daylen, daysec, ihf
14USE comvert_mod, ONLY: preff, pa
15IMPLICIT NONE
16!
17!
18!   Declarations :
19!   --------------
20
21!
22!   local:
23!   ------
24
25! ---------------------------------------------
26! Initialisations de constantes de la dynamique
27! ---------------------------------------------
28! Pi
29pi=2.*asin(1.)
30
31!Reference surface pressure (Pa)
32preff=101325.
33CALL getin('preff', preff)
34! Reference pressure at which hybrid coord. become purely pressure
35! pa=50000.
36pa=preff/2.
37CALL getin('pa', pa)
38! Gravity
39g=9.80665
40CALL getin('g',g)
41! Molar mass of the atmosphere
42molmass = 28.9644
43CALL getin('molmass',molmass)
44! kappa=R/Cp et Cp     
45kappa = 2./7.
46CALL getin('kappa',kappa)
47cpp=8.3145/molmass/kappa*1000.
48CALL getin('cpp',cpp)
49! Radius of the planet
50rad = 6371229.
51CALL getin('radius',rad)
52! Length of a standard day (s)
53daysec=86400.
54CALL getin('daysec',daysec)
55! Rotation rate of the planet:
56! Length of a solar day, in standard days
57daylen = 1.
58CALL getin('daylen',daylen)
59! Number of days (standard) per year:
60year_day = 365.25
61CALL getin('year_day',year_day)
62! Omega
63! omeg=2.*pi/86400.
64omeg=2.*pi/daysec*(1./daylen+1./year_day)
65CALL getin('omeg',omeg)
66
67! Intrinsic heat flux (default: none) (only used if planet_type="giant")
68ihf = 0.
69call getin('ihf',ihf)
70
71END SUBROUTINE conf_planete
Note: See TracBrowser for help on using the repository browser.