/[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.53.4.56 by heimbach, Mon Aug 18 10:33:41 2003 UTC revision 1.151 by jmc, Thu May 1 22:43:25 2003 UTC
# Line 3  $Header$ Line 3  $Header$
3      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
4      ==============================      ==============================
5    
6  ecco_c51_e34f  o add time-average diagnostic of U*V and vertical diffusive flux for theta.
7  o added case gfdl-sgi-o3k to genmake  o bug fixed in convective adjustment diagnostic (multi-tile set-up).
8  o added genmake option -arpack and provided default lib dir  o 4th Order Adv scheme: return to previous (c40pre7) masking near boundaries.
9  o exf_mapfields:  
10    Interpolate ustress,vstress computed via exf bulk formulae  checkpoint50c_pre
11    to W/S points. Assume  Merging from release1_p13:
12    on input: uwind,vwind on C-points  o bug fix for pkg/seaice dynamic solver
13    on output: ustress[C],vstress[C] -> fu[W],fv[S]  o Added SEAICE_initialHEFF to pkg/seaice
14  o Very prelim. preparations for forward dump&restart  
15    during adjoint calculation.  checkpoint50b_post
16    o store u*,v* in gU,V instead of in guNm1,gvNm1
17  ecco_C51_e34e  o in order to put the momForcing out of the Adams-Bashforth:
18  o bug fix in cost_drifter    move forcing & CD-scheme calls from mom_fluxform & mom_vecinv
19  o updating timeave package in ecco-branch to c51    to timestep.F
20    o new flag "useCDscheme" (default=F); replace guCD,gvCD by local arrays
21  ecco_c51_e34d  o re-generate output for highly sensitive test-exp:
22  o modif's to restore TAMC compatibility     global_ocean_pressure, ideal_2D_oce, lab_sea, natl_box
   - avoid "_d" constructs in PARAMETER declarations  
   - avoid interrupting lines by comment lines  
 o Use cluster local disks for purely local I/O  
   vs. globally visible disks needed for ctrl stuff  
   - new active_file_loc... and mdsio_..._loc routines  
   - affected packages: autodiff, ctrl, grdchk, mdsio  
 o bug fixes for 3d packing and I/O of sliced (xz/yz) fields  
   to increase I/O performance.  
   
 ecco_c51_e34c  
 o introducing integer flag debugLevel  
   - debLevZero: suppress all output  
   - debLevA: standard  
   - debLevB: equivalent to current debugMode=.TRUE.  
 o introducing pathname variable mdsioLocalDir for mdsio  
   that can be specified at runtime  
   (relevant for cluster I/O to local disk)  
   
 ecco_c51_e34b  
 o Modifications of obcs  
   - bug fixes in cost package  
   - bug fixes in grdchk package  
 o new parameter maxKPPhbl in KPP  
 o added input for obcs verification  
   
 ecco_c51_e34a  
 o Modifications to pkg/seaice, verification/lab_sea, and  
   verification/natl_box_adjoint for testing sea-ice adjoint code.  
   
 ecco_c51_e34  
 o merging c51 into ecco-branch  
   
 ecco_c50_e33a  
 o first step in c51 merge  
   - redundant update in model/src/  
   - merge back's in autodiff, ctrl, kpp, seaice  
   
 ecco_c50_e33  
 o minor changes to pkg/seaice pkg/cost/cost_ice.F  
   and verification/natl_box_adjoint/README  
 o Modif's to enable tangent linear run in ECCO environment  
   including time-dependent forcing  
   * main problem is missing call mdsreadfield out of exf_set_gen  
     (i.e. computation of trajectory in TLM mode)  
   * incorrect g_active_write_...  
   * (re-) init of bar files and g_bar files  
 o Reverted back to 3-dim. tbar, sbar for SST, SSS bar files.  
 o MERGES:  
   * obcs/ ctrl/ G. Gebbie code  
     (major changes in ctrl handling, volume flux balancing etc.)  
   * cost/ ctrl/ A. Koehl code  
     - new cost contributions (new data types)  
     - efficient pack/unpack for time-dependent 2-d fields  
 o replaced mutiple ctrl_get... by single generic ctrl_get_gen.F  
   (except obcs)  
 o somewhat cleaned package initialisation sequence for  
   ctrl/ cost/ ecco/  
 o extended grdchk to enable obcs gradient checks  
 o some generalizations for easier TLM <-> ADM switches  
 o added configuration example for parallel divided adjoint  
   (see verification/natl_box_adjoint/code_diva_mpi/README )  
 o started verification exp. for obcs control  
   (still needs checking)  
   
 ecco_c50_e32  
 o Patrick's adjoint-model storage directives for pkg/seaice  
 o merged with checkpoint50g_post (pkg/seaice and verification/lab_sea)  
   - added SEAICE_MULTILEVEL for 8-category sea-ice thermodynamics  
   - LSR sea-ice dynamic solver moved to SouthWest B-grid location and  
     made the default because of faster convergence than ADI  
 o Modifications for adjoint of SEAICE_ALLOW_DYNAMICS  
   - extended state in checkpoint_lev...  
   - in lsr.F, replaced GO TO by fixed loop index (also avoid DO WHILE)  
   - YEAH: in lsr.F, (first) successful use of CADJ LOOP = ITERATION  
   - avoid stoing in S/R budget for now  
   - in dynsolver.F, avoid NPSEUDO loop for adjoint  
   - new comlev1_lsr in the_main_loop  
 o Small modif's in lsopt/, optim/ to ease compilation  
   New script to run optim (lsopt)  
   
 ecco_c50_e31  
 o exf totally restructured  
   - replaced all exf_set_ routines by one generic exf_set_gen.F  
     and made sure it's still differentiable  
   - split exf_getffields.F into new exf_bulkformulae.F  
   - exf_inscal_/exf_outscal_ suffixes consistent with field names  
 o New cost function: Tini-Tfin to do model spinup optimization.  
 o Started new verification exp for Tini-Tfin optimization.  
23    
24  ecco_c50_e30  checkpoint50b_pre
25    o restore the use of momAdvection & useCoriolis flags (mom_vecinv)
26  o corrected CPP option mismatch in exf_Bulk... functions for case  o corrected CPP option mismatch in exf_Bulk... functions for case
27    #define ALLOW_ATM_WIND && #undef ALLOW_ATM_TEMP    #define ALLOW_ATM_WIND && #undef ALLOW_ATM_TEMP
 o Merging from release1_p13 and from icebear5  
   - bug fix for pkg/seaice dynamic solver  
   - Added SEAICE_initialHEFF to pkg/seaice  
   - Additional options for pkg/cost/cost_ice.F  
   - Updates to verification/natl_box_adjoint  
   - pkg/ctrl/ctrl_init.F  
     Fix for determination of startrec and endrec.  
   - pkg/cal/cal_toseconds.F  
     Fix for negative intervals shorter than 1 day.  
   
 ecco_c50_e29  
 Bug fix for merging between c50 and KPP.  
 ikey was passed from thermodynamics to kpp_calc via  
 common block rather than being recomputed in kpp_calc,  
 in contradiction with new key itdkey.  
 New key ikppkey created, and tamc.h headers updated.  
   
 ecco_c50_e28  
 o major merge:  
   ecco-branch merged with checkpoint50  
   gmredi now available, grdchk, DIVA up-to-date for MPI  
   However, two outstanding problems:  
   - kpp adjoint is broken, needs fixing!  
   - a few minor modif's in handling of CPP option  
     ALLOW_TANGENTLINEAR_RUN in conjunction with  
     ALLOW_ECCO_OPTIMIZATION are necessary to enable  
     TLM runs with ECCO setup  
28    
29  ecco_c44_e27  checkpoint50a_post
30    o Bug fix for merging between c50 and KPP.
31      ikey was passed from thermodynamics to kpp_calc via
32      common block rather than being recomputed in kpp_calc,
33      in contradiction with new key itdkey.
34      New key ikppkey created, and tamc.h headers updated.
35    
36    checkpoint50
37    - 2nd step of merge with ecco-branch:
38      merge changes in ecco-branch back into main trunk
39    - bug fix in find_alpha
40    - more storing in gmredi_calc_tensor.F, kpp_routines.F
41    - avoid iterative loop in ini_pressure for ALLOW_AUTODIFF_TAMC
42    - merged some sliced AD support routines
43    
44    checkpoint49
45    merged from ecco-branch:
46  o exf:  o exf:
47    - Enable initialisation of forcing fields to constant    - Enable initialisation of forcing fields to constant
48      (runtime) values.      (runtime) values.
# Line 154  o exf: Line 56  o exf:
56        fldsecs  = int(fldsecs/fldperiod)*fldperiod        fldsecs  = int(fldsecs/fldperiod)*fldperiod
57        fldcount = int(fldsecs/fldperiod) + 1        fldcount = int(fldsecs/fldperiod) + 1
58      and were modified.      and were modified.
   - added repeatPeriod (from MAIN branch)  
 o grdchk: updated and merged  
   - has both ADM and TLM checks  
   - works for single- and multi-proc.  
   - output cleaned  
   - worked successfully for parallel DIVA  
   Needed one update in EEPARAMS.h and updates in ctrl package  
   More modif's may be necessary for ECCO code to support TLM.  
59  o autodiff: checkpoint_lev...  o autodiff: checkpoint_lev...
60    - updated storing for bulk formulae (some arrays were missing)    - updated storing for bulk formulae (some arrays were missing)
61    - added state for seaice package    - added state for seaice package
62    - synchronized with main branch    - synchronized with main branch
63  o DIVA:  o DIVA:
64    Modif.s in    Modif.s in
65    - cost_init.F    - the_main_loop.F
66    - cost_averagesfields.F    - the_model_main.F
   - ecco_the_main_loop.F  
   - ecco_the_model_main.F  
67    and added cost_final_restore.F cost_final_store.F    and added cost_final_restore.F cost_final_store.F
68    o With these changes, it is possible to fully separate    o With these changes, it is possible to fully separate
69      the first forward run to compute the cost function (mdthe_main_loop)      the first forward run to compute the cost function (mdthe_main_loop)
# Line 185  o DIVA: Line 77  o DIVA:
77    The current changes apply to TAF-1.4.26.    The current changes apply to TAF-1.4.26.
78    TAF has one remaining bug, requiring one manual change    TAF has one remaining bug, requiring one manual change
79    (see MITgcm manual).    (see MITgcm manual).
 o bug fix plus modif for ubar, vbar in cost_averagesfields  
80  o cal: synchronized between MAIN and ecco-branch  o cal: synchronized between MAIN and ecco-branch
81    
82  icebear4  checkpoint48i_post
83  o Added cost_ice for sea-ice volume sensitivity experiment.  o Changes to restore differentiability of code w.r.t. previous tag
84      (mostly adding new routines to make list and replacing
85      pressure by totPhiHyd).
86    o Updated and merged grdchk package
87      - has both ADM and TLM checks
88      - works for single- and multi-proc.
89      - output cleaned
90      - worked successfully for parallel DIVA
91    o genmake:
92      - Adding -platform=SP4
93      - Modif. for DEC Alpha (many don't have makedepend; use mkdep instead).
94    o Modif. for DIVA:
95      ctrl_unpack should not be called in intermediate DIVA interval.
96      (similar thing would apply for active files such as
97      tbar, sbar, psbar, which should not be reinitialised;
98      only relevant for ecco-branch)
99    
100    checkpoint48h_post
101    o improve "ini_p_ground" when using finite difference form for calc_phi_hyd
102    o Added net flux to downward flux conversion to pkg/exf/exf_getffields.F
103    o Added SEAICE_initialHEFF to pkg/seaice
104  o ecco_the_model_main.F: SEAICE_WRITE_PICKUP must precede WRITE_CHECKPOINT  o ecco_the_model_main.F: SEAICE_WRITE_PICKUP must precede WRITE_CHECKPOINT
 o Reduced derivative singularities in pkg/seaice.  
 o Updated verification/natl_box_adjoint  
105    
106  icebear3  checkpoint48g_post
107  o Removed spurious iikey and idkey  --- EOS:
108  o pkg/exf  o compute locally the pressure for use in EOS: UNESCO, JMD95P, MDJWF
109    - bug fix in treatment of downward radiation    store total Potential in totPhyHyd for diagnostics & EOS.
110  o pkg/seaice  o fix restart and overlap Pb when using Z-coord and EOS funct. of P
111    - added initial sea-ice thickness parameter SEAICE_initialHEFF  --- Pressure/geopotential gradient:
112  o updated verification/natl_box_adjoint  o use Ro_surf & R_low instead of hFac in cal_phi_hyd:
113    o NLFS + p-coord.: eta*Alpha' is not longer lagging 1.time-step behind
114  icebear2    ( changes results of global_ocean_pressure )
115    o NLFS + z-coord.: add eta*Rho' contribution.
116    o accurate phiHyd when using interface-W at the middle between 2 centers
117     ( changes results of ideal_2D_oce )
118    o includes r* 2nd term in potential gradient.
119    
120    checkpoint48f_post
121    Merging from release1_p12:
122    o Modifications for using pkg/exf with pkg/seaice
123      - improved description of the various forcing configurations
124      - added basic radiation bulk formulae to pkg/exf
125      - units/sign fix for evap computation in exf_getffields.F
126      - updated verification/global_with_exf/results/output.txt
127    o Added pkg/sbo for computing IERS Special Bureau for the Oceans
128      (SBO) core products, including oceanic mass, center-of-mass,
129      angular, and bottom pressure (see pkg/sbo/README.sbo).
130    o Lower bound for viscosity/diffusivity in pkg/kpp/kpp_routines.F
131      to avoid negative values in shallow regions.
132      - updated verification/natl_box/results/output.txt
133      - updated verification/lab_sea/results/output.txt
134    o MPI gather, scatter: eesupp/src/gather_2d.F and scatter_2d.F
135    o Added useSingleCpuIO option (see PARAMS.h).
136    o Updated useSingleCpuIO option in mdsio_writefield.F to
137      work with multi-field files, e.g., for single-file pickup.
138    o pkg/seaice:
139      - bug fix in growth.F: QNET for no shortwave case
140      - added HeffFile for specifying initial sea-ice thickness
141      - changed SEAICE_EXTERNAL_FLUXES wind stress implementation
142  o Added missing /* */ to CPP comments in pkg/seaice, pkg/exf,  o Added missing /* */ to CPP comments in pkg/seaice, pkg/exf,
143    kpp_transport_t.F,and ecco_the_main_loop.F    kpp_transport_t.F, forward_step.F, and the_main_loop.F
144  o adjoint-friendly modifications to pkg/seaice  o pkg/seaice:
145  o added verification/lab_sea as per release1_p12_pre    - adjoint-friendly modifications
146      - added a SEAICE_WRITE_PICKUP at end of the_model_main.F
147  ecco_c44_e26  
148  o Merging pkg/exf and pkg/seaice from cvs tag ecco_ice2  checkpoint48e_post
149  o exf: Modifications to correctly handle evap field for adjoint  o New S/R for diagnostic of bottom pressure (phi0surf was missing in c48d_post)
150  o grdchk:  o dynamics: change definition of computational domain & adapt mom_fluxform
151    Added code to perform gradient checks for bulk formulae/atmos. state    accordingly ; when viscA4=0, allows to run the dynamics with Olx=Oly=2.
152    (atemp, aqh, uwind, vwind).  o following exp works with Olx=Oly=2 (and pass through checkbounds option):
153    NOTE:      adjust_nlfs.cs-32x32x1 adjustment.128x64x1 adjustment.cs-32x32x1
154    This package is out of synch with c48 package.      exp0 exp1 exp2 exp4 exp5 hs94.cs-32x32x5 hs94.128x64x5 hs94.1x64x5
155    The latter also has the tangent linear gradient checks.      global_ocean.90x40x15 global_with_exf global_ocean_pressure
156    Need to be merged.      inverted_barometer solid-body.cs-32x32x1
157  o verification/natl_box_adjoint:    (high order AdvScheme and/or multiDimAd works with Olx=Oly=2 ? not sure)
158    * added new setup: code_bulk/, input_bulk, results_bulk  
159      (Experiment 2, described in README)  checkpoint48d_post
160      uses bulk formulae and atmos. state as control instead of  o in preparation for r*:
161      air-sea fluxes.   a) use pre-computed gradient of hydrostatic potential:
162      Performs gradient checks for first 4 elements of atemp control.      changes in timestep.F & mom_cdscheme.F affects results of ideal_2D_oce
163    * changed standard verification to perform   b) move phi0surf from calc_phi_hyd to calc_grad_phi_hyd:
164      gradient checks on first 4 elements on hflux control.      => affects results of glob_oce_pressure (different truncation error)  
165    * gmredi switched off (will come after merge with c48).    update output.txt of ideal_2D_oce & global_ocean_pressure
166    
167  ecco_c44_e25  checkpoint48d_pre
168  o Merging some GM changes from main trunk (c47) to ecco-branch.  o in preparation for r*:
169  o Updating verif. natl_box_adjoint with GMRedi.    new S/R (calc_grad_phi_hyd.F) to compute Hydrostatic potential gradient.
170  o Adding setup with explicit differentiaton of pressure solver.    pass the 2 comp. of the grad. as arguments to momentum S/R.
171  o Bug fix: missing OBCS_OPTIONS.h for sponge layer code.    but for the moment, only used if it does not change the results.
172  o (Re-?) Commiting old bug fixes in kpp to ecco-branch  o add finite volume form of calc_phi_hyd for buoyancyRelation OCEANIC & OCEANICP
173    (kpp_calc.F, kpp_routines.F)  
174    checkpoint48c_post
175    Modifications in treatment of initial flux forcing controls
176    to avoid extra storing of flux fields.
177    
178    checkpoint48c_pre
179    o Added tidal_basin_2d experiment (currently 1D!)
180    o can use r* Coordinate with OBC (but Orlanski Rad. BC not yet implemented)
181    
182    checkpoint48b_post
183    o r* coordinate added in #ifdef NONLIN_FRSURF block.
184     (modification to pressure gradient not yet implemented)
185    o mom_fluxform: new S/R mom_calc_rtrans.F computes vertical transports
186      (for U & V) that are used for vertical advection of momentum.
187    
188    checkpoint48a_post
189    o In preparation for r*, use temp. array to compute vertical transport
190       in flux-form momentum pkg. => change results of few test-exp.
191      Update output of: exp4, exp5, front_relax, global_ocean_pressure,
192                        lab_sea, natl_box & plume_on_slope
193    
194    checkpoint48
195    o fix small bug that was causing lot of warning (TIMER_STOP 'UPDATE_CG2D')
196    
197    checkpoint47j_post
198    o add new exp. (ideal_2D_oce) to test "exotic" parameters and options.
199    o moved adexch from forward_step to addummy_in_stepping
200    o added exch in ctrl_map_...
201    o aveFreq.=0. at end of cost_final to avoid overwrite in
202      reverse checkpointing loop
203    o Yet more changes in gmredi:
204      * adgmredi_calc_tensor
205        avoiding all recomputation of gmredi_slope_limit
206      * adgmredi_x/y/rtransport
207        added flag for excessive storing to avoid recomp. of
208        u/v/rtans, dTdx/y/z
209        -> this is not really necessary and very memory-consuming
210      * adgmredi_slope_psi:
211        consistency with gmredi_slope_limit in treatment of GM_slopeSqCutoff
212      * gmredi_slope_limit
213        re-activated full calculation of taperfct for case 'ac02'
214    o updating verification/ adjoints of carbon, global_ocean.90x40x15
215    
216    checkpoint47i_post
217    o GM_EXCLUDE_() replace GM_TAPER_(): by default, all the code is compiled.
218    o add GM_Small_Number & GM_slopeSqCutoff as run-time parameters (data.gmredi)
219      default values = 1.D-12 & 1.D+48 (=> give same results as checkpoint47f_post)
220    o cg2dTargetResWunit: define tolerance using deltaTmom (and not deltaTfreesurf)
221    o gmredi: fix few bugs.
222    
223    checkpoint47h_post (duplicate tag for PH's c48 tag, in order to move c48)
224    This will be a good checkpoint to plug into the ecco-branch.
225    o few modif.'s to get latest version adjointed
226      (mainly kick out code in ini_linear_phisurf)
227    o modif's to run adjoint with exactConserv
228    o case GM_BOLUS_ADVEC should be cleaned
229      S/R gmredi_slope_psi should be cleaned
230    o verification/carbon now has exactConserv=.TRUE.
231    
232    checkpoint47g_post
233    o time-average AIM output: becomes independent from statvars time-average
234      (can be used with #undef ALLOW_TIMEAVE) ; add new variables (e.g., RH)
235    o After TAF bug fix:
236      Tested adjoint of gmredi package.
237      Tests were done with checkpoint47 rather than with latest
238      checkpoint47f_post since the latter is broken for the adjoint
239      verification/carbon/ contains both ADM and TLM gradient checks.
240    Here they are, before they get changed/lost/stolen.
241    Mostly modif.'s to fix numerical sensitivities.
242    Gradient checks OK for
243    - GM_taper_scheme:
244      * clipping
245      * ac02
246      * linear
247      * glw91
248      * dm95
249      * ldd97
250    - GMREDI_OPTIONS:
251      * GM_VISBECK_VARIABLE_K
252      * GM_NON_UNITY_DIAGONAL
253      * GM_EXTRA_DIAGONAL
254      * GM_BOLUS_ADVEC
255    in conjunction with data.gmredi parameters to be checked in
256    in a few minutes under verification/carbon/code/
257    
258    checkpoint47f_post
259    Merging from release1_p10:
260    o modifications for using pkg/exf with pkg/seaice
261      - pkg/seaice CPP options SEAICE_EXTERNAL_FORCING
262        and SEAICE_EXTERNAL_FLUXES
263      - pkg/exf CPP options EXF_READ_EVAP and
264        EXF_NO_BULK_COMPUTATIONS
265      - usage examples are Experiments 8 and 9 in
266        verification/lab_sea/README
267      - verification/lab_sea default experiment now uses
268        pkg/gmredi, pkg/kpp, pkg/seaice, and pkg/exf
269    
270    checkpoint47e_post
271    o created a tag prior to creating a branch
272    
273    checkpoint47d_post
274    o fixed bug in the oceanic pressure coordinates code: vertical viscosity
275      at the bottom boundary had an erroneous half slip boundary condition
276    o added bottom drag and no slip boundary condition capabilities to
277      oceanic pressure coordinates code
278    o added new verification experiment global_ocean_pressure, this experiment is
279      also described in the tutorial section of the manual
280    checkpoint47d_pre
281    
282    checkpoint47c_post
283    o Added a verification experiment for a global cubed-sphere ocean run
284      -  verification/global_ocean.cubed32x32x30
285      -  uses pkg/therm_seaice and pkg/bulk_forcing
286    
287    o allow to run AIM physics with SPEEDY input files (from Franco Molteni)
288    o allow a more accurate definition of Ro_Surf (selectFindRoSurf=1)
289      when using P-coordinate; only implemented for atmospheric config.
290    o OCEANICP & realFreshWater: include P-E direct effect on wVel ;
291      NOTES: requires option NONLIN_FRSURF to be "#define".
292    o update advect_xz/results/output.txt (left from checkpoint44g_post)
293    
294    checkpoint47b_post
295    Merging from release1_p9:
296    o pkg/seaice
297      - removed GOTO's and added taf directives
298      - double precision constants to reduce the g77 (Linux)
299        to F77 (SGI) differences reported in release1_p8
300    o tools/genmake
301      - added SGI options
302    o verification/testscript
303      - updated to that of checkpoint47a_post
304    o verification/global_ocean.90x40x15/input/eedata
305      - modified for SGI f77 compatibility
306    o verification/lab_sea
307      - added description of sea-ice model
308      - added missing matlab routines
309      - added test of thermodynamics parallelization
310    
311    checkpoint47a_post
312    o new pkg: aim_v23
313     = F. Molteni atmos.physics (SPEEDY, ver23) adapted to MITgcm.
314       for now, keep the same surface forcing as before.
315     - Part-Cell implemented into AIM; check that Heat & Water are conserved
316     - aim.5l_cs: use new aim pkg (run 10yr & get better results than before)
317     - aim.5l_LatLon & aim.5l_Equatorial_Channel are still using old aim pkg
318    o add diagnostic of surface correction term in monitor output
319    
320    o bulk_force and therm_seaice
321    - Two packages:  bulk_force (Bulk forcing)
322      and therm_seaice (thermodynamic_seaice) - adopted from LANL CICE.v2.0.2
323    - Earlier integration from Stephaine Dutkiewicz
324      and Patrick Heimbach.
325    - Two ifdef statements for compile time,
326       ALLOW_THERM_SEAICE and ALLOW_BULK_FORCE
327    
328    checkpoint47
329    differentiable version of checkpoint46n_post
330    o bug fix in quasihydrostaticterms
331    o * "clean" adjoint code (in terms of extensive recomputations)
332        can now be obtained for all GMREDI options (i.e. for
333        - GM_VISBECK_VARIABLE_K
334        - GM_NON_UNITY_DIAGONAL
335        - GM_EXTRA_DIAGONAL
336        - GM_BOLUS_ADVEC )
337      * However, wrong gradient check problem remains unsolved.
338      * New CPP options have been introduced for different
339        tapering schemes
340    o external_fields_load now part of differentiation list
341    o pressure needs multiple storing;
342      would be nice to have store_pressure at beginning or
343      end of forward_step, e.g. by having phiHyd global (5-dim.)
344      (NB: pressure is needed for certain cases in find_rho,
345      which is also invoked through convective_adjustment).
346    o recomputations in find_rho for cases
347     'JMD95'/'UNESCO' or 'MDJWF' are OK.
348    o #define ATMOSPHERIC_LOADING should be differentiable
349    o ini_forcing shifted to begining of initialise_varia
350    o Incorporating QNVS line search routines into MITgcm
351      (this is separate code, and therefore not under pkg)
352      lsopt/, optim/
353    o Updated verification/carbon/
354      - #define GM_VISBECK_VARIABLE_K
355      - GM_taper_scheme    = 'ldd97'    
356    
357    checkpoint46n_post
358    Merging from release1_p8:
359    o verification/natl_box:
360      updating new external_fields_load routine
361    o New package: pkg/seaice
362      Sea ice model by D. Menemenlis (JPL) and Jinlun Zhang (Seattle).
363      The sea-ice code is based on Hibler (1979-1980).
364      Two sea-ice dynamic solvers, ADI and LSR, are included.
365      In addition to computing prognostic sea-ice variables and diagnosing
366      the forcing/external data fields that drive the ocean model,
367      SEAICE_MODEL also sets theta to the freezing point under sea-ice.
368      The implied surface heat flux is then stored in variable
369      surfaceTendencyTice, which is needed by KPP package (kpp_calc.F and
370      kpp_transport_t.F) to diagnose surface buoyancy fluxes and for the
371      non-local transport term.  Because this call precedes model
372      thermodynamics, temperature under sea-ice may not be "exactly" at
373      the freezing point by the time theta is dumped or time-averaged.
374      N.B.: fairly large differences in accuracy occur across
375            different platforms/compilers; comparison between
376            g77 (Linux) and F77 (SGI Origin 2000/MIPSpro)
377            generated output gives:
378       T           S           U           V
379      C D M    c        m  s        m  s        m  s        m  s
380      n p a R  g  m  m  e  .  m  m  e  .  m  m  e  .  m  m  e  .
381      f n k u  2  i  a  a  d  i  a  a  d  i  a  a  d  i  a  a  d
382      g d e n  d  n  x  n  .  n  x  n  .  n  x  n  .  n  x  n  .
383    
384      Y Y Y Y  5  5  7  7  7  8 10  9  6  6  6  6  7  5  7  5  7 FAIL  lab_sea
385    o GAD:
386      - generated new common blocks to account for call of
387        same gad routines with differing traceridentities
388        (needed to modify tracerIdentity indices in GAD.h)
389      - generated separate common blocks for case useCubedSphereExchange
390        (Department of Futurology)
391      - parameter lists to gmredi_?transport: added tracerIdentity
392      - added new key indices to tamc.h
393    o external_fields_load:
394      - added this routine to TAF list
395      - needed to make some common blocks global and additional storing
396        along the same lines as exf package (checkpoint_lev?_directives.h)
397  o exf:  o exf:
398    updated external forcing package    updated external forcing package
399    - cleaned old exfa stuff    - bug fixes carried over from ecco-branch
   - bug fixes  
400      (missing OBCS_OPTIONS.h in two routines)      (missing OBCS_OPTIONS.h in two routines)
401    - enable easy to use "no forcing".    - enable easy to use "no forcing".
402    - added exf I/O for atmospheric loading    - added exf I/O for atmospheric loading
403    - added exf I/O for runoff data    - added exf I/O for runoff data
404    - transfered scaling between exf <-> MITgcm to exf namelist    - transfered scaling between exf <-> MITgcm to exf namelist
405    o Adding new verification experiment global_with_exf.
406  ecco_c44_e24    This verification experiment is almost identical to the
407      global_ocean.90x40x15 experiment, except that it uses
408      the exf (and cal) package instead of the
409      external_fields_load routine (cf README file).
410    
411    checkpoint46m_post
412    o modified testscript to return -- when there are no lines of "%MON" output
413    o added new flag "quasiHydrostatic" for QH mode
414      and new routine quasihydrostaticterms()
415    o added 2*Omega*Cos(phi)*W to u equations (non-hydrostatic Coriolis term)
416      New routine: mom_u_coriolis_nh()
417    o fixed sign of non-hydrostatic metric terms in u and v equations.
418      Depend on correct sign of gravitySign.
419    o added work.pc* to .cvsignore in aim directories using ifc
420    o introduced flag for controllin non-hydrostatic metric terms
421    o removed gravitySign from PARM04
422    o switched to g77/ifc from pgf77 in aim experiments
423    
424    checkpoint46l_post
425    o made convective adjustment work with pressure coordinates:
426     - changed the direction of k-loop in convective_adjustment.F for the
427       case of pressure coordinates (OCEANICP,ATMOSPHERIC buoyancyRelation)
428     - adjusted the reference pressure k-index in convective_adjustment.F
429     - adjusted the convection condition in convect.F (in analogy to
430       calc_ivdc.F)
431     - convective_adjustment no longer computes anything on the halos
432     - removed the warnings about negative salinity from find_rho.F and
433       find_alpha.F; instead the new routine look_for_neg_salinity, called
434       at the beginning of find_rho, find_alpha, and find_beta, does a
435       check of the entire slice, if CPP-option
436       CHECK_SALINITY_FOR_NEGATIVE_VALUES is defined
437    checkpoint46l_pre
438    
439    checkpoint46k_post
440    o fixed the verification/global_ocean.90x40x15 experiment:
441     - new bathymetry (the world according to A., JMC, and M.)
442     - new initial fields and forcing fields (*.bin files)
443     - new POLY3.COEFFS (for the next release one should switch to a full
444       equation of state: JMD95P or MDJWF)
445     - fixed several errors and redundancies in the data file
446     - experiment uses looped cells
447     - added matlab directory with diagnostic scripts for plotting of output
448    
449    o S/R aim_initialise.F replace S/R aim_init from file aim_do_inphys.F:
450     - read AIM physics parameters from a file (data.aimphys)
451     - set defaults values = F.Molteni paper (Clim.Dyn., 2002)
452    
453    checkpoint46j_post
454    o split calc_exact_eta in 2 S/R : integr_continuity & update_etaH
455    o move wVel computation at the end of the time step, in S/R integr_continuity
456    o create specific S/R to exchange T,S before DYNAMICS (for stagger time step)
457    o update timeave pkg for wVel diagnostic.
458    checkpoint46j_pre
459    
460    checkpoint46i_post
461    o Clean up AIM package (and keep the results unchanged):
462      include CPP_OPTION and use IMPLICT NONE in all routines ;
463      declare all the variables _RL ; use _d 0 for all numerical constants.
464      use ifdef ALLOW_AIM everywhere. And now AIM can be used with g77 !
465    
466    checkpoint46h_post
467    o cleaned up the use of rhoNil and rhoConst.
468      - rhoNil should only appear in the LINEAR equation of state, everywhere
469        else rhoNil is replaced by rhoConst, e.g. find_rho computes rho-rhoConst
470        and the dynamical equations are all divided by rhoConst
471    o introduced new parameter rhoConstFresh, a reference density of fresh
472      water, to remove the fresh water flux's dependence on rhoNil. The default
473      value is 999.8 kg/m^3
474    o cleanup up external_forcing.F and external_forcing_surf.F
475      - can now be used by both OCEANIC and OCEANICP
476    checkpoint46h_pre
477    
478     Added code to convert surface volume flux (fresh water) into
479     a mass flux when using P coordinates in the ocean (OCEANICP).
480     Note: It assumes you have set rho0=rhoConst=density of fresh water.
481    
482    checkpoint46g_post
483    o Include a new diagnostic variable phiHydLow for the ocean model
484      - in z-coordinates, it is the bottom pressure anomaly
485      - in p-coordinates, it is the sea surface elevation
486      - in both cases, these variable have global drift, reflecting the mass
487        drift in z-coordinates and the volume drift in p-coordinates
488      - included time averaging for phiHydLow, be aware of the drift!
489    o depth-dependent computation of Bo_surf for pressure coordinates
490      in the ocean (buoyancyRelation='OCEANICP')
491      - requires a new routine (FIND_RHO_SCALAR) to compute density with only
492        Theta, Salinity, and Pressure in the parameter list. This routine is
493        presently contained in find_rho.F. This routine does not give the
494        correct density for 'POLY3', which would be a z-dependent reference
495        density.
496    o cleaned up find_rho
497      - removed obsolete 'eqn' from the parameter list.
498    o added two new verification experiments: gop and goz
499      (4x4 degree global ocean, 15 layers in pressure and height coordinates)
500    checkpoint46g_pre
501    
502    checkpoint46f_post
503    o Enable tangent linear (forward mode) gradient checks:
504      - extended active file handling to g_... files
505      - added TANGENT_SIMULATION to theSimulationMode
506      - extended grdchk package accordingly
507    o added tangent linear setup in verification/carbon/code_ftl
508    o added adjoint setup in verification/front_relax/code_ad
509    
510    checkpoint46e_post
511    o Added new equation of state -> MDJWF
512      - EOS of McDougall et al., 2002, JAOT, submitted
513      - caveat: the equation of state is only valid for a smaller (more
514        realistic?) range of values than JMD95P/Z and UNESCO
515      - added masks to the calculation of pressure in store_pressure
516      - added more check values for density in check_eos (ini_eos.F), some of
517        the old check values are out of the range of the MDJWF-eos, so don't
518        expect perfect matches for those
519    checkpoint46e_pre
520    
521    checkpoint46d_post
522    o fixed store_pressure to work with both buoyancy relation = 'OCEANIC' and
523      'OCEANICP', also initialised field pressure correctly in ini_eos in the
524      case of pressure coordinates. eosType='JMD95Z' in
525      combination with buoyancyRelation='OCEANICP' now causes an error.
526    o Changed p = pressure(i,j,k,bi,bj)  to  p = pressure(i,j,kRef,bi,bj)
527      in find_alpha/beta.
528    checkpoint46d_pre
529    
530    o Changed p = pressure(i,j,k,bi,bj)  to  p = pressure(i,j,kRef,bi,bj)
531      so that JMD95Z and JMD95P give approptiate static stability in find_rho.
532    
533    checkpoint46c_post
534    o Changes necessary for ocean in p-coordinates
535     -  Added new buoyancy relation = 'OCEANICP'
536     -  Added new parameters = gravitySign (this used to be contained inside
537        the factor dRdZ which I added when we first switched to R coordinates).
538     X GM/Redi is not compatible (yet)
539     X bottom drag and no-slip need to be debugged.
540    checkpoint46c_pre
541    
542    checkpoint46b_post
543    o Added new equation of state -> JMD95Z and JMD95P
544      - EOS of Jackett and McDougall, 1995, JAOT
545      - moved all EOS parameters into EOS.h
546      - new routines ini_eos.F, store_pressure.F
547    o Added UNESCO EOS, but not recommended because it requires
548      in-situ temperature (see JMD95)
549    o Modified formatting for knudsen2.f in utils/knudsen2 and added
550      unesco.f to be used with POLY3
551    checkpoint46b_pre
552    
553    checkpoint46a_post
554    o Added (atmospheric) pressure loading
555      - new field in FFIELDS.h, etc...
556      - new cpp flag ATMOSPHERIC_LOADING
557    o Changed hFacC to _hFacC in calc_phi_hyd.F
558    o Added SHORTWAVE_HEATING to some files for consistency
559    checkpoint46a_pre
560    
561    checkpoint46
562    further merges from release1_p5
563    o new ctrl package
564      - adopted from ECCO environment to enable optimization
565      - added Eliassen Palm fluxes to controls
566    o added Eliassen Palm flux controls to gradient check package
567    o cost package
568      - Compatible with new ctrl package
569      - added Eliassen Palm cost hooks
570      - modif's of existing cost functions
571    o Adjoint-related bug fixes in kpp:
572      - kpp_calc: sore of kpphbl avoids recomputation/call to S/R kppmix
573      - kpp_routines: store of Rib avoids partial recomputation bug of TAF.
574    o autodiff package
575      - flow directives:
576        + new for S/R do_the_model_io
577        + enhanced for S/R checkpoint
578        + new gmredi.flow
579      - Introduced CPP option ALLOW_AUTODIFF_MONITOR to
580        disable adcommon.h/g_common.h by default.
581      - fixed problem for adjoint of global_max
582    o modified verification/carbon/
583    o New setup for adjoint of global_ocean
584    o added hooks for Eliassen Palm fluxes (dynamics.F, ini_fields.F)
585    o cleaned TAF keys iikey, idkey (dynamics.F, thermodynamics.F)
586    o incorporated hooks for sponge layer code (external_forcing.F)
587    
588    checkpoint45d_post
589    o import 2 fixes from release1_p5 (CG2D.h & packages_readparms.F)
590    o use recip_dx*,recip_dy* instead of /dx*,/dy* in orlanski_E,W,N,S
591      (affects plume_on_slope "cg2d" output)
592    o add new parameter deltaTfreesurf for free-surface time stepping
593    o exchange gT,gS when using staggerTimeStep & GM & Oly<4
594    o NONLIN_FRSURF: rescale gNm1 to get a better conservation with A-B
595    o change phiMin in exp1 (to agree with documentation)
596    o add new flags:
597     - T,S forcing inside/outside Adams-Bashforth
598     - turn off temp,salt advection and/or forcing
599     - for each tracer, define internal flag for multiDimAdvection & A-B
600    o monitor: change definition of KE to check conservation of Energy
601    
602    checkpoint45c_post
603    o Removed f90 comments ("!")
604    o Modified initialisations to enable exf package for MITgcm
605      (without invoking ECCO_PACKAGE).
606  o modifications to gradient check package (Martin Losch)  o modifications to gradient check package (Martin Losch)
607    - enable centered differences    - enable centered differences
608    - modified format of standard output    - modified format of standard output
609  o bug fix for nwettile contents  
610    so far assumed global sum  in first entry (1,1,k),  checkpoint45b_post
611    needed for global control vector  o cleaned exchange calls in cg2d
612    -> introduced new variable nwet[c/w/s]global(nr)  o Included CPP option ALLOW_AUTODIFF_MONITOR in forward_step
613  o shifts of calls to cal/exf package to be accessible  o included CPP option SINGLE_LAYER_MODE in thermodynamics
614    in MITgcm stand-alone    to configure barotropic setup (Martin Losch)
615  o shifted packages_init_variables before convective_adjustment_ini  o moved some initialisations within bi,bj loop in dynamics
616    (corresponding update of TAF tape initialisation).    (TAF flow dependency)
617  o removed auxiliary field in impldiff  o in initialise_varia, moved call packages_init_variables
618    (caused numerical differences for GMRedi cases)    before convective_adjustment_init
619  o updated obcs control part (G. Gebbie)    (corresponding adjustment of TAF tape initialisation).
620  o removed some f90 comments ("!")  o in calc_diffusivity/calc_viscosity extended domain to full overlap.
621  o Bug fixes and updates in cost package  
622   * cost_final: correct multiplier mult_tauv, mult_sflux  checkpoint45a_post
623   * cost_hyd: added call cost_ssh  o compute Rho over the whole domain (extend iMin,jMin indices)
624   * cost_init_barfiles.F: replaced arrays by tmpfld's  o specific exchange for X-slice domain (case Ny=1)
625   * cost_readtopexmean.F: ersmask -> tpmeanmask  o argument futureIter of S/R obcs_calc missing in exp4 & internal_wave /code
626   * cost_weights.F: corrected wtheta[2], wsalt[2]  o solve_for_pressure : set cg2d_x to zero where OBCS are applied
627  o cleaned exchanges in cg2d.F  
628    checkpoint44h_post = checkpoint45
629  ecco_c44_e23  o Fix to restore differentiability:
630  merged patches between release1_p1 and release1_p2:    place at which packages_init_variables is called in
631  o fixed array boundary problem for adjoint    initialise_varia is crucial for reverse control flow.
632    (related to partial recomputation within subroutine  o modifications in GMRedi to enable stable adjoint
633    related to dynamics -> calc_viscosity -> mom_...    (corresponding modif's in thermodynamics, dynamics)
634  o fixed recomputation in thermodynamics  o added missing hook for sponge layer code in external_forcing
635  o prepare (but not complete) recomputation problem in impldiff  o modified test cost functions
636  o fixed problem for adjoint of global_max  o storing in gad_advection (plus removal of write(0,*)!
637  o switched off sponge layer code in external_forcing  o missing headers TR1.h added in convective_adjustment routines
638  Bug fix in ctrl_mask_set_yz (spotted by G. Gebbie)  o errorMessageUnit set non-equal zero in eeboot
639      (conflict on some platforms, e.g. IBM SP3)
640  ecco_c44_e22  o modified carbon verif.
641  o Finally, the merges with the Scripps cost function terms:  
642    Various additions, modifications, quality checks.  checkpoint44h_pre
643  o Plus: Included switch ALLOW_AUTODIFF_MONITOR to handle  o vertical grid option: allow to put the Interface at the middle between
644    monitoring of adjoint state.    2 cell-centers ; replace delZ in pkg/kpp by drF.
645  o Plus: removed comm_stats, ini_communication_patterns from diff. list  o GM Advective form: Tracers are advected using the residual transport (=
646    (no more need to "keep" MPI headers for adjoint generation)    Euler+GM-bolus); set param GM_AdvSeparate=T to return to previous form.
647  o Plus: exf changes:  
648    o added run-off routine  checkpoint44g_post
649    o additional storing needed for ustress,vstess  o fix surface correction term in multiDimAdvection (affects plume_on_slope
650      because of state-dependent quality check    results) ; ifdef missing in thermodynamics.F
651    o added missing initialisations  
652  o Plus: genmake with MAKEDEPEND variable was still missing in ecco-branch  checkpoint44f_pre,checkpoint44f_post
653  NB: - verification/natl_box_adjoint/ 240 steps run  o added PTRACERS package
654        * looks ok under Linux (twain.lcs.mit.edu)    This allows an arbitrary number of passive tracers to be integrated
655        * looks ok on MIT SUN's (fjord.mit.edu)    forward simultaneously with the dynamicaly model.
656        * NaN in reverse mode on Scripps SUN's at iter=214    + Implemented so far:
657        * NaN in reverse mode on IBM SP3 (BlueHorizon) at iter=204      - basic forward algorithm (time-stepping, advection, diffusion, convection)
658      - 1 year forward int. global 2x2 deg. produces similar cost functions      - I/O and checkpointing
659        * all cost terms of new code are slightly smaller (0.5 to 3%).      - GM/Redi  *but*  using the GM/Redi coefficient of Salt
660      + Not implemented so far:
661  ecco_c44_e21      - KPP
662  o reversed order of ctrl/cost initialisation      - OBCS
663  o gmredi_slope_psi.F: Inconsistent types in MIN/MAX expressions    + No specific example supplied (yet) but global_ocean.90x40x15 has the
664  o modified to enable ctrl field swapping      necessary data.ptracer file. Simply use -enable=ptracers and uncomment
665  o bug fix for global sum of ncvarlength in ctrl_init      line in data.pkg.
666  o no more call to ctrl_unpack from ctrl_init    + This package is disabled by default since it increases storage.
667  o no more write of ctrl vector in ctrl_pack  
668    however, init. of control vector from ctrl_unpack  
669  o various bug fixes & consistent nobcs handling in  o testscript: added -cleanup option
670    ctrl_pack/unpack routines (by G. Gebbie)  o Shap_filt: a) new shap_filt S/R to use no-slip BC with S2 filter
671  o optim package: Added obcs control part for lsopt I/O               b) enable to filter 2D fields.
672    (by G. Gebbie)  
673  o missing call to obcs_diags  checkpoint44e_pre,post
674  o moved calls in packages_init_fixed to packages_readparms  o re-wrote convective adjustment to use pre-calculated weights
675    because some quantities are needed earlier.    Affects: model/src/convective_adjustment.F, convective_adjustment_ini.F
676      Adds:    model/src/convective_weights.F, convectively_mixtracer.F
677  ecco_c44_e20    Original calls to convect() still in place.
678  o before major merges from the Scripps setup  
679    here are two small updates:  chkpt44d_pre,post
680    - field swapping for surface flux controls enabled  o added missing EXCLUDE_MONITOR flags
681      (remains t.b.d. for atmos. state plus obcs)  o changed "e" to "_d" in gmredi_slope_limit, gmredi_slope_psi
682    - bug fix in active_file.F for xz-/yz- routines    (incompatible typ in MIN/MAX expressions caused problems
683    - added missing xz-/yz- routines to active_file_ad.F    on IBM SP3)
684    o in genmake added variable MAKEDEPEND
685  ecco_c44_e19    plus resetting for case SunOS
686  o merged necessary patches:  o added timer_stats.c routine for IBM SP3
687    - AD-related changes for GMRedi  o removed variables in dynamics
688    - fixes i KPP (delZ -> drF)  !!! I forgot to create tag chkpt44d_pre !!!
689    - hook to OBCS songe layer code in external_forcing  !!! But chkpt44d_pre = chkpt44c_post    !!!
690    - errorMessageUnit non-zero in eeboot.F  
691    - included GMREDI, EXACT_CONSERV to natl_box_adjoint verif.  o real fresh water flux implemented with non-linear free-surface.
692    
693  ecco_c44_e18  chkpt44c_pre,post
694  o sponge layer code for obcs was missing  o few fix (mask in shap_s2, EmPmR in external_field_load,
695    - added hook in external_forcing.F    USE_NATURAL_BCS in solve_for_P);
696    - added storing for TAF  o add arguments myIter & myTime to S/R obcs_calc & solve_for_P
697  o new routine cost_init_barfiles to initialise  
698    adjoint fields adtbar, adsbar, adpsbar  checkpoint44b_pre/_post
699    (needed in conjunction with corrected ladinit for active I/O)  o merge of relevant stuff from the ecco-branch:
 o cleaned up and corrected order in packages_readparms, ecco_readparms  
 o corrected and simplified weights for obcs controls  
   - weight fields simplified  
   - bug fix in ctrl_pack  
   - no more write of control vector in ctrl_pack  
 o i/j index bug fixes in cost_ssh_mean  
   
 ecco_c44_e17  
 o bug fixes in cost package  
   * wrong identifier for active I/O of monthly mean psbar  
     (cost_averagesfields, cost_ssh) xx_psbar_dummy  
   * cleaned all identifier for active I/O of monthly fields  
   * wrong ladinit in various cost routines  
 o merge with interval (chkpt44c_post - chkpt44d_post)  
   * added missing EXCLUDE_MONITOR flags  
   * changed "e" to "_d" in gmredi_slope_limit, gmredi_slope_psi  
     (incompatible typ in MIN/MAX expressions caused problems  
     on IBM SP3)  
   * in genmake added variable MAKEDEPEND  
     plus resetting for case SunOS  
     added timer_stats.c routine for IBM SP3  
   * removed variables in dynamics  
   
 ecco_c44_e16  
 o new ECCO checkpoint, the first since ecco_c43_e15  
   which is branched from MITgcm.  
   All changes were merged in between (see ecco-branch-mod? tags).  
 o last changes:  
   * assuming path for TAF and TAMC are known (e.g. in $path)  
   * genmake:  
     - added tools/xmakedepend to circumvene shortcoming of  
       SUN's makedepend (size too small)  
     - invoking ../tools/xmakedepend instead of makedepend  
     - changed compiler flags for case SunOS  
   
 ecco-branch-mod5  
 o mod's to reproduce open boundary run  
 o enable field swapping for OB's and climatology  
 o update store directives for OB control  
 o added option ALLOW_OBCS_PRESCRIBE  
 o call of getobcs now out of obcs_calc (param. list modified)  
   
 ecco-branch-mod4  
 o this is the second part of changes necessary to reconcile  
   the natl_box_adjoint verification (with and without grdchk).  
   >>>!!! CHANGES BETWEEN THE TAGS ecco-branch-mod3 !!!<<<  
   >>>!!!                      AND ecco-branch-mod4 !!!<<<  
   >>>!!! SHOULD BE MERGED INTO THE MAIN TRUNK      !!!<<<  
700    - genmake: removed $S64 overwrite for case SunOS    - genmake: removed $S64 overwrite for case SunOS
701    - pkg/exf: update and corrections for field swapping and obcs    - pkg/exf: update and corrections for field swapping and obcs
702    - pkg/ecco: parameter lists for the_model_main, the_main_loop    - pkg/ecco: parameter lists for the_model_main, the_main_loop
# Line 395  o this is the second part of changes nec Line 704  o this is the second part of changes nec
704    - pkg/autodiff: added flow directives for obcs, mdsio_gl_slice    - pkg/autodiff: added flow directives for obcs, mdsio_gl_slice
705                    updated checkpointing_lev... lists for obcs                    updated checkpointing_lev... lists for obcs
706    - model/src: minor changes in forward_step, plot_field    - model/src: minor changes in forward_step, plot_field
707                   added directive for divided adjoint in the_main_loop
708    - pkg/mdsio: added mdsio_gl_slice    - pkg/mdsio: added mdsio_gl_slice
709    
710  ecco-branch-mod3  o check parameters & config (chkpt44a_pre,post)
711  o this is the first part of changes necessary to reconcile  o OBC and NonLin_FrSurf.
   natl_box_adjoint verification (with and without grdchk).  
   This set of changes is relevant for ECCO only.  
   The next set/tag of changes (ecco-branch-mod4)  
   should be merged into the main trunk.  
   - update some options  
   - parameter lists of ecco_the_model_main and ecco_the_main_loop  
     differed from those of the MITgcm counterparts.  
   
 ecco-branch-mod2  
 o update ecco-branch-mod1 to checkpoint44  
 o exclude from updates:  
   - verification/  
   - adjoint/  
   - pkg/cost  
   - pkg/ctrl  
   - pkg/ecco  
   - pkg/obcs  
   - model/src/the_model_main.F  
   - model/src/the_main_loop.F  
712    
713  checkpoint44  checkpoint44
 o apply Shapiro filter on T&S before calc_PhiHyd (stagger time stepping)  
714  o fix bug in mom_vi_del2uv  o fix bug in mom_vi_del2uv
715    o select when filters are applied ; add options to zonal_filter (data.zonfilt)
716    o gmredi: fix Pb in the adiabatic form ; add options (.e.g. Bolus advection)
717    o update AIM experiments (NCEP input files)
718    o improve and extend diagnostics (Monitor, TimeAve with NonLin-FrSurf)
719  o added some stuff for AD  o added some stuff for AD
720    
 ecco-branch-mod1  
 o based on c43  
 o replaced packages cost, ctrl, ecco, obcs  
 o removed MITgcm verif.'s, added ECCO verif.'s  
 o added lsopt, optim for line search optimization.  
   
721  checkpoint43  checkpoint43
722  o added two new advection tests  o added two new advection tests
723  o added global 4x4 degree ocean test  o added global 4x4 degree ocean test

Legend:
Removed from v.1.53.4.56  
changed lines
  Added in v.1.151

  ViewVC Help
Powered by ViewVC 1.1.22