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

Annotation of /trunk/Sources/dyn3d/test_fxhyp.f

Parent Directory Parent Directory | Revision Log Revision Log


Revision 166 - (hide annotations)
Wed Jul 29 14:32:55 2015 UTC (8 years, 9 months ago) by guez
File size: 949 byte(s)
Split ppm3d.f into files containing a single procedure.

Factorized computations of filtering matrices into a procedure
inifilr_hemisph. Had then to change the matrices from allocatable to
pointer and from customized lower bound to lower bound 1. The change
in lower bounds does not matter because the matrices are only used as
a whole as actual arguments.

Also, in infilr, instead of finding jfilt[ns][uv] from approximately
jjm /2, start at index j1 that corresponds to the equator. This is not
the same if there is a zoom in latitude.

Also, the test (rlamda(modfrst[ns][uv](j)) * cos(rlat[uv](j)) < 1) in
the loops on filtered latitudes is not useful now that we start from
j1: it is necessarily true. See notes.

Just encapsulated lwvn into a module and removed unused argument ktraer.

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 guez 139 use read_serre_m, only: read_serre
7 guez 123 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