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

Last change on this file since 187 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
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
[186]20  USE getin_mod
[32]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.