/[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.1000 by gforget, Wed Jun 11 18:31:22 2008 UTC revision 1.1091 by jmc, Sun Apr 26 02:56:36 2009 UTC
# Line 4  $Name$ Line 4  $Name$
4      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
5      ==============================      ==============================
6    
7    o adjustment.cs-32x32x1: test "blank tiles" when compiling with mpi
8    o eesupp/src: new S/R to stop everyone when 1 (or more) Proc. find an error
9    o pkg/exch2: new S/R EXCH2_CHECK_DEPTHS to check that unconnected tile
10      edges are closed (zero depth); might happen with "blank tiles";
11      for now, disable if using OBCS.
12    o tools/gemake2: delete modifs from March 27-29 2009 (with option "-mpi":
13      link *_mpi files from "MODS" dir) to return to a strait-forward behavior.
14    o move read_write.F to write_utils.F ; new parameter (maxLengthPrt1D)
15      to reduce length of 1.D array print in STDOUT.
16    o adjustment.cs-32x32x1: changed to oceanic test with 1 continent
17      + use 48 tiles (16x8): this will allow to test "blank tiles".
18    o start to test MPI + multi-threaded (change testreport to skip tests
19      where number of tiles do not match).
20    o genmake2: add option "-omp" to turn on Open-MP ;
21      passed to testreport (if using -mth) ; merge few (*+mth) optfiles.
22    o pkg/exch2 (+ utils/exch2): ignore blank tiles (instead of removing them)
23      and keep the original tile numbering. Note: compatible with old
24      W2_EXCH2_TOPOLOGY.h & w2_e2setup.F topology files.
25    
26    checkpoint61l
27    o Overlaps had been forgotten in calculating ijk keys : fixed
28    o fix cyrus makedepend for gcc 4.4
29    o experiment fizhi-cs-aqualev20:
30      switching to new version of S/R GETPWHERE (#define TRY_NEW_GETPWHERE)
31      and generate new output. Note: new output is much closer to the old
32      one (v1.8) generated with pgf77 -O2 on Mar 14 2007 than it is from
33      the previous one (v1.9) generated with pgf77 -O0 on Jun 13 2007.
34    o pkg/fizhi:
35     - change to avoid indices going deliberately over array-bounds
36     - remove unnecessary MPI stuff
37     - perpetual spring equinox with "#define FIZHI_USE_FIXED_DAY"
38     - new version of S/R GETPWHERE available with #define TRY_NEW_GETPWHERE:
39       does not go over array-bounds and therefore returns reliable output.
40     - in fizhi_driver.F & fizhi_turb.F, fix "get_alarm" calls:
41       this S/R has 4 output arguments, which are used inside "get_alarm";
42       when it is called with the same variable "ndum" in place of several
43       output arguments, this produces weird side effect inside "get_alarm"
44      (e.g., datein=dates(n) is overwritten by timein=times(n) since both
45       datein & timein, resp. 2nd and 3rd arg., are pointed to the same "ndum"
46       address) and produces error.
47    o tools/gemake2: (March 27-29 2009)
48      with option "-mpi": link *_mpi files from "MODS" dir (this is what
49      testreport does).
50    o lab_sea experiment: change from 2x1 tiles to 2x2 tiles and regenerate output
51    o pkg/seaice: fix bug in new version of S/R SEAICE_LSR.F:
52       update lab_sea.lsr & lab_sea.salt_plume results
53    o pkg/seaice:
54      Major overhaul of the LSOR solver: change discretization from original
55      finite differences to finite volume. The old routine is still in the
56      file seaice_lsr.F and can be turned on with the CPP-flag
57      SEAICE_OLD_AND_BAD_DISCRETIZATION
58      Along with this change, the discretization of the strain rates is changed.
59      The same flag turns on the old discretization in seaice_calc_strainrates.F
60      All verification experiments with seaice are affected (except for
61      seaice_exf_offline); global_ocean.cs32x15/input.icedyn recovers the old
62      results by setting the above flag.
63      Further changes involve changes in defaults:
64      - ALLOW_FLOODING is defined and turned on by default
65      - SEAICEadvSnow = .true. is now the default
66      - SEAICEadvSalt = .true. is now the default
67      - SEAICEadvAge  = .true. is now the default
68      - SEAICE_clipVelocities  = .false. is now the default
69                                         (as per J. Zhang's recommendation)
70      - B-grid, and thus not tested: SEAICE_TEST_ICE_STRESS_1/EXPLICIT_SSH_SLOPE
71        is defined, SEAICE_TEST_ICE_STRESS_1 is renamed into SEAICE_BICE_STRESS
72      - seaice_growth: replace computation of UG by a simple copy from wspeed
73    
74    checkpoint61k
75    o pkg/seaice: a few bug fixes mostly concerning the metric terms:
76      seaice_lsr: fixed a few metric terms,
77                  affects lab_sea.lsr, lab_sea.salt_plume, seaice_obcs
78                  add more debugging diagnostics
79      seaice_evp: removed superfluous many metric terms of FV discretization
80                  of stress divergence
81                  affects lab_sea, lab_sea.hb87
82      seaice_ocean_stress: removed superfluous many metric terms of FV discretization
83                  of stress divergence for non-evp solver
84                  (this part is not tested anywhere)
85      seaice_calc_strainrates: fixed a bug in no-slip boundary conditions,
86                  they were effectively only half slip
87                  affects lab_sea.hb87
88      updated all relevant verification experiments (except lab_sea//output_adm.txt
89        which was left over, but fails at level 7 & 6)
90    
91    checkpoint61j
92    o pkg/seaice: turn a few hard wired parameters into run time parameters
93    o pkg/seaice/exf/monitor: turn netcdf coordinate variable T into model
94      time in seconds as promised by the usints in seaice_monitor and exf_monitor
95    o pkg/mnc: fix the coordinate variables for multi-tile curvi-linear grids
96      that are not cubed-sphere grids
97    o pkg/seaice: fix a bug for the rare case of SEAICE_VECTORIZE_LSR
98    o pkg/ctrl:
99      fix bug in ctrl_depth_ini.F (EXCH_XYZ was applied to 2-D field xx_r_low)
100    o non-hydrostatic (calc_gw.F):
101      fix viscA4 on CS-grid + fix diagnostics: Wm_Diss & Wm_Advec
102      update output of exp. global_ocean.cs32x15.viscA4
103    o remove backward compatibility for thetaMin,phiMin
104      (replaced by xgOrigin,ygOrigin on Jan 27, 2009)
105    o pkg/exf: fix a serious bug in exf_set_uv.F: do a rotation of the
106      interpolated wind field not only for curvilineargrid=.true., but also
107      for a rotated grid (rotateGrid=.true.)
108    o pkg/exf: add capability to bypass interpolation for runoff
109      by setting yet another CPP-flag
110    o eesupp/src, B-grid exchange:
111      add simple (6 tiles, non MPI) CS exch + standard interface S/R.
112    o autodiff
113      new CPP option ALLOW_TAMC_SINGLEPREC_COMLEV
114      allows to switch to single-prec common blocks;
115      corresponding switch occurs in tamc.h, integer isbyte
116      TAF store directivs "byte" and "kind" seem equivalent;
117      Most tamc.h have been updated correspondingly.
118    
119    checkpoint61i
120    o autodiff:
121      fix autodiff_sore/_restore indices for obcs
122      (spotted by I. Hoteit, M. Mazloff)
123    o pkg/autodiff, pkg/ctrl:
124      Add flag to reduce all tapelev I/O to single-prec.
125      In data.ctrl set
126        doSinglePrecTapelev = .TRUE.
127    o model/src/ini_curvilinear_grid.F, new grid-file format without EXCH2:
128      check in hack for SX8 compiler (for dxV, dyU); remove it afterwards.
129    o pkg/shelfice: fix bug with SHELFICEconserve flag, prepare for
130      real fresh water flux (not yet)
131    o pkg/mdsio & rw:
132     - add argument & option to S/R MDS_WRITEVEC_LOC (avoid Open/Close & write-meta)
133     - add simple interface S/R to MDS_WRITEVEC_LOC in pkg/rw: WRITE_GLVEC_RX.
134    o pkg/flt:
135     - replace FLT_INIT by 3-steps standard initialization S/R: readparms,
136       init_fixed & ini_varia.
137     - clean-up & simplify linear interpolation S/R:
138       move flt_bilinear.F -> flt_interp_linear.F ; add arg. myThid.
139     - fix some indices (mainly vertical index) ; add IMPLICIT NONE.
140     - more efficient I/O, using new options of MDS_WRITEVEC_LOC S/R.
141     - add vert. coordinate & horiz. indices to output files.
142     - change vertical position in initial condition, from index to depth [m].
143     - store horizontal (decimal) indices i,j in common block (instead of coord. x,y)
144       (easier to get to a curvilinear grid) ; delX,delY no longer used.
145    o pkg/ggl90: a few modifications:
146     - fix bottom boundary condition, and interpolation of KappaE
147       near top and bottom
148     - use same set of viscosity/diffusivity for TKE/momnetum/tracer
149     - add smoothing option as in Blanke and Delecluse (1993) (--> ALLOW_GGL90_SMOOTH)
150     - add a few options for mixing-length upper limit (--> mxlMaxFlag)  
151    o rename thetaMin,phiMin -> xgOrigin,ygOrigin
152      (maintain temporary backward compatibility, until next checkpoint)
153    
154    checkpoint61h
155    o pkg/diagnostics with "periodic averaging diagnostics":
156      fix S/R DIAGNOSTICS_COUNT ; hack diagnostics_out.F to get NetCDF output.
157    o pkg/my82: add missing "_d 0" ; update output of exp. vermix.my82
158    o pkg/obcs: fix loop ranges for balancing code
159    o pkg/flt:
160     - improve criteria for finding which tile owns a float.
161     - fix FLT_EXCH for case without MPI ; apply cyclic BC + fix also the MPI case.
162    o Add well-mixed CO2 box in aim_v23 (which includes changes in pkg
163      aim_compon_interf, ocn_compon_interf, gchem, and dic to get coupling
164      with ocean biogeo)
165    o pkg/flt:
166     - do some cleaning and fix the restart ;
167     - add output.txt in flt_example/results to get this pkg tested.
168    o genmake2: add an EXTENDED_SRC option and flag.
169       for now, cannot be turned on on command line (only through genmake_local)
170    o pkg/mdsio: read/write vector:
171     - singleCpuIO: fix RL/RS in MDSWRITEVECTOR by calling new S/R MDS_BUFFERtoRL/RS
172     - singleCpuIO=F: fix arrType & filePrec & _BYTESWAPIO in both MDSREADVECTOR
173       & MDSWRITEVECTOR by calling new S/R MDS_RD/WR_VEC_RL/RS.
174     - affects gradient of test-exp "bottom_ctrl_5x5" (not due to _BYTESWAPIO)
175      (drop to 4 digits matching for grad). But this test does not look very
176      reliable (was already failling on eddy, fc5, at 6.digits). Update output.
177    
178    checkpoint61g
179    o pkg/flt:
180      - move FLT_CPPOPTIONS.h to FLT_OPTIONS.h (standard name)
181       + make OPTIONS file more standard.
182      - rename S/R MDSREADVECTOR_FLT to FLT_MDSREADVECTOR (to mach src file name)
183      - fix _BYTESWAPIO CPP-Flag
184      - standard reading of file-parameter "data.flt"
185      - fix argument list of S/R FLT_RESTART
186      still need to fix _BYTESWAPIO in pkg/mdsio/mdsio_writevector.F
187    o pkg/obcs: fix a CPP-flag, bug reported by Jeff Blundell
188    o gather_2d.F/scatter_2d.F: suppress vectorisation on NEC SX in respond to
189      a still unfixed compiler/optimization bug
190    
191    checkpoint61f
192    o pkg/seaice, LSR solver: fix bug in GLOBAL_MAX calls ;
193      update results of global_ocean.cs32x15.icedyn test.
194    o move getcon.F from model/src to pkg/fizhi.
195    o pkg/exf: add parameter for extrapolation of temp up to the surface
196       and option to prevent negative evap (in bulk_largeyeager04 only).
197    o Prather Advection scheme: fix tracer conservation for divergent flow
198      (e.g., when using NonLin Free-Surface).
199    o kpp: retire KPPmld diagnostics --> use MXLDEPTH instead
200    o gmredi (Visbeck-variable-K):
201     - add more parameters to Visbeck-variable-K.
202     - use a 5 points average for dSigmaR in Visbeck-K calculation
203       (more "homogenous": same stencil and same weights as in dSigmaH)
204     - update output of ideal_2D_oce experiment.
205    o fix CPP-flag REAL4_IS_SLOW, when _RS expands into real*4
206      (#undef REAL4_IS_SLOW)
207     - pkg/mnc/mnc_cw_cvars.F: write _RS-grid variable as netcdf-real*4
208       coordinate variable
209     - pkg/monitor/monitor.F: call mon_printstats_rs for _RS variables
210       (surface forcing fields), removes some of the functionality
211     - model/src/write_grid.F: write _RS variables as RS-netcdf variables
212     - model/src/diags_oceanic_surf_flux.F: fix a few diagnostics, but
213       probably not all: copy _RS surface flux fields to RL field before
214       calling DIAGNOSTICS_FILL
215     - model/inc/CPP_EEMACROS.h: make the _EXCH_*_RS macros expand into
216       real*4 (CALL EXCH_*_RS) versions.
217     - tested in 2 experiments: inverted_barometer & solid-body.cs-32x32x1
218       (generate new output).
219    
220    checkpoint61e
221    o gmredi (Visbeck-variable-K):
222      - fix the "OLD_VISBECK_CALC" option + move the #undef to GMREDI_OPTIONS.h
223      - change computation of Visbeck-K where Slope > Smax :        
224        now: N*min(Slope,Smax) (similar to OLD_VISBECK_CALC with gkw91 taperFct)
225        previously was: M*sqrt(min(Slope,Smax)) ; update ideal_2D_oce output.
226    o bottom-drag in P-coord.: change units of bottomDragLinear & bottomDragQuadratic
227      to be the same as with Z-coord. (i.e., m/s & no-units respectively).
228    o calc_viscosity: fix index range (was wrong for CD-Scheme & explicit Visc)
229    o pkg gmredi, linear & fm07 tapering: put an upper limit on Slope^2 so that Kwz
230      does not become too large (to avoid precision problems in implicit-diffusion
231      solver); update ideal_2D_oce output.
232    o pkg ggl90:
233      - in ggl90_cal.F: compute only one set of visc/diff, used to
234        step forward TKE/T/S/momentum
235      - default value of GGL90viscMax/GGL90diffMax larger and equal (to
236        have Prandlt =1)
237    o gmredi (fm07):
238      - add lower limit on Transition layer thickness + lower limit on 1/lambda
239      - allow to set few more params from data.gmredi
240    o pkg ggl90:
241          - add missing sqrt(2) factor in computation of mixing length
242            scale in ggl90_calc.F
243            (see Eq. (2.35) in Blanke and Delecluse, JPO, 1993)
244    
245    checkpoint61d
246    o adjoint modifs
247      - adapt to new field rhoInSitu
248      - allow to revert from StoreDynVars2D, StoreDynVars3D
249        to individual fields
250        (switches AUTODIFF_USE_OLDSTORE_2D, AUTODIFF_USE_OLDSTORE_3D)
251    o store 3-D (in-situ) density in commom block (DYNVARS.h)
252     - save 1 rho computation.
253     - no longer modify T & S inland when using down_slope pkg.
254    o gad_som : allow to switch to SOM advection for T & S without corresponding
255      pickup files (using new S/R: MDS_CHECK4FILE).
256    o keep removing _BEGIN/_END_MASTER between MDSIO S/R call
257      (for useSingleCpuIO to work in multi-threaded)
258    o pkg/shelfice:
259     - disentangle parameters: retire shelfIceFile in data/PARM05 and substitute
260       with SHELFICEtopoFile in data.shelfice
261     - move code for reading ice shelf topography and modifying hFacC into separate
262       subroutines of the shelfice package
263     - separate initialisation of "fixed" and "variable" variables
264     - kTopC is not (yet) replaced by kSurfC, although this is possible
265    o fix broken diagnostic SIatmFW and changed the units from m/s to kg/m^2/s
266    o add initialisation of deepFac if using Pcoords (new S/R set_grid_factors)
267      and rename:
268      ini_phiref.F          --> set_ref_state.F  (+ set anelastic factors)
269      ini_reference_state.F --> load_ref_files.F (- set anelastic factors)
270    o change FIND_RHO_SCALAR : return rho (instead of rho - rhoConst)
271    
272    checkpoint61c
273    o add option ALLOW_ADDFLUID and run-time selector "selectAddFluid"
274      to add mass source/sink of fluid in the interior (3-D generalisation
275      of oceanic real-fresh water flux).
276    o pkg/monitor: fix recip_dr argument in mon_advcflw2 call (-> advcfl_W_hf_max)
277    o vermix test experiment: switch from LINEAR EOS to MDJWF (not tested before)
278    o new pkg: down_slope
279      density-driven downsloping flow parameterisation (Tellus, 1999, 412-430).
280      - presently, modify T & S inland (dwnslp_calc_rho.F).
281        can be avoided if rho was stored in common bloc and used in calc_phi_hyd.
282      - implementation for ocean in P-coords not tested.
283      - add secondary test to global_ocean.90x40x15 (input.dwnslp) to check this pkg.
284    o move kSurfC,kSurfW & kSurfS from SURFACE.h to GRID.h (next to kLowC)
285      (to facilitate conversion of oceanic pkg to p-coordinate).
286    o change FIND_RHO to FIND_RHO_2D :
287      pass arguments tFld & sFld as local (tile) 2D-arrays ; +fix EOS_CHECK.
288    o new pkg: cheapaml
289    
290    checkpoint61b
291    o exch2:
292      change index bounds in rx2_cube exchanges (new S/R: EXCH2_GET_UV_BOUNDS)
293      - no longer depend on the order sequence (N,S,E,W).
294      - 3rd exchange no longer needed (tested with 24 tiles).
295      - same modif to hand-written adjoint S/R (global_ocean.cs32x15: zero diff)
296      - exch_UV_A-grid readily available (but not yet tested).
297      comment out 3rd exch call and remove copy of V to U at SE & NE face-corners
298      (both are no longer needed) in exch_uv: no effect in forward run; but
299       get truncation error differences with adjoint exch.
300    o exch2:
301     - change index-bounds storage (moved from target to local tile,
302       more intuitive this way) ; rename/remove some variables.
303     - change utils/exch2/matlab-topology-generator & code_mods accordingly.
304    
305    checkpoint61a
306    o exch2/matlab-topology-generator:
307      fix bug in index bounds for N-S internal exchange (hard to tell if this
308      had any effect on the solution).
309    o pkg/mnc: small fix in mnc_handle_err: make sure that length of msgbuf
310      is exceeded
311    o testreport: modify syntax of awk, so that it works also with
312      awk, and not only with nawk and gawk
313    
314    checkpoint61
315    o autodiff
316      Change default for TAF usage by removing argument '-nonew_arg'
317      This changes TAF subroutine argument list (ordering of active var.)
318      Omission of this flag leads to incompatibility w.r.t. TAMC
319      To restore TAMC-compatibility, need following modifs:
320      * use adoptfile tools/adjoint_options/adjoint_tamc_compatibility
321      * use CPP option #define AUTODIFF_TAMC_COMPATIBILITY
322      Tested TAF version is 1.9.22
323      N.B.: exch2 hand-written adjoint code currently not TAMC compatible
324    
325    checkpoint60
326    --->>> This is the last checkpoint for which generating the adjoint
327    --->>> is done via the TAMC-compatibility TAF-flag "-nonewarg"
328    --->>> by default.
329    --->>> Future checkpoints will ommit this flag by default.
330    o pkg/generic_advdiff, os7mp: simplify and avoid division by zero.
331    o pkg/grdchk: Add gradient check to etan0 perturbations.
332    o pkg/exf: read ice area from file (#ifdef ALLOW_ICE_AREAMASK).
333  o pkg/autodiff: more flexibility in zeroadj.  o pkg/autodiff: more flexibility in zeroadj.
334    model/src: modifying do_oceanic_phys.F accordingly.    model/src: modifying do_oceanic_phys.F accordingly.
335  o pkg/ctrl: Bounds for (kapgm,kapredi,eddypsi) parameter estimation.  o pkg/ctrl: Bounds for (kapgm,kapredi,eddypsi) parameter estimation.
# Line 20  o bridging the gap between eddy stress a Line 346  o bridging the gap between eddy stress a
346    -> ALLOW_EDDYPSI_CONTROL for estimation purpose.    -> ALLOW_EDDYPSI_CONTROL for estimation purpose.
347    The key modifications are in model/src/taueddy_external_forcing.F    The key modifications are in model/src/taueddy_external_forcing.F
348    pkg/gmredi/gmredi_calc_*F pkg/gmredi/gmredi_*transport.F    pkg/gmredi/gmredi_calc_*F pkg/gmredi/gmredi_*transport.F
349    
350  checkpoint59r  checkpoint59r
351  o new S/R (pkg/mdsio/mdsio_facef_read.F) reading of sNx+1 x sNy+1 array  o new S/R (pkg/mdsio/mdsio_facef_read.F) reading of sNx+1 x sNy+1 array
352    from per-face file (take S/R READSYMTILE_RS out of ini_curvilinear_grid.F);    from per-face file (take S/R READSYMTILE_RS out of ini_curvilinear_grid.F);
# Line 2144  o add two new packages Line 2471  o add two new packages
2471    - my82 (Mellor and Yamada, 1982) level 2 turbulence closure scheme    - my82 (Mellor and Yamada, 1982) level 2 turbulence closure scheme
2472    - these scheme operate in a similar manner as KPP (without ghat, of    - these scheme operate in a similar manner as KPP (without ghat, of
2473      course), however there is only one diffusivity that is used for both      course), however there is only one diffusivity that is used for both
2474      temperature and salinity (and passive tracers)ç      temperature and salinity (and passive tracers)
2475  o add offline package  o add offline package
2476  o add "verification/dic_example" which contains a dic-gchem-ptracer setup  o add "verification/dic_example" which contains a dic-gchem-ptracer setup
2477    provided by Stephanie Dutkiewicz    provided by Stephanie Dutkiewicz

Legend:
Removed from v.1.1000  
changed lines
  Added in v.1.1091

  ViewVC Help
Powered by ViewVC 1.1.22