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