source: codes/icosagcm/devel/src/base/earth_const.f90 @ 732

Last change on this file since 732 was 732, checked in by dubos, 6 years ago

devel : more cleanup and reorganization in dynamics/
File size: 1.3 KB
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 :: cppv=1860.
12  REAL(rstd),SAVE :: Rv=461.5
13  REAL(rstd),SAVE :: Treff=273.
14  REAL(rstd),SAVE :: preff=101325.
15  REAL(rstd),SAVE :: pa=50000. ! default value set to preff/2 by disvert_std
16  REAL(rstd),SAVE :: scale_height=8000. ! atmospheric scale height (m)
17  REAL(rstd),SAVE :: scale_factor=1.
18  REAL(rstd),SAVE :: gas_constant = 8.3144621 
19  REAL(rstd),SAVE :: mu                 ! molar mass of the atmosphere
20
21CONTAINS
22 
23  SUBROUTINE init_earth_const
24  USE getin_mod
25  IMPLICIT NONE
26  REAL(rstd) :: X=1
27 
28    CALL getin("radius",radius)
29    CALL getin("g",g)
30    CALL getin("scale_factor",scale_factor)
31    CALL getin("omega",omega) 
32    CALL getin("kappa",kappa) 
33    CALL getin("cpp",cpp) 
34    CALL getin("cppv",cppv)
35    CALL getin("Rv",Rv)
36    CALL getin("preff",preff) 
37    CALL getin("Treff",Treff) 
38    CALL getin("scale_height",scale_height)
39   
40    mu=kappa/cpp
41    radius=radius/scale_factor
42    omega=omega*scale_factor
43   
44  END SUBROUTINE init_earth_const
45 
46 
47END MODULE earth_const
Note: See TracBrowser for help on using the repository browser.