trunk/dyn3d/test_fxhyp.f90

program test_fxhyp

  use dynetat0_m, only: fxhyp, xprimm025, rlonv, xprimv, rlonu, xprimu, &
       xprimp025
  use dynetat0_chosen_m, only: read_serre
  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: fxhyp, xprimm025, rlonv, xprimv, rlonu, xprimu, &
4	xprimp025
5	use dynetat0_chosen_m, only: read_serre
6	use unit_nml_m, only: unit_nml, set_unit_nml
7
8	implicit none
9
10	!--------------------------------------------------------
11
12	call set_unit_nml
13	open(unit_nml, file="used_namelists.txt", status="replace", action="write")
14	call read_serre
15	call fxhyp
16	close(unit_nml)
17
18	! 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
28	end program test_fxhyp