1 |
guez |
3 |
program test_disvert |
2 |
|
|
|
3 |
|
|
use dimens_m, only: llm |
4 |
|
|
use comvert, only: pa, disvert, ap, bp, preff |
5 |
|
|
use comconst, only: initialize |
6 |
|
|
use in_out, only: new_unit |
7 |
|
|
|
8 |
|
|
implicit none |
9 |
|
|
|
10 |
|
|
integer unit, l |
11 |
|
|
REAL sigma(llm+1) ! "sigma(l)" : à l'interface des couches "l" et "l-1" |
12 |
|
|
real p(llm + 1) ! pressure at level, in Pa |
13 |
|
|
|
14 |
|
|
!--------------------- |
15 |
|
|
|
16 |
|
|
print *, "llm = ", llm |
17 |
|
|
call initialize |
18 |
|
|
pa = 5e4 |
19 |
|
|
call disvert |
20 |
|
|
sigma = ap / pa + bp |
21 |
|
|
p = ap + bp * preff |
22 |
|
|
|
23 |
|
|
unit = new_unit() |
24 |
|
|
open(unit, file="test_disvert.csv", status="replace", action="write") |
25 |
|
|
write(unit, fmt=*) '"ap (Pa)" "bp" "sigma" "p (Pa)"' ! title line |
26 |
|
|
do l = 1, llm |
27 |
|
|
write(unit, fmt=*) ap(l), bp(l), sigma(l), p(l) |
28 |
|
|
end do |
29 |
|
|
close(unit) |
30 |
|
|
print *, 'The file "test_disvert.csv" has been created.' |
31 |
|
|
|
32 |
|
|
end program test_disvert |