| 4 |
Notes on tags used in MITgcmUV |
Notes on tags used in MITgcmUV |
| 5 |
============================== |
============================== |
| 6 |
|
|
| 7 |
|
o nonHydrostatic: |
| 8 |
|
Fix missing vertical flux of vert. momentum near surface (k=1). This fixes |
| 9 |
|
a spurious source of energy in simple baroclinic adjusment test case. |
| 10 |
|
Update all non-hydrostatic exp output. |
| 11 |
|
o model/src (3-d solver): |
| 12 |
|
Change again 3-D solver free-surface RHS term when exactConserv=T: |
| 13 |
|
simpler, consistent with exactConserv, and works with implicDiv2Dflow < 1 |
| 14 |
|
(which was not the case with previous formulation). |
| 15 |
|
To recover previous version, just reset "oldFreeSurfTerm" to previous value. |
| 16 |
|
Affects results of exp. global_ocean.cs32x15.viscA4 & hs94.cs-32x32x5.impIGW. |
| 17 |
|
o nonHydrostatic: |
| 18 |
|
- compatible with implicDiv2Dflow < 1 ; + implement implicit factor |
| 19 |
|
for Non-Hydrostatic pressure gradient (similar to implicSurfPress). |
| 20 |
|
- change 3-D solver to work with r* (for now, only with linear free-surf.) |
| 21 |
|
- add few debug diagnostics (cg3d_x, residual) |
| 22 |
|
o model/src: |
| 23 |
|
fix missing EmPmR in 3-D solver RHS (was in cg2d_b but missing in cg3d_b). |
| 24 |
|
update output of exp. global_ocean.cs32x15.viscA4 |
| 25 |
|
o update output of exp. global_ocean.cs32x15.icedyn after removing |
| 26 |
|
#define SEAICE_OLD_AND_BAD_DISCRETIZATION from SEAICE_OPTIONS.h |
| 27 |
|
o add entry points for Gianmaria's Nesting pkgs (from Contrib/nesting_sannino) |
| 28 |
|
o mnc: add horizontal grid coordinates (XC,YC) and (XG,YG) to netcdf output |
| 29 |
|
for complex grids, for now only to "state", "adstate", "adexf", "adseaice" |
| 30 |
|
o model: add Christopher Wolfe's implemetation of the single reduction cg2d |
| 31 |
|
(d'Azevedo, Eijkhout, and Romine, 1999) with CPP-flag ALLOW_SRCG and |
| 32 |
|
runtime flag useSRCGSolver. Modify verification experiment global_with_exf |
| 33 |
|
to test this code |
| 34 |
|
o pkg/cd_code: |
| 35 |
|
Change defaut D-grid time-stepping to get the same for D-grid momentum as for |
| 36 |
|
C-grid momentum. D-grid velocity used to be stepped forward in time with |
| 37 |
|
Adams-Bashforth only on surface pressure term. Tests show that using AB on |
| 38 |
|
D-grid coriolis term improves stability (as expected from CD-scheme paper). |
| 39 |
|
To reproduce old results, add 2 options (in CD_CODE_OPTIONS.h). |
| 40 |
|
Update results of exp: global_with_exf (x2), ideal_2D_oce, isomip (x2), |
| 41 |
|
lab_sea (3 Fwd +3 AD), natl_box (x2), tutorial_global_oce_biogeo (Fwd+AD), |
| 42 |
|
tutorial_global_oce_latlon, tutorial_global_oce_optim & tutorial_tracer_adjsens. |
| 43 |
|
Keep old results by adding CD_CODE_OPTIONS.h in exp: cfc_example, exp2, |
| 44 |
|
lobal_ocean.90x40x15(+AD), global_ocean_ebm(+AD) and OpenAD. |
| 45 |
|
o pkg/obcs: |
| 46 |
|
- fix bug in initialisation (broken since Oct 1rst modification) |
| 47 |
|
for useOBCSprescribe without useEXF. |
| 48 |
|
|
| 49 |
|
checkpoint61y |
| 50 |
|
o pkg/seaice |
| 51 |
|
- change (once again) the discretization for the evp solver: rather than |
| 52 |
|
squaring the averages of strain rates, average the squares to give |
| 53 |
|
a better approximation to VP behavior (in simple test); this breaks |
| 54 |
|
verification experiments that use EVP (lab_sea and lab_sea.hb87) |
| 55 |
|
--> update outputs |
| 56 |
|
o add CPP-flag (in mom_calc_visc.F) and TAF directive (in gmredi_slope_limit) |
| 57 |
|
to improve vectorization of adjoint code on TARGET_NEC_SX |
| 58 |
|
o verification/vermix: |
| 59 |
|
- change vermix to a more "classical" 1d column with high vertical |
| 60 |
|
resolution, seasonal heat flux forcing and constant wind forcing |
| 61 |
|
--> update outputs |
| 62 |
|
o pkg/seaice |
| 63 |
|
- fix diagnostics of principal stress components for EVP |
| 64 |
|
- add optional diagnostics of non-linear convergence for LSOR solver |
| 65 |
|
- make maximum zeta a runtime parameter |
| 66 |
|
o model/src (3-d solver): |
| 67 |
|
Undo modif from Dec 12, 2005 (solve_for_pressure.F cvs version v.49), |
| 68 |
|
which consisted in removing the Surf. Non-Hydrostatic pressure from the |
| 69 |
|
whole column and adding it to Eta (inspired from Casualli, 1999): |
| 70 |
|
Only affects 3-D solver initial guess and only if using ExactConserv, |
| 71 |
|
but tests did show a degradation of solver convergence. For this reason, |
| 72 |
|
decide to disable this piece of code. |
| 73 |
|
to recover previous version, just set zeroPsNH = exactConserv (line 78). |
| 74 |
|
Update output of exp. global_ocean.cs32x15.viscA4 & hs94.cs-32x32x5.impIGW |
| 75 |
|
o pkg/ggl90: |
| 76 |
|
- fix previous modif (unfortunatly has not been documented here). |
| 77 |
o pkg/seaice/autodiff: |
o pkg/seaice/autodiff: |
| 78 |
- add CPP flag to include clipping of zeta in EVP code; |
- add CPP flag to include clipping of zeta in EVP code; |
| 79 |
although clipping of zeta is not necessary for EVP and not recommended |
although clipping of zeta is not necessary for EVP and not recommended |
| 80 |
it was done by default, now clipping is disabled by default |
it was done by default, now clipping is disabled by default |
| 81 |
- modify computation and store zetaC/Z in seaice_evp for TAF |
- modify computation and store zetaC/Z in seaice_evp for TAF |
| 82 |
- add CPP-flag SEAICE_ALLOW_CLIPVELS (undefined by default) |
- add CPP-flag SEAICE_ALLOW_CLIPVELS (undefined by default) |
| 83 |
- change seaice_calc_viscosities/strainraties for better |
- change seaice_calc_viscosities/strainrates for better |
| 84 |
vectorization after AD by TAF |
vectorization after AD by TAF |
| 85 |
- rearrange common blocks in SEAICE.h once again for easier dumps of |
- rearrange common blocks in SEAICE.h once again for easier dumps of |
| 86 |
adjoint output, adjust adcommon.h and g_common.h |
adjoint output, adjust adcommon.h and g_common.h |
| 87 |
- change flux form of default advection scheme for better vectorization |
- change flux form of default advection scheme for better vectorization |
| 88 |
of adjoint code |
of adjoint code |
Parent Directory
| 
Revision Log
| 
Revision Graph
| 
Patch