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

Legend:
Removed from v.1.319  
changed lines
  Added in v.1.407

  ViewVC Help
Powered by ViewVC 1.1.22