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

Legend:
Removed from v.1.306  
changed lines
  Added in v.1.527

  ViewVC Help
Powered by ViewVC 1.1.22