trunk/dyn3d/test_fxhyp.f

program test_fxhyp

  USE dimens_m, ONLY: iim
  use fxhyp_m, only: fxhyp
  use jumble, only: new_unit
  use read_serre_m, only: read_serre
  use unit_nml_m, only: unit_nml

  implicit none

  REAL, dimension(iim + 1):: xprimm025, rlonv, xprimv, rlonu, xprimu, xprimp025
  integer unit

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

  call new_unit(unit_nml)
  open(unit_nml, file="used_namelists.txt", status="replace", action="write")
  call read_serre
  call fxhyp(xprimm025, rlonv, xprimv, rlonu, xprimu, xprimp025)
  close(unit_nml)

  unit = unit_nml
  open(unit, file="test_fxhyp_out.txt", status="replace", action="write")
  write(unit, fmt = *) "xprimm025 = ", xprimm025
  write(unit, fmt = *) "rlonv = ", rlonv
  write(unit, fmt = *) "xprimv = ", xprimv
  write(unit, fmt = *) "rlonu = ", rlonu
  write(unit, fmt = *) "xprimu = ", xprimu
  write(unit, fmt = *) "xprimp025 = ", xprimp025
  close(unit)

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