|
Last change
on this file since 5018 was
1673,
checked in by Laurent Fairhead, 12 years ago
|
|
Fin du phasage de la dynamique parallele localisee (petite memoire) avec le tronc LMDZ5 r1671
Il reste quelques routines a verifier (en particulier ce qui touche a l'etude des cas academiques)
et la validation a effectuer
End of the phasing of the localised (low memory) parallel dynamics package with the
LMDZ5 trunk (r1671)
Some routines still need some checking (in particular the academic cases) and some
validation is still required
|
|
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.
