Log of /trunk/Sources/dyn3d/PPM3d/qckxyz.f

---

Revision 176 - (view) (annotate) - [select for diffs]
Modified Tue Feb 23 17:00:39 2016 UTC (8 years, 2 months ago) by guez
File length: 2988 byte(s)
Diff to previous 166

1e-60 was underflowing in simple precision.

---

Revision 166 - (view) (annotate) - [select for diffs]
Added Wed Jul 29 14:32:55 2015 UTC (8 years, 9 months ago) by guez
File length: 2312 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.