source: LMDZ5/branches/LMDZ5V1.0-dev/libf/dyn3dpar/conf_planete.F90 @ 1523

Last change on this file since 1523 was 1437, checked in by Ehouarn Millour, 14 years ago

Implementation of FC's improved Newtonian configuration: more realistic (and tunable) friction & fields towards which to relax and planet parameters (gravity, radius, rotation rate...) can be read in from a planet.def file (this occurs only when running in "newtonian mode"; planet parameters are otherwise read in from the start file).

Checked that these implementation don't change results when running standard Earth bench tests on Vargas.

EM

File size: 1.3 KB
RevLine 
[1437]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
66END SUBROUTINE conf_planete
Note: See TracBrowser for help on using the repository browser.