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

Last change on this file since 182 was 165, checked in by dubos, 11 years ago

boussinesq flag for shallow-water / Ripa equations
File size: 924 bytes
RevLine 
[12]1MODULE earth_const
2  USE prec
3  USE math_const
4 
[32]5  REAL(rstd),SAVE :: radius=6.37122E6
6  REAL(rstd),SAVE :: g=9.80616
[12]7  REAL(rstd),PARAMETER :: daysec=86400
[32]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.
[12]14
[165]15  LOGICAL, SAVE :: boussinesq
16
[32]17CONTAINS
18 
19  SUBROUTINE init_earth_const
20  USE ioipsl
21  IMPLICIT NONE
22  REAL(rstd) :: X=1
23 
24    CALL getin("radius",radius)
25    CALL getin("g",g)
[40]26    CALL getin("scale_factor",scale_factor)
[32]27    CALL getin("omega",omega) 
28    CALL getin("kappa",kappa) 
29    CALL getin("cpp",cpp) 
30    CALL getin("preff",preff) 
31   
[165]32    boussinesq=.FALSE.
33    CALL getin("boussinesq",boussinesq) 
34   
[32]35    radius=radius/scale_factor
[58]36    omega=omega*scale_factor
[32]37   
38  END SUBROUTINE init_earth_const
39 
40 
[12]41END MODULE earth_const
42 
Note: See TracBrowser for help on using the repository browser.