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

Legend:
Removed from v.1.240  
changed lines
  Added in v.1.364

  ViewVC Help
Powered by ViewVC 1.1.22