1 |
guez |
123 |
program test_fxhyp |
2 |
|
|
|
3 |
guez |
313 |
use dynetat0_m, only: fxhyp, xprimm025, rlonv, xprimv, rlonu, xprimu, & |
4 |
|
|
xprimp025 |
5 |
|
|
use dynetat0_chosen_m, only: read_serre |
6 |
guez |
276 |
use unit_nml_m, only: unit_nml, set_unit_nml |
7 |
guez |
123 |
|
8 |
|
|
implicit none |
9 |
|
|
|
10 |
|
|
!-------------------------------------------------------- |
11 |
|
|
|
12 |
guez |
276 |
call set_unit_nml |
13 |
guez |
123 |
open(unit_nml, file="used_namelists.txt", status="replace", action="write") |
14 |
|
|
call read_serre |
15 |
guez |
277 |
call fxhyp |
16 |
guez |
123 |
close(unit_nml) |
17 |
|
|
|
18 |
guez |
276 |
! We can use the same unit number although we are not writing a namelist: |
19 |
|
|
open(unit_nml, file="test_fxhyp_out.txt", status="replace", action="write") |
20 |
|
|
write(unit_nml, fmt = *) "xprimm025 = ", xprimm025 |
21 |
|
|
write(unit_nml, fmt = *) "rlonv = ", rlonv |
22 |
|
|
write(unit_nml, fmt = *) "xprimv = ", xprimv |
23 |
|
|
write(unit_nml, fmt = *) "rlonu = ", rlonu |
24 |
|
|
write(unit_nml, fmt = *) "xprimu = ", xprimu |
25 |
|
|
write(unit_nml, fmt = *) "xprimp025 = ", xprimp025 |
26 |
|
|
close(unit_nml) |
27 |
guez |
123 |
|
28 |
|
|
end program test_fxhyp |