Last change
on this file since 190 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
|
Rev | Line | |
---|
[12] | 1 | MODULE 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] | 17 | CONTAINS |
---|
| 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] | 41 | END MODULE earth_const |
---|
| 42 | |
---|
Note: See
TracBrowser
for help on using the repository browser.