trunk/dyn3d/test_fxhyp.f90

program test_fxhyp

  USE dimensions, ONLY: iim
  use dynetat0_m, only: read_serre
  use fxhyp_m, only: fxhyp
  use unit_nml_m, only: unit_nml, set_unit_nml

  implicit none

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

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

  call set_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)

  ! 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	guez	123	program test_fxhyp
2
3	guez	265	USE dimensions, ONLY: iim
4	guez	276	use dynetat0_m, only: read_serre
5	guez	123	use fxhyp_m, only: fxhyp
6	guez	276	use unit_nml_m, only: unit_nml, set_unit_nml
7	guez	123
8			implicit none
9
10			REAL, dimension(iim + 1):: xprimm025, rlonv, xprimv, rlonu, xprimu, xprimp025
11
12			!--------------------------------------------------------
13
14	guez	276	call set_unit_nml
15	guez	123	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	guez	276	! 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	guez	123
30			end program test_fxhyp