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

Last change on this file since 199 was 186, checked in by ymipsl, 10 years ago

Add new openMP parallelism based on distribution of domains on threads. There is no more limitation of number of threads by MPI process.

YM
File size: 927 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
15  LOGICAL, SAVE :: boussinesq
16
17CONTAINS
18 
19  SUBROUTINE init_earth_const
20  USE getin_mod
21  IMPLICIT NONE
22  REAL(rstd) :: X=1
23 
24    CALL getin("radius",radius)
25    CALL getin("g",g)
26    CALL getin("scale_factor",scale_factor)
27    CALL getin("omega",omega) 
28    CALL getin("kappa",kappa) 
29    CALL getin("cpp",cpp) 
30    CALL getin("preff",preff) 
31   
32    boussinesq=.FALSE.
33    CALL getin("boussinesq",boussinesq) 
34   
35    radius=radius/scale_factor
36    omega=omega*scale_factor
37   
38  END SUBROUTINE init_earth_const
39 
40 
41END MODULE earth_const
42 
Note: See TracBrowser for help on using the repository browser.