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