Directory: | ./ |
---|---|
File: | dyn3d_common/iniconst.f90 |
Date: | 2022-01-11 19:19:34 |
Exec | Total | Coverage | |
---|---|---|---|
Lines: | 22 | 27 | 81.5% |
Branches: | 2 | 6 | 33.3% |
Line | Branch | Exec | Source |
---|---|---|---|
1 | ! | ||
2 | ! $Id: iniconst.F90 2601 2016-07-24 09:51:55Z emillour $ | ||
3 | ! | ||
4 | 2 | SUBROUTINE iniconst | |
5 | |||
6 | USE control_mod | ||
7 | use IOIPSL | ||
8 | USE comconst_mod, ONLY: im, imp1, jm, jmp1, lllm, lllmm1, lllmp1, & | ||
9 | unsim, pi, r, kappa, cpp, dtvr, dtphys | ||
10 | USE comvert_mod, ONLY: disvert_type, pressure_exner | ||
11 | |||
12 | IMPLICIT NONE | ||
13 | ! | ||
14 | ! P. Le Van | ||
15 | ! | ||
16 | ! Declarations: | ||
17 | ! ------------- | ||
18 | ! | ||
19 | include "dimensions.h" | ||
20 | include "paramet.h" | ||
21 | include "iniprint.h" | ||
22 | |||
23 | character(len=*),parameter :: modname="iniconst" | ||
24 | character(len=80) :: abort_message | ||
25 | ! | ||
26 | ! | ||
27 | ! | ||
28 | !----------------------------------------------------------------------- | ||
29 | ! dimension des boucles: | ||
30 | ! ---------------------- | ||
31 | |||
32 | 1 | im = iim | |
33 | 1 | jm = jjm | |
34 | 1 | lllm = llm | |
35 | 1 | imp1 = iim | |
36 | 1 | jmp1 = jjm + 1 | |
37 | 1 | lllmm1 = llm - 1 | |
38 | 1 | lllmp1 = llm + 1 | |
39 | |||
40 | !----------------------------------------------------------------------- | ||
41 | |||
42 | 1 | dtphys = iphysiq * dtvr | |
43 | 1 | unsim = 1./iim | |
44 | 1 | pi = 2.*ASIN( 1. ) | |
45 | |||
46 | !----------------------------------------------------------------------- | ||
47 | ! | ||
48 | |||
49 | 1 | r = cpp * kappa | |
50 | |||
51 | 1 | write(lunout,*) trim(modname),': R CP Kappa ',r,cpp,kappa | |
52 | ! | ||
53 | !----------------------------------------------------------------------- | ||
54 | |||
55 | ! vertical discretization: default behavior depends on planet_type flag | ||
56 |
1/2✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
|
1 | if (planet_type=="earth") then |
57 | 1 | disvert_type=1 | |
58 | else | ||
59 | ✗ | disvert_type=2 | |
60 | endif | ||
61 | ! but user can also specify using one or the other in run.def: | ||
62 | 1 | call getin('disvert_type',disvert_type) | |
63 | 1 | write(lunout,*) trim(modname),': disvert_type=',disvert_type | |
64 | |||
65 | 1 | pressure_exner = disvert_type == 1 ! default value | |
66 | 1 | call getin('pressure_exner', pressure_exner) | |
67 | |||
68 |
1/2✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
|
1 | if (disvert_type==1) then |
69 | ! standard case for Earth (automatic generation of levels) | ||
70 | 1 | call disvert() | |
71 | ✗ | else if (disvert_type==2) then | |
72 | ! standard case for planets (levels generated using z2sig.def file) | ||
73 | ✗ | call disvert_noterre | |
74 | else | ||
75 | ✗ | write(abort_message,*) "Wrong value for disvert_type: ", disvert_type | |
76 | ✗ | call abort_gcm(modname,abort_message,0) | |
77 | endif | ||
78 | |||
79 | 1 | END SUBROUTINE iniconst | |
80 |