/[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.280 by mlosch, Mon Apr 19 22:30:46 2004 UTC revision 1.502 by jmc, Mon Jun 20 23:23:33 2005 UTC
# Line 1  Line 1 
1  $Header$  $Header$
2    $Name$
3    
4      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
5      ==============================      ==============================
6    
7    o add Sadourny like discretization to the 4th.O vorticity advection scheme
8    o fix bug in rStar code (calc_r_star.F) when used with topo/bathy.
9      update output.txt from: aim.5l_cs (& .thSI), cpl_aim+ocn,
10         fizhi-gridalt-hs and global_ocean.cs32x15 (& .thsice).
11    o fix few little things to avoid floating error with -ftrapuv (ifort v9)
12     (ini_spherical_polar_grid.F, calc_r_star.F, monitor/mon_surfcor.F)
13    
14    checkpoint57i_post
15    o Diagnostics package was changed: diagnostic allocation is now deactivated
16      by setting an empty file name in the namelist (used to be by setting 0 frequency)
17    o momentum advection (VI form) options: activate 4th.O vorticity Coriolis ;
18      add 1rst Order upwind vertical Shear.
19    o add a quick fix for A/B checkpointing in OFFLINE mode
20    o add IGNORE_TIME define and -ignore_time option to genmake2 as an
21      end-run around the C/Fortran linking horror in eesup/src/timers.F
22    o added pkg/runclock: allows user to terminate run based on wall clock
23    o added non-hydrostatic scaling parameter: nh_am2
24    o pkg/my82: futile attempt to make results compiler-independent by adding
25      some "_d 0", etc. no effect.
26    o for mnc output, fill the 'T' coordinate var with myTime and create a
27      separate 'iter' variable for iteration count
28    o pkg/diagnostics:
29      add new capability: output of Global/Regional & per-level statistics ;
30      tested against monitor output (including land & thSIce) on various exp.
31    o make seaice-pkg I/O type flags (MNC/MDS) more standard.
32    
33    checkpoint57h_done
34      tag after changing user customized files (DIAGNOSTICS_SIZE.h,
35      data.diagnostics, data.gchem, data.ptracers) is done.
36    
37    o pkg/diagnostics: prepare for Global/Regional statistics.
38      add 2nd namelist in data.diagnostics and update DIAGNOSTICS_SIZE.h
39    o remove "baseTime" (no used) from arg. list of DIFF_BASE_MULTIPLE
40      and rename it: DIFFERENT_MULTIPLE
41    
42    checkpoint57h_post
43    o move tIter0 from data.gchem to PTRACERS_Iter0 in data.ptracers
44    o pkg/diagnostics: change output frequency, now in seconds (instead
45      of Nb of iterations). This requires to change data.diagnostics.
46    checkpoint57h_pre
47    
48    o in mnc, let the user create the output dir if both mnc_outdir_date
49      and mnc_outdir_num are false
50    o add fflags() to pkg/diagnostics which allows one to specify, on a
51      per-file basis, the file precision [default behavior is unchanged]
52    o added diagnostics to pkg/kpp, including computation of mixed layer
53      depth based on a temperature/density criterion
54    o updated verification/natl_box to test the new pkg/kpp diagnostics
55    o bugs fixed (with useAbsVorticity=T & with useJamartMomAdv=T)
56    o test AB-3 code with 4th.Order advection scheme (advect_xy/input.ab3_c4)
57    o fix mnc flags in pkg/seaice so that they are "owned" by seaice,
58      not mnc -- this mistake was found by JMC and is now corrected
59    
60    checkpoint57g_post
61    o add Adams-Bashforth 3rd Order time stepping:
62      - default remains AB-2 (#undef ALLOW_ADAMSBASHFORTH_3)
63      - not yet implemented with NH, CD_Scheme, pTracers, MNC_pickup ...
64    checkpoint57g_pre
65    
66    o make lambda[Theta,Salt]ClimRelax spatially varying
67      (default remains spatially uniform constant based on tau)
68    o added new control variables
69      * init. uVel, vVel, etanN
70      * lambda[Theta,Salt]ClimRelax
71    o separate masks used for ctrl_pack/unpack 'from write_grid' output
72      (suggested by G. Forget)
73    
74    checkpoint57f_post
75    o add parameter "baseTime" = model base time = time at iteration zero
76      new function DIFF_BASE_MULTIPLE (replace DIFFERENT_MULTIPLE) to
77      account for baseTime shift
78    o move Coriolis arrays from PARAMS.h to GRID.h
79    checkpoint57f_pre
80    
81    o mnc-ifying the seaice package
82    o introduced an alternative length scale calculation for use in the limiting
83      maximum viscosity (off by default, see mom_vi_hdissip.F for details)
84    o S/R ptracers_calc_diff: argument list bug for CALL KPP_CALC_DIFF_S
85    o added modified Leith scheme with viscosity proportional to grad(Hdiv)
86      as well as grad(vort3).  Made it compatible with both approximate form
87      for the cubed sphere and full form of these gradients.  Added flag
88      useFullLeith to select between approximate and full form.  
89    o add "d0" in getcon.F to reduce differences between compilers/platforms
90      update output.txt of exp. fizhi-gridalt-hs
91    o change port_rand funtion to allow to set the seed.
92    o fix obscure MNC bug involving conversion between reals and integers
93      - it caused no problems with certain (eg. gcc, ifc) compilers but the
94        new version is now more standards-compliant
95    o adjoint-related modifications
96      * gad_dst3_adv_?: to be able to use advscheme=30 in adjoint
97      * gad_calc_rhs: fool adjoint, revert to advscheme=30 in reverse when using 33
98      * mom_vi_hdissip: to be able to use Leith scheme in adjoint
99    
100    checkpoint57e_post
101    o Adding eddy stress controls a la Ferreira et al.
102    checkpoint57e_pre (actually was tagged: eckpoint57e_pre)
103    
104    o fizhi pkg: fix AtoC and CtoA when called with same input & output array
105      update output.txt of exp. fizhi-gridalt-hs
106    o add more "Z" dimensions and associated coordinate vars to MNC
107    o pkg/exf: added twoDigitYear capability for option useExfYearlyFields
108    o parameter calendarDumps: when set, approximate months (30-31 days) and years
109      (360-372 days) for parameters chkPtFreq, pChkPtFreq, taveFreq, SEAICE_taveFreq,
110      KPP_taveFreq, and freq in pkg/diagnostics are converted to exact calendar
111      months and years.  Requires pkg/cal.
112    o diagnostics: add pickup functionality
113      - off by default for backwards compatibility
114      - current version only uses MDSIO -- MNC will be added
115    o allow to use pkg diagnostics for LAND pkg output & Shapiro Tendencies.
116    o pkg/shap_filt: fix for multi-tiles per face (CS-grid with EXCH2)
117    
118    checkpoint57d_post
119    o pkg/thsice used in the ocean: fix tile Pb related to sIceLoad forcing
120    o modify verification/ideal_2D_oce so that it tests the creation of mnc
121      output from the diagnostics package
122    o temporary fix in EXCH2_UV: add a 3rd call to fix the 24.tiles CS-grid set-up
123    o make genmake2 more portable by using ".$FS" instead of ".f" for more
124      of the test compilations (eg. IBM xlf on Mac OSX w/HFS filesys)
125    o add coordinate variables to the MNC output from the diagnostics package
126    o change MNC-flags default value: true when useMNC=TRUE
127    o add a "readgrid_mnc" flag to MNC and give model/src/ini_curvilinear_grid.F
128      the ability (off by default) to read the grid information using MNC
129    o lab_sea: use LSR_ERROR=1.E-12 for testing ; update output.txt
130    o diagnostics: fix a bug (sub-set of levels); keep full precision when averaging
131    o do snap-shot diagnostics at the middle of the time-interval (instead of
132      at the end).
133    o update exp. aim.5l_cs (use default LAND_OPTIONS.h file) and add
134      new test for AIM with thSIce + slab-ocean
135    o cpl_aim+ocn test: update output (optfile & default GMREDI_OPTIONS.h)
136    o move state variable diagnostics to the beginning of the time step ;
137      (note: still needs some work specially for snap-shot diagnostics)
138      duplicate & replace UVELVVEL diagnostics with UV_VEL_C & UV_VEL_Z
139    
140    checkpoint57c_post
141    o Change testing option (-ieee) to get closer agreement between g77 & ifort
142    o generate new output.txt files on faulks with g77 & the new -ieee option
143    checkpoint57c_pre
144    
145    o add "_d 0" in freeze ; change some "FLOAT" to "DFLOAT" and went back (TAF);
146      update output.txt of exp. global_ocean.90x40x15
147    o fix in ini_eos.F (one "_d 0" was missing): now g77 & ifc give much closer
148      results; update global_ocean_pressure & global_ocean.cs32x15 output.txt
149    o add saltFlux(+File) forcing capability with periodicExternalForcing
150    o allow to use pkg diagnostics for GM-Redi & ThSIce output
151    o diagnostics_fill: fill in up to sNx+1 or sNy+1 (if U,V or Z pt) for mnc.
152    o improvements for the MNC output from DIAGNOSTICS:
153      - correct handling of the unlimitied dimension
154      - all dims now have associated coordinate variables
155      - variables located at U, V, and Corner ("Z") points are now written
156        with the correct horzontal dimensions (eg. 'Xp1', 'Yp1') so that
157        *all* data points are output
158    o added seaice_summary.F and removed obsolete ALLOW_SEAICE's from pkg/seaice
159    o added pkg/seaice/seaice_monitor.F
160    o mdsio_readfield.F can now read files visible only to master MPI process.
161    o pkg/seaice: some test code for ice-modified ocean stress
162    o exf_getffields interpolation defaults to bilinear for all scalar forcing
163      fields (remains bicubic for wind velocity and stress).  This avoids, e.g.,
164      spurious negative numbers for precipitation and humidity.  Will cause
165      some small numerical differences for integrations using pkg/exf/exf_interp.F.
166    o pkg/ctrl/
167      * Add ctrlvec diagnostics in pack/unpack for nondimensional I/O
168        May be enabled via doPackDiag
169      *  introduce z-scaling of
170         + gradient (doZscalePack) and
171         + control (doZscaleUnpack)
172    o pkg/mdsio/
173      * enable useSingleCpuIO for state estimation (ctrl_pack/unpack) mdsio_gl.F
174      * N.B.: leave tiled I/O for tapelev (i.e. mds...vector)
175    
176    checkpoint57b_post
177    o simplify options in diagnostics_fill arguments (remove weird combinations)
178    o fix MNC coordinate variables so that they now contain per-face index
179      values (instead of xC,xG,yC,yG values) when run with EXCH2 since, in
180      general, EXCH2 grids are only aligned to a "local" coordinate system
181    o test diagnostics pkg in dic_example (no output) and ideal_2D_oce (with ouput)
182    o Implement instantaneous output within pkg diagnostics:
183      setting a negative frequency in data.diagnostics produces
184     "snap-shot" output (as opposed to time-average if freq > 0)
185    o move out off EQUIVALENCE all diagnostics that fizhi does not use ;
186      NOTE: change ptracer diagnostics name: use 2 characters to work with
187            more than 9 tracers. Add also diagnostics of advect & diffus flux.
188    o more CF compliance: convert all MNC 'iter' variables to 'T' so that it
189      becomes the coordinate variable for the 'T' dimension
190    o add CF-style coordinate variables to MNC
191      - just a first cut:  numbers are meaningless on the cubesphere and
192        missing along the T axis but otherwise it works!
193    o add surface flux diagnostics: TAUX, TAUY, TFLUX, TRELAX, TICE, SFLUX, SRELAX
194    o pkg/thSIce: read initial state from 6 files ; allow monitor to use MNC.
195    o add diagnostics of advective & diffusive flux for T & S (pTracers coming soon)
196      this is based directly on what the model compute, so that the content of
197      each diagnostics might be affected by switching options (e.g., implicitDiffusion)
198    o small changes for diagnostics (gad_biharm_?.F, arguments of S/R impldiff.F)
199    o clean up build options files for SunOS, in particular
200      - remove -traditional from xmakedepend
201      - change -xarch=v9 to -xarch=native for f77
202      - restore use of make (not gmake)
203      - add CC=gcc if using g77
204    o changes per discussion of AM, JMC, & EH3 today:
205      - move redundant "grid" variables out of "state" -- in preparation of
206        creating proper CF-style coordinate variables and in recognition
207        that any/all grid variables can be easily copied (eg. w/NCO) from
208        the grid files into any of the other NetCDF files
209      - *always* output grid variables with MDSIO (this will eventually be
210        changed)
211      - output MNC grid file if useMNC is true
212    o every instance of _END_MASTER() has been found and replaced with
213      _END_MASTER( myThid ) in order to satisfy certain picky Sun
214      preprocessors
215    o re-arrange diagnostics pkg initialisation:
216     - allow each package to extend the available diagnostics list
217     - add some checking and fix small problems (multi-threaded, ...)
218     - new S/R diagnostics_fill (replace fill_diagnostics):
219       look through the short list of active diag. (instead of the long list)
220       create function DIAGNOSTICS_IS_ON to tell if a diagnostics is active
221    
222    checkpoint57a_post
223    o this set of changes restores TAMC compatibility
224      tested for global_ocean.90x40x15 adjoint on trough
225      (what a pain)
226    
227    checkpoint57a_pre
228    o some modifs in tools/build_options for sunos
229    o call ctrl_pack for ctrl at end for optimcycle=0
230    
231    checkpoint57
232    o add horizontal bi-harmonic viscosity for vertical velocity
233      - new parameter viscA4W (defaults to viscA4)
234      - affects only non-hydrostatic code
235      - to be done: no-slip lateral boundary conditions
236    o switch back to "sflux = sflux - runoff" in exf_getforcing.F
237    o implement Implicit Vertical advection for pTracers
238    o depth convergence accelerator: replace deltaTtracer by dTtracerLev(k)
239      (implemented also for implicit vertical diffusion & advection)
240    o OBCS: extend application of tracer open boundary conditions
241      to a strip beyond the boundary of width Olx/Oly
242    o add call monitor to TAF set of required calls (flow directives)
243      and include in natl_box_adjoint
244    o untangle a few ifdefs related to cost, ctrl, autodiff
245      untangel a few ALLOW_ECCO_OPTIMIZATION
246    o Apply Cooper & Haines representative estimates also at sub-surface (k=Nr)
247      (theta2, salt2 in ecco_cost_weights)
248    
249    checkpoint56c_post
250    o OBCS as control variables
251      - update ctrl_ad_diff.list, obcs_ad_diff.list
252      - remove balance of obcs controls from default
253      - fix index bug nobcs in ctrl_init
254      - fix dummy fields filen in ctrl_pack
255      - add dummy weights for obcs
256    o lsopt: comment out fort.94 output
257    o optim:
258      - remove unnecessary header files
259      - adjusted namelists
260      - add xerbla.F to Makefile
261    o exf:
262      - multiple modifications of exf_check_range
263      - unit change in exf_init_runoff from [m/year] to [m/s]
264    o mdsio:
265      - i/o-vector change debugLevel from debLevA to debLevB
266    
267    checkpoint56b_post
268    o collect all ptracers output in an new routine (PTRACERS_OUTPUT) that
269      is called from DO_THE_MODEL_IO
270      - rename PTRACERS_STATV_WRITE to PTRACERS_WRITE_TIMEAVE
271      - introduce new parameter PTRACERS_dumpFreq (defaults to dumpFreq)
272      - clean up PTRACERS_WRITE_TIMEAVE a little
273    o add another verification experiment that tests both the cfc-pkg and
274      the gchem-pkg with GCHEM_SEPARATE_FORCING undefined, that is with
275      reactive tendencies treated inline with regular timestepping
276    o testreport previously did a lousy job of comparing small numbers, so
277      we (Ed and Martin) changed the  way the cut-off of 1.e-12 is
278      implemented in the little c-progam tmp_cmpnum.c which is embedded in
279      testreport.
280    o ptracers_monitor: bug fix for MPI
281    
282    checkpoint56a_post
283    o GCHEM: finish reorganizating the package
284      - forward_step calls GCHEM_CALC_TENDENDY, which computes gchemTendency
285        (introduces another 3D-array for each passive tracer, but only if
286        GCHEM_SEPARATE_FORCING is undefined. For GCHEM_SEPARATE_FORCING
287        gchemTendency is not needed because the timestep is done separately)
288      - gchemTendency is added to gPtr in GCHEM_ADD_TENDENCY
289      - GCHEM is now---more or less cleanly---separate from PTRACERS
290      - to bo done: handling of surfaceForcingPtr that is needed for KPP not
291        so clear to me how that should be done in a general way.
292    o CFC: improve formatting of output in CFC_ATMOS
293    o PTRACERS:
294      - rename GCHEM_MONITOR to PTRACERS_MONITOR and call it from MONITOR, so
295        that ALL experiments with ptracers enable can be checked. This makes
296        GCHEM_MONITOR obsolete.
297      - include a runtime parameter PTRACERS_monitorFreq that defaults to
298        monitorFreq
299      - set default PTRACERS_write_mdsio to false if PTRACERS_write_mnc is true
300    o add PTRACERS (PTR0[1-5]) to testreport output
301    o exf_check_range:
302      * default changed to .TRUE.
303      * will check values for niter0 if debLevA, every timestep if debLevB
304      * STOP in exf_check_range if 'out of range'
305      * range chosen generous to allow for potential extema
306      TBD: a version with stricter range
307           permitted range would be runtime parameters
308    o GMredi:
309      - fix ldd97 slope limit ; extend valid domain of Psi-Bolus ;
310      - change S/R GRAD_SIGMA to be valid in the overlap on CS-grid.
311      - change the default GMREDI_OPTIONS.h: compile everything except Visbeck
312    o GCHEM:
313      - move call to GCHEM_FORCING_SEP, so that it is now called before
314        the blocking exchanges to avoid computing halos in GCHEM_FORCING_SEP
315      - introduced a useGCHEM runtime flag and applied it everywhere.
316      - renamed PTRACERS_SEPARATE_FORCING to GCHEM_SEPARATE_FORCING;
317        this flag is now only visible within GCHEM-pkg routines.
318      - moved gchem_forcing_int temporarily into ptracers_forcing, this will
319        be replaced by a 3D array that is computed before thermodynamics; for
320        this, the tracer ID iTracer had to be included into the parameter list
321        of ptracers_forcing
322      - set up forward_step to call gchem_forcing, that will replace
323        gchem_forcing_int in ptracers_forcing
324      - cleaning up
325    o small cleaning of exch2: remove files that are also in eesupp/src
326    
327    checkpoint56
328    o AIM: read land-Fraction earlier (from ini-fixed) and monitor land at nIter0
329    o mnc: cleanup "full" files--fixes out-of-memory problem
330    o ECCO-GODAE
331      * model/, pkg/, code/, input/ now complete w.r.t. ECCO-GODAE production setup
332      * fully merged and updated from SIO code
333        windstress sign change affects
334        - exf_inscal_[u/v]stress
335        - sign of xx_tau[u/v] adjustments
336        - signs in cost_scat
337      * ctrl_get_gen: avoid update of tau[u/v] for first few days
338      * natl_box_adjoint updated and compressed
339      * 2-level checkpoint switch added
340        (tested in EXP. 1 of natl_box_adjoint)
341      * new switch doInitXX to avoid init. of xx_... in optimcycle=0
342    o SVD
343      * first complete, non-hacked version for SVD calculations
344      * all interfaces clean w.r.t to handling of xx_..., adxx_..., g_xx_...
345        performed via direct access files admtlm_vector.it????
346        (via ctrl_[un]pack)
347      * optimcycle is used for ARPACK iteration counter
348        (note: adxx_... is iterated, rather than xx_...; latter is constant!)
349      * EV problem ok w.r.t. L2-norm in tangent space; need to check other norms
350        and add linear transforms where necessary
351      * controls are initial/final SST, SSS (can be extended; need to check
352        ctrl_set_[un]pack_xyz.F according to ctrl_set_[un]pack_xy.F
353    o mdsio: filename info to both stderr, stdout
354    o exf: exf_inscal_runoff: separated default conversion [m/year] -> [m/sec]
355      from sign change
356    
357    checkpoint55j_post
358    o isolate dissipation tendency (allow to keep it out off AB) in V.I form.
359      affects (truncation level) output.txt of global_ocean.cs32x15(&.thsice)
360    o add consistent build directories and .cvsignore files to verification
361    o added viscA4GridMin and viscA4GridMax for more control on biharmonic viscosity
362    o make all subroutines 30 or less characters in length for certain older
363      compilers to work (a thorough search only turns up two routines in
364      pkg/ptracers)
365    o allow (vecinv only) a different horiz. viscosity for Divergence and Vorticity;
366      affects (truncation level) output.txt of global_ocean.cs32x15(&.thsice)
367    
368    checkpoint55i_post
369    o change arrays for total vertical diffusivity (2.D if explicit,
370      3.D if implicit) ; each ptracer uses its own background value.
371    o added mnc/NetCDF output capability to ocn_compon_interf
372      - added bits to verification/cpl_aim+ocn so that it can run on
373        ITRDA and other clusters without the mpirun "-wd" option
374    o fixed wayward endif in mom_fluxform.F that messed up cyclindrical
375      configuration
376    o give MNC the ability to create additional files rather than
377      exceeding a specified file size limit "MNC_MAX_FILESIZE"
378      - the default is just less than 2GB which is a practical limit
379        on many systems
380    o Compute the non-local term using either KPP vertical diffusivity (=default)
381       or the total vertical diffusivity (when KPP_ghatUseTotalDiffus=T)
382    
383    checkpoint55h_post
384    o add verification test "vermix"
385    o remove all the *_ioinc flags and replace with the single global
386      outputTypesInclusive flag
387    o allow to set a vertical profile of vertical diffusivity for T & S
388    o call PACKAGE_BOOT & PACKAGES_READPARMS just after INI_PARMS
389    o useMNC moved from data.mnc to data.pkg
390    
391    checkpoint55g_post
392    o Renaming carbon AD experiments to comply with testreport conventions
393    o Updating adjoint of Held-Suarez hs94.1x64x5
394    o Fixing directives for non-ECCO forward_step, after ECCO directives
395      have been fixed to avoid recomp. of do_oceanic_phys
396      Works well for global_ocean.90x40x15,
397      but recomp. issue remains for carbon
398    o Updating adjoint output for natl_box_adjoint after
399      pkg/ecco has been modified (see taga c55d_post)
400    
401    checkpoint55f_post
402    o rename to timeave_init_fixed
403    
404    checkpoint55e_post
405    o some delicate re-shuffle of store directives to avoid one
406      extra call of do_oceanic_physics
407      affects ecco_the_main_loop.F, do_oceanic_physics.F
408      checkpoint_lev1_directives.h
409    o NB: this may break global_ocean adjoint temporarily,
410      but it is clear how to fix it. Will do later, need this now.
411    
412    checkpoint55d_post
413    o ECCO specific cost function terms (up-to-date with 1x1 runs)
414    o ecco_cost_weights is modified to 1x1 runs
415    o modifs to allow observations to be read in as
416      single file or yearly files
417    o enable to read exf forcing fields as either
418      single file or yearly files (flag useExfYearlyFields)
419    checkpoint55d_pre
420    
421    o switch between 2-level and 3-level checkpointing
422    o add ability of MNC to write local and "partial" (eg. 2D slices where
423      the full 3D field is never actually stored) arrays to NetCDF files
424      with the correct (that is, the complete multi-dimensional) set of
425      array indicies
426      - used in mom_vecinv() to write the diagFreq output
427      - tested (demonstrated) in verification/aim.5l_cs
428    o further mnc-ification
429      - add verification/lab_sea to the list with mnc setups
430      - fix undefined mnc v-names in cd_code
431      - add kpp to mnc output (both snapshot and timeave)
432    o non-hydrostatic code
433      - added horizontal viscosity parameter viscAhW for w, that replaces viscAh
434        in calc_gw.F
435      - defaults to viscAh
436      - this is useful, when I want to use viscAh=0 (e.g., only biharmonic
437        friction or Leith) in non-hydrostatic mode, as long as there is
438        no biharmonic friction and/or Leigth or equivalent for vertical velocity.
439    o pkg/cost
440      extend simple cost functions to treat time mean quantities
441      (T,U,V) via cost_accumulate_state from cost_tile
442    o multidimAdvection on CS-grid: implement the full 3 pass ; overlap needs to
443        be 2 times larger than the minimum size the advection scheme requires.
444        now get a perfect conservation of salt in advect_cs test-exp (updated)
445    o GGL90 (Gaspar's TKE mixed layer)
446      - add horizontal diffusion of TKE, requires exchanges in
447        do_fields_blocking_exchanges, horizontal diffusivity is zero by
448        default. In OPA there is no horizontal diffusion of TKE but the
449        mixing coefficients are computed from a horizontal average of TKE of
450        6 points or so. I think that diffusion has a little more physical
451        justification.
452      - clean up ggl90_calc in the hope of reducing memory usage (this hope
453        was in vain)
454      - mask tke-variable in ggl90_init
455    o fix mask Pb in multidimAdvection on CS-grid (still a small Pb with ipass)
456    o fix biharmonic diffusion on CS-grid.
457    o more obcs stuff
458      - added obcs_external_fields_load as an alternative to the EXF-pkg option
459        for reading boundary values. This routine is called from
460        obcs_prescribe_read
461      - updated obcs_check
462      - changed sponge layer relaxation parameter units for days to seconds
463      - turned off sponge layer by default
464      - changed argument list of OBCS_CALC to myTime+deltaTclock
465      - fixed about 35 bugs in obcs_external_fields_load
466    
467    checkpoint55c_post
468    o obcs stuff
469      - new runtime flags useOBCSprescribe, useOBCSbalance
470      - re-arranged obcs_calc
471      - moved exf_getobcs to obcs_prescibe_read
472        obcs_prescibe_read is called after ORLANSKI stuff
473        and independent of it (executable can contain all).
474      - updated obcs_check
475      - added exchanges for U,V,gU,gV for useOBCSprescribe=.TRUE.
476        These may be necessary only if tangential velocities are
477        prescribed at OB, not sure, needs testing.
478    o make the MNC_CW_ADD_VATTR_* subroutines all take scalar arguments
479      which simplifies the calls and passes the ifc check-bounds test
480    
481    checkpoint55b_post
482    o add a 2D scalar, Z location, exchange_2 routine.
483    o EXCH_UV for CS-grid: add one u,v in the corner-halo region, allows
484      to compute vorticity on a wider stencil [e.g., vort3(0,1)&(1,0)]
485      affects advect_cs results (for the wrong reason): update output.txt
486    o use a more standard discretization for biharmonic viscosity ;
487     (original version still available with #define MOM_VI_ORIGINAL_VISCA4 )
488      affects results of dome set-up: update output.txt
489    o fix gad_advection.F (CubedSphere part), back to version 1.23
490    o fix bug in mom_vi_del2uv.F (CubedSphere part).
491    
492    checkpoint55a_post
493    o merged code to
494      * prescribe/read time-dependent open boundaries
495        (works in conjunction with exf, cal)
496      * sponge layer code for open boundaries
497      * each boundary N/S/E/W now has its own CPP option
498        (healthy for the adjoint)
499    
500    checkpoint55
501    o up-to-date adjoint of previous breaks
502      - natl_box_adjoint: problem for exp 12, 13: cost = 0. needs fix
503      - global_ocean_ebm: adjoint still needs fix (TLM ok)
504      - global_ocean.90x40x15: extended GMREDI_OPTIONS.h
505      - carbon: * remove tr1-related adjoint verification
506                * add ptracer as control variable
507    o remove all tr1-related code (ALLOW_PASSIVE_TRACER)
508      (adjoint stuff still has some tr1 'names', but all use ptracer arrays)
509    
510    checkpoint54f_post
511    o in verification/carbin/ remove experiment related to
512      ALLOW_PASSIVE_TRACER (array tr1 etc)
513      code_ad/ input_ad/ results_ad
514    
515    o add two new packages
516      - ggl90, TKE-model of Gaspar et al. (1990), Blanke+Delecluse (1993)
517        + subject to optimization (memory)
518      - opps, OPPS (ocean penetrative plume scheme) of Paluszkiewicz+Romea (1997)
519        + subject to reordering (many GOTOs in the original code)
520        + so far opps is called for every (i,j) point, and then operates in
521          the vertical; one may want to reorder the scheme in order operate in
522          horizontal slabs at the cost of (a lot of) memory, at least ten new 3D
523          fields)
524    o overhaul of IO so that we now have flags for MDSIO and/or MNC
525      - all verification tests compile and run with linux_ia32_g77
526      - defaults are compatible with current input files--nothing
527        should change if you were not previously using MNC
528      - MNC output has been added in numerous places (eg. timeave)
529        but there are still a few writes not yet do-able with MNC
530        (this is in progress)
531      - flags now allow for either/or/both use of MDSIO and MNC and
532        documentation will soon follow
533      - numerous small formatting cleanups for ProTeX
534    o add version info (via BUILD_INFO.h) to MNC and stdout
535    
536    checkpoint54e_post
537    o add two new packages
538      - pp81 (Packanowski and Philander, 1981), Richardson number and
539        stratification dependent mixing
540      - my82 (Mellor and Yamada, 1982) level 2 turbulence closure scheme
541      - these scheme operate in a similar manner as KPP (without ghat, of
542        course), however there is only one diffusivity that is used for both
543        temperature and salinity (and passive tracers)ç
544    o add offline package
545    o add "verification/dic_example" which contains a dic-gchem-ptracer setup
546      provided by Stephanie Dutkiewicz
547    o genmake2: add ability to search for C compilers and cleanup test output
548      - fixes recent SX-6 problems and should help with others (eg. HP-UX)
549    
550    checkpoint54d_post
551    o data.diagnostics: if levels() is unspecified then ALL levels are diagnosed
552    o bug fix in pkg/cal/cal_addtime.F for TheCalendar='model' day count is wrong,
553      leading to additional shift by 5 days each year of integration
554       (i.e. leading to considerable offset over the years)
555    o update global_ocean.cs32x15/thsice: use z* + RealFW + staggerTimeStep
556    o start with no AB and no initial filtering (staggerTimeStep):
557      change the results of test-exp: hs94.128x64x5, hs94.1x64x5 & hs94.cs-32x32x5
558    o move ocean specific S/R calls from forward_step.F to do_oceanic_phys.F
559      (external loading & cpl import/export stay in forward_step.F)
560    o update coupling interface to work with new staggerTimeStep
561    o change MDSIO ptracers_pickup.
562    o fix a bug in Flux deriv. vs Ts (pkg/aim_v23/phy_suflux_sice.F);
563    o genmake2 fixes so netcdf now works on halem ("osf1_alpha_*")
564    o bugs fixed in kpp_calc.F (ustarY) and external_forcing.F (Qsw & Partial-cell)
565    
566    checkpoint54c_post
567    o replace surfaceTendency U,V,S,T,Tice,Ptr by surfaceForcing U,V,S,T,Tice,Ptr
568      this affects (truncation error) output.txt of test-exp:
569       global_ocean_ebm,_pressure,.cs32x15(&.thsice), lab_sea & natl_box
570    o pkg/kpp: replace pMask by maskC to reduce memory size (one 3.D array)
571    o change Qnet to always be the net heat flux, (+upward).
572      affects (truncation error) output of natl_box exp.: updated
573    o define the max number of passive tracers in PTRACERS_SIZE.h
574    o start to test mnc routinely in exp0
575    o add sea-ice loading (sIceLoad) as a new forcing field.
576    
577    checkpoint54b_post
578    o add testreport "-j JOBS" option for parallel builds (tested and pretty
579       impressive on Altix and Opteron systems)
580    o diagnostics can now write NetCDF files through the mnc package
581    o staggerTimeStep & NLFS: update pkg aim_v23 & thsice; update aim.5l_cs output
582    o staggerTimeStep & multi-Dim-Advec: Exch(U,V) added
583    o staggerTimeStep & NLFS: rescale T & S surface tendencies.
584      affects (truncation error) output of lab_sea & natl_box exp.: updated
585    o change pickup _cd & _ptracers (old pickup restart: usePickupBeforeC54=.T.)
586    
587    checkpoint54a_post
588    o re-write staggerTimeStep: step fwd U,V and then T,S (reverse the order).
589    o change pickup file (to restart from old pickup, set usePickupBeforeC54=.T.)
590    o get a perfect restart when using RealFreshWaterFlux + NLFS (not finished)
591    o update output of all exp. that use staggerTimeStep (except hs94.cs & 2D)
592      + internal_wave (shift in writing cg2d_ini_res)
593    checkpoint54a_pre
594    
595    o require forcing_In_AB=.FALSE. and remove surfaceTendencyTice from pkg/seaice
596      - modified verification/lab_sea and natl_box_adjoint accordingly
597    
598    checkpoint54
599    checkpoint53g_post
600    o prepare splitting of thermodynamics: store convect.counter in common block
601    o store d.etaH/dt (instead of Div.hV) in common ; affects (truncation error)
602      output of: global_ocean.cs32x15, global_ocean_pressure & internal_wave
603    o cnh's modifs to gad_advection.F for cube-sphere multi-dim advection
604    o new parameter "nShapS": allow to treat differently S from T (nShapT)
605    o separate Vert.Advec.Scheme from horizontal Advec.Scheme (T & S)
606    o restore adjointability
607    
608    checkpoint53f_post
609    o modifications that affect AIM experiment results (update output.txt):
610      - no advection of S at top interface Nr (pkg/generic_advdiff)
611      - air surf. temp: use lower value between T0 & T1 (phy_suflux_prep.F)
612    o AIM pkg:
613      - include stability function into surf.Flux derivative relative to Tsurf
614      - calculate clear-sky radiation & surface temp. change (for diagnostics)
615    
616    checkpoint53e_post
617    o merged cylindrical coord configuration
618    o added rotating_tank experiment
619    checkpoint53e_pre
620    
621    checkpoint53d_post
622    o faulks upgraded: re-generate output.txt that have changed:
623       global_ocean.cs32x15(.thsice), global_ocean_ebm, lab_sea, natl_box
624    o Added Bryan and Lewis, 1979, vertical diffusivity profile
625    o pkg/land: fix problem in land_impl_temp.F when Ts arround freezing.
626    o added Sadourny discretization of Coriolis
627      - moved some PARAMETERS from mom_*_coriolis.F to PARAMS.h
628      - re-enabled use of omega3 (useAbsVorticity)
629    o added Leith (Phys. Fluids, 1968, 10, 1409-1416) variable viscosity for
630      vector-invariant equations
631      - note: does *NOT* yet use same variable visc. for no-slip BC's.
632    o recoded the lateral boundary conditions for non-hydrostatic w, this time
633      hopefully without new bugs
634      - still to do: biharmonic friction for non-hydrostatic w
635    o bug fix for mom_*_sidedrag: rA -> rAs,rAw
636    o regenerated output for global_ocean_ebm (using -ieee)
637    o fix mom_common/mom_*_sidedrag.F to recover numeric output prior to creation
638      of mom_common/
639      - simple re-arrangement of instructions within computing loop
640    checkpoint53d_pre
641    
642    checkpoint53c_post
643    o new version of coupler for AIM+thSIce and Ocean models
644     - new fields (SSS,h_MixLayer ...) are exchanged through the coupler
645     - land: change Evaporation of snow
646     - aim with thSIce: fix bug in snow-precip.
647    o changes to calc_gw.F
648     - bug fix for w^2 term near the bottom boundary
649     - (hopefully) improve the lateral slip boundary condtions for use with
650        partial/looped cells
651      - because the bug fix changes two verification experiments anyway (exp5
652        and plume on slope), change the lateral boundary condition from half slip
653        to the value of no_slip_sides
654    o ebm package wasnt quite up to the parallel task. Fixed that.
655    o Adjoint of EBM package + verification
656      * TLM OK
657      * ADM not sufficient accuracy yet (needs further testing)
658    
659    checkpoint53b_post
660    o Commiting new energy balance model to repository
661      * package is pkg/ebm
662      * verif. is verification/global_ocean_ebm
663      * references are in ebm_driver.F
664    C      * X. Wang, P. Stone and J. Marotzke, 1999:
665    C        Global thermohaline circulation. Part I:
666    C        Sensitivity to atmospheric moisture transport.
667    C        J. Climate 12(1), 71-82
668    C      * X. Wang, P. Stone and J. Marotzke, 1999:
669    C        Global thermohaline circulation. Part II:
670    C        Sensitivity with interactive transport.
671    C        J. Climate 12(1), 83-91
672    C      * M. Nakamura, P. Stone and J. Marotzke, 1994:
673    C        Destabilization of the thermohaline circulation
674    C        by atmospheric eddy transports.
675    C        J. Climate 7(12), 1870-1882
676    o Will need long integration testing.
677    checkpoint53b_pre
678    
679    o created pkg/mom_common in preparation for variable viscosity code
680    o renamed and split ini_mnc_io
681    o pkg/land: only liquid water can diffuse or run-off
682    o AIM using land: soil wetness availability account for snow
683    o added printMapIncludesZeros to EEPARMS namelist
684    
685    checkpoint53a_post
686    o restore TAMC compatiility for exf_constants.h
687      (no usage of _d in .h)
688    o Additional storing in KPP in an attempt to fix
689      TAMC problem with adjoint of KPP code.
690      (no success...)
691      Identical code works fine with TAF.
692    o Fixed verification/natl_box_adjoint experiments.
693      Initial conditions for Experiments 4 and 5 are now
694      set to "SEAICE_initialHEFF=0", which is more reasonable,
695      and hence gives results closer to that of Experiment 2.
696      Experiment 6 requires "SEAICE_initialHEFF=1" and is working once again.
697    
698    checkpoint53
699    o stable checkpoint with up to date adjoint
700      TLM tests only run for 1st grdchk at the moment. Its unclear if
701      this is a code problem or a problem with resetting code for g_
702      variables that is being autogenerated wrong.
703      Otherwise:
704       o code is testing well on Linux (P4 and AMD), SGI Alitix and Irix.
705       o IBM and HPQ tests could not be run because NCAR (IBM) and NCCS (HPQ)
706         are both still in chaos following their security messes.
707    
708    checkpoint52n_post
709    o various build (mostly genmake2) changes including:
710      - ability to build/run more easily within Cygwin and MacOS X
711      - more complicated Makefiles with a stronger dependence on Gnu Make
712    o pkg/seaice modifications
713      - added SEAICE_deltaTtherm and SEAICE_deltaTdyn in order to reduce cost
714        of pkg/seaice.  For time being, SEAICE_deltaTtherm = deltaTtracer
715        and SEAICE_deltaTdyn must be integer multiple of SEAICE_deltaTtherm
716      - changed proxy for geostrophic velocity and updated
717        verification/lab_sea/results/* accordingly
718      - removed ADI dynamic solver from pkg/seaice
719      - replaced all occurrences of 1/CSTICE or 1/CSUICE
720        they caused a segmentation fault under certain circumstances
721      - removed unused field WATR
722    o bug fix in exf-package:
723      -  u/vwind are defined at tracer/pressure points. The values of u/vwindmask
724         now reflect that.
725    o small fixes in the cal-package:
726      - set reference date to first of January for the 'model'-calendar
727      - failing of cal_CheckDate in cal_FullDate now actually results in a
728        warning
729  o small modifications in exf-package to increase flexibility:  o small modifications in exf-package to increase flexibility:
730    - change climtempfreeze from a constant to a runtime parameter that    - change climtempfreeze from a constant to a runtime parameter that
731      is specified in data.exf_clim (default = -1.9)      is specified in data.exf_clim (default = -1.9)

Legend:
Removed from v.1.280  
changed lines
  Added in v.1.502

  ViewVC Help
Powered by ViewVC 1.1.22