6 |
|
|
7 |
private llm |
private llm |
8 |
|
|
9 |
real ap(llm+1), pa ! in Pa |
real, save:: ap(llm+1), pa ! in Pa |
10 |
real bp(llm+1) |
real, save:: bp(llm+1) |
11 |
real presnivs(llm) ! pressions approximatives des milieux de couches, en Pa |
|
12 |
real, parameter:: preff = 101325. ! in Pa |
real, save:: presnivs(llm) |
13 |
real nivsigs(llm), nivsig(llm+1) |
! pressions approximatives des milieux de couches, en Pa |
14 |
|
|
15 |
save |
real, parameter:: preff = 101325. ! in Pa |
16 |
|
|
17 |
contains |
contains |
18 |
|
|
22 |
! Author: P. Le Van |
! Author: P. Le Van |
23 |
|
|
24 |
! This procedure sets the vertical grid. It defines the host |
! This procedure sets the vertical grid. It defines the host |
25 |
! variables "ap", "bp", "presnivs", "nivsigs" and "nivsig". "pa" |
! variables "ap", "bp", "presnivs". "pa" should be defined before |
26 |
! should be defined before this procedure is called. |
! this procedure is called. |
27 |
|
|
28 |
use jumble, only: new_unit |
use jumble, only: new_unit |
29 |
use nr_util, only: pi, assert |
use nr_util, only: pi, assert |
59 |
|
|
60 |
print *, "Call sequence information: disvert" |
print *, "Call sequence information: disvert" |
61 |
|
|
|
forall (l = 1: llm) nivsigs(l) = REAL(l) |
|
|
forall (l = 1: llm + 1) nivsig(l) = REAL(l) |
|
|
|
|
62 |
print *, "Enter namelist 'disvert_nml'." |
print *, "Enter namelist 'disvert_nml'." |
63 |
read(unit=*, nml=disvert_nml) |
read(unit=*, nml=disvert_nml) |
64 |
write(unit_nml, nml=disvert_nml) |
write(unit_nml, nml=disvert_nml) |