My Project
 All Classes Files Functions Variables Macros
conf_planete.F90
Go to the documentation of this file.
1 !
2 ! $Id$
3 !
4 SUBROUTINE conf_planete
5 !
6 #ifdef CPP_IOIPSL
7 USE ioipsl
8 #else
9 ! if not using IOIPSL, we still need to use (a local version of) getin
10 USE ioipsl_getincom
11 #endif
12 IMPLICIT NONE
13 !
14 !
15 ! Declarations :
16 ! --------------
17 #include "dimensions.h"
18 #include "comconst.h"
19 #include "comvert.h"
20 !
21 ! local:
22 ! ------
23 
24 ! ---------------------------------------------
25 ! Initialisations de constantes de la dynamique
26 ! ---------------------------------------------
27 ! Pi
28 pi=2.*asin(1.)
29 
30 !Reference surface pressure (Pa)
31 preff=101325.
32 CALL getin('preff', preff)
33 ! Reference pressure at which hybrid coord. become purely pressure
34 ! pa=50000.
35 pa=preff/2.
36 CALL getin('pa', pa)
37 ! Gravity
38 g=9.80665
39 CALL getin('g',g)
40 ! Molar mass of the atmosphere
41 molmass = 28.9644
42 CALL getin('molmass',molmass)
43 ! kappa=R/Cp et Cp
44 kappa = 2./7.
45 CALL getin('kappa',kappa)
46 cpp=8.3145/molmass/kappa*1000.
47 CALL getin('cpp',cpp)
48 ! Radius of the planet
49 rad = 6371229.
50 CALL getin('radius',rad)
51 ! Length of a standard day (s)
52 daysec=86400.
53 CALL getin('daysec',daysec)
54 ! Rotation rate of the planet:
55 ! Length of a solar day, in standard days
56 daylen = 1.
57 CALL getin('daylen',daylen)
58 ! Number of days (standard) per year:
59 year_day = 365.25
60 CALL getin('year_day',year_day)
61 ! Omega
62 ! omeg=2.*pi/86400.
63 omeg=2.*pi/daysec*(1./daylen+1./year_day)
64 CALL getin('omeg',omeg)
65 
66 ! Intrinsic heat flux (default: none) (only used if planet_type="giant")
67 ihf = 0.
68 call getin('ihf',ihf)
69 
70 END SUBROUTINE conf_planete