trunk/dyn3d/test_fxhyp.f

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