| 4 |
Notes on tags used in MITgcmUV |
Notes on tags used in MITgcmUV |
| 5 |
============================== |
============================== |
| 6 |
|
|
| 7 |
|
o pkg/seaice: replace erroneously storing e12Csqr on the wrong tape with |
| 8 |
|
the wrong key by inialising it before use -> fixes a recomputation |
| 9 |
|
problem but does not fix the gradient |
| 10 |
|
o pkg/layers: |
| 11 |
|
- move resetting of flux arrays to outside the iLa loop |
| 12 |
|
(this caused layers_diapycnal to not work with more than one layer coord.) |
| 13 |
|
- fix bug in calculation of layers_recip_delta in layers_init_fixed.F |
| 14 |
|
o pkg/layers: |
| 15 |
|
- remove from common block 3D layers fields (used to diagnose transport |
| 16 |
|
in tracer coords) and declare them as local arrays in layers_calc.F |
| 17 |
|
o pkg/seaice: |
| 18 |
|
- Rename SEAICE_VECTORIZE_LSR_ZEBRA to SEAICE_LSR_ZEBRA and define it |
| 19 |
|
in lab_sea forward experiment (changes results -> updated) |
| 20 |
|
- update adjoint experiment lab_sea after cleaning up seaice_lsr_tridiagu/v |
| 21 |
|
o pkg/seaice: modify seaice_lsr.F in order to improve the gradient |
| 22 |
|
computations; for SEAICE_VECTORIZE_LSR |
| 23 |
|
- move the loops over which the tridiagonal solvers (seaice_lsr_tridiagu/v) |
| 24 |
|
sweep to the outside, |
| 25 |
|
- remove store directives and add "CADJ loop sequential" directives |
| 26 |
|
in analogy to model/src/solve_tridiagonal.F |
| 27 |
|
- replace many "#ifdef SEAICE_VECTORIZE_LSR_ZEBRA" by variable loop |
| 28 |
|
boundaries to yield more readable code. This has also the charming side |
| 29 |
|
effect that your can use SEAICE_VECTORIZE_LSR_ZEBRA without |
| 30 |
|
SEAICE_VECTORIZE_LSR (i.e. adjoint without recomputations in |
| 31 |
|
seaice_lsr_tridiagu/v now requires either of these flags, vectorization |
| 32 |
|
still requires SEAICE_VECTORIZE_LSR with SEAICE_VECTORIZE_LSR_ZEBRA as |
| 33 |
|
an additional option) |
| 34 |
|
The tridiagonal solvers are now completely analoguous to solve_tridiagonal.F |
| 35 |
|
but the adjoint gradients (which are affected by this change) still explode. |
| 36 |
|
o coupler (atm & ocn _compon_interf): |
| 37 |
|
- avoid using Exch2 IO mapping if not suitable. |
| 38 |
|
|
| 39 |
|
checkpoint65m (2015/06/15) |
| 40 |
|
o pkg/layers: |
| 41 |
|
- clean-up old parameters setting (from Sep 2012): |
| 42 |
|
+ remove old parameters setting from common block (declared locally); |
| 43 |
|
+ stop if any old parameter is mixed with new parameters; |
| 44 |
|
+ check for valid or consistent parameter setting (e.g., layers_name) |
| 45 |
|
- check that layers_bounds for Pot.Density are in the right range; |
| 46 |
|
this is supposed to catch un-updated setting after recently shifting |
| 47 |
|
potential density field "prho" by -1000 (now in "sigma" scale). |
| 48 |
o pkg/layers: |
o pkg/layers: |
| 49 |
- Added missing call to layers_fill for keeping track of tottend in |
- Added missing call to layers_fill for keeping track of tottend in |
| 50 |
pkg/diagnostics/diagnostics_fill_state.F |
pkg/diagnostics/diagnostics_fill_state.F |
| 62 |
thermodynamics.F, impldiff.F) and gad to save the necessary fields |
thermodynamics.F, impldiff.F) and gad to save the necessary fields |
| 63 |
for layers_thermodynamics calculations |
for layers_thermodynamics calculations |
| 64 |
o pkg/ptracers: |
o pkg/ptracers: |
| 65 |
- add capability of stepping a passive tracer forward |
- add capability of stepping a passive tracer forward without advection |
| 66 |
without advection by setting PTRACERS_advScheme to 0. |
by setting PTRACERS_advScheme to 0. |
| 67 |
o pkg/diagnostics: |
o pkg/diagnostics: |
| 68 |
- add few (threads) _BARRIER around any update of "diag_pkgStatus". |
- add few (threads) _BARRIER around any update of "diag_pkgStatus". |
| 69 |
o pkg/layers: |
o pkg/layers: |
Parent Directory
| 
Revision Log
| 
Revision Graph
| 
Patch