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

Annotation of /trunk/test_fxhyp.f

Parent Directory Parent Directory | Revision Log Revision Log


Revision 123 - (hide 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 guez 123 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