1 |
program test_fxhyp |
program test_fxhyp |
2 |
|
|
3 |
USE dimensions, ONLY: iim |
USE dimensions, ONLY: iim |
4 |
|
use dynetat0_m, only: read_serre |
5 |
use fxhyp_m, only: fxhyp |
use fxhyp_m, only: fxhyp |
6 |
use jumble, only: new_unit |
use unit_nml_m, only: unit_nml, set_unit_nml |
|
use read_serre_m, only: read_serre |
|
|
use unit_nml_m, only: unit_nml |
|
7 |
|
|
8 |
implicit none |
implicit none |
9 |
|
|
10 |
REAL, dimension(iim + 1):: xprimm025, rlonv, xprimv, rlonu, xprimu, xprimp025 |
REAL, dimension(iim + 1):: xprimm025, rlonv, xprimv, rlonu, xprimu, xprimp025 |
|
integer unit |
|
11 |
|
|
12 |
!-------------------------------------------------------- |
!-------------------------------------------------------- |
13 |
|
|
14 |
call new_unit(unit_nml) |
call set_unit_nml |
15 |
open(unit_nml, file="used_namelists.txt", status="replace", action="write") |
open(unit_nml, file="used_namelists.txt", status="replace", action="write") |
16 |
call read_serre |
call read_serre |
17 |
call fxhyp(xprimm025, rlonv, xprimv, rlonu, xprimu, xprimp025) |
call fxhyp(xprimm025, rlonv, xprimv, rlonu, xprimu, xprimp025) |
18 |
close(unit_nml) |
close(unit_nml) |
19 |
|
|
20 |
unit = unit_nml |
! We can use the same unit number although we are not writing a namelist: |
21 |
open(unit, file="test_fxhyp_out.txt", status="replace", action="write") |
open(unit_nml, file="test_fxhyp_out.txt", status="replace", action="write") |
22 |
write(unit, fmt = *) "xprimm025 = ", xprimm025 |
write(unit_nml, fmt = *) "xprimm025 = ", xprimm025 |
23 |
write(unit, fmt = *) "rlonv = ", rlonv |
write(unit_nml, fmt = *) "rlonv = ", rlonv |
24 |
write(unit, fmt = *) "xprimv = ", xprimv |
write(unit_nml, fmt = *) "xprimv = ", xprimv |
25 |
write(unit, fmt = *) "rlonu = ", rlonu |
write(unit_nml, fmt = *) "rlonu = ", rlonu |
26 |
write(unit, fmt = *) "xprimu = ", xprimu |
write(unit_nml, fmt = *) "xprimu = ", xprimu |
27 |
write(unit, fmt = *) "xprimp025 = ", xprimp025 |
write(unit_nml, fmt = *) "xprimp025 = ", xprimp025 |
28 |
close(unit) |
close(unit_nml) |
29 |
|
|
30 |
end program test_fxhyp |
end program test_fxhyp |