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

Legend:
Removed from v.1.167  
changed lines
  Added in v.1.566

  ViewVC Help
Powered by ViewVC 1.1.22