/[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.2040 by jmc, Wed Aug 6 23:19:25 2014 UTC revision 1.2120 by jmc, Sat Jan 31 00:37:39 2015 UTC
# Line 4  $Name$ Line 4  $Name$
4      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
5      ==============================      ==============================
6    
7    o pkg/shelfice:
8      - change setting of kTopC to also be zero where there is no ice-shelf.
9    o pkg/seaice: add the option to prescribe a fixed thickness PDF
10      for thermodynamic calculations
11      - because of divisions of the type 7/7 not exactly equal to 1, this change
12        affects results of some verification experiments by up to 4 digits:
13        global_ocean.cs32x15.seaice, lab_sea, lab_sea.fd, lab_sea.salt_plume,
14        seaice_obcs, seaice_obcs.seaiceSponge, seaice_obcs.tides
15        adjoint: lab_sea, lab_sea.noseaicedyn
16        -> updated
17    
18    checkpoint65i (2015/01/23)
19    o model/src:
20      - change units of frictionHeating field from W to W/m^2
21      - fix bug in frictional heating from stratospheric drag in
22        S/R ATM_PHYS_TENDENCY_APPLY_V
23      - with synchronous time-stepping: move resetting to zero of frictionHeating
24        field from load_fields_driver.F to thermodynamics.F ;
25      - add diagnostics for frictional heating.
26    o model/src & gmredi K3D:
27      - add frictionHeating field to pickup-files (for synchronous time-stepping)
28      - move ALLOW_EDDYPSI block out of DYNVARS.h and merge it into FFIELDS.h
29      - rename uMean,vMean --> uEulerMean,vEulerMean (+ change diagnostics name
30        and name in pickup file).
31      - fix arguments of few DIAGNOSTICS_FILL calls in gmredi_k3d.F
32    o model/src:
33      - move (from ini_forcing.F) initialisation of all arrays in header file
34        FFIELDS.h into new S/R ini_ffields.F ; this fix broken restart
35        (since Nov 2012) with option selectAddFluid=1.
36    o pkg/shelfice:
37      - change units (+ document) of newly added field "shelfIceMassDynTendency"
38      - upgrade from _RS to _RL shelficeMass & shelficeLoadAnomaly arrays.
39    o pkg/obcs: re-organize obcs_calc_stevens.F to prepare for passive tracers
40    o pkg/seaice: two new runtime parameters:
41      - SEAICE_cStar replaces the hard wired "20" in the strength formulation,
42        long overdue
43      - SEAICE_tensilFac: preparation for Koenig-Beatty+Holland (2012)
44        parameterization of tensil stress for fast ice (does not do anything yet)
45    o pkg/seaice/seaice_check.F: left over from restricted addtive Schwarz method
46        for LSR: move checking overlaps out of JFNK-if-block
47    o pkg/atm_phys:
48      - move copy of grid and dynamical fields out of atm_phys_driver.F into
49        new S/R ATM_PHYS_DYN2PHYS.
50      - fix level height for Linear Free-Surf case (to be consistent with
51        fixed-in-time pressure levels).
52    o pkg/mom_common (bottomdrag):
53      - average bottomDragFld (ALLOW_BOTTOMDRAG_CONTROL) from grid-cell center
54        to correct location (U & V point) in S/R MOM_U/V_BOTTOMDRAG;
55        update output (adm & tlm) of experiment global_ocean.90x40x15.bottomdrag.
56    o pkg/shelfice (u/v_drag):
57      - add one argument (the other velocity component) to S/R SHELFICE_U/V_DRAG;
58      - use SHELFICEDragLinear & SHELFICEDragQuadratic in shelfice_v_drag.F
59        (now similar to shelfice_u_drag.F)
60      - remove condition on bottomDragTerms when calling SHELFICE_U/V_DRAG
61        from mom_vecinv.F (similar to mom_fluxform calls).
62    o pkg/mom_common (bottomdrag):
63      - add one argument (the other velocity component) to S/R MOM_U/V_BOTTOMDRAG;
64      - add options for alternative discretisation of velocity norm in quadratic
65        bottom drag: averaging of grid-cell-center KE to U.point and V.point
66        (selectBotDragQuadr=0, original code) can produce noisy patterns; this
67        improves when averaging only the other component (e.g. vVel for U bottom
68        drag) using wet-point method (selectBotDragQuadr=2) or not (=1).
69      - add options (bottomVisc_pCell) to account for partial-cell in bottom
70        friction (no-slip BC) thus getting the true distance from the bottom;
71        off by default since this also reduces model stability (explicit bottom
72        friction). Changes produce different machine truncation with linear-drag
73        + no-slip BC: update output of experiments exp4.stevens (10 digits) and
74        front_relax (bvp & mxl, 10 & 11 digits).
75    o model/src & mom pkgs:
76      - add one more level to vertical viscosity local arrays (Nr+1, previously Nr)
77        since no-slip bottom BC uses viscosity @ k+1 to update velocity @ level k
78      - for now and until vertical mixing scheme are updated to fill up level Nr+1,
79        just copy level Nr value to level Nr+1.
80    o doc:
81      - add notes form Abhisek Chakraborty on how to close Heat & Salt budgets.
82    o pkg/streamice:
83      - shelfice_init_varia.F, shelfice_init_fixed.F, shelfice_readparms.F,
84        SHELFICE.h: move array shelficemass to init_varia, initialise
85        from array shelficemassinit, which is set in init_fixed and
86        has option to be set from file; introduce array
87        shelficeMassDynTend array, set thru file in init_fixed
88      - NEW FILE: shelfice_step_icemass.F: update shelficemass from either
89        shelfIceFreshWaterFlux and shelficeMassDynTend or H_streamice in
90        pkg/streamice
91      - shelfice_thermodynamics.F: call shelfice_step_icemass if
92        SHELFICEMassStepping = .true.
93      - shelfice_ad_diff.list: add new subroutine shelfice_step_icemass.f
94    
95    checkpoint65h (2014/12/17)
96    o verification/global_ocean.90x40x15, global_ocean_ebm:
97      - update to generic controls.
98    o pkg/ctrl:
99      - ctrl_map_ini_genarr.F: test for igen>0
100      - ctrl_init_variables.F: add call to CTRL_MAP_INI_GENTIM2D.
101      - move gentim2d part of ctrl_map_ini_genarr.F to
102        new separate routine : ctrl_map_ini_gentim2d.F
103      - ctrl_cost_driver.F,ctrl_init.F ctrl_map_ini_genarr.F:
104        add time variable weights and cyclic controls (time mean,
105        seasonal cycle, etc) capabilities to gentim2d
106    o pkg/ecco:
107      - move reading of error, conversion to weights from
108        ecco_cost_init_fixed.F to ecco_toolbox.F (ecco_readwei)
109    o global_oce_llc90, cs32:
110      - activte sbo.
111      - CTRL_SIZE.h : add maxCtrlProc.
112      - activate the new genctrl capabilities.
113      - update following reorganization of global_oce_input_fields.tar.gz  
114    o pkg/ctrl/ctrl_ad_diff.list:
115      - add missing ctrl_init_rec.f
116    o tools/OAD_support/keepOriginal:
117      - add ctrl_init_rec.F
118    o model/src/the_main_loop.F, pkg/ctrl/ctrl_map_ini_genarr.F:
119      - avoid recomputations of xx_gentim2d_dummy.
120    o momentum diagnostics:
121      - move filling of diagnostics Um_Diss & Vm_Diss from mom_fluxform.F and
122        mom_vecinv.F to timestep.F (to include Smag-3D contribution).
123    o pkg/seaice/seaice_lsr.F:
124      - fix computation of BU/BV, make sure they are non-zero, only necessary
125        for cubed sphere exchanges and SEAICE_OLx/y>0.
126        This makes calling s/r fill_cs_corner_rl obsolete (removed).
127    o verification/global_ocean.cs32x15/input.seaice
128      - add test for strong implicit coupling and restricted addtive Schwarz
129        methods for LSR
130      - update results/output.seaice.txt
131    o verification/tutorial_global_oce_biogeo:
132      - AD & OAD set-up: switch off ALLOW_OLD_VIRTUALFLUX (like in fwd set-up)
133        and set PTRACERS_EvPrRn(1,2) to zero (for DIC & ALK);
134        update output_adm,_oadm,_tlm.txt
135    o pkg/dic:
136      - change default to "#undef ALLOW_OLD_VIRTUALFLUX" in DIC_OPTIONS.h
137      - add a stop when ALLOW_OLD_VIRTUALFLUX is defined and PTRACERS_EvPrRn
138        for tracer 1 & 2 is set (to avoid double counting of dilution effect).
139    o verification/tutorial_global_oce_biogeo:
140      - AD & OAD set-up: do not set PTRACERS_EvPrRn for tracer 1 & 2 since
141        dilution is already added by pkg/dic ALLOW_OLD_VIRTUALFLUX ;
142        update output_adm,_oadm,_tlm.txt
143    o eesupp/src:
144      - add fill_cs_corner_rl.F (_RL-copy of fill_cs_corner_rs.F)
145        and call it in seaice_lsr when using nonzero SEAICE_OLx/y
146    o pkg/seaice: introduce strongly implicit coupling, intended for stabilizing
147      LSR (following Hutchings et al. 2004)
148      - introduce zetaZ as a global variable (requires adjustments in
149        seaice_evp.F and seaice_jfnk.F) and compute analoguously to etaZ
150        in seaice_calc_viscosities.F
151      - new flag SEAICEuseStrImpCpl (default off)
152      - add new terms zetaZ*du/dy and zetaZ*dv/dx on both sides of the
153        momentum equations
154      - JFNK-solver (preconditioner) adjusted to allow using this method
155        (where it does not have any positive effect)
156    o pkg/seaice
157      - add overlap for restrictive additive Schwartz method also to LSR, but
158        turn it off by default (because it only seems to work for non-cubed
159        sphere/llc topologies)
160    o pkg/autodiff
161      - add extra exchanges in addummy_in_steppging also for mnc output
162    o pkg/seaice:
163      - add extra exchanges before dumping adjoint fields in seaice_ad_dump.F
164        to avoid stripes in the adjoint output (similar to addummy_in_stepping.F)
165    o pkg/obcs:
166      - add option to specify domain-connected piece Id along OB though files
167        (for now, only used to diagnose PhiVel; default =1 -> global correction)
168      - modify obcs_diag_balance.F to process each connected-part separately;
169        exclude OB where Id = 0 but keep the full divergence there (-> constant
170        Psi along this portion of OB and all inflow derived from grad(Phi))
171      - add example of OB-conect files in verification/so_box_biogeo/input.
172    o pkg/diagnostics:
173      - improve treatment of Open-Boundaries in Velocity-Potential calculation:
174      a) use non zero matrix element across OB and solve for Phi 1 grid point
175         outside OB (i.e, where tracer OBCS are applied).
176      b) global imbalance of inFlow/outFlow across OB: instead of modifying
177         velocity across OB, spread it uniformly along the OB in the RHS.
178    
179    checkpoint65g (2014/11/19)
180    o pkg/diagnostcis:
181      - add preconditioner off-diagonal factor (diagCG_pcOffDFac) as run-time param
182    o pkg/my82,pp81,kl10:
183      - fix ALLOW_3D_DIFFKR implementation (previously had no effect): use 3-D
184        field diffKr as background diffusivity to set minimum diffusivity value.
185    o pkg/seaice:
186      - add new parameter SEAICEuseEVPrev (for revisited EVP by
187        Bouillon et al 2013)
188      - add SEAICE_evpAlpha/Beta to namelist
189      - modify seaice_evp code to able to use SEAICEuseEVPrev. These changes
190        involve trivial re-ordering of divisions/multiplications that should not
191        have any effects, but ... they do change the results of lab_sea.hb87 so
192        that only 2 digits of agreement remain, in line with the general
193        instability/fragility of the evp-code. For now, I do not care and just
194        update the verification experiment
195    o pkg/ctrl : add generic pre-processing (xx_genarr2d_preproc etc)
196      - remove xx_genarr2d_numsmooth etc
197      - use xx_genarr2d_preproc etc instead
198      - carry out smooth_correl2D at initialization stage rather than during
199        time stepping.
200      - get smoothOpNb from xx_gentim2d_preproc
201    o global_oce_cs32 :
202      - add test of generic ecco capabilities.
203    o pkg/exf/exf_getsurfacefluxes.F :
204      - add to tmpUE/tmpVN rather than re-initializing them (conventional behavior).
205    o pkg/ctrl,ecco :
206      - omit un-necessary arguments to ctrl_cost_gen.F unless ECCO_CTRL_DEPRECATED
207        is defined.
208    o pkg/seaice: add diagnostics for Delta (deformation parameter)
209    o model/src:
210      - when ALLOW_3D_DIFFKR is defined, add a Warning if any of the vertical
211        diffusivity (background) is specified but ignored.
212    o verification:
213      - in testreport, convert hostname to lower case for default output dir name
214    o model, pkg : add useCTRL, useECCO parameters (on/off run time switch).
215    o pkg/dic: add modifications from Hajoon:
216      - add iron sedimentation flux function of the flux of PO4 out of bottom layer
217      - light attenuation: add self-shading effect by climatological Chlorophyll
218        (read from file).
219    o pkg/seaice/seaice_preconditioner.F : fix cpp flags, so that code also
220        compiles when SEAICE_ALLOW_DYNAMICS is undefined
221    o global_oce_llc90,cs32 : ecco_v4 verification exps now use separate
222        AUTODIFF_OPTIONS.h, CTRL_OPTIONS.h, ECCO_OPTIONS.h rather than
223        ECCO_CPPOPTIONS.h (i.e. replacing the deprecated, circumvoluted,
224        approach where CPP_OPTIONS.h includes an omnipotent
225        ECCO_CPPOPTIONS.h that superseeds pkgs '_OPTIONS.h').
226    o pkg/ctrl/ctrl_map_ini_genarr.F : include GMREDI_OPTIONS.h where
227        ALLOW_3D_KAPGM, ALLOW_3D_KAPREDI (to be added later) belong.
228    o pkg/seaice/SEAICE_SIZE.h : replace ALLOW_AUTODIFF_TAMC with ALLOW_AUTODIFF
229        to avoid needing AUTODIFF_OPTIONS.h anytime SEAICE_SIZE.h is included
230        (it seems that THSICE_SIZE.h, PTRACERS_SIZE.h have the same issue...)
231    o model/src, pkg/smooth, profiles, salt_plume, seaice, exf, ecco, ctrl :
232      add ECCO_OPTIONS.h, AUTODIFF_OPTIONS.h, CTRL_OPTIONS.h when needed
233      - ECCO_OPTIONS.h is needed when including ecco_cost.h, ecco.h
234      - AUTODIFF_OPTIONS.h is needed when including tamc.h, tamc_keys.h
235      - CTRL_OPTIONS.h is needed when including ctrl.h, optim.h, etc
236    o pkg/ctrl,ecco,exf :
237      - remove variables associated with NON-generic costs,
238        unless ECCO_CTRL_DEPRECATED is defined.
239      - remove NON-generic NON-obcs controls,
240        unless ECCO_CTRL_DEPRECATED is defined.
241    o pkg/ctrl : added ctrl_local_params.h
242    o pkg/ecco : added ecco_local_params.h
243    o verification : added ECCO_CTRL_DEPRECATED to most as and oad
244        experiments, which will progressively be removed once the
245        experiments are revised to use only generic controls/costs.
246    
247    checkpoint65f (2014/10/13)
248    o pkg/ctrl, ecco, seaice : minor changes in headers
249    o global_oce_llc90,cs32 : ecco_v4 verification exps now test generic
250        cost function and control capabilities added since checkpoint65c.
251    o pkg/ctrl: regroup all obcs ctrl variables in CTRL_OBCS.h.
252    o pkg/ecco: increase modularity and improve logic
253      - ecco_toolbox.F (new) : basic operations that are often used,
254        including ecco_zero, ecco_cp, ecco_cprsrl, ecco_diffmsk,
255        ecco_addcost, ecco_add, ecco_div, ecco_readbar.
256      - cost_gencal.F (new) : determine filnames and records.
257      - cost_genread.F, cost_generic.F : condense code using ecco_toolbox.
258      - cost_generic.F : use spzeroloc rather than hard-coded -9999,
259        and exploit the above functions leading to more compact code.
260      - cost_genread.F : rm un-necessary arguments, re-arrange args order.
261      - ecco_check.F, ecco_readparms.F : bring logic up to date (using_tpj,
262        using_ers, using_gfo, using_mdt is now only defined within gencost)
263    o pkg/ecco: extend generic cost capabilities
264      - pkg/ecco : added 3D gencost (ALLOW_GENCOST3D), added gencost
265        capabilities via gencost_preproc (+ _c,_i,_r related params),
266        gencost_posproc (+ _c,_i,_r related params), gencost_nrecperiod,
267        gencost_outputlevel arguments to cost_generic
268      - pkg/ecco/cost_genread.F (new) : interface between cost_generic
269        and active_read, making room for gencost_preproc to take effect.
270        Added preproc code ('climmon') to replace GENERIC_BAR_MONTH.
271      - pkg/ecco/cost_generic.F :
272        compute and store misfit maks*(model-data) in localdif that is
273        then used in cost and for output, output misfit if outlev.GT.0
274      - pkg/ecco/cost_averagesfields.F : depend on gencost_barskip
275        rather gencost_barfile for deciding whether to write barfile
276      - pkg/ecco/ecco_cost_driver.F : extract call to cost_profiles from cost_hyd
277      - pkg/ecco/ecco_cost_init_fixed.F : extract gencost weights from
278        ecco_cost_weights, re-activate ecco_summary, determine gencost_barskip
279      - pkg/ecco/ecco_summary.F : report gencost variables
280      - pkg/ecco/cost_gencost_customize.F : treat case of scatterometer costs
281    o pkg/sea/SEAICE_OPTIONS.h : add ALLOW_COST_ICE, ALLOW_SEAICE_COST_SMR_AREA
282    o pkg/ecco: internal re-organization
283      - resolve un-necessary dependancies within pkg/ecco and to pkg/ctrl
284        (rm includes of optim.h and ctrl.h, replace optimcycle with eccoiter,
285        mv relevant variables definition from ecco_cost.h to ecco.h,
286        make sure gencost routines only include ecco.h)
287      - throughout : add missing CCP switches, rm unused variables,
288        and rm un-necessary CPP switches
289      - throughout gencost routines : replace modelstartdate with
290        gencost_startdate when needed
291      - cost_bp_read.F : use arguments, instead of common blocs
292        so that it can be used also by (new) cost_gencost_bpv4.F
293      - ecco_init_varia.F : initialize RHOsumGlob_0, VOLsumGlob_0
294      - ecco_phys.F : compute etanFull
295      - ecco_readparms.F : rename using_topex as using_tpj, add
296        bpv4-grace to gencost list, initialize added parameters
297      - ecco_summary.F : include SIZE.h needed for augmented ecco.h
298    
299    checkpoint65e (2014/09/28)
300    o pkg/obcs:
301      - replaced the just-added CPP options with run-time variables
302        (see pkg/obcs/OBCS_PARAMS.h for details)
303    o pkg/obcs:
304      - add additional CPP options for sponge boundary layers
305        (see pkg/obcs/OBCS_OPTIONS.h for details)
306    o pkg/atm_phys:
307      - implement simple (obliquity only, assuming circular planet orbit)
308        seasonal cycle for incoming insolation ; used when run-time parameter
309        "select_incSW" is set to 1 ; by default (=0), uses old unchanged formula.
310    o pkg/exf, pkg/ecco, pkg/ctrl:
311      - add ctrlUseGen and .NOT.ctrlUseGen exclusive brackets
312      - activate gentim2d for xx_atemp,xx_aqh,xx_swdown,xx_lwdown,xx_precip
313      - implement gentim2d for xx_tauu, xx_tauv
314    o pkg/ctrl, pkg/ecco :
315      - fix and complement the generic controls implementation
316    o pkg/exf/exf_getclim.F, exf_getffields.F, exf_getsurfacefluxes.F,
317        pkg/shelfice/shelfice_thermodynamics.F,
318        verification/hs94.1x64x5/code_ad/ctrl_map_gentim2d.F,
319        verification/hs94.1x64x5/code_oad/ctrl_map_gentim2d.F :
320      - add weight in CTRL_GET_GEN calls
321    o pkg/ctrl :
322      - ctrl_get_gen.F : add weight as input parameter, and apply
323        it directly rather than using smooth_correl2Dw (to be retired).
324      - ctrl_map_gentim2d.F : add wgentim2d in CTRL_GET_GEN call
325    o pkg/ctrl :
326      - introduce run time switch ctrlUseGen (in ctrl.h, ctrl_check.F,
327        ctrl_init_variables.F, ctrl_readparms.F, ctrl_summary.F)
328      - re-activate ctrl_summary (call in ctrl_init.F)
329    o pkg/ctrl/ctrl_cost_gen.F :
330      - fix CCP brackets, remove un-necessary myiter, mytime
331    o pkg/ecco, pkg/ctrl :
332      - remove cost_genctrl.F (now pkg/ctrl/ctrl_cost_driver.F)
333      - remove cost_forcing_gen.F (now pkg/ctrl/ctrl_cost_gen2d.F)
334      - add ctrl_cost_driver.F (replaces pkg/ecco/cost_genctrl.F)
335      - add ctrl_cost_gen2d.F (replaces pkg/ecco/cost_forcing_gen.F)
336    
337    checkpoint65d (2014/09/15)
338    o pkg/ctrl/CTRL_GENARR.h, ctrl_readparms.F :
339      - add xx_genarr2d_bounds,xx_genarr3d_bounds,
340        xx_genarr2d_numsmooth, xx_genarr3d_numsmooth
341    o pkg/ctrl/ctrl_map_ini_genarr.F :
342      - remove compilation failure
343      - replace hard coded example
344      - use ctrl_map_ini_gen2D, 3D
345    o pkg/ctrl/ctrl_pack.F, ctrl_unpack.F :
346      - on/off switch based upon xx_genarr2d_weight,3d
347    o src/packages_init_variables.F :
348      - regroup pkg/estimation calls
349      - remove CTRL_MAP_INI, CTRL_MAP_INI_ECCO
350        (moved to ctrl_init_variables)
351    o cost/cost_init_varia.F :
352      - remove ecco_cost_init_varia
353        (moved to ecco_init_varia)
354    o ctrl/ctrl_init_variables.F :
355      - add CTRL_MAP_INI, CTRL_MAP_INI_ECCO
356      - add CTRL_MAP_INI_GENARR, SEAICE_CTRL_MAP_INI
357    o ctrl/ctrl_map_ini.F :
358      - remove ctrl_map_ini_genarr, seaice_ctrl_map_ini
359        (moved to ctrl_init_variables)
360    o ctrl/ctrl_map_ini_ecco.F :
361      - remove ECCO_PHYS (mv to ecco_init_varia)
362      - remove seaice_ctrl_map_ini (mv to ctrl_init_variables)
363    o ecco/ecco_ad_diff.list :
364      - add ecco_init_varia.f
365    o ecco/ecco_init_varia.F
366      - new routine
367    o pkg/streamice:
368      - move call to STREAMICE_INITIALIZE_PETSC from the_model_main.F
369        to streamice_init_fixed.F ; fix setting of STREAMICEisOn.
370      - from model/src, only call pkg/streamice routines if useStreamIce is TRUE.
371    o pkg/generic_advdiff:
372      - fix bug causing negatives in GAD_SMOLARKIEWICZ_HACK, introduced 8/8/2014.
373    o pkg/exf:
374      - exf_getffields.F : add xxg_atemp, xxg_aqh, xxg_swdown, xxg_lwdown,
375         xxg_precip capability
376      - exf_wind.F : no need for ALLOW_WSPEED_CONTROL cpp switch
377    o pkg/ctrl:
378      - ctrl_pack.F, ctrl_unpack.F : omit (un)packing of undefined xx_gentim2d
379    o model/src/load_fields_driver.F :
380      - move call to CTRL_MAP_GENTIM2D ahead of BULKF_FIELDS_LOAD, EXF_GETFORCING
381        so xx_gentim2d can be applied to e.g. exf fields in a timely manner.
382    o OpenAD verification experiments:
383      - Include explicitly AUTODIFF_OPTIONS.h, COST_OPTIONS.h, and CTRL_OPTIONS.h
384        in src files (to enable to skip the ECCO_CPPOPTIONS.h);
385        for now, only in pkgs used in OpenAD verification experiments.
386      - Change all except one (halfpipe_streamice) OpenAD experiments to use pkg
387        specific option file (AUTODIFF_OPTIONS.h, COST_OPTIONS.h & CTRL_OPTIONS.h)
388        instead of using ECCO_CPPOPTIONS.h (with all relevant options in it).
389        This applies to both TAF (*/code_ad/) and OpenAD (*/code_oad/) tests.
390    o OpenAD verification tests:
391      - finish removing ALLOW_AUTODIFF_OPENAD (replaced by ALLOW_OPENAD);
392      - finish replacing ALLOW_AUTODIFF_TAMC by ALLOW_AUTODIFF (except for
393        tape/storage which are specific to TAF/TAMC) in pkgs used in OpenAD exp.
394      - leave ALLOW_AUTODIFF_TAMC undef in OpenAD verification tests (*/code_oad/)
395        and remove local copy of "tamc.h" (needed only if ALLOW_AUTODIFF_TAMC).
396    o model/src & generic_advdiff:
397      - change hard-coded loop-range in S/R APPLY_FORCING_T/S to set T & S
398        forcing over 0:sNx+1, 0:sNy=1 (instead of just 1:sNx,1:sNy): Model
399        dynamics requires valid tracer forcing only over tile interior but, for
400        some diagnostics, may need to extend it over 1 point in tile halo region.
401      - in gad_calc_rhs.F, update calls to gmredi_x/ytransport to get valid
402        tendency over i,j_Min,Max range. This allows to simplfy i,j_Min,Max
403        setting in temp/salt/ptracer_integrate.F
404    o pkg/dic:
405      - multiply tracer tendencies by maskInC (if ussing OBCS) before updating
406        gchem tracers: this allows to skip the 2nd call to OBCS_APPLY_PTRACER
407        in gchem_forcing_sep.F
408    o pkg/autodiff & obcs:
409      - add storage of current OB values for U & V (was already done for T & S).
410    
411    checkpoint65c (2014/08/30)
412    o pkg/dic:
413      - take out of dic_atmos.F the initialisation part (corresponding to istate=0)
414        and move it to new S/R DIC_INI_ATMOS ; also simplify the dependencies for
415        AD compilers (but still a Pb with multi-threading and Ajoint dependencies).
416    o verification:
417      - new experiment testing OBCS and DIC pkgs, using params & output from
418        tutorial_global_oce_biogeo but limited to a sub-domain around Drake
419        passage: take OBCS from the last year of a 2 yrs simulation and initial
420        conditions from t=1.yr while keeping the same forcing.
421    o model/src:
422      - when allowFreezing=T, only reset below freezing initial Temp to Tfreezing
423        (at any k) if checkIniTemp is not set to FALSE (default checkIniTemp=T).
424    o pkg/diagnostics (Stats-Diags):
425      - re-write S/R DIAGSTATS_CALC (more similar to TARGET_NEC_SX part) to:
426        a) refine region where stats are computed, excluding where arrhFac=0 ;
427        b) fix missing interior mask (use with OBCS) in 3-D fields statistics;
428        c) extend regional stats over the region edges (for U or V pt location).
429    o pkg/gchem:
430      - add a (2nd) call to OBCS_APPLY_PTRACER in gchem_forcing_sep.F (needed since
431        pkg/dic is currently changing/updating ptracers values beyond the OB).
432    o verification (hs94):
433      - convert Held & Suarez external_forcing.F (from hs94.cs, with Sigma-coords
434        bits) to apply_forcing.F and undef USE_OLD_EXTERNAL_FORCING in experiments
435        hs94.cs-32x32x5 and tutorial_held_suarez_cs ; also use the same version of
436        apply_forcing.F in 2D hs94 (fwd, ad & oad) but keep LatLon hs94 unchanged.
437    o pkg/ptracers:
438      - implement Adams-Bashforth on pTracers (instead of on tracer tendency),
439        switched on by setting PTRACERS_doAB_onGpTr=F (default set to doAB_onGtGs)
440      - change exp4 secondary test (input.nlfs) to test AB on T & S & pTracer
441        with AB-2 code (so that it get tested).
442    o model/src:
443      - allow to apply Adams-Bashforth on Temp & Salt (instead of on gT,gS)
444        with AB-2 code (previously only available with AB-3 code).
445        Also move the code (call to ADAMS_BASHFORTH S/R) form thermodynamics.F
446        to temp/salt_integrate.F (since the 2nd part of gt/sNm1 update is done
447        there, by calling CYCLE_AB_TRACER).
448    o pkg/generic_advdiff:
449      - change Tracer argument (drop bi,bj indices) in S/R GAD_CALC_RHS,
450        GAD_C2_ADV_R, GAD_U3_ADV_R, GAD_C4_ADV_R, GAD_DIFF_R, GAD_BIHARM_R
451        + also in GMREDI_X/Y/RTRANSPORT ; and update corresponding calls in
452        S/R temp/salt/ptracers_integrate.F
453    o model/src & pkg/ptracers:
454      - remove common block variables gT,gS (in DYNVARS.h), except when
455         USE_OLD_EXTERNAL_FORCING is defined; update TAF storage accordingly.
456      - remove common block variable gPtr.
457    o model/src & pkgs:
458      - change gTracer (and/or tracer) argument (drop bi,bj indices) in S/R
459         ADAMS_BASHFORTH(2&3), CYCLE_(AB_)TRACER, FREESURF_RESCALE_G,
460         IMPLDIFF, SOLVE_TRIDIAGONAL & _PENTADIAGONAL, GAD_SOM_ADVECT,
461         GAD_ADVECTION, GAD_CALC_RHS, GAD_IMPLICIT_R, TIMEAVE_CUMUL_DIF_1T
462      + change accordingly all calling S/R.
463      - use local temperature/salinity/ptracer tendency array (without bi,bj)
464        in S/R temp/salt/ptracers_integrate.F
465    o pkg/seaice:
466      - move computation of stress contributions to rhsU/V in s/r seaice_lsr
467        to separate s/r seaice_lsr_rhsu/v in order to be able to re-use code
468      - make u/vIceC local variables if SEAICE_CGRID is defined only to be used
469        in seaice_lsr
470      - seaice_lsr.F: make local 4D arrays UXX, VYY, etc 2D
471        in s/r seaice_lsr_calc_coeffs
472    o pkg/mdsio:
473      - in MDS_READ/WRITE_FIELD, stop if file-name (+prefix) is too long
474        (e.g., > MAX_LEN_MBUF -90 ) or, in INI_MODEL_IO, if mdsioLocalDir is
475        too long (e.g., > MAX_LEN_FNAM/2 )
476    
477    checkpoint65b (2014/08/12)
478    o pkg/matrix:
479      - pass, as argument, updated tracer field to S/R MATRIX_STORE_TENDENCY_EXP
480    o model/src:
481      - fix initialisation of geothermalFile; compile geothermal-flux code in
482        natl_box verification experiment.
483    o model/src & pkg/down_slope:
484      - push the k loop inside S/R TIMESTEP_TRACER (now called outside the k loop).
485      - in dwnslp_apply.F, update tracer tendency instead of the future tracer
486        field and call S/R DWNSLP_APPLY before TIMESTEP_TRACER; this affects
487        results at machine truncation level: update output of experiments
488        global_ocean.90x40x15.dwnslp (FWD) and lab_sea.noseaicedyn (AD).
489    o pkg/diagnostics (Stats-Diags):
490      - do not cumulate the full volume when DIAGNOSTICS_FILL is
491        called with bibjFlg < 0 (no increment of the counter for 2D/3D diag);
492        This fix the mean statistics when DIAGNOSTICS_FILL is called multiple
493        times (but Min,Max and StD are still wrong).
494    o model/src:
495      - add geothermal flux forcing: read in ini_forcing.F as time-constant
496        field "geothermalFlux" applied in S/R APPLY_FORCING_T (apply_forcing.F)
497    o model/src:
498      - move CYCLE_TRACER calls from tracers_correction_step.F to temp/salt/ptracer
499        _integrate.F so that, when leaving S/R THERMODYANMICS, theta,salt and pTrac
500        arrays are already updated while adjustments (filters, conv.adjustment) are
501        still applied later, in S/R TRACERS_CORRECTION_STEP.
502  o model/src:  o model/src:
503    - store in common bloc array "rhoInSitu" the virtual potential temperature    - store in common bloc array "rhoInSitu" the virtual potential temperature
504      anomaly that is used to compute geopotential: this make the atmos code      anomaly that is used to compute geopotential: this make the atmos code
# Line 82  o model/src: Line 577  o model/src:
577    - collect T & S forcing tendency into local 2-D array (passed as argument to    - collect T & S forcing tendency into local 2-D array (passed as argument to
578      apply_forcing S/R); add 3-D diagnostics for temp & salt forcing tendency.      apply_forcing S/R); add 3-D diagnostics for temp & salt forcing tendency.
579      affects results at machine truncation level (e.g., ocean: if using      affects results at machine truncation level (e.g., ocean: if using
580      SHORTWAVE_HEATING and atmos:, if using FrictionHeating or (atm_Rq<>0)      SHORTWAVE_HEATING ; atmos: if using FrictionHeating or (atm_Rq<>0) moist
581      moist correction due to vertical advection).      correction due to vertical advection).
582    - update output of experiments: fizhi-cs-32x32x40 (13 digits),    - update output of experiments: fizhi-cs-32x32x40 (13 digits),
583      global_ocean.cs32x15 (fwd: icedyn, seaice, thsice, all 3 @ 11d; adm & tml:      global_ocean.cs32x15 (fwd: icedyn, seaice, thsice, all 3 @ 11d; adm & tml:
584      seaice_dynmix & thsice, fwd grd @ 3 & 8), hs94.1x64x5 (ad grd @ 13),      seaice_dynmix & thsice, fwd grd @ 3 & 8), hs94.1x64x5 (ad grd @ 13),
# Line 112  o pkg/seaice: more clean up Line 607  o pkg/seaice: more clean up
607      nearly a factor of 10 (but keep in mind there is mostly seaice code in      nearly a factor of 10 (but keep in mind there is mostly seaice code in
608      this experiment):      this experiment):
609      > size mitgcmuv_ad      > size mitgcmuv_ad
610      text    data        bss     dec        hex        text    data        bss     dec        hex
611      5100659 1216  220458952  225560827  d71c8fb      5100659 1216  220458952  225560827  d71c8fb
612      5339425 1200 1868783272 1874123897 6fb4dc79 (SEAICE_LSR_ADJOINT_ITER)      5339425 1200 1868783272 1874123897 6fb4dc79 (SEAICE_LSR_ADJOINT_ITER)
613      but only a moderate increase of 36% for global_ocean.cs32x15:      but only a moderate increase of 36% for global_ocean.cs32x15:
614      text    data        bss     dec        hex      text    data        bss     dec        hex
Legend:
Removed from v.1.2040  
changed lines
  Added in v.1.2120
  ViewVC Help
Powered by ViewVC 1.1.22