Log of /trunk/Sources/filtrez/filtreg_hemisph.f

---

Revision 168 - (view) (annotate) - [select for diffs]
Modified Wed Sep 9 10:41:47 2015 UTC (8 years, 8 months ago) by guez
File length: 735 byte(s)
Diff to previous 142

In order to be able to choose finer resolutions, set large memory
model in compiler options and use dynamic libraries.

Variables rlatd, rlond, cuphy and cvphy of module comgeomphy were
never used. (In LMDZ, they are used only for Orchid.)

There is a bug in PGI Fortran 13.10 that does not accept the
combination of forall, pack and spread in regr_pr_av and
regr_pr_int. In order to circumvent this bug, created the function
gr_dyn_phy.

In program test_inifilr, use a single latitude coordinate for north
and south.

---

Revision 142 - (view) (annotate) - [select for diffs]
Modified Fri Jun 5 19:31:37 2015 UTC (8 years, 11 months ago) by guez
File length: 690 byte(s)
Diff to previous 141

This revision changes results.

---

Revision 141 - (view) (annotate) - [select for diffs]
Modified Fri Jun 5 19:21:08 2015 UTC (8 years, 11 months ago) by guez
File length: 766 byte(s)
Diff to previous 136

Procedure filtreg_hemisph is always called with sdd2 = 1 / sdd1. It
would not make sense otherwise. I do not like that this does not
appear explicitly in filtreg_hemisph. So replaced the two arguments
sdd1 and sdd2 by a single argument sdd. The cost is computing a
division instead of multplication for every 3-dimensional
position. Besides, this changes the results.

---

Revision 136 - (view) (annotate) - [select for diffs]
Added Thu Apr 30 18:35:49 2015 UTC (9 years ago) by guez
File length: 785 byte(s)

Clarified the logic in filtreg by creating a procedure
filtreg_hemisph. It was terrible with a loop on hemispheres and tests
on hemisphere inside the loop, plus maddening indirections on latitude
bounds, plus repeated code. Went from 126 lines to much clearer 74 +
32 = 106 lines.

In module inifilr_m, finally made the arrays matrice[uv][ns],
matrinv[ns] dynamic (following LMDZ). Changed the lower bound of
matriceun and matrinvn in the 3rd dimension: 2 instead of 1, the index
1 was not defined (nor used).

In module inifilr_m, changed the bounds of matriceus and matrinvs in
the 3rd dimension: jfiltsu:jjm instead of 1:jjm - jfiltsu + 1. Changed
the bounds of matricevs in the 3rd dimension: jfiltsv:jjm instead of
1:jjm - jfiltsv + 1. It is a little simpler and clearer this way in
procedure inifilr.