27 |
! "pa" should be defined before this procedure is called. |
! "pa" should be defined before this procedure is called. |
28 |
|
|
29 |
use nr_util, only: pi |
use nr_util, only: pi |
30 |
use l_util, only: new_unit |
use jumble, only: new_unit |
31 |
|
|
32 |
REAL s(llm+1) |
REAL s(llm+1) |
33 |
! "s(l)" is the atmospheric hybrid sigma-pressure coordinate at |
! "s(l)" is the atmospheric hybrid sigma-pressure coordinate at |
37 |
! "ds(l)" : épaisseur de la couche "l" dans la coordonnée "s" |
! "ds(l)" : épaisseur de la couche "l" dans la coordonnée "s" |
38 |
|
|
39 |
INTEGER l, unit |
INTEGER l, unit |
40 |
REAL alpha, x(llm), trash |
REAL alpha, x(llm) |
41 |
|
|
42 |
character(len=7):: s_sampling = "LMD5" |
character(len=7):: s_sampling = "LMD5" |
43 |
! (other allowed values are "param", "strato1", "strato2" and "read") |
! (other allowed values are "param", "strato1", "strato2" and "read") |
120 |
position="rewind") |
position="rewind") |
121 |
read(unit, fmt=*) ! skip title line |
read(unit, fmt=*) ! skip title line |
122 |
do l = 1, llm + 1 |
do l = 1, llm + 1 |
123 |
read(unit, fmt=*) trash, s(l) |
read(unit, fmt=*) s(l) |
124 |
end do |
end do |
125 |
close(unit) |
close(unit) |
126 |
! Quick check: |
! Quick check: |