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

Legend:
Removed from v.1.1364  
changed lines
  Added in v.1.1440

  ViewVC Help
Powered by ViewVC 1.1.22