Rev | Line | |
---|
[1632] | 1 | ! |
---|
[1658] | 2 | ! $Id: iniconst.F 1403 2010-07-01 09:02:53Z fairhead $ |
---|
[1632] | 3 | ! |
---|
| 4 | SUBROUTINE iniconst |
---|
| 5 | |
---|
| 6 | USE control_mod |
---|
[1657] | 7 | |
---|
[1632] | 8 | IMPLICIT NONE |
---|
| 9 | c |
---|
| 10 | c P. Le Van |
---|
| 11 | c |
---|
| 12 | c----------------------------------------------------------------------- |
---|
| 13 | c Declarations: |
---|
| 14 | c ------------- |
---|
| 15 | c |
---|
| 16 | #include "dimensions.h" |
---|
| 17 | #include "paramet.h" |
---|
| 18 | #include "comconst.h" |
---|
| 19 | #include "temps.h" |
---|
| 20 | #include "comvert.h" |
---|
[1657] | 21 | #include "iniprint.h" |
---|
[1632] | 22 | |
---|
| 23 | |
---|
| 24 | c |
---|
| 25 | c |
---|
| 26 | c |
---|
| 27 | c----------------------------------------------------------------------- |
---|
| 28 | c dimension des boucles: |
---|
| 29 | c ---------------------- |
---|
| 30 | |
---|
| 31 | im = iim |
---|
| 32 | jm = jjm |
---|
| 33 | lllm = llm |
---|
| 34 | imp1 = iim |
---|
| 35 | jmp1 = jjm + 1 |
---|
| 36 | lllmm1 = llm - 1 |
---|
| 37 | lllmp1 = llm + 1 |
---|
| 38 | |
---|
| 39 | if (planet_type=="earth") then |
---|
| 40 | disvert_type=1 |
---|
| 41 | else |
---|
| 42 | disvert_type=2 |
---|
| 43 | endif |
---|
| 44 | |
---|
| 45 | c----------------------------------------------------------------------- |
---|
| 46 | |
---|
| 47 | dtdiss = idissip * dtvr |
---|
| 48 | dtphys = iphysiq * dtvr |
---|
| 49 | unsim = 1./iim |
---|
| 50 | pi = 2.*ASIN( 1. ) |
---|
| 51 | |
---|
| 52 | c----------------------------------------------------------------------- |
---|
| 53 | c |
---|
| 54 | |
---|
| 55 | r = cpp * kappa |
---|
| 56 | |
---|
[1657] | 57 | write(lunout,*)'iniconst: R CP Kappa ', r , cpp, kappa |
---|
[1632] | 58 | c |
---|
| 59 | c----------------------------------------------------------------------- |
---|
| 60 | |
---|
| 61 | CALL disvert(pa,preff,ap,bp,dpres,presnivs,nivsigs,nivsig) |
---|
| 62 | c |
---|
| 63 | c |
---|
| 64 | RETURN |
---|
| 65 | END |
---|
Note: See
TracBrowser
for help on using the repository browser.