trunk/dyn3d/test_fxhyp.f

program test_fxhyp

  use dynetat0_m, only: read_serre, fxhyp, xprimm025, rlonv, xprimv, rlonu, &
       xprimu, xprimp025
  use unit_nml_m, only: unit_nml, set_unit_nml

  implicit none

  !--------------------------------------------------------

  call set_unit_nml
  open(unit_nml, file="used_namelists.txt", status="replace", action="write")
  call read_serre
  call fxhyp
  close(unit_nml)

  ! We can use the same unit number although we are not writing a namelist:
  open(unit_nml, file="test_fxhyp_out.txt", status="replace", action="write")
  write(unit_nml, fmt = *) "xprimm025 = ", xprimm025
  write(unit_nml, fmt = *) "rlonv = ", rlonv
  write(unit_nml, fmt = *) "xprimv = ", xprimv
  write(unit_nml, fmt = *) "rlonu = ", rlonu
  write(unit_nml, fmt = *) "xprimu = ", xprimu
  write(unit_nml, fmt = *) "xprimp025 = ", xprimp025
  close(unit_nml)

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