trunk/dyn3d/test_fxhyp.f

program test_fxhyp

  use dynetat0_m, only: read_serre, fxhyp, xprimm025, rlonv, xprimv, rlonu, &
       xprimu, xprimp025
  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	277	use dynetat0_m, only: read_serre, fxhyp, xprimm025, rlonv, xprimv, rlonu, &
4			xprimu, xprimp025
5	guez	276	use unit_nml_m, only: unit_nml, set_unit_nml
6	guez	123
7			implicit none
8
9			!--------------------------------------------------------
10
11	guez	276	call set_unit_nml
12	guez	123	open(unit_nml, file="used_namelists.txt", status="replace", action="write")
13			call read_serre
14	guez	277	call fxhyp
15	guez	123	close(unit_nml)
16
17	guez	276	! We can use the same unit number although we are not writing a namelist:
18			open(unit_nml, file="test_fxhyp_out.txt", status="replace", action="write")
19			write(unit_nml, fmt = *) "xprimm025 = ", xprimm025
20			write(unit_nml, fmt = *) "rlonv = ", rlonv
21			write(unit_nml, fmt = *) "xprimv = ", xprimv
22			write(unit_nml, fmt = *) "rlonu = ", rlonu
23			write(unit_nml, fmt = *) "xprimu = ", xprimu
24			write(unit_nml, fmt = *) "xprimp025 = ", xprimp025
25			close(unit_nml)
26	guez	123
27			end program test_fxhyp