/[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.1853 by gforget, Thu Apr 11 14:52:12 2013 UTC revision 1.2010 by jmc, Tue Jun 24 22:03:47 2014 UTC
# Line 4  $Name$ Line 4  $Name$
4      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
5      ==============================      ==============================
6    
7  o verification/offline_exf_seaice : revise data.seaice and  o pkg/seaice (ITD code):
8    update results accordingly in input.thermo and input_ad.    - fix definition of Hlimit for multi-threaded case
9      - use double precision arithmetic to compute Hlimit: this affects results;
10        update output of experiment seaice_itd.
11    o replace content of verification/seaice_itd
12      - instead of a 1D column use the configuration of offline_exf_seaice
13        to test the ice thickness distribution code with ridging
14      - this is a first shot at this experiment
15      - change of some defaults of ridging parameters SEAICEhStar (from
16        100 to 25) and SEAICEmuRidging (from 5 to 3) in seaice_readparms.F
17        and lot of new comments.
18    
19    checkpoint64y (2014/06/20)
20    o pkg/seaice/seaice_reg_ridge.F:
21      - modify the regularization of very small volumes of ice; only check
22        total volume HEFF, and do not apply correction to individual ice
23        classes because that leads to funny patterns in the heat and fw-fluxes
24        into the ocean. Adjust seaice_itd/results/output.txt
25      - replace nITD by SEAICE_multDim in loop-ranges
26    o pkg/exf:
27      - add code to process the case field period = -12. with UV interpolation.
28      - move debug print (formerly active when exf_verbose=T) from EXF_getFFieldRec
29        into S/R that calls it (i.e., exf_set_gen.F, exf_set_uv.F & obcs_exf_load.F)
30    o verification/global_with_exf:
31      - switch all fields period to -12. (was previously untested) in input/data.exf
32        (but keep secondary test ".yearly" unchanged); update the reference output.
33    o pkg/cal:
34      - new calendar type "noLeapYear" for 365 days calendar without any leap year.
35    o pkg/exf:
36      - new integer parameter "exf_debugLev" (default set to main model debugLevel)
37        to replace logical switch "exf_verbose"
38      - report in STDOUT when loading a new record (if exf_debugLev >= 3)
39    o pkg/sbo:
40      - update from Katherine Quinn: evaluate mass directly from density and SSH
41        (budget using Finite Volume method) instead of from bottom pressure.
42      - allow testreport to collect few selected SBO output ; tested in exp.
43        global_ocean.90x40x15 (after updating the results, to include SBO output).
44      - fix some multi-threading issues (in output part).
45    o model/src & most packages:
46      - always call {PKG}_READPARMS even if corresponding use{PKG}=F;
47      - in {PKG}_READPARMS, in case use{PKG}=F: add a call to new S/R
48        PACKAGES_UNUSED_MSG to print a weak warning when parameter file
49        "data.{pkg}" exist ; and leave the S/R (return) just after.
50      - more standard LAND pkg interface (land_init_fixed.F, land_check.F).
51        + call LAND_READPARMS from packages_readparms.F
52    o pkg/seaice:
53      - remove code that resets useThSice in S/R seaice_readparms and
54        replace it by a warning in S/R seaice_check
55      - remove parameter MULTDIM, replace with nITD and runtime
56        parameter SEAICE_multDim (still hard-wired to nITD if
57        SEAICE_ITD is defined)
58    
59    checkpoint64x (2014/05/24)
60    o verification/lab_sea:
61      - update adjoint and TLM output following changes to pkg/salt_plume
62      - un-do the update after putting back local recip_Cp in few places.
63    o pkg/salt_plume: modifications by An Nguyen:
64      - add capability to do brine rejection (salt-plume) as function of
65        volume+tracers (instead of only salinity tracer).
66      - conserve volume by exchanging volume from each kLev to above to
67        accommodate in-coming brine volume
68      - store surface and klev exchanged fraction + volume in 4 new 3-d arrays
69      - switch KPP surf buoyancy due to brine forcing to 3-d
70        (3 more 3-d arrays, 2 for SPforcing[T,S], 1 for boplume)
71      - add salt_plume_volfrac, salt_plume_apply, salt_plume_tendency_apply_t
72      - unrelatedly, also add ad-hoc capability to vary salt_plume scheme
73        in East/West Arctic.
74    o pkg/sbo:
75      - incorporate modifications from Katherine Quinn to enable the use of sbo
76        pkg in ECCO-v4 set-up ; remove SBO-TimeAve output ; add monitor-like
77        output for future testreport verification (unfinished).
78    o pkg/atm_phys:
79      - change diagnostics (AtPh_Trf & AtPh_Qrf) of DARGAN_BETTSMILLER S/R:
80        skip the reseting of t_ref & q_ref to initial t & q ; this preserve the
81        output of S/R CAPECALCNEW (Parcel temp & q) for diagnostics.
82    o pkg/seaice
83      seaice_calc_ice_strength.F: modify expression to avoid divisions
84        (affects truncation)
85      seaice_do_ridging.F: clean up the insectarium
86      - fix a bug that was compensated by areaScaling (-> remove areaScaling code)
87      - added SEAICEsnowFracRidge to redistributed snow (2nd bug)
88      - set closing/openingRate to zero where no extra iteration is need (3rd bug)
89      - add some debugging information and error messages (warnings)
90    o model/src: modifications relatd to atmosphere p*:
91      - Switch to more accurate p* coordinate slope term (calc_grad_phi_hyd.F)
92      - With p* or Sigma-P, use constant reference Pot.Temp (thetaConst) instead
93        of vertical profile tRef in geopotential background and anomaly.
94      - Acount for true p* correction in theta <-> T conversion in pkg/aim_v23
95        and pkg/atm_phys.
96      - Update FWD output of experiments: aim.5l_cs (x2), fizhi-cs-32x32x40,
97        fizhi-cs-aqualev20, fizhi-gridalt-hs, tutorial_held_suarez_cs, cpl_aim+ocn.
98    o pkg/seaice: fix a bug in seaice_prepare_ridging.F, forgot a sqrt
99    o pkg/ptracers: compute ptracer source/sink due to linear free surface.
100        not used right now, but can be used in BGC packages to fix conservation.
101    o pkg/dic, pkg/gchem: move tracer_meanarea to gchem, so it can be used by
102        other packages too.
103    o pkg/seaice/seaice_calc_ice_strength.F: replace SEAICE_rhoSnow by
104        SEAICE_rhoIce in the Rothrock (1975)-type strength calculations (bug fix)
105    o model/src:
106      - since we have several heat capacity (even in model/src), define "recip_Cp"
107        as local variable (no longer in common block) to avoid any confusion.
108      - add thetaConst in PARAMS.h (to be used with pStar or sigma-P coords).
109    o pkg/exf:
110      - from Jeff Scott: add option (#define EXF_CALC_ATMRHO) to calculate local
111        air density as function of air Temp, Humidity and atm pressure.
112    
113    checkpoint64w (2014/05/02)
114    o verification:
115      - update output_adm.txt files for experiment global_ocean.cs32x15 (seaice &
116        seaice_dynmix) and lab_sea (primary test & noseaicedyn): moving/changing
117        code from seaice_growth.F & seaice_advdiff.F to new ridging routines
118        seaice_reg_ridge.F seaice_do_ridging.F affects AD gradient
119        (down to 8-7 matching digits for global_ocean and 5-6 for lab_sea).
120    o model/src: for atmosphere in p* coords:
121      - store factor rStarFacC**atm_kappa (= pStarFacK) in common block
122      - implement more accurate p* coordinate slope term (in calc_grad_phi_hyd.F),
123        but not yet activated.
124    o pkg/seaice: fix restart problem for ITD with new ice strength parameterisation:
125      - move computation of opnWtrFrac = 1-AREA to before computing ice strength
126      - this changes the results (but not for any verification experiment, because
127        none tests this), but avoids having an extra field in the restart files
128    o pkg/seaice: remove global variable TICE (not needed)
129      - replace by TICES or local variable where necessary (cost_seaice_test,
130        seaice_diagnostics_state)
131      - move copying TICE to TICES from seaice_check_pickup to seaice_read_pickup
132      - use specific read/write routines for reading and writing individual leves
133        of a multilevel field
134    o pkg/ctrl & model/inc:
135      - move additional anomaly fields of control vars (related to options:
136         ALLOW_KAPGM_CONTROL, ALLOW_KAPREDI_CONTROL and ALLOW_BOTTOMDRAG_CONTROL)
137        from DYNVARS.h into new header file: CTRL_FIELDS.h
138    o pkg/seaice: some small modifications to the ridging code
139      - make fw2ObyRidge a global variable and add the freshwater flux to
140        seaice_growth
141      - pass actual thickness from seaice_prepare_ridging, so that it is
142        computed only in one place and compute it in the same way as in
143        seaice_growth
144      - initialise more fields in seaice_do_ridging that where cause floating
145        point exceptions over land
146    o pkg/seaice: preparations for ridging code and first real ridging code
147      - move most of "PART 1" of S/R seaice_growth into a new and separate
148        S/R seaice_reg_ridge that will both regularize the state variable after
149        advection call the actual ridging routine;
150      - requires to make 5 fields global fields: d_HEFFbyNEG,d_HSNWbyNEG,
151        d_AREAbyRLX,d_HEFFbyRLX,saltFluxAdjust;
152      - remove ridging code from seaice_advdiff
153      - clean up remaining part 1 of seaice_growth in an attempt to make it
154        more readable (this is probably subjective)
155      - call seaice_itd_redist once in the initialisation phase to have the initial
156        conditions in the correct category
157      - finally: add new ridging routine for ITD code
158    
159    checkpoint64v (2014/04/11)
160    o pkg/ecco:
161      - add runtime parameters mdt[start,end]date[1,2]
162      - rename variable "tpmean" as "mdt" and "topexmeanfile" as "mdtdatfile"
163    o pkg/seaice:
164      - make deltaC a global variable, so that it can be used in ridging code
165      - first step towards ridging: introduce new routines that compute
166        ice strength based on ridging parameters
167      - add diagnostic for HSNOWITD
168      - add opnWtrFrac as a global variable (required for ridging scheme)
169    o pkg/salt_plume + pkg/seaice/seaice_growth.F
170      - add SALT_PLUME_IN_LEADS flag to activate pkg/salt_plume only when ice
171        AREA exceeds a threshold (code contributed by I.Fenty in 2011)
172    o verification/hs94.1x64x5:
173      - leave ALLOW_AUTODIFF_TAMC undefined in code_oad/AUTODIFF_OPTIONS.h
174      - do not use ECCO_CPPOPTIONS.h anymore for TAF & OpenAD Adjoint tests.
175    o autodiff:
176      - Start to include explicitly AUTODIFF_OPTIONS.h, COST_OPTIONS.h,
177        and CTRL_OPTIONS.h in src files (to enable to skip the ECCO_CPPOPTIONS.h)
178        For now, only in pkgs used in verification/hs94.1x64x5.
179      - Replace ALLOW_AUTODIFF_TAMC by ALLOW_AUTODIFF (except for tape/storage
180        which are specific to TAF/TAMC).
181      - Replace ALLOW_AUTODIFF_OPENAD by ALLOW_OPENAD:
182        since ALLOW_OPENAD is defined in PACKAGES_CONFIG.h (any time pkg/openad
183        is compiled), this simplifies/reduces which *_OPTIONS.h file to include.
184      - OpenAD: change Makefile to generate the appropriate AD_CONFIG.h when
185        building mitgcmuv_ad (make adAll).
186    o pkg/ecco:
187      - cost_gencost_seaicev4.F:
188         reformulate proxy cost as function of sst and local param SEAICE_freeze
189         change cost names from siv4-[sst,vol] to siv4-[deconc,exconc]
190      - add stop statement if old names siv4-[sst,vol] are used
191      - add snapshots + cost reports for atmospheric longwave
192      - remove factor 0.1 in front of wtau[u,v]
193      - replace never used flag ALLOW_TMI_COST_CONTRIBUTION with
194         ALLOW_TMI_SST_COST_CONTRIBUTION
195    o pkg/atm_phys:
196      - add some diagnostics for Short-wave and for convection S/R.
197    o pkg/seaice/seaice_growth:
198      - fix a bug in the lateral melt parameterization for SEAICE_ITD
199        in seaice_growth; tempFrz needs to be computed inside the i/j-loop
200      - rearrange loop order for better vector code in the same code block
201      - rewrite the floeDiameter formula to have fewer divisions
202      -  (miraculously) none of this changes the results of seaice_itd
203    o pkg/smooth:
204      - move param smooth2Ddiffnbt to pkg/ecco
205    o pkg/salt_plume:
206      - salt_plume_frac.F:
207         clean up local param names
208         reformulate (1-frac) to frac; original 1-frac was based on swfrac code
209    o pkg/rw:
210      - add RS version of S/R READ_MFLDS_LEV.
211    o pkg/streamice:
212      - fix the restart test by writing/reading pickup_streamice files.
213    o pkg/seaice:
214      - refine checks for jfnk-related parameters
215      - comment out if statement with SEAICEuseIMEX parameter as long as there
216        is no IMEX code to avoid confusion and surprises
217    
218    checkpoint64u (2014/03/07)
219    o verification/MLAdjust:
220      - rename secondary test suffix:
221         1.leith  -> A4FlxF  + switch to flux-form & add side-drag
222         0.leithD -> AhFlxF  + switch to flux-form
223         0.smag   -> AhVrDv  (use Vorticity-Div. visc form) + add side-drag
224         0.leith  -> AhStTn  (use Strain-Tension visc form) + add side-drag
225    o pkg/mom_vecinv:
226      - remove unused arguments from S/R MOM_VI_HDISSIP & MOM_HDISSIP (as it used
227        to be before mom_calc_visc.F)
228      - skip the call to MOM_CALC_TENSION & MOM_CALC_STRAIN if not needed.
229      - add 2nd copy of vort3 & strain that knows about lateral BC (free/no slip):
230        this is simpler for the adjoint and for diagnostics.
231    o model/src + pkg/momentum:
232      - add input files for  2-D additional Western and Southern cell-edge wall
233         (e.g., to add "thin-wall" where it is =1).
234      - fix sideDrag option for thin-walls with Non-Lin Free-Surf
235        using 2nd hFacZ that is computed from initial (fix domain) hFac
236      - known limitations of thin wall implementation: not yet coded for sigma
237        coords; not working with CD-scheme neither.
238    o pkg/seaice:
239      - rename recently introduced parameters from (wrong name) AB to
240        BDF (backward difference time discretisation)
241      - retire runtime parameters SEAICEuseAB2, SEAICE_abEps
242    o tools/adjoint_options:
243      - remove setting of "OPENAD=true" (now set by command line "genmake2 -oad")
244    
245    checkpoint64t (2014/02/01)
246    o pkg/ecco
247      - add compile flag ALLOW_GENCOST_TIMEVARY_WEIGHT in ECCO_OPTIONS.h and
248        runtime logical flag gencost_timevaryweight to read daily sea-ice sigma
249      - replace hard-coded record number and period which only work for monthly
250        with read-in runtime param that works for other periods (month,day,year)
251    o pkg/thsice
252      - corrected store directive fixes thsice adjoint gradient checks
253    o genmake2 & testreport:
254      - add genmake2 option "-oad" to generate a Makefile for OpenAD built and
255        in this case, use tools/adjoint_options/adjoint_oad as default AD optfile.
256        Update testreport accordingly. No need to set "export OPENAD=true"
257        in adjoint_options/adjoint_oad anymore (will be removed later).
258    o pkg/shelfice:
259      - put T & S shelfice forcing into model surfaceForcing_T,_S (for now,
260        only if SHELFICEboundaryLayer=F); Change model/src/external_forcing.F
261        so that surface forcing applies to k=kSurf if using pkg/shelfice.
262    o pkg/ctrl, pkg/grdchk, pkg/ecco:
263      - Changes to make GENARR more modular (i.e. invisible if not used)
264    o model/src/temp_integrate.F and salt_integrate.F
265      - re-include store directives that were removed on Dec 27th.
266    o model/src/external_forcing_surf.F
267      - add storage directives for EmPmR and PmEpR to avoid recomputation of
268        seaice_model in the adjoint that is problematic when seaice_model is
269        turned off in the adjoint. The resolved issue dated back to July when
270        EXTERNAL_FORCING_SURF call was moved.
271        Update global_ocean.cs32x15/*/output_adm.seaice_dynmix.txt accordingly.
272    o pkg/exf/exf_check.F
273      - fix a small bug again (originally found by Wentao Liu, now found again)
274      - streamline error output
275    o pkg/streamice and aim_v23:
276      - replace a few "_EXCH_XY_RL (" with "_EXCH_XY_RL(", so that
277        Apple's /usr/bin/cpp knows what to do:
278        Apple LLVM version 5.0 (clang-500.2.79) (based on LLVM 3.3svn)
279        Target: x86_64-apple-darwin13.0.0
280    
281    checkpoint64s (2014/01/05)
282    o tools/adjoint_options:
283      - switch to fastopt.net TAF server (instead of default fastopt.de) for
284        AD optfiles adjoint_diva & adjoint_f95 (similar to adjoint_default).
285    o model/src:
286      - remove some unnecessary TAF storage directives in forward_step.F and
287        in dynamics.F (note: not always removed by TAF, e.g., double storage
288        of kappaRU,kappaRV in former version of dynamics.F).
289    o model/src & pkg/ptracers:
290      - call directly GAD_IMPLICIT_R & IMPLDIFF + DWNSLP_APPLY from temp_integrate.F
291        salt_integrate.F and ptracers_integrate.F (instead of from thermodynamics.F)
292        and remove ptracers_dwnslp_apply.F & ptracers_implicit.F
293      - call MATRIX_STORE_TENDENCY_EXP for each tracer (in ptracers_integrate.F)
294      - add fZon & fMer as output argument of S/R GAD_CALC_RHS.
295      - update TAF AD-output of experiments global_ocean.90x40x15 (13 for AD-grad).
296    o useSingleCpuInput separates single-CPU input from single-CPU output;
297        for now only implemented for in exf_interp_read.F
298    o verification/advect_xy:
299      - fix missing EXCH call in local version of ini_salt.F and update output.
300    o momentum diagnostics:
301      - fix which deltaT is used in Total-Tendency diagnostics.
302      - fix typo in diagnostics for momentum tendency from AB.
303    o pkg/seaice/seaice_lsr.F: fix an initialisation bug for the
304        SEAICE_VECTORIZE_LSR_ZEBRA option that lead to NaNs on a NEC-SX8
305    
306    checkpoint64r (2013/12/08)
307    o verification/*/results:
308      - Update reference output of experiments that were passing with just 13 digits
309         (make it easier to detect future changes in results)
310    o model/src & pkg/ptracers:
311      - move calls to GAD_ADVECTION (Multi-Dim advection) from thermodynamics.F to
312        inside temp_integrate.F and salt_integrate.F and from ptracers_advection.F
313        to inside ptracers_integrate.F ; remove ptracers_advection.F
314    o coupling interface pkgs:
315        [atm_compon_interf, ocn_compon_interf & atm_ocn_coupler]
316      - change the config check to test the agreement of Ocean Depth
317        with ATM land/sea mask (instead of ATM orography);
318      - send to coupler the number of coupler time-steps to do for this run:
319        this allows to check for consistency and don't need anymore to be
320        specified in parameter file "data";
321      in coupler:
322      - check consistency of RunOff map versus Ocean Depth and ATM land/sea mask.
323      - improve report to log-file of config setting (pass log-file unit to few S/R)
324    o pkg/compon_communic:
325      - add coupling send/receive for non-tiled integer vector;
326      - add S/R to check and stop if any proc has an error (similar to
327        eesupp/src/stop_if_error.F but single thread and across Comm_World).
328    o pkg/exch2:
329      - in exch2*_get_rx[1,2]: use new list "W2_tileIndex" to replace search
330        through list of all tiles my proc owns.
331    o model/src:
332      - AUTODIFF only: in calc_adv_flow.F, recompute rTransKp (instead of copying
333        from rTrans) to simplify dependencies (no need to store rTrans anymore in
334        calling S/R); this affects AD results at machine truncation level:
335        update AD-output (TAF) of experiments: 1D_ocean_ice_column,
336        bottom_ctrl_5x5 & obcs_ctrl (12 digits for AD grad, same FWD grad).
337    o tools/build_options:
338      - simplify Warnings options (-ieee & -devel cases) in linux_amd64_gfortran
339    o coupling interface pkgs:
340        [compon_communic, atm_compon_interf & ocn_compon_interf]
341      - move declaration of local array "ibuf" out of header file "CPLR_SIG.h"
342        to each S/R that needs it; remove un-used variables.
343      - change component tile registration (cpl_register.F in both compon_interf
344        and mitcomponent_tile_register.F) to pass the full list of tiles
345        information (including bi,bj) in 1 vector (instead of in a bi,bj matrix).
346      - with pkg/exch2: adopt the same convention used for I/O to define mapping
347        towards coupler global array (unless W2_useE2ioLayOut=F):
348         a) can now use different tiling between atm & ocn comp.
349         b) no needs to modify the run-Off map when changing tiling (was needed
350            even if using the same tiling between atm & ocn).
351         c) will allow to run with blank-tiles in 1 component (e.g., for ocean
352            full land tiles); not yet tested.
353      - check for valid mapping from comp. tile to coupler global array
354      - improve report of tile registration results to log file.
355    o pkg/exch2:
356      - in eesupp/src/eeboot.F: move call to W2_EEBOOT before Coupler
357        tile-registration (CPL_REGISTER);
358      - activate (un-comment) definition and setting of "W2_tileIndex" list to
359        store the local bi,bj indices of each tile;
360      - also print bi,bj in local report to log file;
361      - fix pkg/debug/debug_cs_corner_uv.F for the case nSy > 1 (left from
362        Jun 2009 modif of pkg/exch2).
363    o testing on baudelaire:
364      - run 2 tests at a time: no significant slow down (since enough cpu) and
365        get the full set of test done earlier
366      - switch the ifort test to "-devel" ; run the OpenAD test everyday.
367    o model/src/thermodynamics.F and routines called from it:
368      - pkg/gmredi: gmredi_residual_flow.F (replacing gmredi_calc_uvflow.F &
369        gmredi_calc_wflow.F) to update eulerian velocity (input) into residual
370        velocity (outp) by adding bolus velocity.
371      - model/src: calc_adv_flow.F (replacing calc_common_factors.F) to compute
372        mass transport from velocity + lateral grid cell area
373      - rename calc_gt.F to temp_integrate.F (includes call to TIMESTEP_TRACER)
374      - rename calc_gs.F to salt_integrate.F (includes call to TIMESTEP_TRACER)
375     in thermodynamics.F
376      - remove k from thermodynamics.F and move it to temp_integrate.F
377        salt_integrate.F and ptracers_integrate.F (now done inside the tracer
378        loop).
379      - compute locally (in thermodynamics.F) 3-D velocity field that is used to
380        advect tracers; pass it as argument to GAD_ADVECTION, GAD_SOM_ADVECT,
381        PTRACERS_ADVECTION, TEMP_INTEGRATE, SALT_INTEGRATE, PTRACERS_INTEGRATE,
382        GAD_IMPLICIT_R and PTRACERS_IMPLICIT
383      - update AD-output (TAF & OAD) of experiments: 1D_ocean_ice_column,
384        global_ocean.90x40x15, tutorial_global_oce_biogeo and
385        tutorial_tracer_adjsens (+ TLM output).
386    o tools/OAD_support and (verification/*/) code_oad:
387      - OAD_support/cb2mGetModules.csh: add PACKAGES_CONFIG.h in header file to
388        module conversion ; remove unneeded CPP_EEOPTIONS.h.
389      - code_oad/CPP_OPTIONS.h: remove PACKAGES_CONFIG.h and use a more standard
390        version (closer to model/inc/CPP_OPTIONS.h).
391    o model/src:
392      - remove parameter "useOldFreezing" from PARAMS.h and push it to the
393        retired list of parameters (routine FREEZE has been replace by surface
394        freezing (10.yrs ago) and by pkg/frazil (interior freezing).
395    o pkg/gmredi:
396      - add option (GM_CALC_TRANSP_FROM_VEL) to compute transport directly from
397        velocity (same scheme but affects machine truncation) for GM_AdvForm.
398      - update output of FWD experiments: front_relax.bvp,
399        global_ocean.90x40x15.dwnslp, global_ocean.cs32x15 and ideal_2D_oce.
400    
401    checkpoint64q (2013/11/17)
402    o verification:
403      - generate new output on baudelaire (gfortran, -ieee) for OpenAD experiments
404        tutorial_global_oce_biogeo and global_ocean.90x40x15.
405      - generate new output (growth & melt) for seaice_itd experiment after
406        updating the input pickup_seaice files (previously wrong).
407    o pkg/cost:
408      - define a "final cost function" for each tile, and use GLOBAL_SUM_TILE_RL
409        (instead of _GLOBAL_SUM_RL) to get the final cost from all procs & threads:
410        o this fixes global-sum value if multi-threads (since GLOBAL_SUM_RL
411          cannot be applied to var in common block).
412        o allows to get results independent of tile-per-proc setting (if using
413          slower GLOBAL_SUM_SEND_RECV option).
414      - update accordingly: seaice_cost_final.F, shelfice_cost_final.F,
415        thsice_cost_final.F and ecco_cost_final.F (this one affects lab_sea
416        forward gradients: update output).
417      - update */code_oad/externalDummies.F by adding dummy GLOBAL_SUM_TILE_RL
418        S/R and, for most exp., also comment out GLOBAL_SUM_R8.
419    o pkg/mom_common:
420      - move to pkg/mom_common and model/src (previously in tutorial_deep_convection
421        code) 2nd version of isotropic 3-D Smagorinsky code interface: strain and
422        viscosity are locally declared in dynmics.F and pass as argument to CALC_GW;
423        ensure that all field value that are used are set.
424    o tutorial_deep_convection (smag3d):
425      - add missing scaling in mom_calc_smag_3d.F;
426      - set smag3D_coeff in input.smag3d such that results do not change.
427    o pkg/dic:
428      - change dic_cost.F to compute cost on each tile first and then apply
429        GLOBAL_SUM_TILE_RL: this affects forward gradient of AD experiments
430        tutorial_dic_adjoffline & tutorial_global_oce_biogeo; update output.
431      - remove kLev argument and use GLOBAL_SUM_TILE_RL in tracer_meanarea.
432    o eesupp/src & pkg/mdsio:
433      - remove gather/scatter_vector.F (from eesupp) and mdsio_read/writevector.F
434        (from pkg/mdsio) since they are no longer used for tape IO.
435      - move simple C subroutine "setdir.c" from pkg/compon_communic to eesupp/src
436        (+ keep a copy in pkg/atm_ocn_coupler).
437    o verification/isomip:
438      - split domain in 8 squared tiles (instead of just 2) for FWD, AD & OAD tests.
439    o tools/build_options:
440      - optfile linux_amd64_ifort11, fast (no ieee) case:
441        to get a working restart for all experiments, add "-fp-model source"
442        to FOPTIM (also remove "-opt-streaming-stores auto" <- this is the default).
443    
444    checkpoint64p (2013/10/20)
445    o pkg/autodiff:
446      - adread_adwrite.F: when type and length do not match ALLOW_AUTODIFF_WHTAPEIO
447        criteria, return to simple (default) tape IO format.
448    o pkg/mdsio (+ eesupp/src):
449      - rewrite/simplify tape-IO S/R (new S/R MDS_READ/WRITE_TAPE replace
450        previous MDSREAD/WRITEVECTOR) with 2 array argument from each type (R4/R8);
451      - fix globalFile and singleCpuIO options using simpler global mapping.
452      - new pairs of scatter/gather vector S/R (one for each type, generated from
453        template) corresponding to simpler global mapping.
454      - adread_write.F: switch to new S/R (but leave old calls commented) and write
455        to globalFile when parameter globalFiles is set to True.
456      - switch on useSingleCpuIO in AD experiments: global_ocean.90x40x15 & isomip.
457    o pkg/autodiff:
458      - adread_adwrite.F: add a hack for RS arrays (when RS=real*4);
459        with ALLOW_AUTODIFF_WHTAPEIO defined, since corresponding (real*4) is
460        missing, just add a STOP if array type do not match.
461    o pkg/exf:
462      - to allow to use seaice-fraction from pkg/exf in other pkgs than pkg/seaice:
463         rename CPP option EXF_ALLOW_SEAICE_RELAX to EXF_SEAICE_FRACTION
464         and rename array obsSIce to exf_iceFraction.
465    o tutorial_deep_convection:
466      - add 1rst version of isotropic 3-D Smagorinsky code (from Louis-Philippe),
467        for now all in this exp. code dir.
468      - add secondary test using Smag-3D viscosity (keep very low backgroud visc for
469        testing) and flux-limter advection scheme for Temp. (no explicit diffusion).
470    o eesupp/src:
471      - added SINGLE_DISK_IO option because reading of parameter files on
472        very large number of processors did not work with existing set-up
473    
474    checkpoint64o (2013/09/25)
475    o pkg/atm_phys & atm_compon_interf:
476      - allows for coupling of gray-atmosphere (atm_phys) with ocean.
477    o tools/build_options:
478      - switch from "-fPIC" to "-mcmodel=medium" in optfiles: linux_amd64_gfortran
479        & linux_amd64_ifort11 and in linux_amd64_g77 & linux_amd64_open64 also
480        remove CPP setting (was identical to genmake2 default).
481    o pkg/gmredi & model/src: Some updating and work done on the residual model.
482      - created new file model/src/calc_eddy_stress.F to calculate the eddy stress.
483      - fixed up the calculation of the eddy stress and how it is added to the
484        external forcing via taueddy_external_forcing.F
485      - at the moment, it is very experimental and only works with the variable eddy
486        diffusivity calculated in gmredi_k3d.F
487    o pkg/seaice/seaice_check.F: add a check to catch the case of specifying
488        lw/swfluxfile instead of lw/swdownfile in data.exf together with pkg/seaice
489    o testreport (& tools/do_tst_2+2):
490      - for MPI test, set default COMMAND to 'mpirun -v -np TR_NPROC ./{executable}'
491    o tools/build_options: linux_amd64_gfortran & linux_amd64_ifort11:
492      - remove CPP='cpp -traditional -P' (this is genmake2 default)
493      - add switch (for big setups) to use 'medium' memory-model (but leave it
494        commented out for now). Note: needed with FC 19 to compile big setups.
495    o pkg/atm_compon_interf:
496      - regroup per field origin (model & pkgs) the coupling storage routines
497         (which accumalate in time each coupling field)
498      - move bi,bj loops inside atm_store_my_data.F and store also wind-stress;
499        update accordingly + simplify aim_do_physics.F
500    
501    checkpoint64n (2013/08/26)
502    o pkg/bbl:
503      - update formulation to reduce dispersion of tracer properties in bbl
504      - replaced horizontal and vertical transport parameters with velocities
505      - changes results for verification/global_with_exf.yearly
506    o verification/halfpipe_streamice:
507      - change to matrix construction; update both Fwd & AD output.
508    o verification/tutorial_tracer_adjsens:
509      - OpenAD test: switch to vector-invariant momentum since fluxform with
510        rStar causes problems with OpenAD; update output.
511    o tools:
512      - genmake2: define a common prefix for all temp-file names ; move them
513        to local /tmp disk (but would be easy to move them back);
514      - xmakedepend: apply the same switch for temp-files.
515    o pkg/offline:
516      - changes for use with rStar (NLFS): assume that we continue to provide
517        UVELMASS/hFacW_ini & VVELMASS/hFacS_ini (hFacWi,S_ini as written in grid
518        output files) and apply r* scaling at each time-step (function of the time
519        varying Eta field); use normal pickup files (needed for restart with NLFS).
520    o adjoint with OpenAD:
521      - move OpenAD-specific files from verification/OpenAD/code_oad_all to a new
522        package openad.  This package has to be compiled for all runs with OpenAD
523        instead of specifying code_oad_all as a modification directory.
524      - link the fortran templates from tools/OAD_support to the build directory
525        for preprocessing to avoid conflicts between different setups compiled
526        simultaneously and simplify cleanup.
527    o pkg/exf:
528      - in exf_getforcing.F: comment out most of the EXCH calls (fields used by
529        main model or other pkgs are exchanged in EXF_MAPFIELDS).
530    o model/src:
531      - do_stagger_fields_exchanges.F: skip EXCH calls when using pkg/offline
532        (similar to do_fields_blocking_exchanges.F).
533    o tools/genmake2:
534      - remove TAF output src code just before calling staf (when not doing a fresh
535        built, will stop if things go wrong with staf).
536      - fix target "adobj" (left from TAF version-2 transition)
537    o pkg/diagnostics:
538      - add a parameter to track the status of pkg/diagnostics activation;
539        this replace/extend the use of logical param "settingDiags" (to check
540        when adding diag to the list) to also check any DIAGNOSTICS_[]FILL* call.
541        New S/R to print error message and stop when pkgStatus is not right.
542      - to try: add a warning when "data.diagnostics" exist but useDiagnostics=F.
543    
544    checkpoint64m (2013/08/13)
545    o model/inc:
546      - move FCORI_MACROS.h from PARAMS(_MACROS).h to GRID_MACROS.h (left from
547        8 yrs ago modif when Coriolis arrays were moved from PARAMS.h to GRID.h)
548    o pkg/kpp & aim_v23:
549      - only call DIAGNOSTICS_FILL when useDiagnostics=TRUE.
550    o pkg/seaice:
551      - restore (last time on 2011/12/24) ability to compile pkg/seaice with
552        undefined pkg/exf ALLOW_ATM_TEMP or ALLOW_RUNOFF CPP-options.
553      - only call DIAGNOSTICS_FILL when useDiagnostics=TRUE.
554    o eesupp/src:
555      - fix check for overlap-size argument (myOLn.GT.MAX_OLY_EXCH)
556        in exch1_rx[_ad].template (thanks to Francois Lefeuvre report).
557    o switch to TAF version-2:
558      - add flow directives to "freeze" naming convention for ad-routines that
559        are used in I/O routines of adjoint state and associated common blocks
560        that are used in adcommon.h (backward compatible with version-1)
561      - update all ad-optfiles (in tools/adjoint_options + few local to some
562        verification exp) to switch to latest TAF version (currently 2.5.5);
563        merge "*_mpi" adoptfiles into the corresponding non-mpi (but commented
564        out except for diva).
565      - the major change with TAF version-2 is that both the differentiated
566        routines and the forward routines (new) are returned as part of the TAF
567        output source code. Update genmake2 accordingly.
568      - new output for AD experiment tutorial_global_oce_biogeo (better AD gradient
569        with version-2 when compared to TLM or FWD gradient); other experiments
570        output show minor changes (machine truncation level); update the few AD &
571        TLM output with 12 matching digits only (some prior to version-2 switch).
572    
573    checkpoint64l (2013/08/05)
574    o model/src:
575      - in dynamics.F, skip the call to CALC_VISCOSITY if momViscosity=F
576      - main_do_loop.F: no longer includes additional (AUTODIFF) header files
577        when using TAF but only if using OpenAD.
578    o pkg/mom_common+fluxform+vecinv:
579      - always set horiz. viscosity arrays to background value before calling
580        MOM_CALC_VISC (in MOM_FLUXFORM & MOM_VECINV) and call S/R MOM_CALC_VISC
581        only when using variable horiz. viscosity (useVariableVisc=T);
582      - simplify mom_vecinv.F (only 1 block for momViscosity).
583    o pkg/mom_common: variable horiz. viscosity in NH code
584      - move horiz-viscosity 3-D arrays viscAh_W & viscA4_W (used in calc_gw.F)
585        from model/inc/NH_VARS.h to pkg/mom_common/MOM_VISC.h
586      - store in common block (in MOM_VISC.h) logical parameters: useHarmonicVisc,
587        useBiharmonicVisc & useVariableVisc.
588      - fix old bug (already in 1rst implementation of variable horizontal viscosity
589        in calc_gw.F, version 1.27, 2006/06/20):
590        horizontal bi-harmonic viscosity in gW was missing when wiscA4W=0
591        even if variable bi-harmonic viscosity (e.g., viscA4Grid) was non-zero.
592      - change background horiz viscosity in vertical momentum (gW) to always be
593        viscAhW,viscA4W (instead of switching to viscAhD, viscA4D when using
594        any type of variable viscosity, e.g., viscAh/4Grid)
595      - change default value of viscAhW,viscA4W to viscAhD,viscA4D (instead of
596        viscAh,viscA4) which default to viscAh,viscA4 :  this avoid missing
597        horiz. viscosity in W eq when just setting viscAhD/Z,viscA4D/Z (since
598        default viscAh,viscA4 is zero).
599    o model/src:
600      - extend length of drC & recip_drC from Nr to Nr+1.
601    o Divided Adjoint (DIVA) with MPI:
602      - genmake2: do not create local ./mpi_headers dir when no file to copy into.
603      - remove option ALLOW_DIVIDED_ADJOINT_MPI (replaced by ALLOW_DIVIDED_ADJOINT
604        and ALLOW_USE_MPI) from model/src code and from ECCO_CPPOPTIONS.h.
605    
606    checkpoint64k (2013/07/23)
607    o tools/OAD_support:
608      - avoid the need for file "oadtempflile" in insertTemplateDir.bash
609        and remove it from genmake2.
610    o tools/build_options:
611      - remove MPI_HEADER_FILES_INC setting (no longer used) and comment out
612        MPI_HEADER_FILES list (let genmake2 find the missing headers)
613    o verification/halfpipe_streamice:
614      - change to driving stress formulation; update both Fwd & AD output.
615    o tools/genmake2 and Divided Adjoint with MPI:
616      - remove default setting of MPIINCLUDEDIR from $MPIHOME (always been broken);
617      - only create local ./mpi_headers dir if using MPI and MPIINCLUDEDIR is set;
618      - add new variable (LOCAL_MPI_HEADERS, replace MPI_HEADER_FILES_INC) in
619        Makefile, which is set by genmake2, and check which mpi-headers are needed.
620        Note: need to update adjoint_options/adjoint_diva_mpi accordingly.
621      - update target "CLEAN" (to also remove ./mpi_headers/ dir).
622    o testreport:
623      - changes for Divided Adjoint Run: when file "run_ADM_DIVA" is present,
624        perform additional runs (number of additional runs is "add_DIVA_runs" taken
625        from file run_ADM_DIVA ; should match "nchklev_3" in file "code_ad/tamc.h").
626    o verification/lab_sea (AD test):
627      - change to test Divided Adjoint Run in this set-up (both standard AD test
628        and 2 secondary AD test: noseaice + noseaicedyn); drop the evp AD test.
629    o pkg atm_ & ocn_compon_interf:
630      - remove cpl_earlyExpImpCall=F option (was mainly a hack added for the old
631        staggerTimeStep implementation) and disable (retired) corresponding
632        parameter.
633    o pkg/offline:
634      - comment out unused field-array: ICEM, ConvectCount, hflx[01], icem[01];
635      - add S/R OFFLINE_CHECK to print offline parameter summary
636        and stop if loading a file that is not used;
637      - enable to calculate GMRedi tensor, KPP mixing and convective index
638        when the corresponding file(s) is/are not specified (but pkg is used).
639      - comment out resetting of monitorFreq (to endTime) since it still useful
640        for monitoring CFL numbers and checking offline fields.
641    o OpenAD:
642      - do not define OAD_DEBUG (in tools/adjoint_options/adjoint_oad) to get
643        smaller standard output) and generate new (smaller) output (on baudelaire
644        using gfortran with -ieee) for exp. tutorial_tracer_adjsens
645    o model/src:
646      - call EXTERNAL_FORCING_SURF outside bi,bj loops (in do_oceanic_phys.F),
647        (where external_forcing_adjust was called) and change list of arguments;
648      - merge external_forcing_adjust.F into external_forcing_surf.F
649      - put surface relaxation of SST & SSS in new S/R (forcing_surf_relax.F)
650        which also contains balancing surface relaxation (ALLOW_BALANCE_RELAX,
651        previously in file "balance_relax.F", now removed); remove corresponding
652        arrays in FFIELDS.h.
653      - EXTERNAL_FORCING_SURF is now called when using pkg/offline; implement
654        a hack for getting back surfaceForcingS (loaded from file) when using
655        ALLOW_OLD_VIRTUALFLUX (pkg/dic, in DIC_OPTIONS.h) in offline set-up.
656      - update output of AD test experiment global_ocean.cs32x15.seaice_dynmix
657        (was getting only 11 digits for AD-grad, due to different truncation
658         error in AD code).
659    
660    checkpoint64j (2013/07/04)
661    o pkg/offline:
662      - to enable the use of other features with pkg offline (e.g., obcs),
663        also call S/R DO_OCEANIC_PHYS in offline run - more modifs to come.
664      - comment out "allowFreezing" (now active) in 2 offline testing exp.
665    o pkg/exf ALLOW_RUNOFTEMP
666      - fix sign and make compatible with temp_EvPrRn .NE. UNSET_RL
667      - ahdded RUNOFTEMP test to verification/global_ocean.cs32x15/input.seaice
668      - changes verification/global_ocean.cs32x15/results/output.seaice.txt
669    o pkg/atm_phys:
670      - allows to add some damping on uVel,vVel (strength only function of k).
671    o verification, ad-tests:
672      - reduce length of the simulation for AD test: global_ocean.cs32x15.seaice,
673        global_ocean_ebm and offline_exf_seaice.thsice.
674    o verification/halfpipe_streamice:
675      - found a bug in s/r pkg/streamice/streamice_advect_thickness which caused
676        flux input boundary conditions to be ignored; a fix for this bug changed
677        verification results for forward test.
678      - add adjoint test.
679    o pkg/ptracers:
680      - allow to restart without SOM pickup (resetting SOM to zero).
681    o pkg/gmredi:
682      - New parameterisation for the calculation of K and (GM_PsiX,GM_PsiY) added.
683      - It uses mixing length theory and suppression by steering level effects to
684        calculate a diffusivity used in the Redi tensor (see Bates et al., submitted
685        to JPO, for more details).
686      - It also uses the diffusivity to calculate an eddy induced streamfunction
687        based on a PV closure (rather than a buoyancy closure as in GM).
688      - New files added to pkg/gmredi:
689        gmredi_k3d.F calculates the diffusivity for the Redi tensor and the eddy
690          induced overturning streamfunction (called from gmredi_calc_tensor)
691        gmredi_calc_eigs.F calculates the eigenvalues and eigenvectors required
692          for use in gmredi_k3d
693        gmredi_calc_urms.F calculates the the rms eddy velocity used to calculate
694          the eddy diffusivity
695    
696    checkpoint64i (2013/06/18)
697    o pkg/cheapaml:
698      - take src files from MITgcm_contrib/verification_other/offline_cheapaml/code
699        that enable the use of seaice (pkg/thsice thermo & pkg/seaice dynamics)
700        with pkg/cheapaml.
701    o utils/python: add subpackage for cubed sphere function, first example is
702        simple cs.pcol to plot 2D fields on a map
703    o pkg/streamice:
704      - This is a new package in MITgcm, which solves dynamic equations for ice
705         sheets and ice streams
706      - At the moment it is independent of atmosphere/ocean dynamics, or any
707        other physical packages in MITgcm, even ShelfIce (although it is anticipat-
708        ed that this will change)
709      - There is a new verification exercise that tests it, halfpipe_streamice,
710        evolves a 40x80 km ice stream/ice shelf (at coarse resolution) in a
711        rectangular embayment
712      - halfpipe_streamice is forward-only for now, will have TAF test soon
713      - for questions, email dgoldber@mit.edu
714    o pkg/thsice:
715      - fix filling of snow-precip diagnostic (SIsnwPrc) (left from the removing
716        of snowPrc from common bloc, on May 2nd, thsice_step_fwd.F revision 1.40);
717      - remove diagnostics 'SI_FrcFx' ( now identical to 'SI_Fract' with new
718         call sequence: ice-Thermo(thsice),ice-Dyn,ice-Advect );
719      - allow precip from pkg/cheapaml to make snow over seaice.
720    o pkg/cfc:
721      - move in separated header file (CFC_ATMOS.h) time-serie arrays of Atmos CFC
722      - remove atmos CFC 2-D arrays from common bloc (now local variables)
723      - change some hard-coded constant to run time parameter for reading and
724        and interpolating atmospheric CFC ; keep same default as previously.
725      - more flexible reading of Atmos CFC (number of records no longer hard-coded)
726      - put arguments in standard order: myTime, myIter, myThid
727      - use standard S/R GET_PERIODIC_INTERVAL for time interpolation of
728        atmospheric CFC ; fix for the case myIter < PTRACERS_iter0.
729      - fix calculation of surf tendencys for partial-cell and/or Non-Lin Free-Surf
730      - new routine to print summary of CFC params and check for consistency.
731    o pkg/seaice: add uice/vice to the parameter list of s/r seaice_advdiff,
732      similar to thsice_advdiff.
733      Note that for the B-grid case the routine parameters uc/vc are overwritten.
734    o pkg/cfc:
735      - change time interpolation for atmospheric CFC-11 (now same as for CFC-12);
736        update output (CFC-11) of experiment cfc_example and tutorial_cfc_offline.
737      - fix interpolation of atmphspheric CFC between Sorthern and Northern
738        hemisphere value.
739    o pkg/seaice:
740      - allow to use seaice dynamics without pkg/exf (e.g., using pkg/cheapaml)
741    o pkg/seaice: add Adams-Bashforth2 time discretization for ice dynamics;
742                  so far only for JFNK-solver, turn on with SEAICEuseAB2 = .TRUE.
743      - requires an additional time level that is stored in the pickup
744      - you can start with AB2 from a pickup without this time level with
745        pickupStrictlyMatch = .false.
746    
747    checkpoint64h (2013/05/26)
748    o pkg/cheapaml:
749      - replace one ()**.5 with SQRT + use temp var "cdu = tau/Du" (coare3_flux);
750        affect machine truncation (but exp cheapAML_box still pass with 13 digits).
751      - consistent treatment of uWind,vWind location: assume everywhere in
752        pkg/cheapaml that they are on C-grid (@ uVel,vVel location)
753        (Note: was already the case for Tair,Qair advection);
754        add option (INCONSISTENT_WIND_LOCATION) to recover old results.
755      - compute diabatic forcing only in the interior (excluding overlap regions)
756      - fill in array "Qsw" (for short-wave heating);
757      - add diagnostic for precip ; fill-in diagnostics for cheapAML state vars
758        before any update (consistent with MITgcm state var diagnostics);
759      - cheapaml_fields_load.F & cheapaml_init_varia.F : only fill outer edges if
760        non periodic ; + adjusted for uWind,vWind (C-grid).
761    o pkg/seaice:
762      - remove dependence of pkg/seaice versus pkg/exf (file: pkg/pkg_depend) and
763        add explicitly pkg/exf in experiments that currently use pkg/seaice.
764      - allows to compile pkg/seaice without pkg/exf.
765    o tools/OAD_support:
766      - use a local copy (ourRWCPinBd) to ensure index ourRWCP stays within bounds
767        of array "ourStepOf" in function rvNextAction (revolve.F90).
768    o pkg/aim_compon_interf renamed to pkg/atm_compon_interf
769      - change pkg name (to be used with other atmospheric physics pkg)
770        but by-pass CVS when making the new pkg so that we keep the CVS history.
771    o pkg/offline:
772      - call OFFLINE_FIELDS_LOAD from pkg/offline/offline_init_varia.F
773        (previously called from pkg/dic_surfforcing_init.F)
774    o model/src:
775      - re-order sequence of calls within packages_init_variables.F, using
776        standard order, except for:
777         DIAGNOSTICS_INIT_VARIA and OFFLINE_INIT_VARIA (called @ the top),
778         CTRL_MAP_INI and OBCS_INIT_VARIABLES (called @ the bottom)
779    o new pkg: atm_phys
780      - add entry point and source code (previously in contrib
781          verification_other/atm_gray/code_gray/) for new pkg "atm_phys"
782        (atmospheric physics pkg from P. O'Gorman and T. Schneider, JCl, 2008).
783      - re-order sequence of calls within packages_readparms.F,
784         packages_init_fixed.F and packages_write_pickup.F
785        (but skip packages_init_variables.F because of TAF issues).
786      - update contrib/*/atm_gray/ experiment.
787    o pkg/seaice (ITD code):
788      - simplify and fix the writing of pickup-files for ITD case
789    o pkg/seaice (ITD code):
790      - introduce floe size dependent lateral melt for ITD case;
791      - cleaning (avoid division by zero, change to f77 syntax, remove tabs)
792      - update seaice_itd experiment main output.
793    o pkg/thsice:
794      - remove snowPrc from common block THSICE_FLUX, in THSICE_VARS.h and
795        declared it locally (without bi,bj) in thsice_main.F & aim_do_physics.F;
796        changed accordingly pkg/atm2d.
797    
798    checkpoint64g (2013/05/01)
799    o pkg/exf ALLOW_RUNOFTEMP
800      - if defined,  river and glacier runoff temperature can be read-in from
801        runoftempfile; runoftempfile needs to have same format as runofffile
802    o pkg/thsice (& pkg/bulk_force):
803      - In case Energy-Reference-Level (ERL) is used (temp_EvPrRn=0), account
804        for energy content of Precip + RunOff & Evap (like with pkg/exf).
805      - change computation domain to always be the interior (without overlap).
806      - fix bug (wrong overlap) when using EXF and stressReduction > 0.
807    o verification/tutorial_global_oce_biogeo:
808      - add preliminary OpenAD setup (still has problems)
809    o pkg/seaice:
810      - turn some parameters for choosing the convergence criterion of
811        the inexact Newton method (JFNK) into runtime parameters for convenience
812      - clean up seaice_calc_viscosities.F: remove/retire two options for
813        SEAICetaZmethod (1,2), introduce area weighting for method 3; this changes
814        the results of lab_sea/input because a division has been replaced by
815        a multiplication with the reciprocal (sumNorm)
816      - modify offline_exf_seaice/input.dyn_jfnk/data.seaice to improve
817        convergence of solver (changes results, of course), and to test some
818        new features
819    o model/src:
820      - remove unused (since Aug 2006) S/R SET_WRITE_GLOBAL_PICKUP and its call.
821      - remove unused forcing field "dQdT" and associated parameter "dQdTFile".
822      - move forcing adjustment (balancing surface forcing) in specific S/R
823        where might also be applied exchanges (if needed) to forcing arrays.
824    o pkg/shelfice:
825      - allow to specify the Ice-Shelf mass per unit area (instead of
826        shelficeLoadAnomaly) and do the conversion from one to the other
827      - move loading from file of shelficeLoadAnomaly from SHELFICE_INIT_VARIA
828        to SHELFICE_INIT_FIXED: (a) avoid some TAF recompution issues and
829        (b) will be loaded in time (before calling INI_PRESSURE) in case it is
830        needed in phi0surf for initial pressure calculation when useDynP_inEos_Zc
831      - update isomip experiment forward output (get machine truncation error
832        from conversion shelficeLoadAnomaly -> shelficeMass -> shelficeLoadAnomaly)
833    o verification/offline_exf_seaice:
834      revise data.seaice in input.thermo and input_ad and update results accordingly
835      HO = 0.2, SEAICE_areaLossFormula = 2,      HO = 0.2, SEAICE_areaLossFormula = 2,
836      SEAICE_growMeltByConv = .TRUE.,      SEAICE_growMeltByConv = .TRUE.,
837      SEAICE_doOpenWaterMelt = .TRUE.,      SEAICE_doOpenWaterMelt = .TRUE.,
838    Improve consistency with input.thsice and    Improve consistency with input.thsice and fix the spurious 500m ice corners
839    fixes the spurious 500m ice corners.  o verification/global_ocean.cs32x15:
840      - set Energy Reference Level to liq. water @ 0 degC: temp_EvPrRn=0.,
841        in icedyn (thsice+seaice-dyn) test and update output.
842    o verification/offline_exf_seaice:
843      - adjust list of tested variables (with testreport) and update
844        standard test output with better slection of monitored fields.
845  o pkg/cd_code:  o pkg/cd_code:
846    - swtich off Adams-Bashforth on the 1rst iteration. Update output of    - swtich off Adams-Bashforth on the 1rst iteration. Update output of
847      experiment lab_sea.hb87      experiment lab_sea.hb87
Legend:
Removed from v.1.1853  
changed lines
  Added in v.1.2010
  ViewVC Help
Powered by ViewVC 1.1.22