source:
LMDZ5/branches/LMDZ5_SPLA/libf/dyn3d_common/conf_planete.F90
Last change on this file was 1945, checked in by , 11 years ago | |
---|---|
|
|
File size: 1.4 KB |
Line | |
---|---|
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 |
Note: See TracBrowser
for help on using the repository browser.