source: codes/icosagcm/trunk/src/earth_const.f90 @ 78

Last change on this file since 78 was 58, checked in by ymipsl, 12 years ago

Correct scaling factor for omega

YM
File size: 854 bytes
Line 
1MODULE earth_const
2  USE prec
3  USE math_const
4 
5  REAL(rstd),SAVE :: radius=6.37122E6
6  REAL(rstd),SAVE :: g=9.80616
7  REAL(rstd),PARAMETER :: daysec=86400
8  REAL(rstd),SAVE :: omega=7.292E-5
9  REAL(rstd),SAVE :: kappa=0.2857143
10  REAL(rstd),SAVE :: cpp=1004.70885
11  REAL(rstd),SAVE :: preff=101325.
12  REAL(rstd),SAVE :: pa=50000.
13  REAL(rstd),SAVE :: scale_factor=1.
14
15CONTAINS
16 
17  SUBROUTINE init_earth_const
18  USE ioipsl
19  IMPLICIT NONE
20  REAL(rstd) :: X=1
21 
22    CALL getin("radius",radius)
23    CALL getin("g",g)
24    CALL getin("scale_factor",scale_factor)
25    CALL getin("omega",omega) 
26    CALL getin("kappa",kappa) 
27    CALL getin("cpp",cpp) 
28    CALL getin("preff",preff) 
29   
30    radius=radius/scale_factor
31    omega=omega*scale_factor
32    PRINT *,"radius = ",radius
33   
34  END SUBROUTINE init_earth_const
35 
36 
37END MODULE earth_const
38 
Note: See TracBrowser for help on using the repository browser.