source: LMDZ5/branches/LMDZ5-DOFOCO/libf/dyn3dmem/conf_planete.F90 @ 3992

Last change on this file since 3992 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!
4SUBROUTINE conf_planete
5!
6#ifdef CPP_IOIPSL
7USE IOIPSL
8#else
9! if not using IOIPSL, we still need to use (a local version of) getin
10USE ioipsl_getincom
11#endif
12IMPLICIT 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
28pi=2.*asin(1.)
29
30!Reference surface pressure (Pa)
31preff=101325.
32CALL getin('preff', preff)
33! Reference pressure at which hybrid coord. become purely pressure
34! pa=50000.
35pa=preff/2.
36CALL getin('pa', pa)
37! Gravity
38g=9.80665
39CALL getin('g',g)
40! Molar mass of the atmosphere
41molmass = 28.9644
42CALL getin('molmass',molmass)
43! kappa=R/Cp et Cp     
44kappa = 2./7.
45CALL getin('kappa',kappa)
46cpp=8.3145/molmass/kappa*1000.
47CALL getin('cpp',cpp)
48! Radius of the planet
49rad = 6371229.
50CALL getin('radius',rad)
51! Length of a standard day (s)
52daysec=86400.
53CALL getin('daysec',daysec)
54! Rotation rate of the planet:
55! Length of a solar day, in standard days
56daylen = 1.
57CALL getin('daylen',daylen)
58! Number of days (standard) per year:
59year_day = 365.25
60CALL getin('year_day',year_day)
61! Omega
62! omeg=2.*pi/86400.
63omeg=2.*pi/daysec*(1./daylen+1./year_day)
64CALL getin('omeg',omeg)
65
66! Intrinsic heat flux (default: none) (only used if planet_type="giant")
67ihf = 0.
68call getin('ihf',ihf)
69
70END SUBROUTINE conf_planete
Note: See TracBrowser for help on using the repository browser.