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

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

  ViewVC Help
Powered by ViewVC 1.1.22