/[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.322 by jmc, Tue Jul 6 02:44:27 2004 UTC revision 1.416 by jmc, Fri Dec 17 00:52:42 2004 UTC
# Line 3  $Header$ Line 3  $Header$
3      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
4      ==============================      ==============================
5    
6    o add diagnostics of advective & diffusive flux for T & S (pTracers coming soon)
7      this is based directly on what the model compute, so that the content of
8      each diagnostics might be affected by switching options (e.g., implicitDiffusion)
9    o small changes for diagnostics (gad_biharm_?.F, arguments of S/R impldiff.F)
10    o clean up build options files for SunOS, in particular
11      - remove -traditional from xmakedepend
12      - change -xarch=v9 to -xarch=native for f77
13      - restore use of make (not gmake)
14      - add CC=gcc if using g77
15    o changes per discussion of AM, JMC, & EH3 today:
16      - move redundant "grid" variables out of "state" -- in preparation of
17        creating proper CF-style coordinate variables and in recognition
18        that any/all grid variables can be easily copied (eg. w/NCO) from
19        the grid files into any of the other NetCDF files
20      - *always* output grid variables with MDSIO (this will eventually be
21        changed)
22      - output MNC grid file if useMNC is true
23    o every instance of _END_MASTER() has been found and replaced with
24      _END_MASTER( myThid ) in order to satisfy certain picky Sun
25      preprocessors
26    o re-arrange diagnostics pkg initialisation:
27     - allow each package to extend the available diagnostics list
28     - add some checking and fix small problems (multi-threaded, ...)
29     - new S/R diagnostics_fill (replace fill_diagnostics):
30       look through the short list of active diag. (instead of the long list)
31       create function DIAGNOSTICS_IS_ON to tell if a diagnostics is active
32    
33    checkpoint57a_post
34    o this set of changes restores TAMC compatibility
35      tested for global_ocean.90x40x15 adjoint on trough
36      (what a pain)
37    
38    checkpoint57a_pre
39    o some modifs in tools/build_options for sunos
40    o call ctrl_pack for ctrl at end for optimcycle=0
41    
42    checkpoint57
43    o add horizontal bi-harmonic viscosity for vertical velocity
44      - new parameter viscA4W (defaults to viscA4)
45      - affects only non-hydrostatic code
46      - to be done: no-slip lateral boundary conditions
47    o switch back to "sflux = sflux - runoff" in exf_getforcing.F
48    o implement Implicit Vertical advection for pTracers
49    o depth convergence accelerator: replace deltaTtracer by dTtracerLev(k)
50      (implemented also for implicit vertical diffusion & advection)
51    o OBCS: extend application of tracer open boundary conditions
52      to a strip beyond the boundary of width Olx/Oly
53    o add call monitor to TAF set of required calls (flow directives)
54      and include in natl_box_adjoint
55    o untangle a few ifdefs related to cost, ctrl, autodiff
56      untangel a few ALLOW_ECCO_OPTIMIZATION
57    o Apply Cooper & Haines representative estimates also at sub-surface (k=Nr)
58      (theta2, salt2 in ecco_cost_weights)
59    
60    checkpoint56c_post
61    o OBCS as control variables
62      - update ctrl_ad_diff.list, obcs_ad_diff.list
63      - remove balance of obcs controls from default
64      - fix index bug nobcs in ctrl_init
65      - fix dummy fields filen in ctrl_pack
66      - add dummy weights for obcs
67    o lsopt: comment out fort.94 output
68    o optim:
69      - remove unnecessary header files
70      - adjusted namelists
71      - add xerbla.F to Makefile
72    o exf:
73      - multiple modifications of exf_check_range
74      - unit change in exf_init_runoff from [m/year] to [m/s]
75    o mdsio:
76      - i/o-vector change debugLevel from debLevA to debLevB
77    
78    checkpoint56b_post
79    o collect all ptracers output in an new routine (PTRACERS_OUTPUT) that
80      is called from DO_THE_MODEL_IO
81      - rename PTRACERS_STATV_WRITE to PTRACERS_WRITE_TIMEAVE
82      - introduce new parameter PTRACERS_dumpFreq (defaults to dumpFreq)
83      - clean up PTRACERS_WRITE_TIMEAVE a little
84    o add another verification experiment that tests both the cfc-pkg and
85      the gchem-pkg with GCHEM_SEPARATE_FORCING undefined, that is with
86      reactive tendencies treated inline with regular timestepping
87    o testreport previously did a lousy job of comparing small numbers, so
88      we (Ed and Martin) changed the  way the cut-off of 1.e-12 is
89      implemented in the little c-progam tmp_cmpnum.c which is embedded in
90      testreport.
91    o ptracers_monitor: bug fix for MPI
92    
93    checkpoint56a_post
94    o GCHEM: finish reorganizating the package
95      - forward_step calls GCHEM_CALC_TENDENDY, which computes gchemTendency
96        (introduces another 3D-array for each passive tracer, but only if
97        GCHEM_SEPARATE_FORCING is undefined. For GCHEM_SEPARATE_FORCING
98        gchemTendency is not needed because the timestep is done separately)
99      - gchemTendency is added to gPtr in GCHEM_ADD_TENDENCY
100      - GCHEM is now---more or less cleanly---separate from PTRACERS
101      - to bo done: handling of surfaceForcingPtr that is needed for KPP not
102        so clear to me how that should be done in a general way.
103    o CFC: improve formatting of output in CFC_ATMOS
104    o PTRACERS:
105      - rename GCHEM_MONITOR to PTRACERS_MONITOR and call it from MONITOR, so
106        that ALL experiments with ptracers enable can be checked. This makes
107        GCHEM_MONITOR obsolete.
108      - include a runtime parameter PTRACERS_monitorFreq that defaults to
109        monitorFreq
110      - set default PTRACERS_write_mdsio to false if PTRACERS_write_mnc is true
111    o add PTRACERS (PTR0[1-5]) to testreport output
112    o exf_check_range:
113      * default changed to .TRUE.
114      * will check values for niter0 if debLevA, every timestep if debLevB
115      * STOP in exf_check_range if 'out of range'
116      * range chosen generous to allow for potential extema
117      TBD: a version with stricter range
118           permitted range would be runtime parameters
119    o GMredi:
120      - fix ldd97 slope limit ; extend valid domain of Psi-Bolus ;
121      - change S/R GRAD_SIGMA to be valid in the overlap on CS-grid.
122      - change the default GMREDI_OPTIONS.h: compile everything except Visbeck
123    o GCHEM:
124      - move call to GCHEM_FORCING_SEP, so that it is now called before
125        the blocking exchanges to avoid computing halos in GCHEM_FORCING_SEP
126      - introduced a useGCHEM runtime flag and applied it everywhere.
127      - renamed PTRACERS_SEPARATE_FORCING to GCHEM_SEPARATE_FORCING;
128        this flag is now only visible within GCHEM-pkg routines.
129      - moved gchem_forcing_int temporarily into ptracers_forcing, this will
130        be replaced by a 3D array that is computed before thermodynamics; for
131        this, the tracer ID iTracer had to be included into the parameter list
132        of ptracers_forcing
133      - set up forward_step to call gchem_forcing, that will replace
134        gchem_forcing_int in ptracers_forcing
135      - cleaning up
136    o small cleaning of exch2: remove files that are also in eesupp/src
137    
138    checkpoint56
139    o AIM: read land-Fraction earlier (from ini-fixed) and monitor land at nIter0
140    o mnc: cleanup "full" files--fixes out-of-memory problem
141    o ECCO-GODAE
142      * model/, pkg/, code/, input/ now complete w.r.t. ECCO-GODAE production setup
143      * fully merged and updated from SIO code
144        windstress sign change affects
145        - exf_inscal_[u/v]stress
146        - sign of xx_tau[u/v] adjustments
147        - signs in cost_scat
148      * ctrl_get_gen: avoid update of tau[u/v] for first few days
149      * natl_box_adjoint updated and compressed
150      * 2-level checkpoint switch added
151        (tested in EXP. 1 of natl_box_adjoint)
152      * new switch doInitXX to avoid init. of xx_... in optimcycle=0
153    o SVD
154      * first complete, non-hacked version for SVD calculations
155      * all interfaces clean w.r.t to handling of xx_..., adxx_..., g_xx_...
156        performed via direct access files admtlm_vector.it????
157        (via ctrl_[un]pack)
158      * optimcycle is used for ARPACK iteration counter
159        (note: adxx_... is iterated, rather than xx_...; latter is constant!)
160      * EV problem ok w.r.t. L2-norm in tangent space; need to check other norms
161        and add linear transforms where necessary
162      * controls are initial/final SST, SSS (can be extended; need to check
163        ctrl_set_[un]pack_xyz.F according to ctrl_set_[un]pack_xy.F
164    o mdsio: filename info to both stderr, stdout
165    o exf: exf_inscal_runoff: separated default conversion [m/year] -> [m/sec]
166      from sign change
167    
168    checkpoint55j_post
169    o isolate dissipation tendency (allow to keep it out off AB) in V.I form.
170      affects (truncation level) output.txt of global_ocean.cs32x15(&.thsice)
171    o add consistent build directories and .cvsignore files to verification
172    o added viscA4GridMin and viscA4GridMax for more control on biharmonic viscosity
173    o make all subroutines 30 or less characters in length for certain older
174      compilers to work (a thorough search only turns up two routines in
175      pkg/ptracers)
176    o allow (vecinv only) a different horiz. viscosity for Divergence and Vorticity;
177      affects (truncation level) output.txt of global_ocean.cs32x15(&.thsice)
178    
179    checkpoint55i_post
180    o change arrays for total vertical diffusivity (2.D if explicit,
181      3.D if implicit) ; each ptracer uses its own background value.
182    o added mnc/NetCDF output capability to ocn_compon_interf
183      - added bits to verification/cpl_aim+ocn so that it can run on
184        ITRDA and other clusters without the mpirun "-wd" option
185    o fixed wayward endif in mom_fluxform.F that messed up cyclindrical
186      configuration
187    o give MNC the ability to create additional files rather than
188      exceeding a specified file size limit "MNC_MAX_FILESIZE"
189      - the default is just less than 2GB which is a practical limit
190        on many systems
191    o Compute the non-local term using either KPP vertical diffusivity (=default)
192       or the total vertical diffusivity (when KPP_ghatUseTotalDiffus=T)
193    
194    checkpoint55h_post
195    o add verification test "vermix"
196    o remove all the *_ioinc flags and replace with the single global
197      outputTypesInclusive flag
198    o allow to set a vertical profile of vertical diffusivity for T & S
199    o call PACKAGE_BOOT & PACKAGES_READPARMS just after INI_PARMS
200    o useMNC moved from data.mnc to data.pkg
201    
202    checkpoint55g_post
203    o Renaming carbon AD experiments to comply with testreport conventions
204    o Updating adjoint of Held-Suarez hs94.1x64x5
205    o Fixing directives for non-ECCO forward_step, after ECCO directives
206      have been fixed to avoid recomp. of do_oceanic_phys
207      Works well for global_ocean.90x40x15,
208      but recomp. issue remains for carbon
209    o Updating adjoint output for natl_box_adjoint after
210      pkg/ecco has been modified (see taga c55d_post)
211    
212    checkpoint55f_post
213    o rename to timeave_init_fixed
214    
215    checkpoint55e_post
216    o some delicate re-shuffle of store directives to avoid one
217      extra call of do_oceanic_physics
218      affects ecco_the_main_loop.F, do_oceanic_physics.F
219      checkpoint_lev1_directives.h
220    o NB: this may break global_ocean adjoint temporarily,
221      but it is clear how to fix it. Will do later, need this now.
222    
223    checkpoint55d_post
224    o ECCO specific cost function terms (up-to-date with 1x1 runs)
225    o ecco_cost_weights is modified to 1x1 runs
226    o modifs to allow observations to be read in as
227      single file or yearly files
228    o enable to read exf forcing fields as either
229      single file or yearly files (flag useExfYearlyFields)
230    checkpoint55d_pre
231    
232    o switch between 2-level and 3-level checkpointing
233    o add ability of MNC to write local and "partial" (eg. 2D slices where
234      the full 3D field is never actually stored) arrays to NetCDF files
235      with the correct (that is, the complete multi-dimensional) set of
236      array indicies
237      - used in mom_vecinv() to write the diagFreq output
238      - tested (demonstrated) in verification/aim.5l_cs
239    o further mnc-ification
240      - add verification/lab_sea to the list with mnc setups
241      - fix undefined mnc v-names in cd_code
242      - add kpp to mnc output (both snapshot and timeave)
243    o non-hydrostatic code
244      - added horizontal viscosity parameter viscAhW for w, that replaces viscAh
245        in calc_gw.F
246      - defaults to viscAh
247      - this is useful, when I want to use viscAh=0 (e.g., only biharmonic
248        friction or Leith) in non-hydrostatic mode, as long as there is
249        no biharmonic friction and/or Leigth or equivalent for vertical velocity.
250    o pkg/cost
251      extend simple cost functions to treat time mean quantities
252      (T,U,V) via cost_accumulate_state from cost_tile
253    o multidimAdvection on CS-grid: implement the full 3 pass ; overlap needs to
254        be 2 times larger than the minimum size the advection scheme requires.
255        now get a perfect conservation of salt in advect_cs test-exp (updated)
256    o GGL90 (Gaspar's TKE mixed layer)
257      - add horizontal diffusion of TKE, requires exchanges in
258        do_fields_blocking_exchanges, horizontal diffusivity is zero by
259        default. In OPA there is no horizontal diffusion of TKE but the
260        mixing coefficients are computed from a horizontal average of TKE of
261        6 points or so. I think that diffusion has a little more physical
262        justification.
263      - clean up ggl90_calc in the hope of reducing memory usage (this hope
264        was in vain)
265      - mask tke-variable in ggl90_init
266    o fix mask Pb in multidimAdvection on CS-grid (still a small Pb with ipass)
267    o fix biharmonic diffusion on CS-grid.
268    o more obcs stuff
269      - added obcs_external_fields_load as an alternative to the EXF-pkg option
270        for reading boundary values. This routine is called from
271        obcs_prescribe_read
272      - updated obcs_check
273      - changed sponge layer relaxation parameter units for days to seconds
274      - turned off sponge layer by default
275      - changed argument list of OBCS_CALC to myTime+deltaTclock
276      - fixed about 35 bugs in obcs_external_fields_load
277    
278    checkpoint55c_post
279    o obcs stuff
280      - new runtime flags useOBCSprescribe, useOBCSbalance
281      - re-arranged obcs_calc
282      - moved exf_getobcs to obcs_prescibe_read
283        obcs_prescibe_read is called after ORLANSKI stuff
284        and independent of it (executable can contain all).
285      - updated obcs_check
286      - added exchanges for U,V,gU,gV for useOBCSprescribe=.TRUE.
287        These may be necessary only if tangential velocities are
288        prescribed at OB, not sure, needs testing.
289    o make the MNC_CW_ADD_VATTR_* subroutines all take scalar arguments
290      which simplifies the calls and passes the ifc check-bounds test
291    
292    checkpoint55b_post
293    o add a 2D scalar, Z location, exchange_2 routine.
294    o EXCH_UV for CS-grid: add one u,v in the corner-halo region, allows
295      to compute vorticity on a wider stencil [e.g., vort3(0,1)&(1,0)]
296      affects advect_cs results (for the wrong reason): update output.txt
297    o use a more standard discretization for biharmonic viscosity ;
298     (original version still available with #define MOM_VI_ORIGINAL_VISCA4 )
299      affects results of dome set-up: update output.txt
300    o fix gad_advection.F (CubedSphere part), back to version 1.23
301    o fix bug in mom_vi_del2uv.F (CubedSphere part).
302    
303    checkpoint55a_post
304    o merged code to
305      * prescribe/read time-dependent open boundaries
306        (works in conjunction with exf, cal)
307      * sponge layer code for open boundaries
308      * each boundary N/S/E/W now has its own CPP option
309        (healthy for the adjoint)
310    
311    checkpoint55
312    o up-to-date adjoint of previous breaks
313      - natl_box_adjoint: problem for exp 12, 13: cost = 0. needs fix
314      - global_ocean_ebm: adjoint still needs fix (TLM ok)
315      - global_ocean.90x40x15: extended GMREDI_OPTIONS.h
316      - carbon: * remove tr1-related adjoint verification
317                * add ptracer as control variable
318    o remove all tr1-related code (ALLOW_PASSIVE_TRACER)
319      (adjoint stuff still has some tr1 'names', but all use ptracer arrays)
320    
321    checkpoint54f_post
322    o in verification/carbin/ remove experiment related to
323      ALLOW_PASSIVE_TRACER (array tr1 etc)
324      code_ad/ input_ad/ results_ad
325    
326    o add two new packages
327      - ggl90, TKE-model of Gaspar et al. (1990), Blanke+Delecluse (1993)
328        + subject to optimization (memory)
329      - opps, OPPS (ocean penetrative plume scheme) of Paluszkiewicz+Romea (1997)
330        + subject to reordering (many GOTOs in the original code)
331        + so far opps is called for every (i,j) point, and then operates in
332          the vertical; one may want to reorder the scheme in order operate in
333          horizontal slabs at the cost of (a lot of) memory, at least ten new 3D
334          fields)
335    o overhaul of IO so that we now have flags for MDSIO and/or MNC
336      - all verification tests compile and run with linux_ia32_g77
337      - defaults are compatible with current input files--nothing
338        should change if you were not previously using MNC
339      - MNC output has been added in numerous places (eg. timeave)
340        but there are still a few writes not yet do-able with MNC
341        (this is in progress)
342      - flags now allow for either/or/both use of MDSIO and MNC and
343        documentation will soon follow
344      - numerous small formatting cleanups for ProTeX
345    o add version info (via BUILD_INFO.h) to MNC and stdout
346    
347    checkpoint54e_post
348    o add two new packages
349      - pp81 (Packanowski and Philander, 1981), Richardson number and
350        stratification dependent mixing
351      - my82 (Mellor and Yamada, 1982) level 2 turbulence closure scheme
352      - these scheme operate in a similar manner as KPP (without ghat, of
353        course), however there is only one diffusivity that is used for both
354        temperature and salinity (and passive tracers)ç
355    o add offline package
356    o add "verification/dic_example" which contains a dic-gchem-ptracer setup
357      provided by Stephanie Dutkiewicz
358    o genmake2: add ability to search for C compilers and cleanup test output
359      - fixes recent SX-6 problems and should help with others (eg. HP-UX)
360    
361    checkpoint54d_post
362    o data.diagnostics: if levels() is unspecified then ALL levels are diagnosed
363    o bug fix in pkg/cal/cal_addtime.F for TheCalendar='model' day count is wrong,
364      leading to additional shift by 5 days each year of integration
365       (i.e. leading to considerable offset over the years)
366    o update global_ocean.cs32x15/thsice: use z* + RealFW + staggerTimeStep
367    o start with no AB and no initial filtering (staggerTimeStep):
368      change the results of test-exp: hs94.128x64x5, hs94.1x64x5 & hs94.cs-32x32x5
369    o move ocean specific S/R calls from forward_step.F to do_oceanic_phys.F
370      (external loading & cpl import/export stay in forward_step.F)
371    o update coupling interface to work with new staggerTimeStep
372    o change MDSIO ptracers_pickup.
373    o fix a bug in Flux deriv. vs Ts (pkg/aim_v23/phy_suflux_sice.F);
374    o genmake2 fixes so netcdf now works on halem ("osf1_alpha_*")
375    o bugs fixed in kpp_calc.F (ustarY) and external_forcing.F (Qsw & Partial-cell)
376    
377    checkpoint54c_post
378    o replace surfaceTendency U,V,S,T,Tice,Ptr by surfaceForcing U,V,S,T,Tice,Ptr
379      this affects (truncation error) output.txt of test-exp:
380       global_ocean_ebm,_pressure,.cs32x15(&.thsice), lab_sea & natl_box
381    o pkg/kpp: replace pMask by maskC to reduce memory size (one 3.D array)
382    o change Qnet to always be the net heat flux, (+upward).
383      affects (truncation error) output of natl_box exp.: updated
384    o define the max number of passive tracers in PTRACERS_SIZE.h
385    o start to test mnc routinely in exp0
386    o add sea-ice loading (sIceLoad) as a new forcing field.
387    
388    checkpoint54b_post
389    o add testreport "-j JOBS" option for parallel builds (tested and pretty
390       impressive on Altix and Opteron systems)
391    o diagnostics can now write NetCDF files through the mnc package
392    o staggerTimeStep & NLFS: update pkg aim_v23 & thsice; update aim.5l_cs output
393    o staggerTimeStep & multi-Dim-Advec: Exch(U,V) added
394    o staggerTimeStep & NLFS: rescale T & S surface tendencies.
395      affects (truncation error) output of lab_sea & natl_box exp.: updated
396    o change pickup _cd & _ptracers (old pickup restart: usePickupBeforeC54=.T.)
397    
398  checkpoint54a_post  checkpoint54a_post
399  o re-write staggerTimeStep: step fwd U,V and then T,S (reverse the order).  o re-write staggerTimeStep: step fwd U,V and then T,S (reverse the order).
400  o change pickup file (to restart from old pickup, set usePickupBeforeC54=.T.)  o change pickup file (to restart from old pickup, set usePickupBeforeC54=.T.)

Legend:
Removed from v.1.322  
changed lines
  Added in v.1.416

  ViewVC Help
Powered by ViewVC 1.1.22