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

Legend:
Removed from v.1.163  
changed lines
  Added in v.1.450

  ViewVC Help
Powered by ViewVC 1.1.22