Last change
on this file since 3946 was
1505,
checked in by Ehouarn Millour, 13 years ago
|
Added possibility to have the base of the atmosphere heated (to mimic a constant uniform heat flux from below), set in .def file (parameter ihf = ... , in W/m2); only active if in "Newtonian mode" (iflag_phys=2) and if planet_type=="giant".
EM
|
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.