/[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.256 by adcroft, Mon Feb 23 18:57:48 2004 UTC revision 1.383 by jmc, Wed Nov 10 19:36:30 2004 UTC
# Line 3  $Header$ Line 3  $Header$
3      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
4      ==============================      ==============================
5    
6    checkpoint55j_post
7    o isolate dissipation tendency (allow to keep it out off AB) in V.I form.
8      affects (truncation level) output.txt of global_ocean.cs32x15(&.thsice)
9    o add consistent build directories and .cvsignore files to verification
10    o added viscA4GridMin and viscA4GridMax for more control on biharmonic viscosity
11    o make all subroutines 30 or less characters in length for certain older
12      compilers to work (a thorough search only turns up two routines in
13      pkg/ptracers)
14    o allow (vecinv only) a different horiz. viscosity for Divergence and Vorticity;
15      affects (truncation level) output.txt of global_ocean.cs32x15(&.thsice)
16    
17    checkpoint55i_post
18    o change arrays for total vertical diffusivity (2.D if explicit,
19      3.D if implicit) ; each ptracer uses its own background value.
20    o added mnc/NetCDF output capability to ocn_compon_interf
21      - added bits to verification/cpl_aim+ocn so that it can run on
22        ITRDA and other clusters without the mpirun "-wd" option
23    o fixed wayward endif in mom_fluxform.F that messed up cyclindrical
24      configuration
25    o give MNC the ability to create additional files rather than
26      exceeding a specified file size limit "MNC_MAX_FILESIZE"
27      - the default is just less than 2GB which is a practical limit
28        on many systems
29    o Compute the non-local term using either KPP vertical diffusivity (=default)
30       or the total vertical diffusivity (when KPP_ghatUseTotalDiffus=T)
31    
32    checkpoint55h_post
33    o add verification test "vermix"
34    o remove all the *_ioinc flags and replace with the single global
35      outputTypesInclusive flag
36    o allow to set a vertical profile of vertical diffusivity for T & S
37    o call PACKAGE_BOOT & PACKAGES_READPARMS just after INI_PARMS
38    o useMNC moved from data.mnc to data.pkg
39    
40    checkpoint55g_post
41    o Renaming carbon AD experiments to comply with testreport conventions
42    o Updating adjoint of Held-Suarez hs94.1x64x5
43    o Fixing directives for non-ECCO forward_step, after ECCO directives
44      have been fixed to avoid recomp. of do_oceanic_phys
45      Works well for global_ocean.90x40x15,
46      but recomp. issue remains for carbon
47    o Updating adjoint output for natl_box_adjoint after
48      pkg/ecco has been modified (see taga c55d_post)
49    
50    checkpoint55f_post
51    o rename to timeave_init_fixed
52    
53    checkpoint55e_post
54    o some delicate re-shuffle of store directives to avoid one
55      extra call of do_oceanic_physics
56      affects ecco_the_main_loop.F, do_oceanic_physics.F
57      checkpoint_lev1_directives.h
58    o NB: this may break global_ocean adjoint temporarily,
59      but it is clear how to fix it. Will do later, need this now.
60    
61    checkpoint55d_post
62    o ECCO specific cost function terms (up-to-date with 1x1 runs)
63    o ecco_cost_weights is modified to 1x1 runs
64    o modifs to allow observations to be read in as
65      single file or yearly files
66    o enable to read exf forcing fields as either
67      single file or yearly files (flag useExfYearlyFields)
68    checkpoint55d_pre
69    
70    o switch between 2-level and 3-level checkpointing
71    o add ability of MNC to write local and "partial" (eg. 2D slices where
72      the full 3D field is never actually stored) arrays to NetCDF files
73      with the correct (that is, the complete multi-dimensional) set of
74      array indicies
75      - used in mom_vecinv() to write the diagFreq output
76      - tested (demonstrated) in verification/aim.5l_cs
77    o further mnc-ification
78      - add verification/lab_sea to the list with mnc setups
79      - fix undefined mnc v-names in cd_code
80      - add kpp to mnc output (both snapshot and timeave)
81    o non-hydrostatic code
82      - added horizontal viscosity parameter viscAhW for w, that replaces viscAh
83        in calc_gw.F
84      - defaults to viscAh
85      - this is useful, when I want to use viscAh=0 (e.g., only biharmonic
86        friction or Leith) in non-hydrostatic mode, as long as there is
87        no biharmonic friction and/or Leigth or equivalent for vertical velocity.
88    o pkg/cost
89      extend simple cost functions to treat time mean quantities
90      (T,U,V) via cost_accumulate_state from cost_tile
91    o multidimAdvection on CS-grid: implement the full 3 pass ; overlap needs to
92        be 2 times larger than the minimum size the advection scheme requires.
93        now get a perfect conservation of salt in advect_cs test-exp (updated)
94    o GGL90 (Gaspar's TKE mixed layer)
95      - add horizontal diffusion of TKE, requires exchanges in
96        do_fields_blocking_exchanges, horizontal diffusivity is zero by
97        default. In OPA there is no horizontal diffusion of TKE but the
98        mixing coefficients are computed from a horizontal average of TKE of
99        6 points or so. I think that diffusion has a little more physical
100        justification.
101      - clean up ggl90_calc in the hope of reducing memory usage (this hope
102        was in vain)
103      - mask tke-variable in ggl90_init
104    o fix mask Pb in multidimAdvection on CS-grid (still a small Pb with ipass)
105    o fix biharmonic diffusion on CS-grid.
106    o more obcs stuff
107      - added obcs_external_fields_load as an alternative to the EXF-pkg option
108        for reading boundary values. This routine is called from
109        obcs_prescribe_read
110      - updated obcs_check
111      - changed sponge layer relaxation parameter units for days to seconds
112      - turned off sponge layer by default
113      - changed argument list of OBCS_CALC to myTime+deltaTclock
114      - fixed about 35 bugs in obcs_external_fields_load
115    
116    checkpoint55c_post
117    o obcs stuff
118      - new runtime flags useOBCSprescribe, useOBCSbalance
119      - re-arranged obcs_calc
120      - moved exf_getobcs to obcs_prescibe_read
121        obcs_prescibe_read is called after ORLANSKI stuff
122        and independent of it (executable can contain all).
123      - updated obcs_check
124      - added exchanges for U,V,gU,gV for useOBCSprescribe=.TRUE.
125        These may be necessary only if tangential velocities are
126        prescribed at OB, not sure, needs testing.
127    o make the MNC_CW_ADD_VATTR_* subroutines all take scalar arguments
128      which simplifies the calls and passes the ifc check-bounds test
129    
130    checkpoint55b_post
131    o add a 2D scalar, Z location, exchange_2 routine.
132    o EXCH_UV for CS-grid: add one u,v in the corner-halo region, allows
133      to compute vorticity on a wider stencil [e.g., vort3(0,1)&(1,0)]
134      affects advect_cs results (for the wrong reason): update output.txt
135    o use a more standard discretization for biharmonic viscosity ;
136     (original version still available with #define MOM_VI_ORIGINAL_VISCA4 )
137      affects results of dome set-up: update output.txt
138    o fix gad_advection.F (CubedSphere part), back to version 1.23
139    o fix bug in mom_vi_del2uv.F (CubedSphere part).
140    
141    checkpoint55a_post
142    o merged code to
143      * prescribe/read time-dependent open boundaries
144        (works in conjunction with exf, cal)
145      * sponge layer code for open boundaries
146      * each boundary N/S/E/W now has its own CPP option
147        (healthy for the adjoint)
148    
149    checkpoint55
150    o up-to-date adjoint of previous breaks
151      - natl_box_adjoint: problem for exp 12, 13: cost = 0. needs fix
152      - global_ocean_ebm: adjoint still needs fix (TLM ok)
153      - global_ocean.90x40x15: extended GMREDI_OPTIONS.h
154      - carbon: * remove tr1-related adjoint verification
155                * add ptracer as control variable
156    o remove all tr1-related code (ALLOW_PASSIVE_TRACER)
157      (adjoint stuff still has some tr1 'names', but all use ptracer arrays)
158    
159    checkpoint54f_post
160    o in verification/carbin/ remove experiment related to
161      ALLOW_PASSIVE_TRACER (array tr1 etc)
162      code_ad/ input_ad/ results_ad
163    
164    o add two new packages
165      - ggl90, TKE-model of Gaspar et al. (1990), Blanke+Delecluse (1993)
166        + subject to optimization (memory)
167      - opps, OPPS (ocean penetrative plume scheme) of Paluszkiewicz+Romea (1997)
168        + subject to reordering (many GOTOs in the original code)
169        + so far opps is called for every (i,j) point, and then operates in
170          the vertical; one may want to reorder the scheme in order operate in
171          horizontal slabs at the cost of (a lot of) memory, at least ten new 3D
172          fields)
173    o overhaul of IO so that we now have flags for MDSIO and/or MNC
174      - all verification tests compile and run with linux_ia32_g77
175      - defaults are compatible with current input files--nothing
176        should change if you were not previously using MNC
177      - MNC output has been added in numerous places (eg. timeave)
178        but there are still a few writes not yet do-able with MNC
179        (this is in progress)
180      - flags now allow for either/or/both use of MDSIO and MNC and
181        documentation will soon follow
182      - numerous small formatting cleanups for ProTeX
183    o add version info (via BUILD_INFO.h) to MNC and stdout
184    
185    checkpoint54e_post
186    o add two new packages
187      - pp81 (Packanowski and Philander, 1981), Richardson number and
188        stratification dependent mixing
189      - my82 (Mellor and Yamada, 1982) level 2 turbulence closure scheme
190      - these scheme operate in a similar manner as KPP (without ghat, of
191        course), however there is only one diffusivity that is used for both
192        temperature and salinity (and passive tracers)ç
193    o add offline package
194    o add "verification/dic_example" which contains a dic-gchem-ptracer setup
195      provided by Stephanie Dutkiewicz
196    o genmake2: add ability to search for C compilers and cleanup test output
197      - fixes recent SX-6 problems and should help with others (eg. HP-UX)
198    
199    checkpoint54d_post
200    o data.diagnostics: if levels() is unspecified then ALL levels are diagnosed
201    o bug fix in pkg/cal/cal_addtime.F for TheCalendar='model' day count is wrong,
202      leading to additional shift by 5 days each year of integration
203       (i.e. leading to considerable offset over the years)
204    o update global_ocean.cs32x15/thsice: use z* + RealFW + staggerTimeStep
205    o start with no AB and no initial filtering (staggerTimeStep):
206      change the results of test-exp: hs94.128x64x5, hs94.1x64x5 & hs94.cs-32x32x5
207    o move ocean specific S/R calls from forward_step.F to do_oceanic_phys.F
208      (external loading & cpl import/export stay in forward_step.F)
209    o update coupling interface to work with new staggerTimeStep
210    o change MDSIO ptracers_pickup.
211    o fix a bug in Flux deriv. vs Ts (pkg/aim_v23/phy_suflux_sice.F);
212    o genmake2 fixes so netcdf now works on halem ("osf1_alpha_*")
213    o bugs fixed in kpp_calc.F (ustarY) and external_forcing.F (Qsw & Partial-cell)
214    
215    checkpoint54c_post
216    o replace surfaceTendency U,V,S,T,Tice,Ptr by surfaceForcing U,V,S,T,Tice,Ptr
217      this affects (truncation error) output.txt of test-exp:
218       global_ocean_ebm,_pressure,.cs32x15(&.thsice), lab_sea & natl_box
219    o pkg/kpp: replace pMask by maskC to reduce memory size (one 3.D array)
220    o change Qnet to always be the net heat flux, (+upward).
221      affects (truncation error) output of natl_box exp.: updated
222    o define the max number of passive tracers in PTRACERS_SIZE.h
223    o start to test mnc routinely in exp0
224    o add sea-ice loading (sIceLoad) as a new forcing field.
225    
226    checkpoint54b_post
227    o add testreport "-j JOBS" option for parallel builds (tested and pretty
228       impressive on Altix and Opteron systems)
229    o diagnostics can now write NetCDF files through the mnc package
230    o staggerTimeStep & NLFS: update pkg aim_v23 & thsice; update aim.5l_cs output
231    o staggerTimeStep & multi-Dim-Advec: Exch(U,V) added
232    o staggerTimeStep & NLFS: rescale T & S surface tendencies.
233      affects (truncation error) output of lab_sea & natl_box exp.: updated
234    o change pickup _cd & _ptracers (old pickup restart: usePickupBeforeC54=.T.)
235    
236    checkpoint54a_post
237    o re-write staggerTimeStep: step fwd U,V and then T,S (reverse the order).
238    o change pickup file (to restart from old pickup, set usePickupBeforeC54=.T.)
239    o get a perfect restart when using RealFreshWaterFlux + NLFS (not finished)
240    o update output of all exp. that use staggerTimeStep (except hs94.cs & 2D)
241      + internal_wave (shift in writing cg2d_ini_res)
242    checkpoint54a_pre
243    
244    o require forcing_In_AB=.FALSE. and remove surfaceTendencyTice from pkg/seaice
245      - modified verification/lab_sea and natl_box_adjoint accordingly
246    
247    checkpoint54
248    checkpoint53g_post
249    o prepare splitting of thermodynamics: store convect.counter in common block
250    o store d.etaH/dt (instead of Div.hV) in common ; affects (truncation error)
251      output of: global_ocean.cs32x15, global_ocean_pressure & internal_wave
252    o cnh's modifs to gad_advection.F for cube-sphere multi-dim advection
253    o new parameter "nShapS": allow to treat differently S from T (nShapT)
254    o separate Vert.Advec.Scheme from horizontal Advec.Scheme (T & S)
255    o restore adjointability
256    
257    checkpoint53f_post
258    o modifications that affect AIM experiment results (update output.txt):
259      - no advection of S at top interface Nr (pkg/generic_advdiff)
260      - air surf. temp: use lower value between T0 & T1 (phy_suflux_prep.F)
261    o AIM pkg:
262      - include stability function into surf.Flux derivative relative to Tsurf
263      - calculate clear-sky radiation & surface temp. change (for diagnostics)
264    
265    checkpoint53e_post
266    o merged cylindrical coord configuration
267    o added rotating_tank experiment
268    checkpoint53e_pre
269    
270    checkpoint53d_post
271    o faulks upgraded: re-generate output.txt that have changed:
272       global_ocean.cs32x15(.thsice), global_ocean_ebm, lab_sea, natl_box
273    o Added Bryan and Lewis, 1979, vertical diffusivity profile
274    o pkg/land: fix problem in land_impl_temp.F when Ts arround freezing.
275    o added Sadourny discretization of Coriolis
276      - moved some PARAMETERS from mom_*_coriolis.F to PARAMS.h
277      - re-enabled use of omega3 (useAbsVorticity)
278    o added Leith (Phys. Fluids, 1968, 10, 1409-1416) variable viscosity for
279      vector-invariant equations
280      - note: does *NOT* yet use same variable visc. for no-slip BC's.
281    o recoded the lateral boundary conditions for non-hydrostatic w, this time
282      hopefully without new bugs
283      - still to do: biharmonic friction for non-hydrostatic w
284    o bug fix for mom_*_sidedrag: rA -> rAs,rAw
285    o regenerated output for global_ocean_ebm (using -ieee)
286    o fix mom_common/mom_*_sidedrag.F to recover numeric output prior to creation
287      of mom_common/
288      - simple re-arrangement of instructions within computing loop
289    checkpoint53d_pre
290    
291    checkpoint53c_post
292    o new version of coupler for AIM+thSIce and Ocean models
293     - new fields (SSS,h_MixLayer ...) are exchanged through the coupler
294     - land: change Evaporation of snow
295     - aim with thSIce: fix bug in snow-precip.
296    o changes to calc_gw.F
297     - bug fix for w^2 term near the bottom boundary
298     - (hopefully) improve the lateral slip boundary condtions for use with
299        partial/looped cells
300      - because the bug fix changes two verification experiments anyway (exp5
301        and plume on slope), change the lateral boundary condition from half slip
302        to the value of no_slip_sides
303    o ebm package wasnt quite up to the parallel task. Fixed that.
304    o Adjoint of EBM package + verification
305      * TLM OK
306      * ADM not sufficient accuracy yet (needs further testing)
307    
308    checkpoint53b_post
309    o Commiting new energy balance model to repository
310      * package is pkg/ebm
311      * verif. is verification/global_ocean_ebm
312      * references are in ebm_driver.F
313    C      * X. Wang, P. Stone and J. Marotzke, 1999:
314    C        Global thermohaline circulation. Part I:
315    C        Sensitivity to atmospheric moisture transport.
316    C        J. Climate 12(1), 71-82
317    C      * X. Wang, P. Stone and J. Marotzke, 1999:
318    C        Global thermohaline circulation. Part II:
319    C        Sensitivity with interactive transport.
320    C        J. Climate 12(1), 83-91
321    C      * M. Nakamura, P. Stone and J. Marotzke, 1994:
322    C        Destabilization of the thermohaline circulation
323    C        by atmospheric eddy transports.
324    C        J. Climate 7(12), 1870-1882
325    o Will need long integration testing.
326    checkpoint53b_pre
327    
328    o created pkg/mom_common in preparation for variable viscosity code
329    o renamed and split ini_mnc_io
330    o pkg/land: only liquid water can diffuse or run-off
331    o AIM using land: soil wetness availability account for snow
332    o added printMapIncludesZeros to EEPARMS namelist
333    
334    checkpoint53a_post
335    o restore TAMC compatiility for exf_constants.h
336      (no usage of _d in .h)
337    o Additional storing in KPP in an attempt to fix
338      TAMC problem with adjoint of KPP code.
339      (no success...)
340      Identical code works fine with TAF.
341    o Fixed verification/natl_box_adjoint experiments.
342      Initial conditions for Experiments 4 and 5 are now
343      set to "SEAICE_initialHEFF=0", which is more reasonable,
344      and hence gives results closer to that of Experiment 2.
345      Experiment 6 requires "SEAICE_initialHEFF=1" and is working once again.
346    
347    checkpoint53
348    o stable checkpoint with up to date adjoint
349      TLM tests only run for 1st grdchk at the moment. Its unclear if
350      this is a code problem or a problem with resetting code for g_
351      variables that is being autogenerated wrong.
352      Otherwise:
353       o code is testing well on Linux (P4 and AMD), SGI Alitix and Irix.
354       o IBM and HPQ tests could not be run because NCAR (IBM) and NCCS (HPQ)
355         are both still in chaos following their security messes.
356    
357    checkpoint52n_post
358    o various build (mostly genmake2) changes including:
359      - ability to build/run more easily within Cygwin and MacOS X
360      - more complicated Makefiles with a stronger dependence on Gnu Make
361    o pkg/seaice modifications
362      - added SEAICE_deltaTtherm and SEAICE_deltaTdyn in order to reduce cost
363        of pkg/seaice.  For time being, SEAICE_deltaTtherm = deltaTtracer
364        and SEAICE_deltaTdyn must be integer multiple of SEAICE_deltaTtherm
365      - changed proxy for geostrophic velocity and updated
366        verification/lab_sea/results/* accordingly
367      - removed ADI dynamic solver from pkg/seaice
368      - replaced all occurrences of 1/CSTICE or 1/CSUICE
369        they caused a segmentation fault under certain circumstances
370      - removed unused field WATR
371    o bug fix in exf-package:
372      -  u/vwind are defined at tracer/pressure points. The values of u/vwindmask
373         now reflect that.
374    o small fixes in the cal-package:
375      - set reference date to first of January for the 'model'-calendar
376      - failing of cal_CheckDate in cal_FullDate now actually results in a
377        warning
378    o small modifications in exf-package to increase flexibility:
379      - change climtempfreeze from a constant to a runtime parameter that
380        is specified in data.exf_clim (default = -1.9)
381      - introduce runtime parameter windstressmax (default = 2.0) that is
382        specified in data.exf
383    o fix in pkg/thsice:
384      - update Evap (Tsf changes) to be consistent with Latent heat flux
385      - maintain snow precipitation unchanged (even if Tsf >=0).
386     both affect global_ocean.cs32x15/(.thsice) results ; update output.txt.thsice
387    
388    checkpoint52m_post
389    o commented out call to cal_GetMonthsRec in swfrac.F and added some comments;
390      this does not change any of the results, because myIter=0 in swfrac always
391      lead to jwtype=2 (the default anyway), but the call makes little sense
392      as a default
393    o no SHORTWAVE HEATING lost at the bottom; affects lab_sea output; updated
394    o re-work pkg/thsice interface with ocean & bulk_force pkg to allow
395      atmospheric model (AIM) to use it ; add a "slab ocean" component.
396    o pkg/aim_v23 : add interface to call thermodynamic sea-ice (thsice)
397    o update global_ocean.cs32x15/output (.thsice) after changes in pkg/thsice
398    o extending useSingleCpuIO option to work with new exch2 I/O format
399      - old-style, missing-tile I/O is still accessible by defining CPP
400        option MISSING_TILE_IO in pkg/mdsio/MDSIO_OPTIONS.h
401    o monitor package can write to NetCDF files and/or STDOUT
402    o fix re-definition bug in MNC_GRID_INIT_ALL()
403    o doc/api_reference added
404      - a framework for building an API Reference Manual using all
405        the "protex" comments embedded in the code
406      - see pkg/generic_advdiff, pkg/mnc for examples
407    o MNC:
408      - make myThid the last argument for all subroutines
409      - new default for the grid files: NO halos
410      - add useMNC and other run-time flags
411      - separate per-file-group handling of the unlimited dimension
412      - early initialization of MNC w/ flags read at the eeboot stage
413        and lookup tables set in initialize_fixed
414      - write all MNC output to a freshly-created directory
415    o thermodynamics sea-ice (thsice) is now tested with testreport.
416    o test exp. advect_cs: use the standard cs32 grid ; update output.
417    o Added functionality to grdchk:
418      pick global i,j,k position (or nearest wet) where to perform check.
419    
420    checkpoint52l_post
421    o new land formulation:
422      a) use ground enthalpy as prognostic variable to ensure exact
423         energy conservation.
424      b) account for water temperature and for latent heat of freezing
425         in all processes (rain, run-off, ground storage)
426      c) compute surface and ground temperature implicitly.
427    o aim_v23 modification for new land model:
428      a) Treat separately land / ocean / sea-ice surface fluxes
429         to allow implicit computation of land & sea-ice surface temp.
430      b) add snow precipitation.
431    checkpoint52l_pre
432    
433    o fix problem (when 1 tile is empty) in min/max monitor.
434    o fix bug in timeAve (u,vZetatave in common block)
435    o no restoring terms under ice for pkg/seaice
436    o moved SEAICE_DO_DIAGS to S/R do_the_model_io
437    o bug fix for pkg/seaice pickups
438    o change MNC call locations and add coordinate attributes
439    o mask vorticity at the corner ; affects global_ocean.cs output
440    o fix ini_curvilinear_grid ; test CS experiments with mpi
441    o aim.5l_cs: modify sea-ice fraction input file ; update output
442    o exf: Adding simple range check routine for exf fields after read
443    
444    checkpoint52k_post
445  o removed top-level directory "diags"  o removed top-level directory "diags"
446    - we migrated these functions to pkg/timeave a long time ago    - we migrated these functions to pkg/timeave a long time ago
447  o fixed bug introduced in pkg/monitor  o fixed bug introduced in pkg/monitor
448    o accuracy ctrlprec = 32 insuffient for gradient checks using
449      averaged fields (I/O via cost_averages)
450      -> use ctrl.h in active_file*.F to control I/O precision.
451    o paramter list of CALL ADTHE_MAIN_LOOP
452    o added initial do_the_model_io to ecco
453    o set all ...Freq to zero after cost_final
454      taveFreq, dumpFreq, pChkptFreq
455    o MAIN and ECCO differentiability restored
456    
457  checkpoint52j_post  checkpoint52j_post
458  o fix nchklev in tamc.h for global_ocean...  o fix nchklev in tamc.h for global_ocean...

Legend:
Removed from v.1.256  
changed lines
  Added in v.1.383

  ViewVC Help
Powered by ViewVC 1.1.22