60 |
real:: vert_z0hig = 80. ! height at which resolution reaches dz |
real:: vert_z0hig = 80. ! height at which resolution reaches dz |
61 |
real:: vert_h_hig = 20. ! width of the transition |
real:: vert_h_hig = 20. ! width of the transition |
62 |
|
|
63 |
real, pointer:: p(:) ! (llm + 1) pressure (in hPa) |
real, allocatable:: p(:) ! (llm + 1) pressure (in hPa) |
64 |
|
|
65 |
namelist /disvert_nml/vert_sampling, vert_scale, vert_dzmin, vert_dzlow, & |
namelist /disvert_nml/vert_sampling, vert_scale, vert_dzmin, vert_dzlow, & |
66 |
vert_z0low, vert_dzmid, vert_z0mid, vert_h_mid, vert_dzhig, & |
vert_z0low, vert_dzmid, vert_z0mid, vert_h_mid, vert_dzhig, & |
125 |
p = preff * EXP(- zz(2:) / vert_scale) |
p = preff * EXP(- zz(2:) / vert_scale) |
126 |
ya = pa / preff |
ya = pa / preff |
127 |
s(2: llm) = hybrid(p) |
s(2: llm) = hybrid(p) |
128 |
deallocate(p) ! pointer |
|
129 |
call compute_ab |
call compute_ab |
130 |
|
|
131 |
case("read_hybrid") |
case("read_hybrid") |