1 | ! |
---|
2 | ! $Id$ |
---|
3 | ! |
---|
4 | SUBROUTINE conf_planete |
---|
5 | ! |
---|
6 | USE IOIPSL |
---|
7 | USE comconst_mod, ONLY: pi, g, molmass, kappa, cpp, omeg, rad, & |
---|
8 | year_day, daylen, daysec, ihf |
---|
9 | USE comvert_mod, ONLY: preff, pa |
---|
10 | IMPLICIT NONE |
---|
11 | ! |
---|
12 | ! |
---|
13 | ! Declarations : |
---|
14 | ! -------------- |
---|
15 | |
---|
16 | ! |
---|
17 | ! local: |
---|
18 | ! ------ |
---|
19 | |
---|
20 | ! --------------------------------------------- |
---|
21 | ! Initialisations de constantes de la dynamique |
---|
22 | ! --------------------------------------------- |
---|
23 | ! Pi |
---|
24 | pi=2.*asin(1.) |
---|
25 | |
---|
26 | !Reference surface pressure (Pa) |
---|
27 | ! 101080 : specific value for CMIP5 aqua/terra planets |
---|
28 | ! "Specify the initial dry mass to be equivalent to |
---|
29 | ! a global mean surface pressure (101325 minus 245) Pa." |
---|
30 | preff=101080. |
---|
31 | CALL getin('preff', preff) |
---|
32 | |
---|
33 | ! Reference pressure at which hybrid coord. become purely pressure |
---|
34 | ! pa=50000. |
---|
35 | pa=preff/2. |
---|
36 | CALL getin('pa', pa) |
---|
37 | |
---|
38 | ! Gravity |
---|
39 | g=9.80665 |
---|
40 | |
---|
41 | CALL getin('g',g) |
---|
42 | ! Molar mass of the atmosphere |
---|
43 | |
---|
44 | molmass = 28.9644 |
---|
45 | CALL getin('molmass',molmass) |
---|
46 | ! kappa=R/Cp et Cp |
---|
47 | |
---|
48 | kappa = 2./7. |
---|
49 | CALL getin('kappa',kappa) |
---|
50 | |
---|
51 | cpp=8.3145/molmass/kappa*1000. |
---|
52 | CALL getin('cpp',cpp) |
---|
53 | ! Radius of the planet |
---|
54 | |
---|
55 | rad = 6371229. |
---|
56 | CALL getin('radius',rad) |
---|
57 | |
---|
58 | ! Length of a standard day (s) |
---|
59 | daysec=86400. |
---|
60 | CALL getin('daysec',daysec) |
---|
61 | |
---|
62 | ! Rotation rate of the planet: |
---|
63 | ! Length of a solar day, in standard days |
---|
64 | daylen = 1. |
---|
65 | |
---|
66 | CALL getin('daylen',daylen) |
---|
67 | ! Number of days (standard) per year: |
---|
68 | |
---|
69 | year_day = 365.25 |
---|
70 | CALL getin('year_day',year_day) |
---|
71 | ! Omega |
---|
72 | ! omeg=2.*pi/86400. |
---|
73 | |
---|
74 | omeg=2.*pi/daysec*(1./daylen+1./year_day) |
---|
75 | CALL getin('omeg',omeg) |
---|
76 | |
---|
77 | ! Intrinsic heat flux (default: none) (only used if planet_type="giant") |
---|
78 | ihf = 0. |
---|
79 | call getin('ihf',ihf) |
---|
80 | |
---|
81 | END SUBROUTINE conf_planete |
---|