/[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.158 by heimbach, Fri May 23 18:24:17 2003 UTC revision 1.526 by heimbach, Sat Jul 30 11:47:39 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    checkpoint57o_post
8    o fix MNC/MDSIO flags in PTRACERS so that (1) it works, and
9      (2) its ready to be documented during the 20050801 DocFest
10    o MNC-ify the gmredi package -- tested and works w/ and w/o mnc
11      for ideal_2D_oce and lab_sea on linux_ia32_g77
12    o set to #undef CTRL_SET_OLD_MAXCVARS_30
13      in global1x1_tot/code_seaice_adj to access indices 32, 33
14      (worth another tag)
15    
16    checkpoint57n_post
17    o adding precip control xx_precip 32
18    o adding swflux control xx_swflux 33
19    o Replacing 10 forcing penalty cost routines by 1 generic.
20    o Implementing partitioning in time-mean and variable adjustment
21      for atmos. state.
22    o Time-mean constant error is now read in data.ecco (instead of hard-coded)
23    o in natl_box_adjoint:
24      * updated exp 1 to 10
25        (but does not have time-mean vs. variable split results yet)
26    o updated global1x1/*_seaice_adj/
27      to incorporate latest changes
28    o enable useSingleCpuIO for mdsio vector routines
29      (applies in particular to tapelev I/O)
30    o removed inefficient K-loop in all active_file*control*.F
31    
32    checkpoint57m_post
33    o re-arrange exp. adjustment.cs-32x32x1 (adjust_nlfs.cs-32x32x1 merge into)
34      and advect_cs (use EXCH-2 and standard cs-32 grid files).
35    o pkg/exch2: filling of face-corner halo regions is optional
36      (ifdef W2_FILL_NULL_REGIONS) and is turned off by default.
37    
38    checkpoint57l_post
39     tag before updating tutorials (pre-release-2)
40    o add cos & sin of grid-direction orientation angle (curvilinear grid only,
41       for output and interpolation of input fields onto model grid)
42    o fix radius in mom_v_metric_cylinder.F (wrong by dy/2).
43       Affects results of exp. rotating_tank ; update output.txt.
44    o pkg/diagnostics: 2 new S/R to fill a diagnostics using a scaling factor
45          and taking the square quantity and/or using a fraction-weight field.
46    o fix aim_diagnostics_init.F mate numbers.
47    o add mnc output capability to diagnostics/diagstat and update
48      our cvsignore files for the ACSII output generated
49    
50    checkpoint57k_post
51    o lab_sea now uses diagnostics package and provides a suite of
52      new EXF diagnostics output.
53    
54    checkpoint57j_post
55    o exf
56      * bug fix in exf_bulkformulae for psimh (spotted by Eliyahu Biton):
57        Bracket missing in eqn. for psimh (eqn. 7 of Large & Pond, 1982);
58        terms -2*ATAN(x) + pi/2 were added to both stable and unstable case.
59      * extended diagnostics output
60        - now includes all EXF variables
61        - hs, hl moved to common block
62        - separate S/R for diagnostics, called after ctrl vector is added
63        - name changed for EXFsw -> EXFswnet, EXFlw -> EXFlwnet
64        - EXFqnet contains true Qnet (i.e. includes qsw)
65      * updated verif. lab_sea, global_with_exf
66      * Rearranging hflux (exf_mapfields -> exf_getforcing)
67        to ensure that diagnostics sees a true EXFqnet
68    o add matlab scripts for CS-grid vectors (in utils/cs_grid)
69    o diagnostics pkg: change pointers so that 1 diagnostic can be used
70      several times with different frequencies.
71    o mnc-ify the thsice package as requested by Daniel Enderton
72      - the monitor--MNC output needs work
73      - many attributes need to be added (most are currently blank)
74    o mnc-ify theaim_v23 package as requested by Daniel Enderton
75      - many attributes need to be added (most are currently blank)
76    o remove "groundAtK1" & "rkFac" (replaced by usingPCoords & -rkSign)
77    o add Sadourny like discretization to the 4th.O vorticity advection scheme
78    o fix bug in rStar code (calc_r_star.F) when used with topo/bathy.
79      update output.txt from: aim.5l_cs (& .thSI), cpl_aim+ocn,
80         fizhi-gridalt-hs and global_ocean.cs32x15 (& .thsice).
81    o fix few little things to avoid floating error with -ftrapuv (ifort v9)
82     (ini_spherical_polar_grid.F, calc_r_star.F, monitor/mon_surfcor.F)
83    
84    checkpoint57i_post
85    o Diagnostics package was changed: diagnostic allocation is now deactivated
86      by setting an empty file name in the namelist (used to be by setting 0 frequency)
87    o momentum advection (VI form) options: activate 4th.O vorticity Coriolis ;
88      add 1rst Order upwind vertical Shear.
89    o add a quick fix for A/B checkpointing in OFFLINE mode
90    o add IGNORE_TIME define and -ignore_time option to genmake2 as an
91      end-run around the C/Fortran linking horror in eesup/src/timers.F
92    o added pkg/runclock: allows user to terminate run based on wall clock
93    o added non-hydrostatic scaling parameter: nh_am2
94    o pkg/my82: futile attempt to make results compiler-independent by adding
95      some "_d 0", etc. no effect.
96    o for mnc output, fill the 'T' coordinate var with myTime and create a
97      separate 'iter' variable for iteration count
98    o pkg/diagnostics:
99      add new capability: output of Global/Regional & per-level statistics ;
100      tested against monitor output (including land & thSIce) on various exp.
101    o make seaice-pkg I/O type flags (MNC/MDS) more standard.
102    
103    checkpoint57h_done
104      tag after changing user customized files (DIAGNOSTICS_SIZE.h,
105      data.diagnostics, data.gchem, data.ptracers) is done.
106    
107    o pkg/diagnostics: prepare for Global/Regional statistics.
108      add 2nd namelist in data.diagnostics and update DIAGNOSTICS_SIZE.h
109    o remove "baseTime" (no used) from arg. list of DIFF_BASE_MULTIPLE
110      and rename it: DIFFERENT_MULTIPLE
111    
112    checkpoint57h_post
113    o move tIter0 from data.gchem to PTRACERS_Iter0 in data.ptracers
114    o pkg/diagnostics: change output frequency, now in seconds (instead
115      of Nb of iterations). This requires to change data.diagnostics.
116    checkpoint57h_pre
117    
118    o in mnc, let the user create the output dir if both mnc_outdir_date
119      and mnc_outdir_num are false
120    o add fflags() to pkg/diagnostics which allows one to specify, on a
121      per-file basis, the file precision [default behavior is unchanged]
122    o added diagnostics to pkg/kpp, including computation of mixed layer
123      depth based on a temperature/density criterion
124    o updated verification/natl_box to test the new pkg/kpp diagnostics
125    o bugs fixed (with useAbsVorticity=T & with useJamartMomAdv=T)
126    o test AB-3 code with 4th.Order advection scheme (advect_xy/input.ab3_c4)
127    o fix mnc flags in pkg/seaice so that they are "owned" by seaice,
128      not mnc -- this mistake was found by JMC and is now corrected
129    
130    checkpoint57g_post
131    o add Adams-Bashforth 3rd Order time stepping:
132      - default remains AB-2 (#undef ALLOW_ADAMSBASHFORTH_3)
133      - not yet implemented with NH, CD_Scheme, pTracers, MNC_pickup ...
134    checkpoint57g_pre
135    
136    o make lambda[Theta,Salt]ClimRelax spatially varying
137      (default remains spatially uniform constant based on tau)
138    o added new control variables
139      * init. uVel, vVel, etanN
140      * lambda[Theta,Salt]ClimRelax
141    o separate masks used for ctrl_pack/unpack 'from write_grid' output
142      (suggested by G. Forget)
143    
144    checkpoint57f_post
145    o add parameter "baseTime" = model base time = time at iteration zero
146      new function DIFF_BASE_MULTIPLE (replace DIFFERENT_MULTIPLE) to
147      account for baseTime shift
148    o move Coriolis arrays from PARAMS.h to GRID.h
149    checkpoint57f_pre
150    
151    o mnc-ifying the seaice package
152    o introduced an alternative length scale calculation for use in the limiting
153      maximum viscosity (off by default, see mom_vi_hdissip.F for details)
154    o S/R ptracers_calc_diff: argument list bug for CALL KPP_CALC_DIFF_S
155    o added modified Leith scheme with viscosity proportional to grad(Hdiv)
156      as well as grad(vort3).  Made it compatible with both approximate form
157      for the cubed sphere and full form of these gradients.  Added flag
158      useFullLeith to select between approximate and full form.  
159    o add "d0" in getcon.F to reduce differences between compilers/platforms
160      update output.txt of exp. fizhi-gridalt-hs
161    o change port_rand funtion to allow to set the seed.
162    o fix obscure MNC bug involving conversion between reals and integers
163      - it caused no problems with certain (eg. gcc, ifc) compilers but the
164        new version is now more standards-compliant
165    o adjoint-related modifications
166      * gad_dst3_adv_?: to be able to use advscheme=30 in adjoint
167      * gad_calc_rhs: fool adjoint, revert to advscheme=30 in reverse when using 33
168      * mom_vi_hdissip: to be able to use Leith scheme in adjoint
169    
170    checkpoint57e_post
171    o Adding eddy stress controls a la Ferreira et al.
172    checkpoint57e_pre (actually was tagged: eckpoint57e_pre)
173    
174    o fizhi pkg: fix AtoC and CtoA when called with same input & output array
175      update output.txt of exp. fizhi-gridalt-hs
176    o add more "Z" dimensions and associated coordinate vars to MNC
177    o pkg/exf: added twoDigitYear capability for option useExfYearlyFields
178    o parameter calendarDumps: when set, approximate months (30-31 days) and years
179      (360-372 days) for parameters chkPtFreq, pChkPtFreq, taveFreq, SEAICE_taveFreq,
180      KPP_taveFreq, and freq in pkg/diagnostics are converted to exact calendar
181      months and years.  Requires pkg/cal.
182    o diagnostics: add pickup functionality
183      - off by default for backwards compatibility
184      - current version only uses MDSIO -- MNC will be added
185    o allow to use pkg diagnostics for LAND pkg output & Shapiro Tendencies.
186    o pkg/shap_filt: fix for multi-tiles per face (CS-grid with EXCH2)
187    
188    checkpoint57d_post
189    o pkg/thsice used in the ocean: fix tile Pb related to sIceLoad forcing
190    o modify verification/ideal_2D_oce so that it tests the creation of mnc
191      output from the diagnostics package
192    o temporary fix in EXCH2_UV: add a 3rd call to fix the 24.tiles CS-grid set-up
193    o make genmake2 more portable by using ".$FS" instead of ".f" for more
194      of the test compilations (eg. IBM xlf on Mac OSX w/HFS filesys)
195    o add coordinate variables to the MNC output from the diagnostics package
196    o change MNC-flags default value: true when useMNC=TRUE
197    o add a "readgrid_mnc" flag to MNC and give model/src/ini_curvilinear_grid.F
198      the ability (off by default) to read the grid information using MNC
199    o lab_sea: use LSR_ERROR=1.E-12 for testing ; update output.txt
200    o diagnostics: fix a bug (sub-set of levels); keep full precision when averaging
201    o do snap-shot diagnostics at the middle of the time-interval (instead of
202      at the end).
203    o update exp. aim.5l_cs (use default LAND_OPTIONS.h file) and add
204      new test for AIM with thSIce + slab-ocean
205    o cpl_aim+ocn test: update output (optfile & default GMREDI_OPTIONS.h)
206    o move state variable diagnostics to the beginning of the time step ;
207      (note: still needs some work specially for snap-shot diagnostics)
208      duplicate & replace UVELVVEL diagnostics with UV_VEL_C & UV_VEL_Z
209    
210    checkpoint57c_post
211    o Change testing option (-ieee) to get closer agreement between g77 & ifort
212    o generate new output.txt files on faulks with g77 & the new -ieee option
213    checkpoint57c_pre
214    
215    o add "_d 0" in freeze ; change some "FLOAT" to "DFLOAT" and went back (TAF);
216      update output.txt of exp. global_ocean.90x40x15
217    o fix in ini_eos.F (one "_d 0" was missing): now g77 & ifc give much closer
218      results; update global_ocean_pressure & global_ocean.cs32x15 output.txt
219    o add saltFlux(+File) forcing capability with periodicExternalForcing
220    o allow to use pkg diagnostics for GM-Redi & ThSIce output
221    o diagnostics_fill: fill in up to sNx+1 or sNy+1 (if U,V or Z pt) for mnc.
222    o improvements for the MNC output from DIAGNOSTICS:
223      - correct handling of the unlimitied dimension
224      - all dims now have associated coordinate variables
225      - variables located at U, V, and Corner ("Z") points are now written
226        with the correct horzontal dimensions (eg. 'Xp1', 'Yp1') so that
227        *all* data points are output
228    o added seaice_summary.F and removed obsolete ALLOW_SEAICE's from pkg/seaice
229    o added pkg/seaice/seaice_monitor.F
230    o mdsio_readfield.F can now read files visible only to master MPI process.
231    o pkg/seaice: some test code for ice-modified ocean stress
232    o exf_getffields interpolation defaults to bilinear for all scalar forcing
233      fields (remains bicubic for wind velocity and stress).  This avoids, e.g.,
234      spurious negative numbers for precipitation and humidity.  Will cause
235      some small numerical differences for integrations using pkg/exf/exf_interp.F.
236    o pkg/ctrl/
237      * Add ctrlvec diagnostics in pack/unpack for nondimensional I/O
238        May be enabled via doPackDiag
239      *  introduce z-scaling of
240         + gradient (doZscalePack) and
241         + control (doZscaleUnpack)
242    o pkg/mdsio/
243      * enable useSingleCpuIO for state estimation (ctrl_pack/unpack) mdsio_gl.F
244      * N.B.: leave tiled I/O for tapelev (i.e. mds...vector)
245    
246    checkpoint57b_post
247    o simplify options in diagnostics_fill arguments (remove weird combinations)
248    o fix MNC coordinate variables so that they now contain per-face index
249      values (instead of xC,xG,yC,yG values) when run with EXCH2 since, in
250      general, EXCH2 grids are only aligned to a "local" coordinate system
251    o test diagnostics pkg in dic_example (no output) and ideal_2D_oce (with ouput)
252    o Implement instantaneous output within pkg diagnostics:
253      setting a negative frequency in data.diagnostics produces
254     "snap-shot" output (as opposed to time-average if freq > 0)
255    o move out off EQUIVALENCE all diagnostics that fizhi does not use ;
256      NOTE: change ptracer diagnostics name: use 2 characters to work with
257            more than 9 tracers. Add also diagnostics of advect & diffus flux.
258    o more CF compliance: convert all MNC 'iter' variables to 'T' so that it
259      becomes the coordinate variable for the 'T' dimension
260    o add CF-style coordinate variables to MNC
261      - just a first cut:  numbers are meaningless on the cubesphere and
262        missing along the T axis but otherwise it works!
263    o add surface flux diagnostics: TAUX, TAUY, TFLUX, TRELAX, TICE, SFLUX, SRELAX
264    o pkg/thSIce: read initial state from 6 files ; allow monitor to use MNC.
265    o add diagnostics of advective & diffusive flux for T & S (pTracers coming soon)
266      this is based directly on what the model compute, so that the content of
267      each diagnostics might be affected by switching options (e.g., implicitDiffusion)
268    o small changes for diagnostics (gad_biharm_?.F, arguments of S/R impldiff.F)
269    o clean up build options files for SunOS, in particular
270      - remove -traditional from xmakedepend
271      - change -xarch=v9 to -xarch=native for f77
272      - restore use of make (not gmake)
273      - add CC=gcc if using g77
274    o changes per discussion of AM, JMC, & EH3 today:
275      - move redundant "grid" variables out of "state" -- in preparation of
276        creating proper CF-style coordinate variables and in recognition
277        that any/all grid variables can be easily copied (eg. w/NCO) from
278        the grid files into any of the other NetCDF files
279      - *always* output grid variables with MDSIO (this will eventually be
280        changed)
281      - output MNC grid file if useMNC is true
282    o every instance of _END_MASTER() has been found and replaced with
283      _END_MASTER( myThid ) in order to satisfy certain picky Sun
284      preprocessors
285    o re-arrange diagnostics pkg initialisation:
286     - allow each package to extend the available diagnostics list
287     - add some checking and fix small problems (multi-threaded, ...)
288     - new S/R diagnostics_fill (replace fill_diagnostics):
289       look through the short list of active diag. (instead of the long list)
290       create function DIAGNOSTICS_IS_ON to tell if a diagnostics is active
291    
292    checkpoint57a_post
293    o this set of changes restores TAMC compatibility
294      tested for global_ocean.90x40x15 adjoint on trough
295      (what a pain)
296    
297    checkpoint57a_pre
298    o some modifs in tools/build_options for sunos
299    o call ctrl_pack for ctrl at end for optimcycle=0
300    
301    checkpoint57
302    o add horizontal bi-harmonic viscosity for vertical velocity
303      - new parameter viscA4W (defaults to viscA4)
304      - affects only non-hydrostatic code
305      - to be done: no-slip lateral boundary conditions
306    o switch back to "sflux = sflux - runoff" in exf_getforcing.F
307    o implement Implicit Vertical advection for pTracers
308    o depth convergence accelerator: replace deltaTtracer by dTtracerLev(k)
309      (implemented also for implicit vertical diffusion & advection)
310    o OBCS: extend application of tracer open boundary conditions
311      to a strip beyond the boundary of width Olx/Oly
312    o add call monitor to TAF set of required calls (flow directives)
313      and include in natl_box_adjoint
314    o untangle a few ifdefs related to cost, ctrl, autodiff
315      untangel a few ALLOW_ECCO_OPTIMIZATION
316    o Apply Cooper & Haines representative estimates also at sub-surface (k=Nr)
317      (theta2, salt2 in ecco_cost_weights)
318    
319    checkpoint56c_post
320    o OBCS as control variables
321      - update ctrl_ad_diff.list, obcs_ad_diff.list
322      - remove balance of obcs controls from default
323      - fix index bug nobcs in ctrl_init
324      - fix dummy fields filen in ctrl_pack
325      - add dummy weights for obcs
326    o lsopt: comment out fort.94 output
327    o optim:
328      - remove unnecessary header files
329      - adjusted namelists
330      - add xerbla.F to Makefile
331    o exf:
332      - multiple modifications of exf_check_range
333      - unit change in exf_init_runoff from [m/year] to [m/s]
334    o mdsio:
335      - i/o-vector change debugLevel from debLevA to debLevB
336    
337    checkpoint56b_post
338    o collect all ptracers output in an new routine (PTRACERS_OUTPUT) that
339      is called from DO_THE_MODEL_IO
340      - rename PTRACERS_STATV_WRITE to PTRACERS_WRITE_TIMEAVE
341      - introduce new parameter PTRACERS_dumpFreq (defaults to dumpFreq)
342      - clean up PTRACERS_WRITE_TIMEAVE a little
343    o add another verification experiment that tests both the cfc-pkg and
344      the gchem-pkg with GCHEM_SEPARATE_FORCING undefined, that is with
345      reactive tendencies treated inline with regular timestepping
346    o testreport previously did a lousy job of comparing small numbers, so
347      we (Ed and Martin) changed the  way the cut-off of 1.e-12 is
348      implemented in the little c-progam tmp_cmpnum.c which is embedded in
349      testreport.
350    o ptracers_monitor: bug fix for MPI
351    
352    checkpoint56a_post
353    o GCHEM: finish reorganizating the package
354      - forward_step calls GCHEM_CALC_TENDENDY, which computes gchemTendency
355        (introduces another 3D-array for each passive tracer, but only if
356        GCHEM_SEPARATE_FORCING is undefined. For GCHEM_SEPARATE_FORCING
357        gchemTendency is not needed because the timestep is done separately)
358      - gchemTendency is added to gPtr in GCHEM_ADD_TENDENCY
359      - GCHEM is now---more or less cleanly---separate from PTRACERS
360      - to bo done: handling of surfaceForcingPtr that is needed for KPP not
361        so clear to me how that should be done in a general way.
362    o CFC: improve formatting of output in CFC_ATMOS
363    o PTRACERS:
364      - rename GCHEM_MONITOR to PTRACERS_MONITOR and call it from MONITOR, so
365        that ALL experiments with ptracers enable can be checked. This makes
366        GCHEM_MONITOR obsolete.
367      - include a runtime parameter PTRACERS_monitorFreq that defaults to
368        monitorFreq
369      - set default PTRACERS_write_mdsio to false if PTRACERS_write_mnc is true
370    o add PTRACERS (PTR0[1-5]) to testreport output
371    o exf_check_range:
372      * default changed to .TRUE.
373      * will check values for niter0 if debLevA, every timestep if debLevB
374      * STOP in exf_check_range if 'out of range'
375      * range chosen generous to allow for potential extema
376      TBD: a version with stricter range
377           permitted range would be runtime parameters
378    o GMredi:
379      - fix ldd97 slope limit ; extend valid domain of Psi-Bolus ;
380      - change S/R GRAD_SIGMA to be valid in the overlap on CS-grid.
381      - change the default GMREDI_OPTIONS.h: compile everything except Visbeck
382    o GCHEM:
383      - move call to GCHEM_FORCING_SEP, so that it is now called before
384        the blocking exchanges to avoid computing halos in GCHEM_FORCING_SEP
385      - introduced a useGCHEM runtime flag and applied it everywhere.
386      - renamed PTRACERS_SEPARATE_FORCING to GCHEM_SEPARATE_FORCING;
387        this flag is now only visible within GCHEM-pkg routines.
388      - moved gchem_forcing_int temporarily into ptracers_forcing, this will
389        be replaced by a 3D array that is computed before thermodynamics; for
390        this, the tracer ID iTracer had to be included into the parameter list
391        of ptracers_forcing
392      - set up forward_step to call gchem_forcing, that will replace
393        gchem_forcing_int in ptracers_forcing
394      - cleaning up
395    o small cleaning of exch2: remove files that are also in eesupp/src
396    
397    checkpoint56
398    o AIM: read land-Fraction earlier (from ini-fixed) and monitor land at nIter0
399    o mnc: cleanup "full" files--fixes out-of-memory problem
400    o ECCO-GODAE
401      * model/, pkg/, code/, input/ now complete w.r.t. ECCO-GODAE production setup
402      * fully merged and updated from SIO code
403        windstress sign change affects
404        - exf_inscal_[u/v]stress
405        - sign of xx_tau[u/v] adjustments
406        - signs in cost_scat
407      * ctrl_get_gen: avoid update of tau[u/v] for first few days
408      * natl_box_adjoint updated and compressed
409      * 2-level checkpoint switch added
410        (tested in EXP. 1 of natl_box_adjoint)
411      * new switch doInitXX to avoid init. of xx_... in optimcycle=0
412    o SVD
413      * first complete, non-hacked version for SVD calculations
414      * all interfaces clean w.r.t to handling of xx_..., adxx_..., g_xx_...
415        performed via direct access files admtlm_vector.it????
416        (via ctrl_[un]pack)
417      * optimcycle is used for ARPACK iteration counter
418        (note: adxx_... is iterated, rather than xx_...; latter is constant!)
419      * EV problem ok w.r.t. L2-norm in tangent space; need to check other norms
420        and add linear transforms where necessary
421      * controls are initial/final SST, SSS (can be extended; need to check
422        ctrl_set_[un]pack_xyz.F according to ctrl_set_[un]pack_xy.F
423    o mdsio: filename info to both stderr, stdout
424    o exf: exf_inscal_runoff: separated default conversion [m/year] -> [m/sec]
425      from sign change
426    
427    checkpoint55j_post
428    o isolate dissipation tendency (allow to keep it out off AB) in V.I form.
429      affects (truncation level) output.txt of global_ocean.cs32x15(&.thsice)
430    o add consistent build directories and .cvsignore files to verification
431    o added viscA4GridMin and viscA4GridMax for more control on biharmonic viscosity
432    o make all subroutines 30 or less characters in length for certain older
433      compilers to work (a thorough search only turns up two routines in
434      pkg/ptracers)
435    o allow (vecinv only) a different horiz. viscosity for Divergence and Vorticity;
436      affects (truncation level) output.txt of global_ocean.cs32x15(&.thsice)
437    
438    checkpoint55i_post
439    o change arrays for total vertical diffusivity (2.D if explicit,
440      3.D if implicit) ; each ptracer uses its own background value.
441    o added mnc/NetCDF output capability to ocn_compon_interf
442      - added bits to verification/cpl_aim+ocn so that it can run on
443        ITRDA and other clusters without the mpirun "-wd" option
444    o fixed wayward endif in mom_fluxform.F that messed up cyclindrical
445      configuration
446    o give MNC the ability to create additional files rather than
447      exceeding a specified file size limit "MNC_MAX_FILESIZE"
448      - the default is just less than 2GB which is a practical limit
449        on many systems
450    o Compute the non-local term using either KPP vertical diffusivity (=default)
451       or the total vertical diffusivity (when KPP_ghatUseTotalDiffus=T)
452    
453    checkpoint55h_post
454    o add verification test "vermix"
455    o remove all the *_ioinc flags and replace with the single global
456      outputTypesInclusive flag
457    o allow to set a vertical profile of vertical diffusivity for T & S
458    o call PACKAGE_BOOT & PACKAGES_READPARMS just after INI_PARMS
459    o useMNC moved from data.mnc to data.pkg
460    
461    checkpoint55g_post
462    o Renaming carbon AD experiments to comply with testreport conventions
463    o Updating adjoint of Held-Suarez hs94.1x64x5
464    o Fixing directives for non-ECCO forward_step, after ECCO directives
465      have been fixed to avoid recomp. of do_oceanic_phys
466      Works well for global_ocean.90x40x15,
467      but recomp. issue remains for carbon
468    o Updating adjoint output for natl_box_adjoint after
469      pkg/ecco has been modified (see taga c55d_post)
470    
471    checkpoint55f_post
472    o rename to timeave_init_fixed
473    
474    checkpoint55e_post
475    o some delicate re-shuffle of store directives to avoid one
476      extra call of do_oceanic_physics
477      affects ecco_the_main_loop.F, do_oceanic_physics.F
478      checkpoint_lev1_directives.h
479    o NB: this may break global_ocean adjoint temporarily,
480      but it is clear how to fix it. Will do later, need this now.
481    
482    checkpoint55d_post
483    o ECCO specific cost function terms (up-to-date with 1x1 runs)
484    o ecco_cost_weights is modified to 1x1 runs
485    o modifs to allow observations to be read in as
486      single file or yearly files
487    o enable to read exf forcing fields as either
488      single file or yearly files (flag useExfYearlyFields)
489    checkpoint55d_pre
490    
491    o switch between 2-level and 3-level checkpointing
492    o add ability of MNC to write local and "partial" (eg. 2D slices where
493      the full 3D field is never actually stored) arrays to NetCDF files
494      with the correct (that is, the complete multi-dimensional) set of
495      array indicies
496      - used in mom_vecinv() to write the diagFreq output
497      - tested (demonstrated) in verification/aim.5l_cs
498    o further mnc-ification
499      - add verification/lab_sea to the list with mnc setups
500      - fix undefined mnc v-names in cd_code
501      - add kpp to mnc output (both snapshot and timeave)
502    o non-hydrostatic code
503      - added horizontal viscosity parameter viscAhW for w, that replaces viscAh
504        in calc_gw.F
505      - defaults to viscAh
506      - this is useful, when I want to use viscAh=0 (e.g., only biharmonic
507        friction or Leith) in non-hydrostatic mode, as long as there is
508        no biharmonic friction and/or Leigth or equivalent for vertical velocity.
509    o pkg/cost
510      extend simple cost functions to treat time mean quantities
511      (T,U,V) via cost_accumulate_state from cost_tile
512    o multidimAdvection on CS-grid: implement the full 3 pass ; overlap needs to
513        be 2 times larger than the minimum size the advection scheme requires.
514        now get a perfect conservation of salt in advect_cs test-exp (updated)
515    o GGL90 (Gaspar's TKE mixed layer)
516      - add horizontal diffusion of TKE, requires exchanges in
517        do_fields_blocking_exchanges, horizontal diffusivity is zero by
518        default. In OPA there is no horizontal diffusion of TKE but the
519        mixing coefficients are computed from a horizontal average of TKE of
520        6 points or so. I think that diffusion has a little more physical
521        justification.
522      - clean up ggl90_calc in the hope of reducing memory usage (this hope
523        was in vain)
524      - mask tke-variable in ggl90_init
525    o fix mask Pb in multidimAdvection on CS-grid (still a small Pb with ipass)
526    o fix biharmonic diffusion on CS-grid.
527    o more obcs stuff
528      - added obcs_external_fields_load as an alternative to the EXF-pkg option
529        for reading boundary values. This routine is called from
530        obcs_prescribe_read
531      - updated obcs_check
532      - changed sponge layer relaxation parameter units for days to seconds
533      - turned off sponge layer by default
534      - changed argument list of OBCS_CALC to myTime+deltaTclock
535      - fixed about 35 bugs in obcs_external_fields_load
536    
537    checkpoint55c_post
538    o obcs stuff
539      - new runtime flags useOBCSprescribe, useOBCSbalance
540      - re-arranged obcs_calc
541      - moved exf_getobcs to obcs_prescibe_read
542        obcs_prescibe_read is called after ORLANSKI stuff
543        and independent of it (executable can contain all).
544      - updated obcs_check
545      - added exchanges for U,V,gU,gV for useOBCSprescribe=.TRUE.
546        These may be necessary only if tangential velocities are
547        prescribed at OB, not sure, needs testing.
548    o make the MNC_CW_ADD_VATTR_* subroutines all take scalar arguments
549      which simplifies the calls and passes the ifc check-bounds test
550    
551    checkpoint55b_post
552    o add a 2D scalar, Z location, exchange_2 routine.
553    o EXCH_UV for CS-grid: add one u,v in the corner-halo region, allows
554      to compute vorticity on a wider stencil [e.g., vort3(0,1)&(1,0)]
555      affects advect_cs results (for the wrong reason): update output.txt
556    o use a more standard discretization for biharmonic viscosity ;
557     (original version still available with #define MOM_VI_ORIGINAL_VISCA4 )
558      affects results of dome set-up: update output.txt
559    o fix gad_advection.F (CubedSphere part), back to version 1.23
560    o fix bug in mom_vi_del2uv.F (CubedSphere part).
561    
562    checkpoint55a_post
563    o merged code to
564      * prescribe/read time-dependent open boundaries
565        (works in conjunction with exf, cal)
566      * sponge layer code for open boundaries
567      * each boundary N/S/E/W now has its own CPP option
568        (healthy for the adjoint)
569    
570    checkpoint55
571    o up-to-date adjoint of previous breaks
572      - natl_box_adjoint: problem for exp 12, 13: cost = 0. needs fix
573      - global_ocean_ebm: adjoint still needs fix (TLM ok)
574      - global_ocean.90x40x15: extended GMREDI_OPTIONS.h
575      - carbon: * remove tr1-related adjoint verification
576                * add ptracer as control variable
577    o remove all tr1-related code (ALLOW_PASSIVE_TRACER)
578      (adjoint stuff still has some tr1 'names', but all use ptracer arrays)
579    
580    checkpoint54f_post
581    o in verification/carbin/ remove experiment related to
582      ALLOW_PASSIVE_TRACER (array tr1 etc)
583      code_ad/ input_ad/ results_ad
584    
585    o add two new packages
586      - ggl90, TKE-model of Gaspar et al. (1990), Blanke+Delecluse (1993)
587        + subject to optimization (memory)
588      - opps, OPPS (ocean penetrative plume scheme) of Paluszkiewicz+Romea (1997)
589        + subject to reordering (many GOTOs in the original code)
590        + so far opps is called for every (i,j) point, and then operates in
591          the vertical; one may want to reorder the scheme in order operate in
592          horizontal slabs at the cost of (a lot of) memory, at least ten new 3D
593          fields)
594    o overhaul of IO so that we now have flags for MDSIO and/or MNC
595      - all verification tests compile and run with linux_ia32_g77
596      - defaults are compatible with current input files--nothing
597        should change if you were not previously using MNC
598      - MNC output has been added in numerous places (eg. timeave)
599        but there are still a few writes not yet do-able with MNC
600        (this is in progress)
601      - flags now allow for either/or/both use of MDSIO and MNC and
602        documentation will soon follow
603      - numerous small formatting cleanups for ProTeX
604    o add version info (via BUILD_INFO.h) to MNC and stdout
605    
606    checkpoint54e_post
607    o add two new packages
608      - pp81 (Packanowski and Philander, 1981), Richardson number and
609        stratification dependent mixing
610      - my82 (Mellor and Yamada, 1982) level 2 turbulence closure scheme
611      - these scheme operate in a similar manner as KPP (without ghat, of
612        course), however there is only one diffusivity that is used for both
613        temperature and salinity (and passive tracers)ç
614    o add offline package
615    o add "verification/dic_example" which contains a dic-gchem-ptracer setup
616      provided by Stephanie Dutkiewicz
617    o genmake2: add ability to search for C compilers and cleanup test output
618      - fixes recent SX-6 problems and should help with others (eg. HP-UX)
619    
620    checkpoint54d_post
621    o data.diagnostics: if levels() is unspecified then ALL levels are diagnosed
622    o bug fix in pkg/cal/cal_addtime.F for TheCalendar='model' day count is wrong,
623      leading to additional shift by 5 days each year of integration
624       (i.e. leading to considerable offset over the years)
625    o update global_ocean.cs32x15/thsice: use z* + RealFW + staggerTimeStep
626    o start with no AB and no initial filtering (staggerTimeStep):
627      change the results of test-exp: hs94.128x64x5, hs94.1x64x5 & hs94.cs-32x32x5
628    o move ocean specific S/R calls from forward_step.F to do_oceanic_phys.F
629      (external loading & cpl import/export stay in forward_step.F)
630    o update coupling interface to work with new staggerTimeStep
631    o change MDSIO ptracers_pickup.
632    o fix a bug in Flux deriv. vs Ts (pkg/aim_v23/phy_suflux_sice.F);
633    o genmake2 fixes so netcdf now works on halem ("osf1_alpha_*")
634    o bugs fixed in kpp_calc.F (ustarY) and external_forcing.F (Qsw & Partial-cell)
635    
636    checkpoint54c_post
637    o replace surfaceTendency U,V,S,T,Tice,Ptr by surfaceForcing U,V,S,T,Tice,Ptr
638      this affects (truncation error) output.txt of test-exp:
639       global_ocean_ebm,_pressure,.cs32x15(&.thsice), lab_sea & natl_box
640    o pkg/kpp: replace pMask by maskC to reduce memory size (one 3.D array)
641    o change Qnet to always be the net heat flux, (+upward).
642      affects (truncation error) output of natl_box exp.: updated
643    o define the max number of passive tracers in PTRACERS_SIZE.h
644    o start to test mnc routinely in exp0
645    o add sea-ice loading (sIceLoad) as a new forcing field.
646    
647    checkpoint54b_post
648    o add testreport "-j JOBS" option for parallel builds (tested and pretty
649       impressive on Altix and Opteron systems)
650    o diagnostics can now write NetCDF files through the mnc package
651    o staggerTimeStep & NLFS: update pkg aim_v23 & thsice; update aim.5l_cs output
652    o staggerTimeStep & multi-Dim-Advec: Exch(U,V) added
653    o staggerTimeStep & NLFS: rescale T & S surface tendencies.
654      affects (truncation error) output of lab_sea & natl_box exp.: updated
655    o change pickup _cd & _ptracers (old pickup restart: usePickupBeforeC54=.T.)
656    
657    checkpoint54a_post
658    o re-write staggerTimeStep: step fwd U,V and then T,S (reverse the order).
659    o change pickup file (to restart from old pickup, set usePickupBeforeC54=.T.)
660    o get a perfect restart when using RealFreshWaterFlux + NLFS (not finished)
661    o update output of all exp. that use staggerTimeStep (except hs94.cs & 2D)
662      + internal_wave (shift in writing cg2d_ini_res)
663    checkpoint54a_pre
664    
665    o require forcing_In_AB=.FALSE. and remove surfaceTendencyTice from pkg/seaice
666      - modified verification/lab_sea and natl_box_adjoint accordingly
667    
668    checkpoint54
669    checkpoint53g_post
670    o prepare splitting of thermodynamics: store convect.counter in common block
671    o store d.etaH/dt (instead of Div.hV) in common ; affects (truncation error)
672      output of: global_ocean.cs32x15, global_ocean_pressure & internal_wave
673    o cnh's modifs to gad_advection.F for cube-sphere multi-dim advection
674    o new parameter "nShapS": allow to treat differently S from T (nShapT)
675    o separate Vert.Advec.Scheme from horizontal Advec.Scheme (T & S)
676    o restore adjointability
677    
678    checkpoint53f_post
679    o modifications that affect AIM experiment results (update output.txt):
680      - no advection of S at top interface Nr (pkg/generic_advdiff)
681      - air surf. temp: use lower value between T0 & T1 (phy_suflux_prep.F)
682    o AIM pkg:
683      - include stability function into surf.Flux derivative relative to Tsurf
684      - calculate clear-sky radiation & surface temp. change (for diagnostics)
685    
686    checkpoint53e_post
687    o merged cylindrical coord configuration
688    o added rotating_tank experiment
689    checkpoint53e_pre
690    
691    checkpoint53d_post
692    o faulks upgraded: re-generate output.txt that have changed:
693       global_ocean.cs32x15(.thsice), global_ocean_ebm, lab_sea, natl_box
694    o Added Bryan and Lewis, 1979, vertical diffusivity profile
695    o pkg/land: fix problem in land_impl_temp.F when Ts arround freezing.
696    o added Sadourny discretization of Coriolis
697      - moved some PARAMETERS from mom_*_coriolis.F to PARAMS.h
698      - re-enabled use of omega3 (useAbsVorticity)
699    o added Leith (Phys. Fluids, 1968, 10, 1409-1416) variable viscosity for
700      vector-invariant equations
701      - note: does *NOT* yet use same variable visc. for no-slip BC's.
702    o recoded the lateral boundary conditions for non-hydrostatic w, this time
703      hopefully without new bugs
704      - still to do: biharmonic friction for non-hydrostatic w
705    o bug fix for mom_*_sidedrag: rA -> rAs,rAw
706    o regenerated output for global_ocean_ebm (using -ieee)
707    o fix mom_common/mom_*_sidedrag.F to recover numeric output prior to creation
708      of mom_common/
709      - simple re-arrangement of instructions within computing loop
710    checkpoint53d_pre
711    
712    checkpoint53c_post
713    o new version of coupler for AIM+thSIce and Ocean models
714     - new fields (SSS,h_MixLayer ...) are exchanged through the coupler
715     - land: change Evaporation of snow
716     - aim with thSIce: fix bug in snow-precip.
717    o changes to calc_gw.F
718     - bug fix for w^2 term near the bottom boundary
719     - (hopefully) improve the lateral slip boundary condtions for use with
720        partial/looped cells
721      - because the bug fix changes two verification experiments anyway (exp5
722        and plume on slope), change the lateral boundary condition from half slip
723        to the value of no_slip_sides
724    o ebm package wasnt quite up to the parallel task. Fixed that.
725    o Adjoint of EBM package + verification
726      * TLM OK
727      * ADM not sufficient accuracy yet (needs further testing)
728    
729    checkpoint53b_post
730    o Commiting new energy balance model to repository
731      * package is pkg/ebm
732      * verif. is verification/global_ocean_ebm
733      * references are in ebm_driver.F
734    C      * X. Wang, P. Stone and J. Marotzke, 1999:
735    C        Global thermohaline circulation. Part I:
736    C        Sensitivity to atmospheric moisture transport.
737    C        J. Climate 12(1), 71-82
738    C      * X. Wang, P. Stone and J. Marotzke, 1999:
739    C        Global thermohaline circulation. Part II:
740    C        Sensitivity with interactive transport.
741    C        J. Climate 12(1), 83-91
742    C      * M. Nakamura, P. Stone and J. Marotzke, 1994:
743    C        Destabilization of the thermohaline circulation
744    C        by atmospheric eddy transports.
745    C        J. Climate 7(12), 1870-1882
746    o Will need long integration testing.
747    checkpoint53b_pre
748    
749    o created pkg/mom_common in preparation for variable viscosity code
750    o renamed and split ini_mnc_io
751    o pkg/land: only liquid water can diffuse or run-off
752    o AIM using land: soil wetness availability account for snow
753    o added printMapIncludesZeros to EEPARMS namelist
754    
755    checkpoint53a_post
756    o restore TAMC compatiility for exf_constants.h
757      (no usage of _d in .h)
758    o Additional storing in KPP in an attempt to fix
759      TAMC problem with adjoint of KPP code.
760      (no success...)
761      Identical code works fine with TAF.
762    o Fixed verification/natl_box_adjoint experiments.
763      Initial conditions for Experiments 4 and 5 are now
764      set to "SEAICE_initialHEFF=0", which is more reasonable,
765      and hence gives results closer to that of Experiment 2.
766      Experiment 6 requires "SEAICE_initialHEFF=1" and is working once again.
767    
768    checkpoint53
769    o stable checkpoint with up to date adjoint
770      TLM tests only run for 1st grdchk at the moment. Its unclear if
771      this is a code problem or a problem with resetting code for g_
772      variables that is being autogenerated wrong.
773      Otherwise:
774       o code is testing well on Linux (P4 and AMD), SGI Alitix and Irix.
775       o IBM and HPQ tests could not be run because NCAR (IBM) and NCCS (HPQ)
776         are both still in chaos following their security messes.
777    
778    checkpoint52n_post
779    o various build (mostly genmake2) changes including:
780      - ability to build/run more easily within Cygwin and MacOS X
781      - more complicated Makefiles with a stronger dependence on Gnu Make
782    o pkg/seaice modifications
783      - added SEAICE_deltaTtherm and SEAICE_deltaTdyn in order to reduce cost
784        of pkg/seaice.  For time being, SEAICE_deltaTtherm = deltaTtracer
785        and SEAICE_deltaTdyn must be integer multiple of SEAICE_deltaTtherm
786      - changed proxy for geostrophic velocity and updated
787        verification/lab_sea/results/* accordingly
788      - removed ADI dynamic solver from pkg/seaice
789      - replaced all occurrences of 1/CSTICE or 1/CSUICE
790        they caused a segmentation fault under certain circumstances
791      - removed unused field WATR
792    o bug fix in exf-package:
793      -  u/vwind are defined at tracer/pressure points. The values of u/vwindmask
794         now reflect that.
795    o small fixes in the cal-package:
796      - set reference date to first of January for the 'model'-calendar
797      - failing of cal_CheckDate in cal_FullDate now actually results in a
798        warning
799    o small modifications in exf-package to increase flexibility:
800      - change climtempfreeze from a constant to a runtime parameter that
801        is specified in data.exf_clim (default = -1.9)
802      - introduce runtime parameter windstressmax (default = 2.0) that is
803        specified in data.exf
804    o fix in pkg/thsice:
805      - update Evap (Tsf changes) to be consistent with Latent heat flux
806      - maintain snow precipitation unchanged (even if Tsf >=0).
807     both affect global_ocean.cs32x15/(.thsice) results ; update output.txt.thsice
808    
809    checkpoint52m_post
810    o commented out call to cal_GetMonthsRec in swfrac.F and added some comments;
811      this does not change any of the results, because myIter=0 in swfrac always
812      lead to jwtype=2 (the default anyway), but the call makes little sense
813      as a default
814    o no SHORTWAVE HEATING lost at the bottom; affects lab_sea output; updated
815    o re-work pkg/thsice interface with ocean & bulk_force pkg to allow
816      atmospheric model (AIM) to use it ; add a "slab ocean" component.
817    o pkg/aim_v23 : add interface to call thermodynamic sea-ice (thsice)
818    o update global_ocean.cs32x15/output (.thsice) after changes in pkg/thsice
819    o extending useSingleCpuIO option to work with new exch2 I/O format
820      - old-style, missing-tile I/O is still accessible by defining CPP
821        option MISSING_TILE_IO in pkg/mdsio/MDSIO_OPTIONS.h
822    o monitor package can write to NetCDF files and/or STDOUT
823    o fix re-definition bug in MNC_GRID_INIT_ALL()
824    o doc/api_reference added
825      - a framework for building an API Reference Manual using all
826        the "protex" comments embedded in the code
827      - see pkg/generic_advdiff, pkg/mnc for examples
828    o MNC:
829      - make myThid the last argument for all subroutines
830      - new default for the grid files: NO halos
831      - add useMNC and other run-time flags
832      - separate per-file-group handling of the unlimited dimension
833      - early initialization of MNC w/ flags read at the eeboot stage
834        and lookup tables set in initialize_fixed
835      - write all MNC output to a freshly-created directory
836    o thermodynamics sea-ice (thsice) is now tested with testreport.
837    o test exp. advect_cs: use the standard cs32 grid ; update output.
838    o Added functionality to grdchk:
839      pick global i,j,k position (or nearest wet) where to perform check.
840    
841    checkpoint52l_post
842    o new land formulation:
843      a) use ground enthalpy as prognostic variable to ensure exact
844         energy conservation.
845      b) account for water temperature and for latent heat of freezing
846         in all processes (rain, run-off, ground storage)
847      c) compute surface and ground temperature implicitly.
848    o aim_v23 modification for new land model:
849      a) Treat separately land / ocean / sea-ice surface fluxes
850         to allow implicit computation of land & sea-ice surface temp.
851      b) add snow precipitation.
852    checkpoint52l_pre
853    
854    o fix problem (when 1 tile is empty) in min/max monitor.
855    o fix bug in timeAve (u,vZetatave in common block)
856    o no restoring terms under ice for pkg/seaice
857    o moved SEAICE_DO_DIAGS to S/R do_the_model_io
858    o bug fix for pkg/seaice pickups
859    o change MNC call locations and add coordinate attributes
860    o mask vorticity at the corner ; affects global_ocean.cs output
861    o fix ini_curvilinear_grid ; test CS experiments with mpi
862    o aim.5l_cs: modify sea-ice fraction input file ; update output
863    o exf: Adding simple range check routine for exf fields after read
864    
865    checkpoint52k_post
866    o removed top-level directory "diags"
867      - we migrated these functions to pkg/timeave a long time ago
868    o fixed bug introduced in pkg/monitor
869    o accuracy ctrlprec = 32 insuffient for gradient checks using
870      averaged fields (I/O via cost_averages)
871      -> use ctrl.h in active_file*.F to control I/O precision.
872    o paramter list of CALL ADTHE_MAIN_LOOP
873    o added initial do_the_model_io to ecco
874    o set all ...Freq to zero after cost_final
875      taveFreq, dumpFreq, pChkptFreq
876    o MAIN and ECCO differentiability restored
877    
878    checkpoint52j_post
879    o fix nchklev in tamc.h for global_ocean...
880    o Bringing flow and store directives up-to-date
881    o Really weird thing!
882      call do_the_model_io from the_main_loop with parameters
883      starttime, niter0, instead of mytime, myiter, screws up the adjoint:
884      TAF places a faulty re-call of initialise_varia
885      after cost_final call which should not be there!
886    checkpoint52j_pre
887    
888    o tweaked pkg/mdsio/mdsio_readfield.F to have more useful logic and
889      error messages
890    o fix errorMessageUnit problem. remove old version of mdsio from eesup/src
891    o modified pkg/mdsio/mdsio_readfield.F -- commented out diskspace-eating
892      debug output
893    o modfied mdsio_writefield.F to handle non-global files under exch2
894    o update AIM Cubed-Sphere exp (aim.5l_cs):
895      use land, p*, and Franco.M forcing fields with fractional land-sea mask.
896    o update AIM Equatorial Channel exp:
897      use standard aim_v23 pkg and change the SST fields to be more realistic.
898    o globalFile now works with exch2
899    o merged mom_vi_calc_relvort3.F, mom_vi_del2uv.F and
900      mon_vort3.F to make exch2 friendly with >1 tile per cube face
901    o limit timeave output for hi-res integrations
902    
903    checkpoint52i_post
904    o removed CPP_EEOPTIONS.h from verification/*/code
905    o updated the tiling for global_ocean_pressure
906    checkpoint52i_pre
907    
908    o update aim.5l_LatLon using the standard aim_v23 pkg (instead of
909      the old pkg/aim).
910    
911    checkpoint52h_pre
912    o shortened offending lines in exch2_send_rx?.template
913    o modified W2_EXCH2_TOPOLOGY.h to reflect current topology generator
914    o genmake2: add build & execute test for NetCDF code
915    o modded aim.5l_cs and global_ocean.cs32x15 to use exch2, removed
916      old grid files, moved CPP_EEMACROS.h to pkg/exch2
917    o added useSingleCPUIO flag to monitor.F and mon_out.F
918    o added temporary hack for sea-ice dynamics on the cube:
919      metric terms are ignored and it is assumed that
920      there is no sea-ice on cube faces 1, 2, 4, and 5.
921    
922    checkpoint52f_post
923    o added exch2 cube routines to pkg/exch2
924    o modified hs94.cs-32x32x5 to use it
925    o added new grid files to hs94.cs-32x32x5
926    checkpoint52f_pre
927    
928    checkpoint52e_post
929    o first attempt to solve vertical advection (momemtum & T,S) Implicitly
930      for now, only implemented with Lax-Wendroff (tracer) and VectInv (mom)
931    o step.1 : rewrite (as in MultiDimAdvec) explicit tracer stepping (gad_calc_rhs.F)
932      to work with implicit vertical advection and AB ;
933    o Implic. vert. advect. implemented with flux-limiter, c2, u3 and c4 scheme.
934    o update output (gad_calc_rhs.F changes affect truncation error) of exp:
935     front_relax, global_ocean .cs32x15 & _pressure, internal_wave, lab_sea, natl_box
936    checkpoint52e_pre
937    
938    o fix bug in impldiff when Nr=1 or Nr=2.
939    o change the albedo in pkg/thsice (as GISS model); update output_thsice.txt
940    o little changes (fftpack, aim/phy_suflux ...) to pass the "checkbounds" test.
941    o fix deltaTfreesurf in update_cg2d.F (left from c45d_post): affects 2
942      ocean exp. that use NLFS: global_ocean.cs32x15 & global_ocean_pressure
943    o added exf_offset_atemp for Celsius to Kelvin conversion of atemp
944    o eesupp/inc/EEPARAMS.h: MAX_NO_PROCS=1024
945    o pkg/cal/cal_toseconds.F fix for long (>60 years) integrations
946    
947    checkpoint52d_post
948    o re-tiled many expts to allow comparison in multi-processor mode
949      - updated output accordingly
950      - changed hs94.128x64x5 to use intial T from file rather than use the
951        random number generator
952    o add infrastructure in model routines for fizhi and gridalt packages
953    o change global_ocean.cs32x15 because of Pb with GM on CS-grid.
954    o changes to permit hi-res, cubed-sphere, configuration
955      - added useSingleCpuIO capability to mdsio_readfield.F and exf_interp.F
956      - added "#undef ALLOW_USE_MPI" support to eesupp/src/scatter_2d.F
957      - added pkg/exf/exf_set_uv.F for on-the-fly interpolation and rotation
958        of surface winds for the cube
959      (A verification experiment, based on global_ocean.cs32x15, but with
960       pkg/seaice turned on and with on-the-fly interpolation from the NCEP
961       Gaussian grid is described in MITgcm_contrib/high_res_cube/README_ice,
962       complete with example surface forcing files and matlab scripts to look
963       at the output and compare it with that of global_ocean.cs32x15.)
964    o packages pickup: $PKG_write_pickup called from packages_write_pickup ;
965       but call to $PKG_read_pickup remains in $PKG_ini_vars: change CD_code
966       to conform to the common rule.
967    o check-in the first version of the coupler + interfaces for multi-components
968       set-up (Chris's version, adapted to work on cubed-sphere grid); use MPI.
969       - imply modification in eesupp (MPI_COMM_MODEL < MPI_COMM_WORLD).
970       - atmosphere: use pkg aim_v23 ; no sea-ice for now.
971       - ocean and atmos. share the same horizontal grid (no interpolation).
972       - exchanges between component and coupler is done in forward_step
973    o new verification experiment, cpl_aim+ocn:
974       coupled ocean - AIM atmosphere on cs32 grid.
975    
976    checkpoint52d_pre
977    o surface flux time-av diagnostic modified for NLFS & RealFreshWater
978    o add call to AIM pendant S/R in external_forcing & external_fields_load.F
979    
980    checkpoint52c_post
981    o Fixing OPTIONS files for use of exf:
982      - ALLOW_ECCO not needed
983      - ECCO_CPPOPTIONS.h should not be included from
984        EXF_OPTIONS.h/CAL_OPTIONS.h, but from CPP_OPTIONS.h
985        as has been the case in the past
986      - All cal/exf related options are currently set in
987        single header file. Should remain so until we
988        agree on a consistent new policy.
989    o Removed all CPP options for cal package.
990      cal has no more CPP options.
991    o modified STOREs in GAD_ADVECTION
992    o corrected key comp. for passkey
993    
994    checkpoint52b_post
995    o new version: "thsice" of thermodynamic sea-ice pkg:
996     - only 1 call from "forward_step"; conserve energy, fresh-water & salt.
997     - modify fluxes but do not change theta directly;
998     - light cleaning of bulk_force pkg (still needed with thsice).
999     - add forcing fields + data files in global_ocean.cs32x15/inp_thsice
1000       + output_thsice.txt to test thsice & bulk_force pkgs.
1001    o put back grid & coriolis monitor output.
1002    checkpoint52b_pre
1003    
1004    o call "do_the_model_io" at the beginning (=> dump initial state)
1005    o use rhoConstFresh to define convertEmP2rUnit (was done only in ocean-P)
1006    o sort out interplay between tamc.h and PTRACERS_OPTIONS.h
1007      (need info on NUMBER_OF_PTRACERS)
1008    o missing PTRACERS headers in some routines
1009    o default tamc.h in pkg/autodiff supposed to crash
1010      to ensure proper customization
1011    
1012    checkpoint52a_post
1013    o modifications to make FREEZE (_SURFACE) flux visible to pkg/kpp
1014      - moved surfaceTendencyTice from pkg/seaice to main code
1015      - subroutine FREEZE_SURFACE now limits only surface temperature
1016        and is called in FORWARD_STEP (if allowFreezing=.TRUE.)
1017        => new output.txt for global_ocean.cs32x15 & global_with_exf
1018      - keep the old version of FREEZE for backward compatibility;
1019        global_ocean.90x40x15 left unchanged using useOldFreezing=.TRUE.
1020      - corresponding modifs for TAF directives
1021    o added surface flux output variables to TIMEAVE_STATVARS
1022    o time-averaged output for pkg/ptracers
1023    o fix problem with ECCO-related divided adjoint
1024      snapshot file was missing some ad-variables,
1025      but no warning issued by TAF
1026    
1027    checkpoint52a_pre
1028    o some bug fixes for #undef REAL4_IS_SLOW
1029    
1030    checkpoint52 = ecco_c52_e35
1031    o after testing adjoints of checkpoint51u_post
1032    
1033    checkpoint51u_post
1034    o this is the post-merge tag of the ecco-branch
1035      (needs another clean checkout and re-testing)
1036    o various natl_box_adjoint configs and results
1037      but no input fields
1038    o obcs adjoint not yet merged!
1039    o adding ECCO global 1x1 and 2x2 degree configs
1040    o updated/extended some store directives for seaice
1041    o pkg/ecco now containes ecco-specific part of cost function
1042    o top level routines the_main_loop, forward_step
1043      supersede those in model/src/
1044      previous input data.cost now in data.ecco
1045      (new namelist ecco_cost_nml)
1046    o in ctrl/ cleaned some cross-dependencies and updated CPP options
1047    
1048    checkpoint51t_post
1049    o this is the pre-merge tag of the ecco-branch
1050    
1051    checkpoint51s_post
1052    o cleanup: removed various obsolete #define-s
1053    o pkg/grdchk
1054      - renamed 'ph-grd 3' to 'grad-res'
1055      - merged from ecco-branch
1056        (remaining bug fixes for obcs gradient checks)
1057      - additional high-precision output for testreport
1058        (grep for 'precision_grdchk_result')
1059    
1060    checkpoint51r_post
1061    o cleaned up or removed CPP_OPTIONS.h files
1062    o updated results for some of the tests to reflect the new
1063      optfiles (-ieee: -O0 instead of -O3)
1064    
1065    checkpoint51q_post
1066    o cd_code:
1067      - moved cd_scheme.F -> cd_code_scheme.F
1068      - separate read_checkpoint from cd_code_ini_vars.F
1069      - separated cd_code part from write_checkpoint
1070      - updated AD_SOURCE, generated .flow
1071      - added CD_CODE_VARS.h to the_main_loop
1072      - Moving cd_code specific store dir.s to pkg/cd_code
1073    o ctrl package totally restructured
1074      - pack/unpack now optional and decoupled from
1075        xx_/adxx_ I/O
1076      - ctrl_pack/unpack cleaned
1077        (new routines ctrl_init_ctrlvar.F, pkg/ctrl/ctrl_init_wet.F)
1078      - confined inclusion of AD_CONFIG.h to where necessary.
1079    o genmake2: changed AD dependencies from SRCFILES to AD_FILES
1080    
1081    checkpoint51p_post
1082    o create initial cd_code package
1083      - WARNING: the initialization of fields within packages
1084        is broken because INI_FIELDS (where READ_CHECKPOINT is
1085        called) is called *before* PACKAGES_INIT_VARIABLES.
1086        This seems to be backwards since
1087        1) its important to initially zero the entire field including
1088           the "halo" region
1089        2) zeroing needs to be done *before* reading pickup info
1090           and the pickup call might or might not happen in any
1091           particular run
1092      - this needs an overhaul for R2
1093    o cleanup verification/*/code/GMREDI_OPTIONS.h
1094    o cleanup testreport
1095    
1096    checkpoint51o_post
1097    o bringing AD-related verif's up to date with genmake2
1098      - carbon/
1099      - global_ocean.90x40x15/
1100      - hs.94.1x64x5/
1101      - (front_relax/ )
1102    o (more) consistent directory structure for AD-related verif.'s
1103    o cleaning ALLOW_GRADIENT_CHECK -> ALLOW_GRDCHK
1104    o cleaning some ALLOW_TANGENTLINEAR_RUN -> ALLOW_AUTODIFF
1105    o bug fix in find_alpha.F for MDJWF:
1106      - modif. to alpha = 1/D*( dN/dT - rho*dD/Dt) to account for
1107        change rho -> rho-rhoConst
1108      - replace call find_rho to find_rhonum
1109    o initialisation of rFlx extended to full array (required by TAF)
1110      and shifted to thermodynamics
1111    o removed PTRACERS.h in ptracers routine
1112    o added surfacetendencyPtr to S/R parameter list pracers_forcing
1113    o Correcting initialisations for fVerU, fVerV.
1114    o helping TAF to recognize dependency on kUp for
1115      fVerU, fVerV
1116    o Replacing delZ by delR in pkg/ctrl/
1117    checkpoint51o_pre
1118    
1119    o undid all of cp51 pending some cleanups and discussion
1120    o added the [#include "AD_CONFIG.h"] statement to all files that need
1121       it for adjoint/tl #defines
1122     o re-worked the build logic in genmake2 to support AD_CONFIG.h
1123     o removed tools/genmake since it no longer works
1124    
1125    checkpoint51n_post
1126    o modifications to make FREEZE flux visible to pkg/kpp
1127      - moved surfaceTendencyTice from pkg/seaice to main code
1128      - FREEZE & EXTERNAL_FORCING_SURF moved to FORWARD_STEP
1129      - subroutine FREEZE now limits only surface temperature
1130        (this means new output.txt for global_ocean.90x40x15,
1131         global_ocean.cs32x15, and global_with_exf)
1132    o added surface flux output variables to TIMEAVE_STATVARS
1133    
1134    checkpoint51n_pre
1135    o added "#undef ALLOW_USE_MPI" support to eesupp/src/gather_2d.F
1136    o added [#include "AD_CONFIG.h"] to all files that use the following
1137      defines:
1138      - ALLOW_ADJOINT_RUN
1139      - ALLOW_TANGENTLINEAR_RUN
1140      - ALLOW_ECCO_OPTIMIZATION
1141    
1142    checkpoint51m_post
1143    o incorporating changes to get TAMC/TAF working in genmake2
1144      - added new TAMC/TAF targets to the Makefile
1145      - various cleanups related to package handling
1146      - these changes should have no impact on the forward model
1147    
1148    checkpoint51l_post
1149    o changes to pkg/cal and pkg/exf that allow and speed-up long integrations
1150     - faster algorithm for cal_addtime.F
1151     - higher precision arithmetic, where needed
1152     - pkg/exf *startdate variables changed from calendar date to
1153       model integration time
1154    
1155    checkpoint51l_pre
1156    o add salt flux to forcing fields (used when salty sea-ice forms or melts)
1157    o enable to apply restoring on S & T only in a latitude band.
1158    o update output (external_forcing_surf modif. change truncation errors)
1159       global_ocean.cs32x15, global_ocean_pressure, lab_sea, natl_box
1160    
1161    checkpoint51k_post
1162    o change default to: useNHMTerms=.FALSE.
1163     - new output.txt for atmos. exp (NHMTerms were wrong but now turned off):
1164       adjustment.128x64x1, aim.5l_Eq_Ch, aim.5l_LatLon, hs94.128x64x5, hs94.1x64x5
1165     - new output.txt for idealized oceanic exp (NHMTerms are now turned off):
1166        exp1, ideal_2D_oce
1167     - change data file of oceanic exp (leave the results unchanged for now)
1168       exp2, global_ocean[90x40x15,pressure,with_exf], lab_sea, natl_box
1169    o set gravitySign automatically according to the vertical coordinate used
1170      (and change sign: back to "logical" value: +1 when R increases downward)
1171    o added an include barrier to the default CPP_OPTIONS.h and
1172      removed the #include "PACKAGES_CONFIG.h" from it
1173    o fix a bug in mom_fluxform (from chk51j_post)
1174    
1175    checkpoint51j_post
1176    o Added some AD-related initialisations in mom_vecinv/ mom_fluxform/
1177    o adjusted some flow directives (dynamics, thermodynamics)
1178    o added cubed-sphere case to test cost function
1179    o heckpoint_lev?_directives.h now include
1180      pkg-specific directives, whenever those are enabled.
1181    o Moving pkg-specific flow directives to corresponding pkg
1182      Moving pkg-specific differentiation list to corresponding pkg
1183      Moving pkg-specific checkpoint-lev storing to corresponding pkg
1184      ---> conventions are
1185           - pkg_ad.flow
1186           - pkg_ad_diff.list
1187           - pkg_ad_check_lev?_dir.h
1188    
1189    checkpoint51i_post
1190    o merge of the branch-genmake2
1191      - knowingly breaks the older genmake
1192      - automatic generation of all ALLOW_${PKG_NAME} defines
1193      - new testreport supersedes older testscript
1194    o large number of package inclusion and header cleanups
1195    checkpoint51i_pre
1196    
1197    o new output for adjoin of global ocean vector invariant
1198    o fix for checkpoint of thermodynamic seaice and bulk_force
1199    o Steph gchem various cleanups
1200    
1201    checkpoint 51h_pre
1202    o fixed some single quotes
1203    o modify gchem pkg to be more versatile
1204    o modifications to clean up dic pkg
1205    o add pkg/cfc
1206    
1207    checkpoint51g_post
1208    o provide links for ARPACK and ADM*TLM
1209    o start removing verification/vero/
1210    o update AD-related stuff
1211      (in particular some IF-statements in model/src/ routines)
1212    o adding to verification/global_ocean.90x40x15/ adjoint of
1213      vector invariant setup
1214    
1215    checkpoint51f_post
1216    o update advect_xz/output.txt (left from chekpoint50f_pre, ini_vel modif)
1217    o Mods and bug fixes to pkg/cal, pkg/exf, etc., needed for computation
1218      of tracer Green's fucntions for ocean inversion project.
1219    
1220    checkpoint51f_pre
1221    o Added on-the-fly spatial interpolation capability
1222        "USE_EXF_INTERPOLATION" to pkg/exf.
1223      - This is a temporary Cartesian-grid hack until
1224        the super-duper ESMF coupler becomes available.
1225      - See verification/global_with_exf/README for usage example.
1226      - Removed obsolete EXFwindOnBgrid and SEAICEwindOnCgrid
1227        flags and modified pkg/seaice accordingly.
1228    o Bug fix to pkg/ptracers, pkg/generic_advdiff/gad_calc_rhs.F,
1229        and pkg/kpp/kpp_transport_ptr.F for dealing with tracer
1230        non-local transport term.
1231    
1232    checkpoint51e_post
1233    o pkg/mom_vecinv:
1234      - changes related to hFac: use S/R from aim.5l_cs/code (coriolis & vort.advec)
1235      - add time-ave diagnostic of vorticity advection.
1236    o add diagnostic (instantaneous) of Shapiro Filter effect for T,S & UV.
1237    o atmospheric geopotential: - include water vapor effect;
1238                                - use p^kappa in p* scaling.
1239    o land pkg: read initial state from files.
1240    o AIM : add stratospheric drag in the upper level
1241    
1242    checkpoint51d_post
1243    o change function ILNBLNK (=> last non-blank char.) to return zero
1244      (instead of the string length) when the string. char is empty
1245    o fix bug in dxV,dyU definition of cubed-sphere grid (ini_curvilinear_grid)
1246    o modif's to restore TAMC compatibility
1247      - avoid "_d" constructs in PARAMETER declarations
1248      - avoid interrupting lines by comment lines
1249    o Use cluster local disks for purely local I/O
1250      vs. globally visible disks needed for ctrl stuff
1251      (merged from ecco-branch)
1252      - new active_file_loc... and mdsio_..._loc routines
1253      - affected packages: autodiff, ctrl, grdchk, mdsio
1254    
1255    checkpoint51c_post
1256    o introducing integer flag debugLevel
1257      - debLevZero: suppress all output
1258      - debLevA: standard
1259      - debLevB: equivalent to current debugMode=.TRUE.
1260    o introducing pathname variable mdsioLocalDir for mdsio
1261      that can be specified at runtime
1262      (relevant for cluster I/O to local disk)
1263    `
1264    >>>checkpoint51b_post
1265    adjoint of ptracers
1266    o disentangled ALLOW_PTRACERS using new ALLOW_GCHEM
1267    o extended state in checkpoint_lev?, etc.
1268    o alternative ctrl/cost with tr1 or ptracer
1269    o new comlev1_bibj_ptracers
1270    o new ptracers.flow
1271    o made ptracers_forcing same as external_forcing_s
1272    o added verification for adjoint with ptracers
1273      - adm. grdchk fail so far
1274      - tlm. grdchk OK
1275      - checked tapelev/comlev recomputations: OK
1276    <<<checkpoint51b_pre
1277    
1278    checkpoint51a_post
1279    o added code for biogeochemistry package
1280      - new pkgs: dic, gchem
1281      - new subroutine in ptracers
1282    o modify code for biogeo packages
1283      - in model/src and pkg/ptracers
1284    
1285    checkpoint51
1286    o merged improved packages from ecco-branch (ecco_c51_e34) to MAIN
1287      - autodiff, ctrl, exf, grdchk, seaice
1288      kept separate packages
1289      - cost, ecco, obcs
1290    o updating verification/internal_wave, verification/lab_sea
1291      (all verif's pass on faulks.lcs.mit.edu RH7.3, g77)
1292    
1293    checkpoint50i_post
1294    o Preparing next differentiable checkpoint and sync
1295      of MAIN vs. ecco-branch
1296      - 1st step: this one is differentiable
1297      - 2nd step: next(! not this) checkpoint will have updated packages
1298                  from ecco-branch)
1299    
1300    checkpoint50h_post
1301    o new package "land": provide surface BC for AIM.
1302    o pkg/ptracers: apply exchange to all tracers after reading a pickup.
1303    o AIM: use true log-P extrapolation for near surface temp (= as in
1304        F.Molteni paper) ; update aim.5l_cs output.
1305    o New variable in PARM03 for more flexible pickups
1306    o somewhat cleaned package initialisation sequence for
1307      ctrl/ cost/ ecco/ in packages_... routines
1308    
1309    checkpoint50g_post
1310    o merged with release1_p17 (pkg/seaice and verification/lab_sea)
1311      - added SEAICE_MULTILEVEL for 8-category sea-ice thermodynamics
1312      - LSR sea-ice dynamic solver moved to SouthWest B-grid location and
1313        made the default because of faster convergence than ADI
1314    
1315    checkpoint50f_post
1316    o exf totally restructured
1317      - replaced all exf_set_ routines by one generic exf_set_gen.F
1318        and made sure it's still differentiable
1319      - split exf_getffields.F into new exf_bulkformulae.F
1320      - exf_inscal_/exf_outscal_ suffixes consistent with field names
1321      - updated adjoint/makefile and relevant verif. accordingly
1322    
1323  checkpoint50f_pre  checkpoint50f_pre
1324   o AIM: near surface temp evaluated from Pot.Temp: fixed for Part Cell ;   o AIM: near surface temp evaluated from Pot.Temp: fixed for Part Cell ;
1325        + change Evap over land. (update aim.5l_cs output)        + change Evap over land. (update aim.5l_cs output)

Legend:
Removed from v.1.158  
changed lines
  Added in v.1.526

  ViewVC Help
Powered by ViewVC 1.1.22