/[MITgcm]/MITgcm/doc/tag-index
ViewVC logotype

Diff of /MITgcm/doc/tag-index

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1789 by gforget, Sat Dec 22 18:35:27 2012 UTC revision 1.1844 by heimbach, Fri Mar 29 23:21:39 2013 UTC
# Line 4  $Name$ Line 4  $Name$
4      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
5      ==============================      ==============================
6    
7    o verification/isomip, pkg/shelfic and OpenAD:
8      - add OpenAD setup to this exp and make shelfice package work with OpenAD
9    o verification/offline_exf_seaice: replace adjoint experiment.
10    o verification/tutorial_tracer_adjsens, pkg/ptracers and OpenAD:
11      - add OpenAD setup to this exp and make ptracers package work with OpenAD
12    o pkg/seaice/seaice_evp.F: mask the asymmetric contribution of the ice-ocean
13      stress (when the turningAngle is non-zero) over open water to avoid
14      stripes in the EVP-solutions after many EVP-substeps. This does not change
15      the verification experiment (because turningAngle=0 in that experiment), but
16      will affect other EVP solutions with turningAngle.ne.0
17    o verification, AD & TLM output:
18      - new TAF version 2.3.10 (previously 2.3.8) add more parenthesis in multiple
19        additions and multiplications which affect results at machine truncation
20        level. Update few AD output (1D_ocean_ice_column, bottom_ctrl_5x5,
21        lab_sea.noseaicedyn & offline_exf_seaice) and TLM output
22        (1D_ocean_ice_column & tutorial_tracer_adjsens.som81) which had only 12
23        digits of matching.
24    o pkg/seaice:
25       modify lsr-algorithm so that it vectorizes AND converges by using a
26       "zebra" method: split vector loop into two with step size 2; this
27       solves a general convergence problem for the JFNK solver on a vector
28       computer. It has the potential to even improve convergence for
29       non-vector machines, but it does change the results so it is turned
30       off by default. To turn on define CPP-flags SEAICE_VECTORIZE_LSR and
31       SEAICE_VECTORIZE_LSR_ZEBRA in SEAICE_OPTIONS.h
32    
33    checkpoint64e (2013/03/05)
34    o pkg/ecco and several verification experiments:
35       Remove special versions of the_model_main.F, the_main_loop.F, forward_step.F
36       from ecco package, and revert to unified routine in model/src/
37       Adjust verification experiments accordingly, especially gloal_oce_llc90
38       Should conclude long-standing work by Gael to reconsile the_main_loop with
39       ECCO and will enable OpenAD's revolve loop in ECCO setup.
40    o pkg/generic_advdiff:
41      - implement compressible flow method for multi-dim advection (similar to
42        gad_som_advect.F); add new option "GAD_MULTIDIM_COMPRESSIBLE" (since
43        TAF generates many recomputations) to use it; no yet coded with implicit
44        vertical advection.
45      - change experiment advect_cs (new output) and advect_xy (same output)
46        to test compressible flow multi-dim advection.
47    o pkg/generic_advdiff:
48      - fix SOM advection for unusual cases (linear Free-Surf in r*, other
49        cases where linear Free-Surf is not @ k=1, e.g., useShelfIce)
50    o pkg/seaice:
51      - fix bug in the no-slip implementation in seaice_calc_strainrates.F, this
52        bug affects two verification experiments (lab_sea and lab_sea.hb87),
53        recover old (wrong) results by setting noSlipFac = 2. in
54        seaice_calc_strainrates.F;
55      - add another averaging method for etaZ (etaZmethod=3), a mix of 2 and 0
56      - add runtime parameter SEAICE_JFNKepsilon for that finite difference
57        gradient in the Jacobian times vector operation (s/r seaice_jacvec)
58      - fix line search in seaice_jfnk.F (s/r seaice_jfnk_update)
59    o pkg/seaice: reorganize seaice_lsr.F to be able to reuse more code:
60      - new subroutines seaice_lsr_tridiagu/v
61      - introduce loop bounds i/jmin/max to be able to use a restricted
62        additive Schwarz method, but not used for now.
63      - rename SEAICE_CALC_LSR_COEFFS into SEAICE_LSR_CALC_COEFFS
64      - use new subroutines seaice_lsr_tridiagu/v in s/r seaice_preconditioner
65    o OpenAD support
66       Reorganize S/R call for overall time-stepping loop to facilitate
67       OpenAD revolve loop replacement
68       (new S/R main_do_loop.F)
69    o model/src:
70      - add 3-D array 3rd dimension to argument list of S/R ADAMS_BASHFORTH2/3;
71      - pkg/cheapaml: remove local S/R ADAMS2d and use instead standard S/R
72        ADAMS_BASHFORTH2.
73    o pkg/cheapaml:
74      - add new parameter (default=F, same as before) to disable full domain
75        (including land) periodicity in X & Y dir.
76      - new S/R CHEAPAML_COPY_EDGES to copy relaxation field into edges of
77        prognostics field ; merge the 2 S/R cycle_2d_tracer and timestep_2d_tracer
78        into single S/R CHEAPAML_TIMESTEP.
79    
80    checkpoint64d (2013/02/17)
81    o pkg/monitor (angular momentum):
82      - do not add Eta contribution to the total AM if using rigid-lid;
83      - correct for Adams-Bashforth in zonal-wind part (instead of in AM_eta).
84      - horizontal discretisation: compute Zonal wind contribution separately for
85        each component (instead of from cell centered averaged) and use the 4
86        grid-cell corner averaged value of Omega*(R*cos_Phi)^2 in AM_eta part.
87      - Check that AM is conserved with linear dynamics (momAdvection=F and linear
88        FreeSurf) if using vectorInvariant momentum.
89    o model/src (horizontal grid):
90      - move initialisation of all horizontal grid arrays to null value
91        from specific horizontal grid S/R to main S/R INI_GRID;
92      - add 2 new arrays to project model C-grid velocity on to zonal direction
93        (similar to AngleCosC & -AngleSinC but at U and V point location) and
94        compute them in S/R CALC_GRID_ANGLES.
95    o pkg/seaice: replace a hard coded parameter by a runtime parameter
96        SEAICE_JFNK_tolIter
97    o pkg/seaice: introduce a restricted additive Schwarz (RAS) method in
98        preconditioner for JFNK (run LSR into the overlaps as far as possible,
99        i.e OLx/y-2, but overwrite the solution in the overlaps)
100      - extend loop range in seaice_calc_viscosity
101      - add new parameters SEAICE_OLx/y (default = OLx/y-2)
102      - add flexible loop ranges to seaice_preconditioner, seaice_calc_lsr_coeffs
103      - slightly modify loops in tridiagonal solver lsr-preconditioner
104      - by default remove the extra exchange in seaice_preconditioner with a
105        replacement of an undocumented CPP flag by SEAICE_PRECOND_EXTRA_EXCHANGE
106      - modify offline_exf_seaice so that old results without RAS are reproduced
107        for output.dyn_jfnk.txt
108      in seaice_fgmres.F (unrelated to RAS)
109      - add a cpp-flag option for a re-ordered mapping (potentially more
110        efficient), but because this changes the results, turn it off by
111        default
112      - replace the save statement by local common blocks
113    o utils/matlab:
114      - fix tile-offset in rdmnc.m when output was produced using EXCH2.
115      - update cs_grid/bk_line scripts to use load_grid.m (in utils/matlab)
116        instead of local script "load_cs.m".
117    o pkg/seaice:
118      - fix a serious bug in seaice_fgmres.F, s/r seaice_scalprod, connected
119        to the scalarProductMetric and the unrolling. By chance this does not
120        change the verification experiment, because of the particular domain
121        decomposition
122    o pkg/diagnostics:
123      - change default missing_value from UNSET_RL to -999.
124    o pkg/layers: improve vectorization of layers_fluxcalc.F
125      - step 1: split loops, make TatU/V 2D arrays. This alone reduces the cpu
126        time spent in this routine by more than a factor of 3 on an NEC-SX8
127      - step 2: replace search algorithm by a bisection algorithm found in
128        Press et al. Numerical Recipes in Fortran. This algorithm can be
129        vectorized and finds the correct bin in approx. log2(Nlayers) steps.
130        For now, use this algorithm only for TARGET_NEC_SX defined.
131      - moved the search algorithm into a separate subroutine, that is kept
132        within the same file for automatic inlining.
133      - modified catching potential errors
134    o pkg/ecco:
135      - add compile switches ALLOW_SHALLOW_ALTIMETRY and ALLOW_HIGHLAT_ALTIMETRY.
136    o model/src, pkg/autodiff, pkg/ctrl, pkg/grdchk :
137      - modifications for ALLOW_GENTIM2D_CONTROL
138      -- forward_step.F, the_main_loop.F : pass options and arrays for adjoint tapes.
139      -- load_fields_driver.F : add call to CTRL_MAP_GENTIM2D( myTime, myIter, myThid )
140      -- add xx_gentim2d0/1 in ctrl_ad_check_lev1_dir.h checkpoint_lev1_directives.h etc
141      -- use ctrl_swapffields.F (new) instead of exf_SwapFFields.F in ctrl_get_gen.F
142    o model/src:
143      - account for moist thermodynamics in (dry) potential temperature equation;
144        update results of experiments: aim.5l_cs (x2), fizhi-cs (std + aqualev)
145        and cpl_aim+ocn.
146      - also switch to standard compiler (gfortran, was using ifort previously)
147        for fizhi-cs experiments reference output.
148    o pkg/ctrl:
149      - 1st step toward generic time-varying control variable (2D) via CPP option
150        ALLOW_GENTIM2D_CONTROL
151      - note new namelist ctrl_nml_genarr in data.ctrl
152    o verification/global_ocean.cs32x15:
153      - switch off OLD_THSICE_CALL_SEQUENCE (both forward and AD code) to use new
154        thsice calling sequence and update "output.icedyn.txt" output file.
155    o pkg/ptracers:
156      - from Ryan A.: code to allow to restart a tracer at regular time interval
157        specifying a reset frequency and a reset-phase (for each tracers).
158    o pkg/thsice:
159      - implement new sequence of calls for thsice+seaice:
160          previously:   ice-Dyn,ice-Advect,ice-Thermo(thsice)
161          new sequence: ice-Thermo(thsice),ice-Dyn,ice-Advect
162      - allows (with temporary CPP option "#define OLD_THSICE_CALL_SEQUENCE"
163        in CPP_OPTIONS.h) to recover old sequence;
164      - update diagnostics for Intermediate State thsice diagnostics: now filled
165        just before advection (instead of just after advection in old sequence);
166      - update output of experiment offline_exf_seaice (standard test);
167      - in global_ocean.cs32x15 exp., postpone output.icedyn.txt update
168        by setting #define OLD_THSICE_CALL_SEQUENCE in CPP_OPTIONS.h
169    o pkg/aim_v23:
170      - use tiled array for aim_CO2 and surface wind speed for export to coupler;
171      - use standard pkg/mdsio routines to read & write CO2 pickup file;
172      - in-line S/R AIM_AIM2DYN_EXCHANGES in aim_do_physics.F;
173      - Implement new sequence of calls for seaice advection (with thsice and
174        seaice diffusion): do ice advection after ice-Thermodynamic (previously
175        was done before); change aim_do_physics.F (split the bi,bj loops block);
176      - use same CPP option (OLD_THSICE_CALL_SEQUENCE) to recover old sequence;
177    o pkg/thsice:
178      - add a missing EXCH of sIceLoad in thsice_main.F (ocean set-up) for the case
179        where just seaice diffusion is used without seaice-dynamics.
180    
181    checkpoint64c (2013/01/18)
182    o pkg/seaice:
183      - add a parameter JFNKres_tFac that can replace JFNKres_t
184        (= JFNKres_tFac*JFNKresidual before JFNK solver)
185      - add a line search option to the JFNK solver
186      - move update of u/vIce and computation of initial residual into a separate
187        routine (same file); makes the code a little easier to read (?)
188      - the line search requires computation of the residual and for some reason
189        I cannot do the update of u/vIce at the beginning of the Newton loop
190        without changing the results, so we move computation of the new residual
191        also to the end of the Newton loop, even without the line search
192    o pkg/mdsio:
193      - write missing value (corresponding to undefined array value) to meta files;
194        skip this if missing value is one (UNSET_RL would be a better value for the
195        case of undefined missing value but it's current the default missing-value);
196        Also update mdsio_read_meta.F (reading timeInterval & missing-value).
197    o pkg/exf & thsice:
198      - for diagnostics purpose, fill-up and update Qsw array (+ fill the
199        diagnostics) even when SHORTWAVE_HEATING is #undef;
200    o verification/cfc_example:
201      - turn on pkg/layers with pkg/diagnostics (to test pkg/layers and to
202        illustrate how to use it).
203    o verification/offline_exf_seaice:
204      in input: add SEAICErestoreUnderIce=.TRUE., and reduce size of forcing
205        files to only one time-record.
206      in input.thermo: narrow the 2 bands of low initial ice-concentration near
207        N & S edges; switch on SEAICErestoreUnderIce=T; use pkg/diagnostics for
208        snapshot output; same mcPheePiston value as in thsice with
209        surf-current=0.2m/s.
210      in input.thsice: similar set-up as in input.thermo (channel with same seaice
211        initial conditions).
212    o verification/cfc_example:
213      - update file "data" with more standard choice (default rotationPeriod,
214        exactConserv=T, reduce ivdc_kappa to 10, remove multiDimAdvection=F,
215        set deltaTFreeSurf to deltaTtracer) and update the output.
216    o pkg/layers:
217      - fix bolus part so that it compiles with #undef GM_BOLUS_ADVEC;
218      - ensure that "layers_bolus" is switched to False when useGMRedi=False.
219    o verification/offline_exf_seaice:
220      - set all {field}period to zero for constant-in-time forcing (not tested
221        anywhere else).
222    o verification/offline_exf_seaice:
223      - dyn_lsr  : use pkg/seaice S/R (instead of pkg/thsice) to advect sea-ice;
224      - dyn_jfnk : use pkg/thsice S/R (instead of pkg/seaice) to advect sea-ice.
225    o model/src:
226      - implement Finite-Volume method for (hydrostatic) presure gradient
227        from S.-J. Lin (QJRMS 1997), for atmosphere using sigma-coordinate.
228    o pkg/seaice: implement EVP* by Lemieux et al. (2012).
229      - In order to be compatible with the default EVP code, this requires changing
230        the EVP code around a little; and because any small change immediately
231        breaks the lab_sea.hb87 verification experiment, the reference output is
232        updated two times (each time the remaining aggrement was only 2 digits):
233        - after replacing 1/dt by 1*recip_dt with a precomputed recip_dt
234        - after replacing 1/(1+0.5*dt/tau) by a/(a+1), with a=2*tau/dt
235        The adjoint test lab_sea.evp also needed to be updated.
236      - Turn on EVP* by setting SEAICE_deltaTevp as before, and setting
237        SEAICEuseEVPstar=.true.; for EVP* SEAICE_deltaTevp does not need to be
238        a factor of SEAICE_deltaTdyn, and you can set the number of EVP steps
239        explicitly with SEAICEnEVPstarSteps (thereby overriding the default of
240        INT(SEAICE_deltaTdyn/SEAICE_deltaTevp)). Controlling EVP* is deliberately
241        kept obscure, because I need to test it a little more.
242    o pkg/cheapaml:
243      - fix few incorrect type (_RL/_RS) in subroutine calls.
244    o pkg/seaice:
245      - new S/R SEAICE_MAP_RS2VEC (similar to SEAICE_MAP2VEC) to map _RS arrays
246        to _RL vector; this fixes JFNK code when compiling with "-use_r4" option.
247    o model/src:
248      - add a warning if using select_rStar=2 with nonlinFreeSurf <> 4;
249        change (in input_ad*/data) select_rStar from 2 to 1 for AD verification
250        experiments global_ocean.cs32x15 & tutorial_tracer_adjsens (+ update output)
251    o pkg/cheapaml:
252      - replace "273.1 _d 0" by celsius2K in pkg/cheapaml/cheapaml.F and
253        update output of experiment cheapAML_box.
254      - store in CHEAPAML.h coeffs ssq[0:2] for saturation specific humidity,
255        and set them only in one place (cheapaml_readparms.F).
256      - replace "273.16 _d 0" by celsius2K & replace "3.14159 _d 0" by constant "PI"
257      - remove re-setting of main-model parameter "celsius2K" (in cheapaml_readparms)
258        and re-setting of cheapAML params "lath" & "p0" (in cheapaml_init_varia).
259    o model/src/the_main_loop.F:
260      - fix revision 129 : need both ALLOW_ECCO and ALLOW_ECCO_EVOLUTION.
261    o pkg/seaice:
262      - remove legacy branch code.
263      - retire SEAICE_GROWTH_LEGACY and SEAICE_CAP_HEFF accordingly.
264    o verification/*/code*/SEAICE_OPTIONS.h:
265      - remove documentation for retired CPP flags (SEAICE_GROWTH_LEGACY,
266        SEAICE_CAP_HEFF).
267    o pkg/cheapaml:
268      - use standard routine (S/R GET_PERIODIC_INTERVAL) to decide when it is time
269        to load new forcing fields.
270      - add missing "_d 0" where it's needed; + some editing (improve indentation).
271      - remove switch to "-r8" compiler option (no longer needed) in verification
272        cheapAML_box, keep same output (small trucation error differences).
273    o verification/offline_exf_seaice:
274      - add thermodynamics only test.
275    o verification/cheapAML_box:
276      - split the domain in 8 (2x4) tiles; add MPI and multi-threaded test;
277        keep same output (same 14 digits for cg2d).
278      - switch on "-r8" compiler option (some numerical values are missing
279        the "_d 0") and update output (only 6 digits for cg2d were preserved).
280    
281    checkpoint64b (2012/12/24)
282    o pkg/cheapaml:
283      - remove un-used variables (Qnet,EmPmR) from pickup files.
284      - store tendency of advected variables into cheapaml pickup file
285        and fix the AB-2 restart.
286  o verification/offline_exf_seaice and lab_sea adjoint:  o verification/offline_exf_seaice and lab_sea adjoint:
287   - switch to new code (#undef SEAICE_GROWTH_LEGACY).    - switch to new code (#undef SEAICE_GROWTH_LEGACY).
288   - switch to 4 time steps in offline, consistent with lab_sea.    - switch to 4 time steps in offline, consistent with lab_sea.
289   - match paramaters from SEAICE_GROWTH_LEGACY branch in data.seaice.    - match parameters from SEAICE_GROWTH_LEGACY branch in data.seaice.
290   - update results after change of experiment seaice.    - update results after change of experiment seaice.
291    o pkg/seaice:
292      - remove option SEAICE_MULTICATEGORY; has been replaced with run-time
293        parameter SEAICE_multDim. Issue a warning if SEAICE_MULTICATEGORY is
294        still defined since SEAICE_multDim needs to be set to MULTDIM value (from
295        SEAICE_SIZE.h) to recover previous #define SEAICE_MULTICATEGORY default.
296  o model/src:  o model/src:
297    - add run-time parameter to select CG2D solver minimum-residual solution    - add run-time parameter to select CG2D solver minimum-residual solution
298      (in case of a poor convergence).      (in case of a poor convergence).
# Line 18  o pkg/seaice/seaice_lsr.F Line 302  o pkg/seaice/seaice_lsr.F
302      file; the main purpose of this change is to be able to reuse the code      file; the main purpose of this change is to be able to reuse the code
303      in S/R SEAICE_PRECONDITIONER      in S/R SEAICE_PRECONDITIONER
304    - modify the computation of the residual so that the mass of sea-ice    - modify the computation of the residual so that the mass of sea-ice
305      is not used twice; this simplyfies the calculations, but since      is not used twice; this simplifies the calculations, but since
306      the residual is never used except when LSR_mixIniGuess.ge.4---and this is      the residual is never used except when LSR_mixIniGuess.ge.4---and this is
307      not tested in testreport---this does not affect the testreport results.      not tested in testreport---this does not affect the testreport results.
308  o verification/seaice_obcs:  o verification/seaice_obcs:
# Line 52  o pkg/seaice: Line 336  o pkg/seaice:
336    - add a verification experiment for JFNK in offline_exf_seaice    - add a verification experiment for JFNK in offline_exf_seaice
337  o verification:  o verification:
338    - change (for now, only fwd test) offline_exf_seaice experiment:    - change (for now, only fwd test) offline_exf_seaice experiment:
339      use a re-entrant channel with idealized costline (trapezoidal shape)      use a re-entrant channel with idealized coastline (trapezoidal shape)
340      at southern boundary to explore seaice dynamics and thermodynamics      at southern boundary to explore seaice dynamics and thermodynamics
341      with simple forcing. Main set-up has dynamics+thermo and SST relaxation;      with simple forcing. Main set-up has dynamics+thermo and SST relaxation;
342      secondary test (input.thsice) use thermodynamics only with SST relaxation.      secondary test (input.thsice) use thermodynamics only with SST relaxation.

Legend:
Removed from v.1.1789  
changed lines
  Added in v.1.1844

  ViewVC Help
Powered by ViewVC 1.1.22