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

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

Parent Directory Parent Directory | Revision Log Revision Log


Revision 166 - (show 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 program test_fxhyp
2
3 USE dimens_m, ONLY: iim
4 use fxhyp_m, only: fxhyp
5 use jumble, only: new_unit
6 use read_serre_m, 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