/[lmdze]/trunk/test_fxhyp.f
ViewVC logotype

Contents of /trunk/test_fxhyp.f

Parent Directory Parent Directory | Revision Log Revision Log


Revision 123 - (show annotations)
Thu Feb 5 12:41:08 2015 UTC (9 years, 3 months ago) by guez
File size: 942 byte(s)
Added some test programs.

In fxhyp_loop_ik, changed precision from 1e-3 to 1e-6. Reset initial
value of xo1 to first guess of xi instead of final value of xi for
previous i (better logic).

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 serre, 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

  ViewVC Help
Powered by ViewVC 1.1.21