/[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.10 by cnh, Wed Jul 8 15:40:09 1998 UTC revision 1.169 by heimbach, Fri Jun 27 02:15:54 2003 UTC
# Line 3  $Header$ Line 3  $Header$
3      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
4      ==============================      ==============================
5    
6    >>>checkpoint51b_post
7    adjoint of ptracers
8    o disentangled ALLOW_PTRACERS using new ALLOW_GCHEM
9    o extended state in checkpoint_lev?, etc.
10    o alternative ctrl/cost with tr1 or ptracer
11    o new comlev1_bibj_ptracers
12    o new ptracers.flow
13    o made ptracers_forcing same as external_forcing_s
14    o added verification for adjoint with ptracers
15      - adm. grdchk fail so far
16      - tlm. grdchk OK
17      - checked tapelev/comlev recomputations: OK
18    <<<checkpoint51b_pre
19    
20    checkpoint51a_post
21    o added code for biogeochemistry package
22      - new pkgs: dic, gchem
23      - new subroutine in ptracers
24    o modify code for biogeo packages
25      - in model/src and pkg/ptracers
26    
27    checkpoint51
28    o merged improved packages from ecco-branch (ecco_c51_e34) to MAIN
29      - autodiff, ctrl, exf, grdchk, seaice
30      kept separate packages
31      - cost, ecco, obcs
32    o updating verification/internal_wave, verification/lab_sea
33      (all verif's pass on faulks.lcs.mit.edu RH7.3, g77)
34    
35    checkpoint50i_post
36    o Preparing next differentiable checkpoint and sync
37      of MAIN vs. ecco-branch
38      - 1st step: this one is differentiable
39      - 2nd step: next(! not this) checkpoint will have updated packages
40                  from ecco-branch)
41    
42    checkpoint50h_post
43    o new package "land": provide surface BC for AIM.
44    o pkg/ptracers: apply exchange to all tracers after reading a pickup.
45    o AIM: use true log-P extrapolation for near surface temp (= as in
46        F.Molteni paper) ; update aim.5l_cs output.
47    o New variable in PARM03 for more flexible pickups
48    o somewhat cleaned package initialisation sequence for
49      ctrl/ cost/ ecco/ in packages_... routines
50    
51    checkpoint50g_post
52    o merged with release1_p17 (pkg/seaice and verification/lab_sea)
53      - added SEAICE_MULTILEVEL for 8-category sea-ice thermodynamics
54      - LSR sea-ice dynamic solver moved to SouthWest B-grid location and
55        made the default because of faster convergence than ADI
56    
57    checkpoint50f_post
58    o exf totally restructured
59      - replaced all exf_set_ routines by one generic exf_set_gen.F
60        and made sure it's still differentiable
61      - split exf_getffields.F into new exf_bulkformulae.F
62      - exf_inscal_/exf_outscal_ suffixes consistent with field names
63      - updated adjoint/makefile and relevant verif. accordingly
64    
65    checkpoint50f_pre
66     o AIM: near surface temp evaluated from Pot.Temp: fixed for Part Cell ;
67          + change Evap over land. (update aim.5l_cs output)
68     o added missing .cvsignores and made uniform across expts
69     o changed configuration of tidal_basin_2d expt
70     o added pkg/debug functions for tracing code in fatal crashes
71     o tidied pkg/monitor to use MONITOR_OPTIONS.h
72     o new fn in pkg/monitor to stop model if solution grows too big
73     o removed lots of single quotes from comments
74    
75    checkpoint50e_post
76    o make KPP work with PTRACERS
77     - fix gad_calc_rhs to call new routine kpp_transport_ptr, which is
78       nearly a copy of kpp_transport_s
79     - there is no analogue to SurfaceTendencyS, so I have to use
80       gPtr(of the surface layer) instead
81    o add a new platform SunFire+mpi (SunFire 15000) to genmake
82    checkpoint50e_pre
83    
84    checkpoint50d_post
85    o change kpp output from multiple-record state files to single-record state
86      files analogous to write_state.F
87    o reduce the output frequency of cg3d-related stuff to the monitor frequency,
88      analogous to the cg2d-related output.
89    o fix small problem with in ptracers_write_checkpoint.F: len(suff)=512,
90      so that writing to internal file fn (with length 512) fails.
91    checkpoint50d_pre
92    
93    checkpoint50c_post
94    o fix problem in ini_procs.F when using MPI (was there since checkpoint48f)
95    o add monitor diagnostic of vorticity (mon_vort3.F) and potential energy.
96    o add time-average diagnostic of U*V and vertical diffusive flux for theta.
97    o bug fixed in convective adjustment diagnostic (multi-tile set-up).
98    o 4th Order Adv scheme: return to previous (c40pre7) masking near boundaries.
99    
100    checkpoint50c_pre
101    Merging from release1_p13:
102    o bug fix for pkg/seaice dynamic solver
103    o Added SEAICE_initialHEFF to pkg/seaice
104    
105    checkpoint50b_post
106    o store u*,v* in gU,V instead of in guNm1,gvNm1
107    o in order to put the momForcing out of the Adams-Bashforth:
108      move forcing & CD-scheme calls from mom_fluxform & mom_vecinv
109      to timestep.F
110    o new flag "useCDscheme" (default=F); replace guCD,gvCD by local arrays
111    o re-generate output for highly sensitive test-exp:
112       global_ocean_pressure, ideal_2D_oce, lab_sea, natl_box
113    
114    checkpoint50b_pre
115    o restore the use of momAdvection & useCoriolis flags (mom_vecinv)
116    o corrected CPP option mismatch in exf_Bulk... functions for case
117      #define ALLOW_ATM_WIND && #undef ALLOW_ATM_TEMP
118    
119    checkpoint50a_post
120    o Bug fix for merging between c50 and KPP.
121      ikey was passed from thermodynamics to kpp_calc via
122      common block rather than being recomputed in kpp_calc,
123      in contradiction with new key itdkey.
124      New key ikppkey created, and tamc.h headers updated.
125    
126    checkpoint50
127    - 2nd step of merge with ecco-branch:
128      merge changes in ecco-branch back into main trunk
129    - bug fix in find_alpha
130    - more storing in gmredi_calc_tensor.F, kpp_routines.F
131    - avoid iterative loop in ini_pressure for ALLOW_AUTODIFF_TAMC
132    - merged some sliced AD support routines
133    
134    checkpoint49
135    merged from ecco-branch:
136    o exf:
137      - Enable initialisation of forcing fields to constant
138        (runtime) values.
139      - in exf_getffields.F
140        Reduce i-/j-loop to interior domain, discarding overlaps.
141        That also fixes wrong TAF-key computations for key_1, key_2
142        with bulf formulae.
143      - exf_init.F modify #ifdef for exf_init_evap
144      - exf_getffieldrec.F, ctrl_getrec.F
145        The following INT-usages are not safe:
146          fldsecs  = int(fldsecs/fldperiod)*fldperiod
147          fldcount = int(fldsecs/fldperiod) + 1
148        and were modified.
149    o autodiff: checkpoint_lev...
150      - updated storing for bulk formulae (some arrays were missing)
151      - added state for seaice package
152      - synchronized with main branch
153    o DIVA:
154      Modif.s in
155      - the_main_loop.F
156      - the_model_main.F
157      and added cost_final_restore.F cost_final_store.F
158      o With these changes, it is possible to fully separate
159        the first forward run to compute the cost function (mdthe_main_loop)
160        and adjoint runs over individual DIVA intervals (adthe_main_loop)
161        into separate model executions.
162        This enables better tuning of DIVA to given queue sizes.
163      o Loss of cost function value in consecutive DIVA runs is
164        fixed by storing to 'costfinal' file.
165      o Appropriate call of ctrl_unpack/_pack,grdchk_main at
166        beginning/end of adjoint integration is accounted for.
167      The current changes apply to TAF-1.4.26.
168      TAF has one remaining bug, requiring one manual change
169      (see MITgcm manual).
170    o cal: synchronized between MAIN and ecco-branch
171    
172    checkpoint48i_post
173    o Changes to restore differentiability of code w.r.t. previous tag
174      (mostly adding new routines to make list and replacing
175      pressure by totPhiHyd).
176    o Updated and merged grdchk package
177      - has both ADM and TLM checks
178      - works for single- and multi-proc.
179      - output cleaned
180      - worked successfully for parallel DIVA
181    o genmake:
182      - Adding -platform=SP4
183      - Modif. for DEC Alpha (many don't have makedepend; use mkdep instead).
184    o Modif. for DIVA:
185      ctrl_unpack should not be called in intermediate DIVA interval.
186      (similar thing would apply for active files such as
187      tbar, sbar, psbar, which should not be reinitialised;
188      only relevant for ecco-branch)
189    
190    checkpoint48h_post
191    o improve "ini_p_ground" when using finite difference form for calc_phi_hyd
192    o Added net flux to downward flux conversion to pkg/exf/exf_getffields.F
193    o Added SEAICE_initialHEFF to pkg/seaice
194    o ecco_the_model_main.F: SEAICE_WRITE_PICKUP must precede WRITE_CHECKPOINT
195    
196    checkpoint48g_post
197    --- EOS:
198    o compute locally the pressure for use in EOS: UNESCO, JMD95P, MDJWF
199      store total Potential in totPhyHyd for diagnostics & EOS.
200    o fix restart and overlap Pb when using Z-coord and EOS funct. of P
201    --- Pressure/geopotential gradient:
202    o use Ro_surf & R_low instead of hFac in cal_phi_hyd:
203    o NLFS + p-coord.: eta*Alpha' is not longer lagging 1.time-step behind
204      ( changes results of global_ocean_pressure )
205    o NLFS + z-coord.: add eta*Rho' contribution.
206    o accurate phiHyd when using interface-W at the middle between 2 centers
207     ( changes results of ideal_2D_oce )
208    o includes r* 2nd term in potential gradient.
209    
210    checkpoint48f_post
211    Merging from release1_p12:
212    o Modifications for using pkg/exf with pkg/seaice
213      - improved description of the various forcing configurations
214      - added basic radiation bulk formulae to pkg/exf
215      - units/sign fix for evap computation in exf_getffields.F
216      - updated verification/global_with_exf/results/output.txt
217    o Added pkg/sbo for computing IERS Special Bureau for the Oceans
218      (SBO) core products, including oceanic mass, center-of-mass,
219      angular, and bottom pressure (see pkg/sbo/README.sbo).
220    o Lower bound for viscosity/diffusivity in pkg/kpp/kpp_routines.F
221      to avoid negative values in shallow regions.
222      - updated verification/natl_box/results/output.txt
223      - updated verification/lab_sea/results/output.txt
224    o MPI gather, scatter: eesupp/src/gather_2d.F and scatter_2d.F
225    o Added useSingleCpuIO option (see PARAMS.h).
226    o Updated useSingleCpuIO option in mdsio_writefield.F to
227      work with multi-field files, e.g., for single-file pickup.
228    o pkg/seaice:
229      - bug fix in growth.F: QNET for no shortwave case
230      - added HeffFile for specifying initial sea-ice thickness
231      - changed SEAICE_EXTERNAL_FLUXES wind stress implementation
232    o Added missing /* */ to CPP comments in pkg/seaice, pkg/exf,
233      kpp_transport_t.F, forward_step.F, and the_main_loop.F
234    o pkg/seaice:
235      - adjoint-friendly modifications
236      - added a SEAICE_WRITE_PICKUP at end of the_model_main.F
237    
238    checkpoint48e_post
239    o New S/R for diagnostic of bottom pressure (phi0surf was missing in c48d_post)
240    o dynamics: change definition of computational domain & adapt mom_fluxform
241      accordingly ; when viscA4=0, allows to run the dynamics with Olx=Oly=2.
242    o following exp works with Olx=Oly=2 (and pass through checkbounds option):
243        adjust_nlfs.cs-32x32x1 adjustment.128x64x1 adjustment.cs-32x32x1
244        exp0 exp1 exp2 exp4 exp5 hs94.cs-32x32x5 hs94.128x64x5 hs94.1x64x5
245        global_ocean.90x40x15 global_with_exf global_ocean_pressure
246        inverted_barometer solid-body.cs-32x32x1
247      (high order AdvScheme and/or multiDimAd works with Olx=Oly=2 ? not sure)
248    
249    checkpoint48d_post
250    o in preparation for r*:
251     a) use pre-computed gradient of hydrostatic potential:
252        changes in timestep.F & mom_cdscheme.F affects results of ideal_2D_oce
253     b) move phi0surf from calc_phi_hyd to calc_grad_phi_hyd:
254        => affects results of glob_oce_pressure (different truncation error)  
255      update output.txt of ideal_2D_oce & global_ocean_pressure
256    
257    checkpoint48d_pre
258    o in preparation for r*:
259      new S/R (calc_grad_phi_hyd.F) to compute Hydrostatic potential gradient.
260      pass the 2 comp. of the grad. as arguments to momentum S/R.
261      but for the moment, only used if it does not change the results.
262    o add finite volume form of calc_phi_hyd for buoyancyRelation OCEANIC & OCEANICP
263    
264    checkpoint48c_post
265    Modifications in treatment of initial flux forcing controls
266    to avoid extra storing of flux fields.
267    
268    checkpoint48c_pre
269    o Added tidal_basin_2d experiment (currently 1D!)
270    o can use r* Coordinate with OBC (but Orlanski Rad. BC not yet implemented)
271    
272    checkpoint48b_post
273    o r* coordinate added in #ifdef NONLIN_FRSURF block.
274     (modification to pressure gradient not yet implemented)
275    o mom_fluxform: new S/R mom_calc_rtrans.F computes vertical transports
276      (for U & V) that are used for vertical advection of momentum.
277    
278    checkpoint48a_post
279    o In preparation for r*, use temp. array to compute vertical transport
280       in flux-form momentum pkg. => change results of few test-exp.
281      Update output of: exp4, exp5, front_relax, global_ocean_pressure,
282                        lab_sea, natl_box & plume_on_slope
283    
284    checkpoint48
285    o fix small bug that was causing lot of warning (TIMER_STOP 'UPDATE_CG2D')
286    
287    checkpoint47j_post
288    o add new exp. (ideal_2D_oce) to test "exotic" parameters and options.
289    o moved adexch from forward_step to addummy_in_stepping
290    o added exch in ctrl_map_...
291    o aveFreq.=0. at end of cost_final to avoid overwrite in
292      reverse checkpointing loop
293    o Yet more changes in gmredi:
294      * adgmredi_calc_tensor
295        avoiding all recomputation of gmredi_slope_limit
296      * adgmredi_x/y/rtransport
297        added flag for excessive storing to avoid recomp. of
298        u/v/rtans, dTdx/y/z
299        -> this is not really necessary and very memory-consuming
300      * adgmredi_slope_psi:
301        consistency with gmredi_slope_limit in treatment of GM_slopeSqCutoff
302      * gmredi_slope_limit
303        re-activated full calculation of taperfct for case 'ac02'
304    o updating verification/ adjoints of carbon, global_ocean.90x40x15
305    
306    checkpoint47i_post
307    o GM_EXCLUDE_() replace GM_TAPER_(): by default, all the code is compiled.
308    o add GM_Small_Number & GM_slopeSqCutoff as run-time parameters (data.gmredi)
309      default values = 1.D-12 & 1.D+48 (=> give same results as checkpoint47f_post)
310    o cg2dTargetResWunit: define tolerance using deltaTmom (and not deltaTfreesurf)
311    o gmredi: fix few bugs.
312    
313    checkpoint47h_post (duplicate tag for PH's c48 tag, in order to move c48)
314    This will be a good checkpoint to plug into the ecco-branch.
315    o few modif.'s to get latest version adjointed
316      (mainly kick out code in ini_linear_phisurf)
317    o modif's to run adjoint with exactConserv
318    o case GM_BOLUS_ADVEC should be cleaned
319      S/R gmredi_slope_psi should be cleaned
320    o verification/carbon now has exactConserv=.TRUE.
321    
322    checkpoint47g_post
323    o time-average AIM output: becomes independent from statvars time-average
324      (can be used with #undef ALLOW_TIMEAVE) ; add new variables (e.g., RH)
325    o After TAF bug fix:
326      Tested adjoint of gmredi package.
327      Tests were done with checkpoint47 rather than with latest
328      checkpoint47f_post since the latter is broken for the adjoint
329      verification/carbon/ contains both ADM and TLM gradient checks.
330    Here they are, before they get changed/lost/stolen.
331    Mostly modif.'s to fix numerical sensitivities.
332    Gradient checks OK for
333    - GM_taper_scheme:
334      * clipping
335      * ac02
336      * linear
337      * glw91
338      * dm95
339      * ldd97
340    - GMREDI_OPTIONS:
341      * GM_VISBECK_VARIABLE_K
342      * GM_NON_UNITY_DIAGONAL
343      * GM_EXTRA_DIAGONAL
344      * GM_BOLUS_ADVEC
345    in conjunction with data.gmredi parameters to be checked in
346    in a few minutes under verification/carbon/code/
347    
348    checkpoint47f_post
349    Merging from release1_p10:
350    o modifications for using pkg/exf with pkg/seaice
351      - pkg/seaice CPP options SEAICE_EXTERNAL_FORCING
352        and SEAICE_EXTERNAL_FLUXES
353      - pkg/exf CPP options EXF_READ_EVAP and
354        EXF_NO_BULK_COMPUTATIONS
355      - usage examples are Experiments 8 and 9 in
356        verification/lab_sea/README
357      - verification/lab_sea default experiment now uses
358        pkg/gmredi, pkg/kpp, pkg/seaice, and pkg/exf
359    
360    checkpoint47e_post
361    o created a tag prior to creating a branch
362    
363    checkpoint47d_post
364    o fixed bug in the oceanic pressure coordinates code: vertical viscosity
365      at the bottom boundary had an erroneous half slip boundary condition
366    o added bottom drag and no slip boundary condition capabilities to
367      oceanic pressure coordinates code
368    o added new verification experiment global_ocean_pressure, this experiment is
369      also described in the tutorial section of the manual
370    checkpoint47d_pre
371    
372    checkpoint47c_post
373    o Added a verification experiment for a global cubed-sphere ocean run
374      -  verification/global_ocean.cubed32x32x30
375      -  uses pkg/therm_seaice and pkg/bulk_forcing
376    
377    o allow to run AIM physics with SPEEDY input files (from Franco Molteni)
378    o allow a more accurate definition of Ro_Surf (selectFindRoSurf=1)
379      when using P-coordinate; only implemented for atmospheric config.
380    o OCEANICP & realFreshWater: include P-E direct effect on wVel ;
381      NOTES: requires option NONLIN_FRSURF to be "#define".
382    o update advect_xz/results/output.txt (left from checkpoint44g_post)
383    
384    checkpoint47b_post
385    Merging from release1_p9:
386    o pkg/seaice
387      - removed GOTO's and added taf directives
388      - double precision constants to reduce the g77 (Linux)
389        to F77 (SGI) differences reported in release1_p8
390    o tools/genmake
391      - added SGI options
392    o verification/testscript
393      - updated to that of checkpoint47a_post
394    o verification/global_ocean.90x40x15/input/eedata
395      - modified for SGI f77 compatibility
396    o verification/lab_sea
397      - added description of sea-ice model
398      - added missing matlab routines
399      - added test of thermodynamics parallelization
400    
401    checkpoint47a_post
402    o new pkg: aim_v23
403     = F. Molteni atmos.physics (SPEEDY, ver23) adapted to MITgcm.
404       for now, keep the same surface forcing as before.
405     - Part-Cell implemented into AIM; check that Heat & Water are conserved
406     - aim.5l_cs: use new aim pkg (run 10yr & get better results than before)
407     - aim.5l_LatLon & aim.5l_Equatorial_Channel are still using old aim pkg
408    o add diagnostic of surface correction term in monitor output
409    
410    o bulk_force and therm_seaice
411    - Two packages:  bulk_force (Bulk forcing)
412      and therm_seaice (thermodynamic_seaice) - adopted from LANL CICE.v2.0.2
413    - Earlier integration from Stephaine Dutkiewicz
414      and Patrick Heimbach.
415    - Two ifdef statements for compile time,
416       ALLOW_THERM_SEAICE and ALLOW_BULK_FORCE
417    
418    checkpoint47
419    differentiable version of checkpoint46n_post
420    o bug fix in quasihydrostaticterms
421    o * "clean" adjoint code (in terms of extensive recomputations)
422        can now be obtained for all GMREDI options (i.e. for
423        - GM_VISBECK_VARIABLE_K
424        - GM_NON_UNITY_DIAGONAL
425        - GM_EXTRA_DIAGONAL
426        - GM_BOLUS_ADVEC )
427      * However, wrong gradient check problem remains unsolved.
428      * New CPP options have been introduced for different
429        tapering schemes
430    o external_fields_load now part of differentiation list
431    o pressure needs multiple storing;
432      would be nice to have store_pressure at beginning or
433      end of forward_step, e.g. by having phiHyd global (5-dim.)
434      (NB: pressure is needed for certain cases in find_rho,
435      which is also invoked through convective_adjustment).
436    o recomputations in find_rho for cases
437     'JMD95'/'UNESCO' or 'MDJWF' are OK.
438    o #define ATMOSPHERIC_LOADING should be differentiable
439    o ini_forcing shifted to begining of initialise_varia
440    o Incorporating QNVS line search routines into MITgcm
441      (this is separate code, and therefore not under pkg)
442      lsopt/, optim/
443    o Updated verification/carbon/
444      - #define GM_VISBECK_VARIABLE_K
445      - GM_taper_scheme    = 'ldd97'    
446    
447    checkpoint46n_post
448    Merging from release1_p8:
449    o verification/natl_box:
450      updating new external_fields_load routine
451    o New package: pkg/seaice
452      Sea ice model by D. Menemenlis (JPL) and Jinlun Zhang (Seattle).
453      The sea-ice code is based on Hibler (1979-1980).
454      Two sea-ice dynamic solvers, ADI and LSR, are included.
455      In addition to computing prognostic sea-ice variables and diagnosing
456      the forcing/external data fields that drive the ocean model,
457      SEAICE_MODEL also sets theta to the freezing point under sea-ice.
458      The implied surface heat flux is then stored in variable
459      surfaceTendencyTice, which is needed by KPP package (kpp_calc.F and
460      kpp_transport_t.F) to diagnose surface buoyancy fluxes and for the
461      non-local transport term.  Because this call precedes model
462      thermodynamics, temperature under sea-ice may not be "exactly" at
463      the freezing point by the time theta is dumped or time-averaged.
464      N.B.: fairly large differences in accuracy occur across
465            different platforms/compilers; comparison between
466            g77 (Linux) and F77 (SGI Origin 2000/MIPSpro)
467            generated output gives:
468       T           S           U           V
469      C D M    c        m  s        m  s        m  s        m  s
470      n p a R  g  m  m  e  .  m  m  e  .  m  m  e  .  m  m  e  .
471      f n k u  2  i  a  a  d  i  a  a  d  i  a  a  d  i  a  a  d
472      g d e n  d  n  x  n  .  n  x  n  .  n  x  n  .  n  x  n  .
473    
474      Y Y Y Y  5  5  7  7  7  8 10  9  6  6  6  6  7  5  7  5  7 FAIL  lab_sea
475    o GAD:
476      - generated new common blocks to account for call of
477        same gad routines with differing traceridentities
478        (needed to modify tracerIdentity indices in GAD.h)
479      - generated separate common blocks for case useCubedSphereExchange
480        (Department of Futurology)
481      - parameter lists to gmredi_?transport: added tracerIdentity
482      - added new key indices to tamc.h
483    o external_fields_load:
484      - added this routine to TAF list
485      - needed to make some common blocks global and additional storing
486        along the same lines as exf package (checkpoint_lev?_directives.h)
487    o exf:
488      updated external forcing package
489      - bug fixes carried over from ecco-branch
490        (missing OBCS_OPTIONS.h in two routines)
491      - enable easy to use "no forcing".
492      - added exf I/O for atmospheric loading
493      - added exf I/O for runoff data
494      - transfered scaling between exf <-> MITgcm to exf namelist
495    o Adding new verification experiment global_with_exf.
496      This verification experiment is almost identical to the
497      global_ocean.90x40x15 experiment, except that it uses
498      the exf (and cal) package instead of the
499      external_fields_load routine (cf README file).
500    
501    checkpoint46m_post
502    o modified testscript to return -- when there are no lines of "%MON" output
503    o added new flag "quasiHydrostatic" for QH mode
504      and new routine quasihydrostaticterms()
505    o added 2*Omega*Cos(phi)*W to u equations (non-hydrostatic Coriolis term)
506      New routine: mom_u_coriolis_nh()
507    o fixed sign of non-hydrostatic metric terms in u and v equations.
508      Depend on correct sign of gravitySign.
509    o added work.pc* to .cvsignore in aim directories using ifc
510    o introduced flag for controllin non-hydrostatic metric terms
511    o removed gravitySign from PARM04
512    o switched to g77/ifc from pgf77 in aim experiments
513    
514    checkpoint46l_post
515    o made convective adjustment work with pressure coordinates:
516     - changed the direction of k-loop in convective_adjustment.F for the
517       case of pressure coordinates (OCEANICP,ATMOSPHERIC buoyancyRelation)
518     - adjusted the reference pressure k-index in convective_adjustment.F
519     - adjusted the convection condition in convect.F (in analogy to
520       calc_ivdc.F)
521     - convective_adjustment no longer computes anything on the halos
522     - removed the warnings about negative salinity from find_rho.F and
523       find_alpha.F; instead the new routine look_for_neg_salinity, called
524       at the beginning of find_rho, find_alpha, and find_beta, does a
525       check of the entire slice, if CPP-option
526       CHECK_SALINITY_FOR_NEGATIVE_VALUES is defined
527    checkpoint46l_pre
528    
529    checkpoint46k_post
530    o fixed the verification/global_ocean.90x40x15 experiment:
531     - new bathymetry (the world according to A., JMC, and M.)
532     - new initial fields and forcing fields (*.bin files)
533     - new POLY3.COEFFS (for the next release one should switch to a full
534       equation of state: JMD95P or MDJWF)
535     - fixed several errors and redundancies in the data file
536     - experiment uses looped cells
537     - added matlab directory with diagnostic scripts for plotting of output
538    
539    o S/R aim_initialise.F replace S/R aim_init from file aim_do_inphys.F:
540     - read AIM physics parameters from a file (data.aimphys)
541     - set defaults values = F.Molteni paper (Clim.Dyn., 2002)
542    
543    checkpoint46j_post
544    o split calc_exact_eta in 2 S/R : integr_continuity & update_etaH
545    o move wVel computation at the end of the time step, in S/R integr_continuity
546    o create specific S/R to exchange T,S before DYNAMICS (for stagger time step)
547    o update timeave pkg for wVel diagnostic.
548    checkpoint46j_pre
549    
550    checkpoint46i_post
551    o Clean up AIM package (and keep the results unchanged):
552      include CPP_OPTION and use IMPLICT NONE in all routines ;
553      declare all the variables _RL ; use _d 0 for all numerical constants.
554      use ifdef ALLOW_AIM everywhere. And now AIM can be used with g77 !
555    
556    checkpoint46h_post
557    o cleaned up the use of rhoNil and rhoConst.
558      - rhoNil should only appear in the LINEAR equation of state, everywhere
559        else rhoNil is replaced by rhoConst, e.g. find_rho computes rho-rhoConst
560        and the dynamical equations are all divided by rhoConst
561    o introduced new parameter rhoConstFresh, a reference density of fresh
562      water, to remove the fresh water flux's dependence on rhoNil. The default
563      value is 999.8 kg/m^3
564    o cleanup up external_forcing.F and external_forcing_surf.F
565      - can now be used by both OCEANIC and OCEANICP
566    checkpoint46h_pre
567    
568     Added code to convert surface volume flux (fresh water) into
569     a mass flux when using P coordinates in the ocean (OCEANICP).
570     Note: It assumes you have set rho0=rhoConst=density of fresh water.
571    
572    checkpoint46g_post
573    o Include a new diagnostic variable phiHydLow for the ocean model
574      - in z-coordinates, it is the bottom pressure anomaly
575      - in p-coordinates, it is the sea surface elevation
576      - in both cases, these variable have global drift, reflecting the mass
577        drift in z-coordinates and the volume drift in p-coordinates
578      - included time averaging for phiHydLow, be aware of the drift!
579    o depth-dependent computation of Bo_surf for pressure coordinates
580      in the ocean (buoyancyRelation='OCEANICP')
581      - requires a new routine (FIND_RHO_SCALAR) to compute density with only
582        Theta, Salinity, and Pressure in the parameter list. This routine is
583        presently contained in find_rho.F. This routine does not give the
584        correct density for 'POLY3', which would be a z-dependent reference
585        density.
586    o cleaned up find_rho
587      - removed obsolete 'eqn' from the parameter list.
588    o added two new verification experiments: gop and goz
589      (4x4 degree global ocean, 15 layers in pressure and height coordinates)
590    checkpoint46g_pre
591    
592    checkpoint46f_post
593    o Enable tangent linear (forward mode) gradient checks:
594      - extended active file handling to g_... files
595      - added TANGENT_SIMULATION to theSimulationMode
596      - extended grdchk package accordingly
597    o added tangent linear setup in verification/carbon/code_ftl
598    o added adjoint setup in verification/front_relax/code_ad
599    
600    checkpoint46e_post
601    o Added new equation of state -> MDJWF
602      - EOS of McDougall et al., 2002, JAOT, submitted
603      - caveat: the equation of state is only valid for a smaller (more
604        realistic?) range of values than JMD95P/Z and UNESCO
605      - added masks to the calculation of pressure in store_pressure
606      - added more check values for density in check_eos (ini_eos.F), some of
607        the old check values are out of the range of the MDJWF-eos, so don't
608        expect perfect matches for those
609    checkpoint46e_pre
610    
611    checkpoint46d_post
612    o fixed store_pressure to work with both buoyancy relation = 'OCEANIC' and
613      'OCEANICP', also initialised field pressure correctly in ini_eos in the
614      case of pressure coordinates. eosType='JMD95Z' in
615      combination with buoyancyRelation='OCEANICP' now causes an error.
616    o Changed p = pressure(i,j,k,bi,bj)  to  p = pressure(i,j,kRef,bi,bj)
617      in find_alpha/beta.
618    checkpoint46d_pre
619    
620    o Changed p = pressure(i,j,k,bi,bj)  to  p = pressure(i,j,kRef,bi,bj)
621      so that JMD95Z and JMD95P give approptiate static stability in find_rho.
622    
623    checkpoint46c_post
624    o Changes necessary for ocean in p-coordinates
625     -  Added new buoyancy relation = 'OCEANICP'
626     -  Added new parameters = gravitySign (this used to be contained inside
627        the factor dRdZ which I added when we first switched to R coordinates).
628     X GM/Redi is not compatible (yet)
629     X bottom drag and no-slip need to be debugged.
630    checkpoint46c_pre
631    
632    checkpoint46b_post
633    o Added new equation of state -> JMD95Z and JMD95P
634      - EOS of Jackett and McDougall, 1995, JAOT
635      - moved all EOS parameters into EOS.h
636      - new routines ini_eos.F, store_pressure.F
637    o Added UNESCO EOS, but not recommended because it requires
638      in-situ temperature (see JMD95)
639    o Modified formatting for knudsen2.f in utils/knudsen2 and added
640      unesco.f to be used with POLY3
641    checkpoint46b_pre
642    
643    checkpoint46a_post
644    o Added (atmospheric) pressure loading
645      - new field in FFIELDS.h, etc...
646      - new cpp flag ATMOSPHERIC_LOADING
647    o Changed hFacC to _hFacC in calc_phi_hyd.F
648    o Added SHORTWAVE_HEATING to some files for consistency
649    checkpoint46a_pre
650    
651    checkpoint46
652    further merges from release1_p5
653    o new ctrl package
654      - adopted from ECCO environment to enable optimization
655      - added Eliassen Palm fluxes to controls
656    o added Eliassen Palm flux controls to gradient check package
657    o cost package
658      - Compatible with new ctrl package
659      - added Eliassen Palm cost hooks
660      - modif's of existing cost functions
661    o Adjoint-related bug fixes in kpp:
662      - kpp_calc: sore of kpphbl avoids recomputation/call to S/R kppmix
663      - kpp_routines: store of Rib avoids partial recomputation bug of TAF.
664    o autodiff package
665      - flow directives:
666        + new for S/R do_the_model_io
667        + enhanced for S/R checkpoint
668        + new gmredi.flow
669      - Introduced CPP option ALLOW_AUTODIFF_MONITOR to
670        disable adcommon.h/g_common.h by default.
671      - fixed problem for adjoint of global_max
672    o modified verification/carbon/
673    o New setup for adjoint of global_ocean
674    o added hooks for Eliassen Palm fluxes (dynamics.F, ini_fields.F)
675    o cleaned TAF keys iikey, idkey (dynamics.F, thermodynamics.F)
676    o incorporated hooks for sponge layer code (external_forcing.F)
677    
678    checkpoint45d_post
679    o import 2 fixes from release1_p5 (CG2D.h & packages_readparms.F)
680    o use recip_dx*,recip_dy* instead of /dx*,/dy* in orlanski_E,W,N,S
681      (affects plume_on_slope "cg2d" output)
682    o add new parameter deltaTfreesurf for free-surface time stepping
683    o exchange gT,gS when using staggerTimeStep & GM & Oly<4
684    o NONLIN_FRSURF: rescale gNm1 to get a better conservation with A-B
685    o change phiMin in exp1 (to agree with documentation)
686    o add new flags:
687     - T,S forcing inside/outside Adams-Bashforth
688     - turn off temp,salt advection and/or forcing
689     - for each tracer, define internal flag for multiDimAdvection & A-B
690    o monitor: change definition of KE to check conservation of Energy
691    
692    checkpoint45c_post
693    o Removed f90 comments ("!")
694    o Modified initialisations to enable exf package for MITgcm
695      (without invoking ECCO_PACKAGE).
696    o modifications to gradient check package (Martin Losch)
697      - enable centered differences
698      - modified format of standard output
699    
700    checkpoint45b_post
701    o cleaned exchange calls in cg2d
702    o Included CPP option ALLOW_AUTODIFF_MONITOR in forward_step
703    o included CPP option SINGLE_LAYER_MODE in thermodynamics
704      to configure barotropic setup (Martin Losch)
705    o moved some initialisations within bi,bj loop in dynamics
706      (TAF flow dependency)
707    o in initialise_varia, moved call packages_init_variables
708      before convective_adjustment_init
709      (corresponding adjustment of TAF tape initialisation).
710    o in calc_diffusivity/calc_viscosity extended domain to full overlap.
711    
712    checkpoint45a_post
713    o compute Rho over the whole domain (extend iMin,jMin indices)
714    o specific exchange for X-slice domain (case Ny=1)
715    o argument futureIter of S/R obcs_calc missing in exp4 & internal_wave /code
716    o solve_for_pressure : set cg2d_x to zero where OBCS are applied
717    
718    checkpoint44h_post = checkpoint45
719    o Fix to restore differentiability:
720      place at which packages_init_variables is called in
721      initialise_varia is crucial for reverse control flow.
722    o modifications in GMRedi to enable stable adjoint
723      (corresponding modif's in thermodynamics, dynamics)
724    o added missing hook for sponge layer code in external_forcing
725    o modified test cost functions
726    o storing in gad_advection (plus removal of write(0,*)!
727    o missing headers TR1.h added in convective_adjustment routines
728    o errorMessageUnit set non-equal zero in eeboot
729      (conflict on some platforms, e.g. IBM SP3)
730    o modified carbon verif.
731    
732    checkpoint44h_pre
733    o vertical grid option: allow to put the Interface at the middle between
734      2 cell-centers ; replace delZ in pkg/kpp by drF.
735    o GM Advective form: Tracers are advected using the residual transport (=
736      Euler+GM-bolus); set param GM_AdvSeparate=T to return to previous form.
737    
738    checkpoint44g_post
739    o fix surface correction term in multiDimAdvection (affects plume_on_slope
740      results) ; ifdef missing in thermodynamics.F
741    
742    checkpoint44f_pre,checkpoint44f_post
743    o added PTRACERS package
744      This allows an arbitrary number of passive tracers to be integrated
745      forward simultaneously with the dynamicaly model.
746      + Implemented so far:
747        - basic forward algorithm (time-stepping, advection, diffusion, convection)
748        - I/O and checkpointing
749        - GM/Redi  *but*  using the GM/Redi coefficient of Salt
750      + Not implemented so far:
751        - KPP
752        - OBCS
753      + No specific example supplied (yet) but global_ocean.90x40x15 has the
754        necessary data.ptracer file. Simply use -enable=ptracers and uncomment
755        line in data.pkg.
756      + This package is disabled by default since it increases storage.
757    
758    
759    o testscript: added -cleanup option
760    o Shap_filt: a) new shap_filt S/R to use no-slip BC with S2 filter
761                 b) enable to filter 2D fields.
762    
763    checkpoint44e_pre,post
764    o re-wrote convective adjustment to use pre-calculated weights
765      Affects: model/src/convective_adjustment.F, convective_adjustment_ini.F
766      Adds:    model/src/convective_weights.F, convectively_mixtracer.F
767      Original calls to convect() still in place.
768    
769    chkpt44d_pre,post
770    o added missing EXCLUDE_MONITOR flags
771    o changed "e" to "_d" in gmredi_slope_limit, gmredi_slope_psi
772      (incompatible typ in MIN/MAX expressions caused problems
773      on IBM SP3)
774    o in genmake added variable MAKEDEPEND
775      plus resetting for case SunOS
776    o added timer_stats.c routine for IBM SP3
777    o removed variables in dynamics
778    !!! I forgot to create tag chkpt44d_pre !!!
779    !!! But chkpt44d_pre = chkpt44c_post    !!!
780    
781    o real fresh water flux implemented with non-linear free-surface.
782    
783    chkpt44c_pre,post
784    o few fix (mask in shap_s2, EmPmR in external_field_load,
785      USE_NATURAL_BCS in solve_for_P);
786    o add arguments myIter & myTime to S/R obcs_calc & solve_for_P
787    
788    checkpoint44b_pre/_post
789    o merge of relevant stuff from the ecco-branch:
790      - genmake: removed $S64 overwrite for case SunOS
791      - pkg/exf: update and corrections for field swapping and obcs
792      - pkg/ecco: parameter lists for the_model_main, the_main_loop
793                  harmonized between ECCO and MITgcm
794      - pkg/autodiff: added flow directives for obcs, mdsio_gl_slice
795                      updated checkpointing_lev... lists for obcs
796      - model/src: minor changes in forward_step, plot_field
797                   added directive for divided adjoint in the_main_loop
798      - pkg/mdsio: added mdsio_gl_slice
799    
800    o check parameters & config (chkpt44a_pre,post)
801    o OBC and NonLin_FrSurf.
802    
803    checkpoint44
804    o fix bug in mom_vi_del2uv
805    o select when filters are applied ; add options to zonal_filter (data.zonfilt)
806    o gmredi: fix Pb in the adiabatic form ; add options (.e.g. Bolus advection)
807    o update AIM experiments (NCEP input files)
808    o improve and extend diagnostics (Monitor, TimeAve with NonLin-FrSurf)
809    o added some stuff for AD
810    
811    checkpoint43
812    o added two new advection tests
813    o added global 4x4 degree ocean test
814    o added test of H&S94 with free surface on cube
815    o some final mods for gradient check, adjoint,...
816    
817    checkpoint42
818      Fixed AD-related problems:
819    o Store directives up-to-date with re-arranged Adams-Bashforth
820      (mainly thermodynamics.F)
821    o New store directives for multi-dim. advection schemes
822      * new CPP flag ALLOW_MULTI_DIM_ADVECTION
823      * new common block and key passkey
824      (mainly gad_advection.F)
825    o Modified store directives for split of dynamics/thermodynamics
826      for the case ALLOW_KPP
827    o Cleaned argument list for timestep_tracer.F
828    
829    checkpoint41
830    o re-formatted/added comments in prologues of key subroutines
831      to work with protex
832    
833    checkpoint40
834    o moved calc_gs, calc_gt and calc_tr1 back to model/src
835    o added FLT (floats) package from Arne along with demo. experiment
836    o re-arranged Adams-Bashforth scheme for tracers to allow easier
837      implementation of forward-in-time forcing terms
838    o more mods for non-linear free-surface including new variable
839    o modified multi-dim method to work on cube
840    o bug fix in KPP (from DM)
841    
842    checkpoint40pre9
843    o added DST advection schemes
844    o fix MPI bug on SGI's
845    o modifications for non-linear free-surface
846    o fixed conditionals for fixed flow/passive tracer only mode
847    o added tension/strain form of stress tensor
848    o changes to commenting of arguments for TAF/TAMC
849    o added/updated many matlab scripts
850    o added multi-dimensional advection method for tracers
851      
852    checkpoint40pre8
853    o new package rw/
854      Added method for dumping intermediate local arrays:
855      mdsio_writetile - same as mdsio_writefield except works from inside
856                        bi,bj loop
857      mdsio_writelocal - same as mdsio_writetile except works for local
858                         arrays
859      write_local_r? - higher-level wrapper for mdsio_writelocal
860      Controlled by diagFreq. Defaults to zero (ie. no dumps)
861      Example given at end of mom_vecinv.F that dumps some local arrays.
862    o replaced calc_mom_rhs.F by mom_fluxform.F/mom_vecinc.F to
863      control flux-form versus vector invariant code.
864      Switch is runtime flag vectorInvariantMomentum (defaults to false)
865      or genmake -DISABLE ...
866    o TAMC-related update to cope with mom_fluxform.F/mom_vecinc.F
867      (additional storing is now needed for case vecinv;
868      currently nothing implemented to switch this storing off
869      in case fluxform is used since no CPP flag available).
870    
871    checkpoint40pre7
872    o bug fix in dynamics.F
873      calc_viscosity called before k=1,Nr loop.
874    o added some #ifdef's ALLOW_PASSIVE_TRACER
875    o cleaned initialisations in ctrl_map_ini
876      related to TAF
877    
878    checkpoint40pre6
879    o split dynamics into thermodynamics & dynamics
880    o corresponding split of calc_diffusivity into
881      calc_diffusivity, calc_viscosity
882      (kpp_calc_diff, kpp_calc_visc).
883    o Added diffkr, kapgm to control variables
884    o bug fix in gmredi_slope_limit for ldd97
885    
886    checkpoint40pre5
887    o proper initialization of Adams-Bashforth time stepping
888      + uses forward step for very first step of model (iter=0)
889    o re-generated *all* output since almost all output is modified by this
890      change (sorry but we've been putting this one off for about 2 years)
891      + natl_box, aim.5l_Equatorial_Channel and aim.5l_zon-ave were
892        *not* regenerated since they are already out of date and I don't
893        want to create the impression that they are suddenly working.
894    o removed KAP stuff for DEC's
895    o modified/added many .cvsignore files in */code and */input
896    o added new expt "solid-body.cs-32x32x1"
897      + this uses an alternative form of ini_curvilinear_grid.F which
898        is not necessarily any better than the current one...
899    
900    checkpoint40pre4
901    o added re-initialisations of utrans,vtrans,rtrans in
902      calc_common_factors which are crucial for TAF.
903    o extended iMin,jMin range for calc_common_factors,
904      calc_diffusivity
905    o added option ALLOW_PASSIVE_TRACER
906    o a few minor updates and fixes
907    
908    checkpoint40pre3
909    o created calc_gtr1.F for gad package
910    o split the_main_loop and re-introduced forward_step
911    o some recomputations of nIter0 for TAMC
912    
913    checkpoint40pre2
914    o Updated to latest branch of AJA, JMC
915    o Added gradient check package
916    o Added routines for passive tracer handling
917    o Added carbon sequestration sensitivity experiment
918    
919    checkpoint39
920     o Modifications and fixes to support TAMC
921       (mainly missing or wrong directives)
922     o Corrected common blocks when fooling TAMC to produce adjoint I/O
923       in addummy_in_stepping
924     o Updated verification/vero/ (bug fixes in code/, correct pickups in input/)
925     o Enabled "sliced" (xz-,yz-) exchanges, needed for obcs.
926     o Included following packages:
927       cal: calendar
928       exf: external forcing
929       ecco: essential routines to control use of cal, exf
930             (NB: this is not the full ECCO environment)
931     o disabled some packages by default in genmake
932    
933    checkpoint38
934     o see notes_c37_adj.txt
935     o corrected wrong layer index in external_forcing_t
936     o replaced the_main_loop/forward_step by the_main_loop
937     o removed forceoutput from write_state
938       (replaced criteria for initial/final state)
939     o added runtime flag for pre-C35 pickups: usePickupBeforeC35
940    
941    checkpoint37
942     o mdsio.F routine split and packaged in pkg/mdsio
943       (mdsio.F and mdsio_gl.F still exist in eesupp/src but are not
944        used unless genmake -disable=mdsio)
945     o allow a more precise surface-geopotential gradient for the atmosphere
946       taking into account the change in space of Alpha (=1/rho) at the surface.
947     o variable "etaN" (replace cg2d_x) is the surface R-anomaly [r unit]
948       (ocean: sea surface height [m] / atmos: surface pressure anomaly [Pa])
949     o Change units: Phi_Hydrostatic (=phiHyd), Phi_Surface(=cg2d_x) and
950       Phi_NonHydrostatic (=cg3d_x), all have now the dimension of a potential
951       [m2/s2] (pressure/rhoConst in the ocean ; Geo-Potential in the atmosphere).
952     o Time-Average diagnostics defined as a pseudo package "timeave".
953     o pickup files : old pickup (before c35) did not work with new code (c35
954       and after) unless activating 2 commented lines in initialise-varia.
955     o Non-Hydrostatic pickup file Pb fixed.
956     o bug with Rigid-lid version (since c35) fixed.  
957    
958    checkpoint36
959     o implement Crank-Nickelson barotropic time-stepping
960     o recover Time-Average diagnostics
961     o fix Pb in checkpoint (PR#19)
962     o AIM become a standard package ; clean unused CPP_kees
963    
964    checkpoint35
965     o subtantial rewrite of dynamics.f
966     o allows staggered time-stepping
967     o packaged obcs, aim (Atmospheric Physics), zonal_filt, shap_filt
968     o added bottom drag (linear + quadratic)
969    
970    checkpoint34
971     o Modified exchange routines for correct reverse simulation.
972     o Added verification output to comply with bug fix in impldiff.F
973       as of c33 (exp2, natl_box).
974     o Corrected store directives in dynamics.F
975     o Corrected array boundaries in impldiff.F
976     o Corrected array dimensions in gmredi.
977       Added temporary scalars to avoid storage by TAMC (gmredi, kpp).
978     o Added routine GLOBAL_SUM_INT to global_sum.F
979     o Added mdsio_gl.F for ECCO purposes.
980    
981    checkpoint33
982     o Fixed bug in lower BC in impldiff.F
983     o Fixed confusion about units of forcing arrays in FFIELDS.h
984       namely Fu,Fv,Qnet,Qsw,EmPmR:
985       - Removed verification/natl_box/code/external_fields_scale.F
986         (did not differ from that in model/src)
987       - Changed units of fu,fv,Qnet,Qsw,EmPmR back to proper units
988         (see FFIELDS.h for description)
989       - Scale fu,fv,Qnet,Qsw,EmPmR when used in external_forcing_surf.F,
990         kpp_calc.F and kpp_transport_t.F
991       - Removed model/src/external_fields_scale.F and calls to it
992       - verification/natl_box uses flux data with "atmospheric" sign so
993         a special version of external_fields_load.F is used to
994         change the data as it's read in. This way, the arrays
995         have the right units and signs at all times tha a user could
996         possibly use them.
997     o Corrected genmake to delete a temporary file during configuration
998    
999    checkpoint32
1000     o Scaling of forcing fields moved from external_forcing.F to
1001       external_fields_scale.F, called directly after loading fields.
1002     o Surface relaxation terms added to surface tendency fields in
1003       extra subroutine.
1004     o impldiff.F changed to avoid storing by TAMC of huge intermediate
1005       fields.
1006     o Changes in GMRedi to comply with packaging structure.
1007     o Modifications and bug fixes in KPP package.
1008     o Verification experiment for KPP added (verification/nat_box/)
1009    
1010    checkpoint31
1011     o Packages interface brought in line with defined standard
1012       -> new routines packages_...F
1013       -> changes in initialise_fixed.F, initialise_varia.F
1014       -> enables initialisation interface for ECCO package
1015     o Rescaling of external forcing fields immediately after read
1016       (removed from external_forcing.F)
1017       -> new: surface_fields_scale.F
1018       -> affects ini_foorcing.F, external_fields_load.F
1019     o Computation of surface tendencies before tendency updates.
1020       -> new: external_forcing_surf.F
1021     o Change of variable names useGMRedi, useKPP, useECCO
1022     o Further changes in dynamics.F, impldiff.F for TAMC
1023     o Tested for experiments 0, 2, 4.
1024       - not yet tested for GMRedi.
1025       - not yet tested for KPP.
1026    
1027    checkpoint30
1028     o Updates for OBCs
1029     o New experiment for test OBCs (internal_wave)
1030     o New output.txt in verification/exp4/results/
1031     o Updates for KPP: calculation of viscosity corrected, new variables
1032       kpp_dumpFreq and kpp_taveFreq
1033     o CPP'd calls to system routines for TARGET_CRAY_VECTOR (P.H.)
1034       and entry for T90 in genmake
1035    
1036    checkpoint29
1037     o packaged KPP and GM/Redi
1038    
1039    checkpoint28
1040     o Corrected initialization of hFacMin for difference vertical coordinates
1041     o Modified calculation of partial cell thickness to more robust algorithm
1042     o config_summary.F: replaced write_1D_... by write_0D_...
1043       read_write.F: added subroutines write_0D_...
1044     o routines calc_common_factors, calc_gs, calc_gt, calc_phi_hyd, find_rho:
1045       included some initialisations required for TAMC.
1046     o routine calc_isoslopes: modified sqrt computation to make it
1047       differentiable for zero argument.
1048     o routines dynamics, impldiff:
1049       included store directives and tape key computations required for TAMC.
1050     o ALL TAMC related changes are between ifdef ALLOW_AUTODIFF_TAMC
1051    
1052    checkpoint27
1053     o fixed calc_isoslopes() to so that cancellation of terms in tensor
1054       formulation works properly with variable resolution.
1055     o restructured the calling tree between the_model_main() and dynamics()
1056     o split initialise() into initialise_fixed() and initialise_varia()
1057     o introduced initialization file pSurfInitFile for free surface (cg2d_x)
1058       in new routine ini_psurf()
1059    
1060    checkpoint26
1061     o moved some #ifdef in optional routines to encompass everything
1062       except the argument declarations. This was using up some memory.
1063     o fixed macros like _ADT to work with both versions of Linux cpp
1064     o cleaned up some unused and uninitialized variables
1065       (helps when debugging with strict compile options)
1066     o split up CPP_EEOPTIONS.h into options and macros
1067       (macros are now in eesupp/inc/CPP_EEMACROS.h)
1068     o patch for 2 processor JAM mode
1069    
1070    checkpoint25
1071     o updates for OBCs and NH
1072     o dramatically reduced memory usage
1073     o added "JAM" routines for use on Artic network
1074     o parameterization of convection by implicit vertical diffusion
1075     o added a diagnostic of convective index
1076     o moved call to calc_divg_hat() from dynamics() to solve_for_pressure()
1077       (this is was partly for OBCs and NH but also in anticipation of
1078        an explicit free-surface option)
1079     o initial state input files for U and V
1080    
1081    checkpoint24
1082     o Removed some spurious "REAL"s
1083     o Updated KPP to latest version from JPL
1084     o Added a COS(latitude) in horizontal viscosity
1085    
1086    checkpoint23
1087     o Added del^4 diffusion for salt and temperature
1088     o Multiple minor fixes (implicit none's, arguments to global_sum, etc.)
1089    
1090    checkpoint22
1091     o Made many small changes: implicit none, format(i5), etc...
1092     o Introduced KPP mixing scheme (compiles and runs but *not* tested)
1093    
1094    checkpoint21
1095     o Introduced new I/O package that supports direct-access records (mdsio.F)
1096     o Split-up mid-level I/O routines
1097     o Updated input data for experiments 0,1,2,4 and 5
1098     o Finally removed the genmake.* scripts (everyone seems happy with genmake)
1099    
1100    checkpoint20
1101     o Tidied up open-boundaries some more
1102     o Added non-hydrostatic code
1103     o Fixed some input data
1104     o Modified MPI start-up to work in "coupled" context.
1105    
1106    checkpoint19
1107     o Tidied up directory structure for verification experiments.
1108     o Added new code for open boundaries.
1109     o Corrected advective terms near boundaries.
1110     o Added CPP flag for "old UV" geometry and corrected spherical geometry.
1111     o Implemented free-slip/no-slip BCs and del^4 dissipation.
1112     o Split default parameters out of ini_parms.F into set_defaults.F.
1113    
1114    checkpoint18
1115     o Further changes for optional code generation with
1116       CPP - ugh!
1117    
1118    checkpoint17
1119     o Beginning to incorporating modularity
1120       for supporting atmos. and ocean configs.
1121     o Small changes for TAMC complinace
1122     o New ultra simple test case. Barotropic box
1123       configuration "exp0".
1124    
1125    checkpoint16
1126     o Minor changes to exchange routines to support
1127       accumulation in reverse mode. Adjoint of an
1128       assignment is an addition.
1129     o Changes to support compilation by TAMC (LR(1)
1130       syntax) and by g77 under Linux ( 72 character
1131       limit ).
1132    
1133    checkpoint15
1134     Rewrote exchange routines to allow
1135     o General tile <-> tile connectivity.
1136     o TAMC reverse mode flag ( is this really needed? )
1137     o Variable width overlap and exchange region widths.
1138     o Hooks for sharing data via DMA, Arctic, Memory Channel,
1139       shmput/shmget, VIA, SYSV shmem and every other cool
1140       communication method.
1141    
1142    checkpoint14
1143     Consistent isomorphism chages made. These allow dynamical core
1144     to support both pressure and meters as a vertical coordinate.
1145    
1146    checkpoint13
1147     Release which corrects global_max problem and a couple
1148     of KGM time-averaging diagnostic errors. Also contains
1149     a double-gyre experiment setup as well as the baseline
1150     4 degree global template. The double gyre template is
1151     organised to allow it to be applied automatically to
1152     the baseline configuration.
1153    
1154    branch-point-rdot
1155     A few tidy-ups have been made bt the real point of this
1156     check-point (excuse the pun) is to allow a branch to be
1157     made for the purposes of adding the rDot alterations.
1158     I don't want to commit the rDot to the main-trunk until
1159     we've actually tested it...
1160    
1161  checkpoint12  checkpoint12
1162   Tidy up and added generic genmake script (formerly   Tidy up and added generic genmake script (formerly
1163   configure script).   configure script).

Legend:
Removed from v.1.10  
changed lines
  Added in v.1.169

  ViewVC Help
Powered by ViewVC 1.1.22