/[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.1353 by jmc, Sun Dec 19 23:20:01 2010 UTC revision 1.1431 by jmc, Fri May 13 14:18:22 2011 UTC
# Line 4  $Name$ Line 4  $Name$
4      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
5      ==============================      ==============================
6    
7    checkpoint62x (2011/05/13)
8    o pkg/mom_vecinv:
9       high-order vorticity scheme: use interior masks (maskInW & S)
10       to zero-out (outside OB) gradient of vorticity.
11    o verification exp.:
12       turn on "useEnergyConservingCoriolis" (better discretisation) in exp.
13       hs94.128x64x5 and seaice_obcs ; update results.
14    o pkg/obcs:
15      fix bug in OBeta for YZ sections in obcs_exf_load.F
16    o pkg/ctrl, pkg/shelfice, and optim:
17      - add new control variable xx_shifwflx (fresh water flux underneath ice
18        shelves). This is almost as tedious as obcs-ctrl, because the
19        variables needs its own mask.
20      - add a cost function term that penalizes the deviation of xx_shifwflx
21        from the first guess (zero) in analogy to other control variables
22    o pkg/shelfice:
23      - start basic adjoint
24    o pkg/ctrl:
25      - remove extra 2D field "for good measure" when initializing in
26        ctrl_set_globfld_[xy/xz/yz].F
27      - replace hard coded precision 32 (a bug) in ctrl_map_ini_gen by ctrlprec
28      - improve error messages in ctrl_set_unpack_*
29    o OBC implementation in momentum:
30       use interior masks (maskInC,W & S) in momentum S/R:
31       - multiply divergence, tension and some del^i.u or del^j.v by maskInC
32       - multiply del2u,del2v by maskInW/S (biharmonic viscosity)
33       - multiply horiz. grad of vorticity by maskInW/S (in Leith scheme)
34       - multiply gradient of wVel by maskInW/S in del2w calculation (NonHyd)
35       Motivation: fix/clean implementation of OBC with biharmonic-viscosity
36        (and also with Leith/Smagorinsky harmonic-viscosity).
37       Update output of experiment dome and exp4 (+nlfs).
38    o exp4 test-exp:
39      - turn on: nonHydrostatic, exactConserv, viscA4 (with momDissip_In_AB=F)
40        and prescribe time-varying OB (without pkg/exf) for zonal-velocity
41        at Western & Eastern Boundaries (+/- 1 per-cent); Also prescribe
42        time-varying OB at Western-side for wVel (NonHyd); update output.
43      - add secondary set-up (.nlfs) using NonLin FreeSurf (z*) + mom_vecinv
44        and biharmonic-viscosity ; also use prescribed time varying OB.
45    o pkg/obcs:
46      fix bug in OBWeta,OBWeta indices in obcs_apply_r_star.F
47    o model/src:
48      - fix startTime setting when baseTime<>0 & nIter0=0 are specified ;
49      - initialise startTime & nIter0 with unlikely values: allows to set one
50        from the other without ambiguity.
51    o seaice_obcs test-exp.:
52      add small diffK4T,S to test biharmonic diffusivity with OBCS ; update output.
53    o pkg/obcs:
54      - remove calls to OBCS_COPY_TRACER after EXCH for T,S & pTracer:
55        no longer needed with maskInC,W,S in pkg/generic_advdiff.
56    o pkg/seaice:
57      - retired SEAICE_SALINITY and SEAICE_salinity in favor of
58        SEAICE_VARIABLE_SALINITY and SIsalFRAC, respectively.
59      - fixed bug in the initialization of sea ice salinity variable HSALT
60        in seaice_init_varia.F  (removed dependence on ICE2WATR)
61    o fizhi-cs test exp.:
62       switch reference output convention for these 2 non-standard test-exp,
63        from: pgi  6.1-5 on faulks (FC.6, -of linux_ia32_pgf77+authors_fc5)
64          to: ifort 11.1 on baudelaire (FC.13, -of linux_amd64_ifort11)
65       reason: faulks might get retired soon and pkg/fizhi has still Pb with gfortran
66    o model/inc:
67      - change type of PmEpR array (from _RS to _RL) to fix restart (with "-use_r4")
68        of exp. global_ocean.90x40x15 (Suncronous time-stepping + Real-FW + NonLin-FS)
69    o pkg/seaice:
70      - fix calculation of saltFlux and saltPlumeFlux when using SEAICE_SALINITY
71    o verification/lab_sea, 1D_ocean_ice_column, seaice_obcs, global_ocean.cs32x15
72      - updated output*.txt files due to bug fixes in salt, saltPlumeFlux,
73        and initial HSALT
74    
75    checkpoint62w (2011/04/26)
76    o pkg/mom_fluxform:
77      - fix calculation of del2u & del2v (no-slip BC) when using Non-Lin Free-Surf.
78      - use parameter "sideDragFactor" in no-slip BC side-drag part.
79      - update output of exp. global_ocean.90x40x15
80    o pkg/mom_vecinv:
81      mask del2u & del2v (no effect on results, but cleaner & closer to mom_fluxform)
82    o update seaice exp.:
83      use more consistent parameter values (remove "backward compatability only":
84       ICE2WATR, SEAICE_lhSublim, SEAICE_cpAir, SEAICE_rhoAir + set rhoConstFresh)
85      related to thermodynamics seaice. Update output (fwd+AD) of exp.:
86       1D_ocean_ice_column (x2), lab_sea (x8), ioffline_exf_seaice (x2) & seaice_obcs
87    o update test-exp global_ocean.90x40x15:
88      o use default CD_CODE_OPTIONS file (Adams-Bashforth on CD Mom & Coriolis).
89      o increase Overlap size from 2 to 3 (to allow to test bi-harmonic Visc.)
90      o new pickup-files from a 100.yr run with pkg/down_slope (input.dwnslp params)
91      o turn on zStar & viscA4 in standard set-up (dir input/)
92        and update output files: output.txt & output.dwnslp.txt
93    o new AD-test experiment: obcs_ctrl
94      (From Matt MITgcm_Contrib/SOSE/BoxAdj) to test OBCS as control vector.
95    o pkg/dic:
96      fix two bugs in handling water vapour pressure in carbonate chemistry
97      and air-sea fluxes. Found by Val Bennington and Galen McKinley
98      (add DIC Option = WATERVAP_BUG to recover old results)
99      update output of tutorial_global_oce_biogeo (fwd+AD) & tutorial_dic_adjoffline
100    o adding allowInteriorFreezing option to check for water below freezing point
101      at depth and bring the negative heat anomaly to the surface level
102    o pkgs _fields_load:
103      - extend GET_PERIODIC_INTERVAL to deal with non-periodic case (from pkg/rbcs)
104      - use new S/R GET_PERIODIC_INTERVAL in pkgs:
105        bulk_force, cheapaml, cfc, dic, obcs (+ ctrl/ctrl_get_gen_rec), offline and
106        rbcs to get interpolation weights and time reccord number.
107      - fix multi-threading in aim_fields_load.F
108      - apply EXCH to temp arrays [0,1] after reading in dic_fields_load.F
109    o external_fields_load:
110      new S/R GET_PERIODIC_INTERVAL to calculate which time-record (and linear
111       interpolation weights) to use in a periodic time-series ; do not divide by
112       time-step to convert to integer:
113       - this fixes case where externForcingPeriod is not multiple of time-step
114          (Thanks to Christopher L. Wolfe)
115       - this fixes case where externForcingPeriod is an odd multiple of time-step
116         (was wrong by half of time-step); update results of tutorial_global_oce_in_p.
117      store index of time-record which is currently loaded:
118       - works without any assumption on sequence of calls to EXTERNAL_FIELDS_LOAD ;
119       - needed for the adjoint but not yet used because of extensive recomputations;
120         (to be fixed soon)
121    
122    checkpoint62v (2011/04/02)
123    o verification/1D_ocean_ice_column added to test Ian Fenty's
124      pkg/seaice_growth.F modifications for improving sea ice code differentiability
125    o update output of exp4, internal_wave, seaice_obcs, tutorial_plume_on_slope
126      (new velocity statistics after changing maskInW & maskInS definition)
127    o pkg/generic_advdiff with OBCS:
128      use maksInW & maskInS in advective flux and in biharmonic diffusion calculation:
129      equivalent to assumption of no tracer gradient outside the OB region.
130      This will affect the results if using obcs+biharmonic diffusion.
131      update output_adm.txt of experiment bottom_ctrl_5x5 (not expected since no obc)
132    o pkg/obcs:
133      change definition of maskInW & maskInS: leave OB edges inside (maskIn=1)
134      (e.g., maskInW(OB_Ie)=1 ) so that velocity normal-component at OB is still
135      in Interior region. No impact on solution but affects Velocity monitor stats.
136    o pkg/ptracers (adjoint)
137      New option AUTODIFF_PTRACERS_SPLIT_FILES for writing separate
138      tape file for each tracer
139    o model/inc/src and pkg/autodiff: add parameters that allow to turn
140      of parts of the adjoint code at run time. For now, there are only
141      useKPPinAdMode, useGMrediInAdMode, and useSEAICEinAdMode. These
142      parameters are set in data.pkg.
143      This is the lazy coward's version of doing this, but still useful for
144      debugging. In order to do it right, everything needs to be moved to
145      pkg/autodiff.
146    o pkg/autodiff: fix global_sum_tile_ad.F by resetting adsumPhi at
147      the end; this fixes the adjoint balancing code for open boundaries.
148    o pkg/ctrl:
149      - change behavior of xx_${varname}period to be analogous to exf-conventions:
150        xx_${varname}period =   0.,  means time-independent control part
151        xx_${varname}period = -12., means one-year-periodic monthly control part
152      - introduce helper routine ctrl_init_rec (part of ctrl_init.F) to
153        clean up ctrl_init a little, remove unused variables
154      - ctrl_init: fixed the handling of the start/endrec for obcs-variables
155      - move "call cal_GetMonthsRec" from s/r ctrl_get_gen to s/r ctrl_get_gen_rec
156    
157    checkpoint62u (2011/03/15)
158    o model/src:
159      - for the adjoint, change NonLin-FreeSurf initialisation sequence: fill-in
160        h0Fac from INITIALISE_FIXED and reset hFac to h0Fac in INITIALISE_VARIA ;
161      - in initialise_varia.F (except for ALLOW_DEPTH_CONTROL) remove 2nd call to
162        INI_DEPTHS, INI_MASKS_ETC & INI_LINEAR_PHISURF.
163    o pkg/ctrl/obcs/ecco:
164      - remove obsolete, unused, and partially broken code:
165        ctrl_obcsvol/ctrl_obcsbal/ctrl_volflux and code within
166        ALLOW_CTRL_OBCS_BALANCE
167      - put a stop statement into ecco/cost_obcsvol.F, but also add
168        cost_obcsvol to ecco_ad_diff.list
169    o pkg/ecco: fix counting cost function contributions for obcs (code
170      by Matt M.) and diffkr for diagnostics
171    o pkg/ecco, pkg/ctrl, pkg/grdchk:
172      add cost penalty terms for initial uvel, vvel, eta controls
173      and I/O for prior error files
174    o pkg/mdsio:
175       move mds_check4file.F from eesupp/src to mdsio_check4file.F ; add option
176       (new argument) to check for tiled file in mdsioLocalDir : this fix the
177       restart from pickup_som , pickup_dic (& dic_co2atm) with mdsioLocalDir.
178    o pkg/ctrl:
179     - ctrl_obcsbal.F, ctrl_obcsvol.F: replace last calls of ctrl_getrec by
180       ctrl_get_gen_rec, put a stop statement into ctrl_getrec. Remove
181       unused and superfluous variables
182     - ctrl_get_gen_rec.F: catch the case of xx_genperiod=0.
183    o pkg/ecco
184     - enable adjoint of Stevens OB
185    o pkg/ctrl/ecco/grdchk:
186      - remove global declaration of tmpfldx/yz and tmpfldx/yz2.
187        Fortunately, we do not have any tests for the numerous cpp-flagged
188        option of "obcs as control parameters" so we will never know
189        if this is an improvement (but at least now things compile).
190      - remove an obsolete initialization from ecco_cost_init_barfiles.F
191      - fix some problems within never used ifdefs to make the code compile
192    o pkg/obcs:
193     - Refine Interior mask definition to work with more complex OB (e.g.
194       if OB_Ie < OB_Iw); In case default setting is not sufficient, allows to
195       read in Interior mask from file. Check if Interior mask is consistent
196       with OB list of indices.
197    o pkg/seaice
198     - beginning code for a seaice tracer package,
199       but currently only used for seaice age: SEAICE_AGE
200       SEAICE_num = 4
201       iTracer = 1: area-weighted ice age
202               = 2: volume-weighted ice age
203               = 3, 4: template for one-year vs. multi-year ice
204       data.seaice has new namelist SEAICE_PARM03 for arrays related to
205       SEAICE_num
206    
207    checkpoint62t (2011/03/05)
208    o pkg/obcs:
209     - obcs_apply_uv.F: change to a 2 steps method, first tangential component;
210       then normal component. This ensures that the normal component is set
211       correctly even when it conficts with tangential setting from an other OB.
212       Affects results of test-experiment seaice_obcs: update the output.txt
213    o pkg/obcs:
214     - add code to monitor statistics of dynamical fields at OB.
215     - take OBCS_BALANCE code out of obcs_calc.F into new S/R OBCS_BALANCE_FLOW :
216       - fix to OBCS_BALANCE code (right hFac ; deal with simple OB intersection).
217       - change default: balance the Net inflow through all OB (instead of balancing
218         each OB separately); keep the old method if setting OBC_balanceFacX < 0.
219     - switch on useOBCSbalance in dome experiment: only correct Western OB outflow
220       to match Nortern OB inflow (OBCS_balanceFacN=0).
221    o pkg/seaice:
222     - redefined actual sea ice age to be ICEAGE/HEFF ifdef SEAICE_AGE_VOL
223    o pkg/seaice (seaice_growth v 1.113):
224     - reverting seaice_growth to v 1.111 in order to have code with exact bit
225       compatibility with MITgcm_contrib/ifenty/Fenty_seaice_thermo_code_merge
226       CPP option SEAICE_OCN_MELT_ACT_ON_AREA removed in v 1.112 is reintroduced
227       runtime flag SEAICEareaFormula is kept but is inactive
228    o model/src:
229     - add new param to scale curvilinear grid distance and area when using
230       a different sphere radius from the grid-files.
231     - test/check multi-threads EXCH applied to local variable (cg2d_x)
232    o adjoint:
233     - enable adjoint of Stevens OB
234     - add dummy_in_dynamics to enable extraction of ADJrhoinsitu
235       (a bit delicate because of potentially modified adj. common block structure;
236       current switch is ALLOW_AUTODIFF_MONITOR_DIAG, could be changed later)
237    
238    checkpoint62s (2011/02/15)
239    o pkg/seaice (seaice_growth v 1.112):
240     - introduce SEAICEareaFormula run time parameter, which selects the recipe for
241       computing the ice cover tendency from effective thickness tendencies. Include
242       Ian Fenty's formula.
243     - remove 'undef SEAICE_OCN_MELT_ACT_ON_AREA' which may now be achieved using
244       SEAICEareaFormula.
245     - simplify notations for the various d_HEFF/d_AREA terms.
246    o pkg/gmredi:
247     - add Boundary-Value-Problem bolus transport (Ferrari etal, OM 2010) code
248       from Christopher Wolfe (UCSD); add partial cell factor in vertical derivative;
249       add a verification test-run (front_relax.bvp) to test GM_useBVP.
250    o pkg/seaice (seaice_growth v 1.111):
251     - for SEAICE_DO_OPEN_WATER_GROWTH, replace Gael's implementation with Ian's.
252     - the new SEAICE_DO_OPEN_WATER_MELT allows open-water air-sea heat fluxes to
253       also melt ice.
254    o pkg/seaice:
255     - McPhee formula for Ice-Ocean fluxes (provided by Ian Fenty).
256     - can be used by #define MCPHEE_OCEAN_ICE_HEAT_FLUX in SEAICE_OPTIONS.h
257    o pkg/cheapaml:
258     - from Nico Wienders: New version of this package.
259    o pkg/flt:
260     - add CPP options to prevent floats to re-enter the opposite side of
261       a periodic domain (default: keep periodicity)
262    o pkg/gmredi:
263     - apply exchange to GM_PsiX,GM_PsiY to get valid values near CS-corners
264       (needed for multi-dim advection).
265     - switch on GM Advective form in global_ocean.cs32x15 and update output.txt
266    o model/src:
267     - fix wrong re-setting of addMass array (ALLOW_ADDFLUID) for atmospheric run.
268    o tutorial_tracer_adjsens
269     - extend to NLFS and replace maskC by hFacC
270       (change reference output accordingly)
271    
272    checkpoint62r (2011/01/24)
273    o testreport:
274     - add option -deloutp (-do) to remove output file after successful run
275     - add option -mfile (-mf) for MPI machile-file and for each test-exp,
276       extract exact list of procs from this file ; adjust COMMAND by
277       replacing TR_MFILE & TR_NPROC with local machile-file & Nb of MPI procs.
278    o pkg/ecco: replace old ineffective way of penalizing deviations from
279      first guess (of obcs) with code by Matt Mazloff
280    o pkg/ctrl: add code that allows the use of ALLOW_OBCS?_CONTROL without
281      exf/cal
282      - ctrl_init: fix initialisation (computation of diffrec) for obcs-ctrl
283        and handle case of obcs?period == 0
284        Not affected by these changes: There is still an inconsistency in
285        the computation of start/endrec for the obcs variables for ALLOW_CAL
286        defined (see comment by Menemenlis, 7-Mar-2003)
287      - make S/R ctrl_get_gen_rec return something meaningful even if
288        ALLOW_CAL is not defined
289      - ctrl_getobcs?: simplify code so that it does not need to use
290                       S/R exf_swapfields
291    o testreport:
292       allow to test MPI with different number of processors for each experiment:
293       create file "SIZE.h.mpi" in each build dir.
294       Note: update of "-command" argument not yet implemented.
295    o model/src, pkg/autodiff, pkg/ecco, pkg/seaice:
296      alleviate the need for additional 3D tapes when using
297      the non-inear free surface in the adjoint. This is done
298      by adding an 'update' of hfacc etc. to their current
299      value at the beginning of forward_step.F.
300    o pkg/seaice
301      seaice_lsr: move stress-divergence computations out of outer pseudo-
302      timestepping loop
303    o pkg/obcs/obcs_prescribe_read.F: move 'CALL OBCS_EXTERNAL_FIELDS_LOAD'
304      up, so that the calls of the ctrl_* subroutines are really at the end
305      of the routine.
306    o pkg/gmredi:
307      allow to specify a 2-D horizontal map and a 1-D vertical profil
308      as scaling factor for Isopycnal diffusivity (Redi) and GM diffusivity.
309    o pkg/seaice:
310      move lsr iterations from seaice_dynsolver to seaice_lsr.
311    o pkg/seaice:
312      reduce memory footprint of seaice mutldim advection tapes,
313      by using a 2D tape independent of the generic_advdiff ones,
314      and setting GAD_AREA=1 etc. rather than GAD_AREA=101 etc.
315    
316    checkpoint62q (2010/12/29)
317    o model/src:
318      - check horizontal grid-spacing for unset or zero/negative values.
319      - change AD test hs94.1x64x5 to be truly a 2-D set-up.
320    o pkg/flt:
321      - use (recent) mdsio S/R MDS_READVEC_LOC instead of local FLT_MDSREADVECTOR
322      - precision of input and output files set according to readBinaryPrec
323        and writeBinaryPrec (previously always 64.b), pickup remains at 64 bits.
324      - initial float position from tiled file was not converted to index space
325        (looks like a bug): fixed to make it similar to initial pos. global file.
326    o global_ocean.cs32x15:
327       change SIZE.h of AD set-up to get same tiling as with MPI.
328       Update output_adm, output_adm.thsice (pass@14) & output_adm.seaice (fail@11)
329    o model/src:
330       add flag to control printing of domain fields (bathymetry, hFac ...) in
331       ini_mask_etc (call to plot_fieds_xy*).
332    o pkg/thsice:
333       if useEXF, some forcing fields are not valid in overlap: fixed by computing
334       over the interior only + apply EXCH at the end (like with SEAICE advection).
335    o pkg/mdsio:
336      - add 2 arguments to S/R MDS_PASS_R4/8toRL/S for overlap size of buffer
337        to copy to/from, so that it can be used to replace MDS_PASS_WH_R8toRL/RS.
338      - change barrier call for safe multi-threads access to 3-D shared buffer.
339      - MDS_READ/WRITEVEC_LOC use buffers from MDSIO_BUFF_3D.h (instead of local ones)
340      - replace call to MDS_PASS_WH_R4/8toRL with call to MDS_PASS_R4/8toRL and
341        remove the 4 S/R MDS_PASS_WH_R4/8toRL/S.
342    o pkg/flt:
343      - add fourth-order Runge-Kutta integration
344      - add basic exch2 support (only 1 facet supported so far!)
345    o pkg/mdsio+rw:
346      new S/R to read local (or tile independent) vector (read counterpart of
347       mds_writevec_loc, and front-end write_glvec_rl/rs).
348    o global_ocean.cs32x15: additional input_ad.seaice experiment
349        that uses pkg/seaice rather than pkg/thsice.
350        (using multidim advection for seaice)
351  o towards more F95 syntax for array specification in namelist:  o towards more F95 syntax for array specification in namelist:
352    - change "data.diagnostics" (data.ptracers was changed before checkpoint62p)    - change "data.diagnostics" (data.ptracers was changed before checkpoint62p)
353      of all verification experiments.      of all verification experiments.

Legend:
Removed from v.1.1353  
changed lines
  Added in v.1.1431

  ViewVC Help
Powered by ViewVC 1.1.22