5 |
============================== |
============================== |
6 |
|
|
7 |
o pkg/exf (interpolation): |
o pkg/exf (interpolation): |
8 |
|
- extended input field (2 rows) near the N & S poles: |
9 |
|
* fill in with the symetric value (when even Nb of data in longitude) |
10 |
|
* add average value at the poles only for scalar quantities ; for vector |
11 |
|
component interpolation, skip the averaging and keep duplicated values, |
12 |
|
unless the 2 components have same location (uvInterp=T) in which case |
13 |
|
interpolate the 2 components together using S/R EXF_INTERP_UV. |
14 |
|
- new S/R EXF_INTERP_UV to interpolate the 2 components of a vector field |
15 |
|
together: this allows to account properly for local orientation when |
16 |
|
averaging at N & S pole. |
17 |
|
- add CPP option EXF_USE_OLD_INTERP_POLE to recover old results regarding |
18 |
|
exf-interpolation near N & S poles. Undef by default. |
19 |
|
- add debug check/print on 1rst iter if debugLevel >=2. |
20 |
|
o pkg/exf: |
21 |
|
- change masking flags definition (no longer "parameter"), put them in |
22 |
|
common block and set them in S/R EXF_INIT_FIXED. |
23 |
|
- change recognized masking flag from s,u,v to c,w,s since the mask |
24 |
|
to use (maskC,W,S) depends on position on model grid and has little to |
25 |
|
do with scalar/vector difference; fix masking if using P-coordinates. |
26 |
|
- fix masking of wind-stress (was wrong if interpolated or read on A-grid) |
27 |
|
and update results of exp. global_with_exf (x2). |
28 |
|
- fix setting of zenith-angle table when SHORTWAVE_HEATING is undef. |
29 |
|
o pkg/obcs |
30 |
|
- S/R obcs_calc_stevens: fix a small bug (gammat -> gammas where appropriate) |
31 |
|
o pkg/exf (interpolation): |
32 |
- assume periodicity in X only if input field cover full longitude range. |
- assume periodicity in X only if input field cover full longitude range. |
33 |
- improve search for lat. index (supposed to be faster, in ~log2(ny) steps, |
- improve search for lat. index (supposed to be faster, in ~log2(ny) steps, |
34 |
and should vectorise). |
and should vectorise). |