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

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

  ViewVC Help
Powered by ViewVC 1.1.22