source: codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/dyn3d_common/conf_planete.F90 @ 245

Last change on this file since 245 was 222, checked in by ymipsl, 10 years ago

Creating temporary dynamico/lmdz/saturn branche

YM

File size: 1.5 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
24real :: year_day_dyn
25
26! ---------------------------------------------
27! Initialisations de constantes de la dynamique
28! ---------------------------------------------
29! Pi
30pi=2.*asin(1.)
31
32!Reference surface pressure (Pa)
33preff=101325.
34CALL getin('preff', preff)
35! Reference pressure at which hybrid coord. become purely pressure
36! pa=50000.
37pa=preff/2.
38CALL getin('pa', pa)
39! Gravity
40g=9.80665
41CALL getin('g',g)
42! Molar mass of the atmosphere
43molmass = 28.9644
44CALL getin('molmass',molmass)
45! kappa=R/Cp et Cp     
46kappa = 2./7.
47CALL getin('kappa',kappa)
48cpp=8.3145/molmass/kappa*1000.
49CALL getin('cpp',cpp)
50! Radius of the planet
51rad = 6371229. 
52CALL getin('radius',rad)
53! Length of a standard day (s)
54daysec=86400.
55CALL getin('daysec',daysec)
56! Rotation rate of the planet:
57! Length of a solar day, in standard days
58daylen = 1.
59CALL getin('daylen',daylen)
60! Number of days (standard) per year:
61year_day_dyn = 365.25
62CALL getin('year_day',year_day_dyn)
63! Omega
64! omeg=2.*pi/86400.
65omeg=2.*pi/daysec*(1./daylen+1./year_day_dyn)
66CALL getin('omeg',omeg)
67
68! Intrinsic heat flux [default is none]
69! Aymeric -- for giant planets
70! [matters only if planet_type="giant"]
71ihf = 0.
72CALL getin('ihf',ihf)
73
74
75
76END SUBROUTINE conf_planete
Note: See TracBrowser for help on using the repository browser.