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

Legend:
Removed from v.1.193  
changed lines
  Added in v.1.505

  ViewVC Help
Powered by ViewVC 1.1.22