/[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.460 by heimbach, Mon Feb 28 17:19:47 2005 UTC revision 1.581 by heimbach, Thu Nov 24 00:07:29 2005 UTC
# Line 1  Line 1 
1  $Header$  $Header$
2    $Name$
3    
4      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
5      ==============================      ==============================
6    
7    o Non-Hydrostatic:
8     - put all NH variables (formerly in DYNVARS.h & GW.h) in NH_VARS.h
9     - fix exchange calls for CS-grid using the new EXCH_S3D_RL.
10    o Allow to apply AB on T,S rather than on AB(gT,gS):
11     - implemented within #ifdef ALLOW_ADAMSBASHFORTH_3
12     - use the same arrays (gtNm,gsNm) to hold tracer field at previous
13       time-steps (if AB(T,S)) and tendencies (if AB(gT,gS)).
14     - need a "clever" read_pickup to restart from the other AB formulation
15    o move local commom bloc /GlobalLo/ (used for singleCpuIO) in EESUPPORT.h
16    o remove some unused variables (reduces number of compiler warnings)
17    o adjoint
18      * Disable useVariableViscosity for adjoint for time being
19        (there are a number of extensive recomputations)
20      * Initialize certain variables needed for adjoint
21    
22  eckpoint57e_pre  checkpoint57w_post
23    o pkg/mdsio: use MAX_LEN_FNAM (instead of hard coded 80) for file-names
24                 fix the singleCpuIO RS version of MDSWRITEFIELD
25       and add a WARNING if using globalFiles=T with multi-processors.
26    o write grid-file(s) in only one format (MDS or MNC): default=MNC (if
27      using MNC); can force MDS output by setting writegrid_mnc=F in data.mnc
28    o tracer advection scheme:
29      add DST-2 (=Lax-Wendroff) & 1rst order upwind schemes (the way to get
30       those 2 schemes before, changing GAD_FLUX_LIMITER.h + to recompile,
31       was not very convenient).
32      add implicit version for DST-3 & DST3_Flux_Limit vertical advection;  
33       (further tests needed for the implicit DST & Flux-Limiter schemes).
34    o PTRACERS: added namelist parameter PTRACERS_ref in analogy to tRef
35      and sRef for convenience, handy, when you want to initialize
36      ptracers in a large domain with a constant non-zero value. Also,
37      untangle ptracers_init.F a little.
38    
39    checkpoint57v_post
40    o change dome exp. to use variable viscosity in side-drag terms.
41    o remove all the files from the cfc_offline test since they are now
42      incorporated into tutorial_cfc_offline which has been producing
43      identical results for a few days of testing
44    o add new parameter: sideDragFactor to enable half-slip-side BC
45      and change side-drag S/R to use viscosity from mom_calc_visc
46      (can switch back by setting sideDragFactor=0.); affects results at
47      truncation level; regenerate output for exp.: global_ocean.cs32x10(x3),
48       global_ocean_ebm, global_ocean_pressure, matrix_example & MLAdjust(x5).
49    o vorticity & strain are now computed without masking, to account for
50      no-slip-side BC in variable viscosity calculation ; mask is applied
51      anyway later on, before computing advection/coriolis terms.
52    o modify exp4 to test new obcs support for reading from external files
53      and for passive tracers
54    o fix behavior of ptracers_monitor when monitor output is written to a
55      netcdf-file; still not independent of monitor.F (writes in the same
56      file).
57    o change default value of tRef (put a uniform, meaningful value for ocean
58      & atmos) and initialize sRef (was missing, as Martin mentioned).
59    o Redefine viscC?Leith*.  Change the coefficients of Leith viscosity
60      to be proportional to enstrophy dissipation length.  Thus,
61       viscC?leith?_new=pi*(viscC?leith?_old)^(1/3)
62      This is now consistent with Smagorinsky and with the manual.
63    
64    checkpoint57u_post
65    o OBCS and PTRACERS: add open boundary support for passive tracers
66      - either use homogenous (pseudo) v.Neumann conditions or prescribe
67        OB-values from file; this is not different from the way theta and salinity
68        are treated
69      - however, Orlanski-radiation conditions are not supported, and the model
70        will stop if you use pTracers and Orlanski at the same time.
71      - beefed up the rountine obcs_external_fields_load: now only those open
72        boundary values are overwritten with values from files for which there
73        are really files, otherwise the OB-fields remain untouched. This makes
74        it possible to use different OBs at different ends of the domain (as
75        with EXF)
76    o reduce number of exch calls in shap_filt (s2).
77    o CS-grid: fix divergence contribution to Leith viscosity. Affects (@
78      truncation level) 4/5 MLAdjust test results (11&12 digits match).
79    
80    checkpoint57t_post
81    o test biharmonic viscosity + no-slip with exp: global_ocean.cs32x15.viscA4
82    o finish adding brackets (coriolis S/R, vecinv) in multi-terms sum to get
83      same truncation, independent of the face number: affects affects several
84      test. exp.: Update global_ocean.cs32x15 (+.thsice) & MLAdjust (all 5).
85    o fix bug in mom_v_sidedrag (biharmonic viscosity was counted 2 times);
86      note that biharmonic viscosity + no-slip was not tested in any experiment.
87    o fix bug in non-linear free-surf side-drag ; generate new output for
88      global_ocean.cs32x15(+.thsice) & global_ocean_pressure.
89    o make mnc honor the writeBinaryPrec flag for most kinds of output
90    o eesupp: split EEPARAMS into EESIZE/EEPARAMS and add a vector global sum
91    o minor change in averaging W (mom_vi_?_vertshear.F) that affects truncation
92      error; update output from global_ocean.cs32x15(+.thsice) & MLAdjust (all 5).
93    o add brackets in multi-terms sum to get same truncation, independent of
94      the face number (manage to get exactly same value of gU,gV at the cs-edge)
95      affects several test. exp. Update global_ocean.cs32x15 (+.thsice).
96    o Clean, new version of mom_calc_visc.  Adds a new feature, viscAhRemax
97      and viscA4Remax, which allow one to limit the minimum viscosity by specifying
98      a maximum grid-scale Reynolds number.  Also, changes all length calculations
99      in the case of grid anisotropy to geometric mean of dx**2 and dy**2.
100      Suggested values of viscosity parameters found in mom_calc_visc.
101      Many new viscosity diagnostics available.
102    
103    checkpoint57s_post
104    o Add clearing-house for viscosity calculations, mom_calc_visc and
105      correct numerous bugs in LeithD and Smagorinsky discretizations.  
106      mom_calc_visc now allows for Smagorinsky, Leith and LeithD to be
107      used simultaneously, and to be used in mom_vecinv and mom_fluxform.
108      Mom_vecinv can be used with standard viscous calculation (mom_vi_hdissip.F)
109      and with Strain-Tension viscous term calculation (mom_hdissip.F).  
110      The Strain-Tension form is now called with useStrainTensonVisc=.true.
111      ViscAhStrain and viscAhTension are now retired.
112    o diagnostics "VISCAHD" & "VISCA4D" replace "VISCAH" & "VISCA4".
113    o In response to JMC's bug report, genmake2 now ignores any soft-links
114      that it finds in the *current* directory when searching for files.
115      Soft-links found in any other source directories will still work as
116      they did before.
117    o add HAVE_SETRLSTK define and useSETRLSTK flag (eedata) which calls
118      a C routine to unlimit the stack size
119    o mnc-ify pkg/aim snapshot vars
120    o mnc-ify pkg/land as requested by Daniel
121    o various changes to mnc including:
122      - all files use the new "BASENAME[[.ITER].{t|f}NUM].nc" format
123      - output can now be grouped so that all files within a group
124          change the ITER portion of their names in lock-step together
125      - can now read ("global") PER-FACE (in addition to PER-TILE) files
126          and works with both EXCH1 and EXCH2 (but needs more testing)
127      - writing works for all verification test cases w/ g77 on Linux
128    
129    checkpoint57r_post
130    o make MNC variable initialization happen earlier
131    o add diagnostics for (almost) each term of the momentum Eq.
132    o remove SST,SSS,UVEL_k2,VVEL_k2 diagnostics.
133    o move pkg timeave out of pkg-groups "oceanic" & "atmospheric" (pkg/pkg_groups)
134    o remove MINIMAL_TAVE_OUTPUT option ;
135      add "dumpInitAndLast" to replace "nodumps" (and MINIMAL_TAVE_OUTPUT).
136    o move some basic S/R (used elsewhere in the code) from pkg/mdsio to eesupp/src
137      allows to compile without pkg/mdsio (& test successfully hs94.1x64x5)
138    o add MNC output for adjoint of EXF (adexf), SEAICE (adseaice)
139    o add time-dependent SST, SSS control
140    o cost terms in ecco/ modularized (new cost_generic.F)
141    o cost_averages routines are modularized (new cost_averagesgeneric.F)
142    o new infrastructure for seaice model vs. obs misfits in seaice/
143      * driver from the_main_loop.F is seaice_cost_driver.F
144      * first term is AREA vs. SMR observation misfit (ALLOW_SEAICE_COST_SMR_AREA)
145      * gradient check in Exp. 6a of natl_box_adjoint/
146    o natl_box_adjoint input_seaice for Lab. Sea now available online via wget
147      (see README file)
148    
149    checkpoint57q_post
150    o fix initialisation Pb in calc_gw (get NANs in exp5 with g77)
151    o fix annoying months-old bug in mnc/mdsio output for seaice pkg
152      - bug was reported by Jean-Michel and Ian Fenty
153    o fix initialization Pb in AIM-Ocean coupled model (aim_sice2aim.F)
154      update exp. cpl_aim+ocn output.
155    o Implement 4-level checkpointing: AUTODIFF_4_LEVEL_CHECKPOINT
156      Needed for high-res. runs with large global files
157      to be able to limit tapelev size
158    o Changed some tests for adjoint:
159      * carbon/ now tests AUTODIFF_2_LEVEL_CHECKPOINT
160      * natl_box_adjoint/code_bulk/ tests AUTODIFF_4_LEVEL_CHECKPOINT
161      * global_ocean.90x40x15/ back to default 3-level checkpointing
162    o Modifications to enable scripted COMMON to MODULE conversion
163      * verification/OpenAD/code/ has necessary scripts and description;
164      * converts, compiles and runs successfully;
165    
166    checkpoint57p_post
167    o reading of delXfile & delYfile moved from ini_parms.F to ini_grid.F
168    o move ${PKG}_MNC_INIT from ${PKG}_READ_PARAMS to ${PKG}_INIT_FIXED
169      (already the case for some pkgs, including recent MNC_init, e.g. thsice)
170    o call CALC_GW from DYNAMICS (instead of from FORWARD_STEP)
171    o Dissipation & phiHyd gradient are always added to gU,gV in timestep.F
172      (was already the case for dissipation with mom_vecinv,
173       and also the case for grad.PhiHyd when staggered-timeStep)
174      This will allow to put dissipation out-off the AB time-stepping.
175     Unfortunately, affects truncation error. Update output of exp. that "fails"
176      (cg2d matching number) and leave the others for later: front_relax(10),
177      global_ocean_pressure(10),ideal_2D_oce(10),internal_wave(9),lab_sea(12),
178      matrix_example(7!),natl_box(12),plume_on_slope(11),vermix(12x3,11x2).
179    
180    checkpoint57o_post
181    o fix MNC/MDSIO flags in PTRACERS so that (1) it works, and
182      (2) its ready to be documented during the 20050801 DocFest
183    o MNC-ify the gmredi package -- tested and works w/ and w/o mnc
184      for ideal_2D_oce and lab_sea on linux_ia32_g77
185    o set to #undef CTRL_SET_OLD_MAXCVARS_30
186      in global1x1_tot/code_seaice_adj to access indices 32, 33
187      (worth another tag)
188    
189    checkpoint57n_post
190    o adding precip control xx_precip 32
191    o adding swflux control xx_swflux 33
192    o Replacing 10 forcing penalty cost routines by 1 generic.
193    o Implementing partitioning in time-mean and variable adjustment
194      for atmos. state.
195    o Time-mean constant error is now read in data.ecco (instead of hard-coded)
196    o in natl_box_adjoint:
197      * updated exp 1 to 10
198        (but does not have time-mean vs. variable split results yet)
199    o updated global1x1/*_seaice_adj/
200      to incorporate latest changes
201    o enable useSingleCpuIO for mdsio vector routines
202      (applies in particular to tapelev I/O)
203    o removed inefficient K-loop in all active_file*control*.F
204    
205    checkpoint57m_post
206    o re-arrange exp. adjustment.cs-32x32x1 (adjust_nlfs.cs-32x32x1 merge into)
207      and advect_cs (use EXCH-2 and standard cs-32 grid files).
208    o pkg/exch2: filling of face-corner halo regions is optional
209      (ifdef W2_FILL_NULL_REGIONS) and is turned off by default.
210    
211    checkpoint57l_post
212     tag before updating tutorials (pre-release-2)
213    o add cos & sin of grid-direction orientation angle (curvilinear grid only,
214       for output and interpolation of input fields onto model grid)
215    o fix radius in mom_v_metric_cylinder.F (wrong by dy/2).
216       Affects results of exp. rotating_tank ; update output.txt.
217    o pkg/diagnostics: 2 new S/R to fill a diagnostics using a scaling factor
218          and taking the square quantity and/or using a fraction-weight field.
219    o fix aim_diagnostics_init.F mate numbers.
220    o add mnc output capability to diagnostics/diagstat and update
221      our cvsignore files for the ACSII output generated
222    
223    checkpoint57k_post
224    o lab_sea now uses diagnostics package and provides a suite of
225      new EXF diagnostics output.
226    
227    checkpoint57j_post
228    o exf
229      * bug fix in exf_bulkformulae for psimh (spotted by Eliyahu Biton):
230        Bracket missing in eqn. for psimh (eqn. 7 of Large & Pond, 1982);
231        terms -2*ATAN(x) + pi/2 were added to both stable and unstable case.
232      * extended diagnostics output
233        - now includes all EXF variables
234        - hs, hl moved to common block
235        - separate S/R for diagnostics, called after ctrl vector is added
236        - name changed for EXFsw -> EXFswnet, EXFlw -> EXFlwnet
237        - EXFqnet contains true Qnet (i.e. includes qsw)
238      * updated verif. lab_sea, global_with_exf
239      * Rearranging hflux (exf_mapfields -> exf_getforcing)
240        to ensure that diagnostics sees a true EXFqnet
241    o add matlab scripts for CS-grid vectors (in utils/cs_grid)
242    o diagnostics pkg: change pointers so that 1 diagnostic can be used
243      several times with different frequencies.
244    o mnc-ify the thsice package as requested by Daniel Enderton
245      - the monitor--MNC output needs work
246      - many attributes need to be added (most are currently blank)
247    o mnc-ify theaim_v23 package as requested by Daniel Enderton
248      - many attributes need to be added (most are currently blank)
249    o remove "groundAtK1" & "rkFac" (replaced by usingPCoords & -rkSign)
250    o add Sadourny like discretization to the 4th.O vorticity advection scheme
251    o fix bug in rStar code (calc_r_star.F) when used with topo/bathy.
252      update output.txt from: aim.5l_cs (& .thSI), cpl_aim+ocn,
253         fizhi-gridalt-hs and global_ocean.cs32x15 (& .thsice).
254    o fix few little things to avoid floating error with -ftrapuv (ifort v9)
255     (ini_spherical_polar_grid.F, calc_r_star.F, monitor/mon_surfcor.F)
256    
257    checkpoint57i_post
258    o Diagnostics package was changed: diagnostic allocation is now deactivated
259      by setting an empty file name in the namelist (used to be by setting 0 frequency)
260    o momentum advection (VI form) options: activate 4th.O vorticity Coriolis ;
261      add 1rst Order upwind vertical Shear.
262    o add a quick fix for A/B checkpointing in OFFLINE mode
263    o add IGNORE_TIME define and -ignore_time option to genmake2 as an
264      end-run around the C/Fortran linking horror in eesup/src/timers.F
265    o added pkg/runclock: allows user to terminate run based on wall clock
266    o added non-hydrostatic scaling parameter: nh_am2
267    o pkg/my82: futile attempt to make results compiler-independent by adding
268      some "_d 0", etc. no effect.
269    o for mnc output, fill the 'T' coordinate var with myTime and create a
270      separate 'iter' variable for iteration count
271    o pkg/diagnostics:
272      add new capability: output of Global/Regional & per-level statistics ;
273      tested against monitor output (including land & thSIce) on various exp.
274    o make seaice-pkg I/O type flags (MNC/MDS) more standard.
275    
276    checkpoint57h_done
277      tag after changing user customized files (DIAGNOSTICS_SIZE.h,
278      data.diagnostics, data.gchem, data.ptracers) is done.
279    
280    o pkg/diagnostics: prepare for Global/Regional statistics.
281      add 2nd namelist in data.diagnostics and update DIAGNOSTICS_SIZE.h
282    o remove "baseTime" (no used) from arg. list of DIFF_BASE_MULTIPLE
283      and rename it: DIFFERENT_MULTIPLE
284    
285    checkpoint57h_post
286    o move tIter0 from data.gchem to PTRACERS_Iter0 in data.ptracers
287    o pkg/diagnostics: change output frequency, now in seconds (instead
288      of Nb of iterations). This requires to change data.diagnostics.
289    checkpoint57h_pre
290    
291    o in mnc, let the user create the output dir if both mnc_outdir_date
292      and mnc_outdir_num are false
293    o add fflags() to pkg/diagnostics which allows one to specify, on a
294      per-file basis, the file precision [default behavior is unchanged]
295    o added diagnostics to pkg/kpp, including computation of mixed layer
296      depth based on a temperature/density criterion
297    o updated verification/natl_box to test the new pkg/kpp diagnostics
298    o bugs fixed (with useAbsVorticity=T & with useJamartMomAdv=T)
299    o test AB-3 code with 4th.Order advection scheme (advect_xy/input.ab3_c4)
300    o fix mnc flags in pkg/seaice so that they are "owned" by seaice,
301      not mnc -- this mistake was found by JMC and is now corrected
302    
303    checkpoint57g_post
304    o add Adams-Bashforth 3rd Order time stepping:
305      - default remains AB-2 (#undef ALLOW_ADAMSBASHFORTH_3)
306      - not yet implemented with NH, CD_Scheme, pTracers, MNC_pickup ...
307    checkpoint57g_pre
308    
309    o make lambda[Theta,Salt]ClimRelax spatially varying
310      (default remains spatially uniform constant based on tau)
311    o added new control variables
312      * init. uVel, vVel, etanN
313      * lambda[Theta,Salt]ClimRelax
314    o separate masks used for ctrl_pack/unpack 'from write_grid' output
315      (suggested by G. Forget)
316    
317    checkpoint57f_post
318    o add parameter "baseTime" = model base time = time at iteration zero
319      new function DIFF_BASE_MULTIPLE (replace DIFFERENT_MULTIPLE) to
320      account for baseTime shift
321    o move Coriolis arrays from PARAMS.h to GRID.h
322    checkpoint57f_pre
323    
324    o mnc-ifying the seaice package
325    o introduced an alternative length scale calculation for use in the limiting
326      maximum viscosity (off by default, see mom_vi_hdissip.F for details)
327    o S/R ptracers_calc_diff: argument list bug for CALL KPP_CALC_DIFF_S
328    o added modified Leith scheme with viscosity proportional to grad(Hdiv)
329      as well as grad(vort3).  Made it compatible with both approximate form
330      for the cubed sphere and full form of these gradients.  Added flag
331      useFullLeith to select between approximate and full form.  
332    o add "d0" in getcon.F to reduce differences between compilers/platforms
333      update output.txt of exp. fizhi-gridalt-hs
334    o change port_rand funtion to allow to set the seed.
335    o fix obscure MNC bug involving conversion between reals and integers
336      - it caused no problems with certain (eg. gcc, ifc) compilers but the
337        new version is now more standards-compliant
338    o adjoint-related modifications
339      * gad_dst3_adv_?: to be able to use advscheme=30 in adjoint
340      * gad_calc_rhs: fool adjoint, revert to advscheme=30 in reverse when using 33
341      * mom_vi_hdissip: to be able to use Leith scheme in adjoint
342    
343    checkpoint57e_post
344    o Adding eddy stress controls a la Ferreira et al.
345    checkpoint57e_pre (actually was tagged: eckpoint57e_pre)
346    
347  o fizhi pkg: fix AtoC and CtoA when called with same input & output array  o fizhi pkg: fix AtoC and CtoA when called with same input & output array
348    update output.txt of exp. fizhi-gridalt-hs    update output.txt of exp. fizhi-gridalt-hs

Legend:
Removed from v.1.460  
changed lines
  Added in v.1.581

  ViewVC Help
Powered by ViewVC 1.1.22