Log of /trunk/src_gcm

---

Revision 136 - (view) (annotate) - [select for diffs]
Modified Thu Apr 30 18:35:49 2015 UTC (9 years ago) by guez
File length: 3640 byte(s)
Diff to previous 134

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.

---

Revision 134 - (view) (annotate) - [select for diffs]
Modified Wed Apr 29 15:47:56 2015 UTC (9 years ago) by guez
File length: 3622 byte(s)
Copied from: trunk/Compilation/src_gcm revision 133
Diff to previous 132

Sources inside, compilation outside.

---

Revision 132 - (view) (annotate) - [select for diffs]
Modified Fri Mar 20 16:31:06 2015 UTC (9 years, 1 month ago) by guez
Original Path: trunk/Compilation/src_gcm
File length: 3622 byte(s)
Diff to previous 131

Removed procedure jacobi, which was a copy of the file from Numerical
Recipes in Fortran 77. Refer to the Numer_Rec_95 library instead.

There was a strange line in procedure coordij: j cannot be equal to 0
after the loop on j.

---

Revision 131 - (view) (annotate) - [select for diffs]
Modified Fri Feb 27 16:44:07 2015 UTC (9 years, 2 months ago) by guez
Original Path: trunk/Compilation/src_gcm
File length: 3631 byte(s)
Diff to previous 124

Renamed procedure fxhyp_loop_ik to invert_zoom_x.

Bug fix. clat is now in rad so there should be no conversion in
fyhyp. (This bug had an effect only if clat was /= 0.)

No need for heavyside to be double precision.

Removed variable tnom of module iniadvtrac_m. Was redundant with tname.

---

Revision 124 - (view) (annotate) - [select for diffs]
Modified Thu Feb 5 15:19:37 2015 UTC (9 years, 2 months ago) by guez
Original Path: trunk/Compilation/src_gcm
File length: 3631 byte(s)
Diff to previous 123

Moved some processing from fxhyp_loop_ik to fxhyp. Now fxhyp_loop_ik
does not necessarily give longitudes near [-pi, pi]. In fxhyp, we look
in rlonm025 whether we need to move the array toward [-pi, pi]. If so,
we apply the same move to all grids: rlonm025, rlonv, rlonp025, rlonu
and the corresponding derivatives. The move itself is done by the new
procedure principal_cshift. This revision makes the logic
clearer. (For example, we do not have a saved variable is2 in
fxhyp_loop_ik any longer and we remove a test on ik in fxhyp_loop_ik.)

Fixed a bad error message in fxhyp_loop_ik: talked about rlonu when
xvrai is not always rlonu.

---

Revision 123 - (view) (annotate) - [select for diffs]
Added Thu Feb 5 12:41:08 2015 UTC (9 years, 2 months ago) by guez
Original Path: trunk/Compilation/src_gcm
File length: 3612 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).