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

Legend:
Removed from v.1.1338  
changed lines
  Added in v.1.1418

  ViewVC Help
Powered by ViewVC 1.1.22