trunk/dyn3d/test_fxhyp.f

program test_fxhyp

  USE dimensions, 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	guez	123	program test_fxhyp
2
3	guez	265	USE dimensions, ONLY: iim
4	guez	123	use fxhyp_m, only: fxhyp
5			use jumble, only: new_unit
6	guez	139	use read_serre_m, only: read_serre
7	guez	123	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