/[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.1045 by mlosch, Wed Nov 26 15:27:02 2008 UTC revision 1.1132 by jmc, Thu Jun 25 15:55:40 2009 UTC
# Line 4  $Name$ Line 4  $Name$
4      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
5      ==============================      ==============================
6    
7    checkpoint61r
8    o pkg/seaice:
9     - clean up the 3-time levels of UICE,VICE,HEFF,AREA in three steps:
10      1. introduce fields for time step (N-1)
11      2. do not use levels 2 and 3 anymore (but keep them)
12      3. replace 3D versions of UICE,VICE,HEFF,AREA by 2D versions. This
13         includes changes to pkg/thsice/thsice_get_velocity.F,
14         model/src/external_forcing_surf.F, pkg/obcs/obcs_apply_seaice.F, and
15         pkg/obcs/obcs_apply_uvice.F.
16     - completely remove code for leap-frog time stepping for second order
17       advection
18     - adjusted adjoint accordingly;
19       gradient checks have improved by an order of magnitude!
20    
21    checkpoint61q
22    o pkg/seaice: fix bug in free slip boundary conditons for LSOR solver
23     - the boundary conditions along open boundaries are unclear, because the
24       associated masks are set to zero and are thus treated as closed when
25       applying the slip-boundary conditions.
26     - move initialisation of grid variables from seaice_init_varia to _fixed
27    o adding doc/diags_changes.txt
28    o model/src:
29      do_the_model_io: rename few pkg S/R which are called from do_the_model_io,
30       to avoid "_diags" naming if no relation to pkg/diagnostics.
31      ini_theta,ini_salt: fix (at least compile the code) reading from NetCDF
32       file and only call PLOT_FIELD_XYZ if debugMode=T.
33      remove few PLOT_FIELD_XY calls outside initialisation.
34      remove few unnecessary BARRIER.
35    o pkg/salt_plume: add SaltPlumeSouthernOcean to allow disabling
36      of salt_plume package in Southern Ocean
37    o pkg/thsice:
38     - fix hOceMxL in overlap (needed for sea-ice advection);
39       update output of experiment global_ocean.cs32x15.icedyn.
40    o pkg/seaice:
41     - make the (placeholder-)tracer IceAge a little more meaningful
42      by avoiding ridging of age due to convergent ice flow
43     - add diagnostics for surface temperature, atmospheric and oceanic
44       heat flux, and atmospheric and oceanic ice growth
45    
46    checkpoint61p
47    o eesupp (global_max & global_sum):
48     - add one more element to buffer (start at index 0) for GL_MAX/SUM output;
49     - remove starting & endding barrier (no longer needed)
50    o pkg/mdsio, writelocal:
51       works also in multi-threaded when LOCBIN_IO_THREAD_SAFE is undefined
52       (remove the STOP): uses shared buffer IO to store data from all threads
53       and then let master-thread write nThreads tiles.
54    o pkg/mdsio (read/write field):
55     - tiled IO done in 1 piece (all levels at a time);
56     - new header file "MDSIO_BUFF_3D.h" with 3-D buffers.
57     - multi-threaded: allow to read/write local (non-shared) array
58       (was already working with singleCpuIO ; now works also without);
59       => no longer needs barrier call after reading a file (added inside
60       mds_read/write field).
61     - move barrier calls outside gather/scatter_2d to mds_read/write_field
62    o pkg/obcs/exf
63      - exf_getffieldrec.F: add a few ". _d 0" and
64                            add code to catch the case fldperiod = 0.  
65      - obcs_precribe_read.F: change if statements so that they match
66                              the variable type (_RL)
67      - add a separate exf_iprec_obcs to exf_set_obcs (by default equal to
68        exf_iprec) as part of namelist EXF_NML_01
69    o pkg/mdsio:
70      - read/write tiled (local) files: read/write 1-level tile chunk at a time
71        (instead of segment of length sNx); expected to speed up tiled IO.
72    o pkg/exch2:
73      - take buffer copy from/to array out of S/R exch2_send/recv into new
74        S/R exch2_put/get ; Exch of local variable now works with MPI+MTH
75        (tested by removing commom block statement in SOLVE_FOR_PRESSURE.h).
76      - ad version of send/recv no longer needed (but needs exch2_ad_put/ad_get)
77      - simplify argument list of S/R exch2_get_uv_bounds and exch2_get_scal_bounds
78        (which replaces exch2_get_recv_bounds & exch2_get_send_bounds) by
79        including Topology header file (and common blocks).
80      - implement EXCH_IGNORE_CORNERS in scalar exchange (rx1);
81        remove 2nd exch call in exch2_s3d_rx (no longer needed).
82      - document S/R arg. list (+ update some comments)
83    o pkg/seaice:
84      - added SEAICE_CLIM_AIR code, which permits modulation of surface air
85        temperature and humidity over sea ice based on climatological values.
86      - clean up computation of Hibler+Bryan (1987) stress coupling for the case
87        of LSR (this change is expected to change the results slightly because now
88        it uses slightly difference moduli from the second last LSR solution, but
89        that is more consistent with the stress computations; this part of the code
90        is not tested in the verification experiments)
91      - fix diagnostics for flux of ice in case of B-grid
92      - turn strain rates eij and press into global fields within global
93        common block in SEAICE.h for more straighforward diagnostics
94      - fix diagnostics for principal stress components, press, zeta, eta
95        to give something meaning full
96    o change tiling of global_ocean.cs32x15 (12 tiles, 32x16) and generate
97      new output (all 4 tests) using default rotationPeriod.
98    o pkg KPP:
99      - add flag (LimitHblStable) to avoid limiting BL depth under stable conditions
100      - fix vertical indices in RI_IWMIX (match between KPP grid from 0 to Nr+1
101        and normal grid from 1 to Nr)
102    o Fix computations of L4rdt in mom_calc_visc.F
103    o eesupp:
104      new S/R ALL_PROC_DIE (like EEDIE, without irrelevant threads err msg):
105      should be called before a "stop" if we know that all process are going
106      to stop; if one Proc does call this S/R, will hang in MPI_FINALISE.
107    
108    checkpoint61o
109    o pkg/exch2: change check_depths to issue just a warning (instead of error+stop)
110    o model/src/ini_curvilinear_grid.F:
111       compute angleSin/CosC when they are not read from a file
112    o pkg/seaice:
113       add pseudo time stepping for LSOR
114    o pkg/exch2:
115       temporary fix for exch2_s3d_rx; update results of exp. hs94.cs-32x32x5.impIGW
116    o pkg/seaice:
117       Adding separate sea ice albedo, drag coefficients, and HO in Southern
118       Ocean from those in Arctic Ocean.  See *_south in SEAICE_PARAMS.h
119    o SingleCpuIO:
120       fix bug put before chkp61n (May 11), when more than 2 MPI cpu are used.
121    o pkg/seaice:
122       put exchanged variable in (local) common block for MPI+Multi-Threads
123       to work (necessary with exch1, and for now with exch2 too).
124    o verification: to get more experiments tested with mpi, mth and mpi+mth,
125      change tiling of:
126       aim.5l_Equatorial_Channel, aim.5l_LatLon, dome, exp4, global_ocean.90x40x15,
127       hs94.cs-32x32x5, MLAdjust, offline_exf_seaice & tutorial_plume_on_slope
128      update output of exp:
129       exp4, global_ocean.90x40x15(+.dwslp), MLAdjust(all 5),
130       tutorial_plume_on_slope and hs94.cs-32x32x5 (use default rotationPeriod)
131      + ideal_2D_ocean: use default rotation period + test Global_Sum_SINGLECPU
132                       (note: output.txt generated on faulks using only 1 tile)
133    o eesupp & mdsio (SingleCpuIO):
134     - remove MISSING_TILE_IO (from MDSIO_OPTIONS.h), replaced by run-time
135       parameter "W2_useE2ioLayOut" (read from file "data.exch2").
136       Note: was not effective for globalFile & useSingleCpuIO=F; now fixed.
137     - move MDSIO_SCPU.h (pkg/mdsio) to EEBUFF_SCPU.h (eesupp/inc).
138     - fix GLOBAL_SUM_SINGLECPU when using Exch2; re-use same buffers and same
139       gather/scatter S/R as with SingleCpuIO (=> 1 less 2D global RL array).
140    o pkg/exch2:
141     - new code to set-up W2-Exch2 topology (replace matlab-topology-generator):
142       read parameter file "data.exch2" if it exists; otherwise try default regular
143       cube without blank-tile; new header file "W2_EXCH2_SIZE.h";
144     - update matlab-topology-generator (exch2_setup_cs6_print.m) to be compatible
145       with new code (not to be used, just as a backup solution).
146     - add examples (e.g., for lat-lon-cap_120) of parameter file "data.exch2"
147       in utils/exch2/input.
148    
149    checkpoint61n
150    o eesupp & mdsio (SingleCpuIO):
151      - move mapping to global io-buffer inside gather_2d/scater_2d: this save
152        memory (1 less 2D global RL array); skip 1 buffer copy; and only send/
153        receive real*4 array (instead of real*8) to read/write 32.bit files.
154        (new template for gather_2d/scater_2d, new set of S/R MDS_PASS_R4,8toRL,S
155         and remove new mapping S/R)
156    o pkg/mdsio: in preparation for topology-generator within pkg/exch2:
157      - Change declaration of SingleCpuIO buffer, reduced to only 1 dim.:
158        xyBuffer_size, which is set to W2_ioBufferSize from W2_EXCH2_PARAMS.h
159        when using exch2. Add S/R to map global io-buffer to global model array.
160      - uses exch2_global_Nx,y instead of exch2_domain_nxt,nyt.
161    
162    checkpoint61m
163    o write_utils : change WRITE_1D_R8 to WRITE_1D_RL
164      and add 3 other S/R for other types (R4,R8,RS).
165      Change calls to WRITE_0D/1D_R[4,8,S,L] to match the type of argument.
166    o eesupp/inc/CPP_EEMACROS.h (EXCH & GLOBAL_SUM/MAX macros):
167     - ensure suffix (_R4,_R8,_RS,_RL) correspondance between macro and S/R ;
168     - reduce number of macro definitions (keep only those which are used);
169       can always call directly the corrresponding S/R if needed.
170     - everywhere in the code: change macros (EXCH & GLOBAL_SUM/MAX) suffix
171       _R4/_R8 to _RS/_RL when applied to _RS/_RL variable.
172     Fixed wrong type of arg in:
173       atm2d/put_ocnvars.F cheapaml/cheapaml.F ocn_compon_interf/ocn_apply_import.F
174       ctrl/ctrl_map_ini_ecco.F seaice/seaice_do_diags.F seaice/seaice_model.F
175     To fix: global_sum of var in common block is wrong (multi-threaded) in:
176       dic/dic_atmos.F ebm/ebm_area_t.F ebm/ebm_zonalmean.F sbo/sbo_calc.F
177     To fix: aim_v23/aim_do_co2.F (BARRIER & GLOBAL_SUM within bi,bj loops)
178    o adjustment.cs-32x32x1: test "blank tiles" when compiling with mpi; 2nd
179      "nlfs" mpi test is going to stop when checking for unconnected tile egdes.
180    o eesupp/src: new S/R to stop everyone when 1 (or more) Proc. find an error
181    o pkg/exch2: new S/R EXCH2_CHECK_DEPTHS to check that unconnected tile
182      edges are closed (zero depth); might happen with "blank tiles";
183      for now, disabled if using OBCS.
184    o tools/gemake2: delete modifs from March 27-29 2009 (with option "-mpi":
185      link *_mpi files from "MODS" dir) to return to a strait-forward behavior.
186    o move read_write.F to write_utils.F ; new parameter (maxLengthPrt1D)
187      to reduce length of 1.D array print in STDOUT.
188    o adjustment.cs-32x32x1: changed to oceanic test with 1 continent
189      + use 48 tiles (16x8): this will allow to test "blank tiles".
190    o start to test MPI + multi-threaded (change testreport to skip tests
191      where number of tiles do not match).
192    o genmake2: add option "-omp" to turn on Open-MP ;
193      passed to testreport (if using -mth) ; merge few (*+mth) optfiles.
194    o pkg/exch2 (+ utils/exch2): ignore blank tiles (instead of removing them)
195      and keep the original tile numbering. Note: compatible with old
196      W2_EXCH2_TOPOLOGY.h & w2_e2setup.F topology files.
197    
198    checkpoint61l
199    o Overlaps had been forgotten in calculating ijk keys : fixed
200    o fix cyrus makedepend for gcc 4.4
201    o experiment fizhi-cs-aqualev20:
202      switching to new version of S/R GETPWHERE (#define TRY_NEW_GETPWHERE)
203      and generate new output. Note: new output is much closer to the old
204      one (v1.8) generated with pgf77 -O2 on Mar 14 2007 than it is from
205      the previous one (v1.9) generated with pgf77 -O0 on Jun 13 2007.
206    o pkg/fizhi:
207     - change to avoid indices going deliberately over array-bounds
208     - remove unnecessary MPI stuff
209     - perpetual spring equinox with "#define FIZHI_USE_FIXED_DAY"
210     - new version of S/R GETPWHERE available with #define TRY_NEW_GETPWHERE:
211       does not go over array-bounds and therefore returns reliable output.
212     - in fizhi_driver.F & fizhi_turb.F, fix "get_alarm" calls:
213       this S/R has 4 output arguments, which are used inside "get_alarm";
214       when it is called with the same variable "ndum" in place of several
215       output arguments, this produces weird side effect inside "get_alarm"
216      (e.g., datein=dates(n) is overwritten by timein=times(n) since both
217       datein & timein, resp. 2nd and 3rd arg., are pointed to the same "ndum"
218       address) and produces error.
219    o tools/gemake2: (March 27-29 2009)
220      with option "-mpi": link *_mpi files from "MODS" dir (this is what
221      testreport does).
222    o lab_sea experiment: change from 2x1 tiles to 2x2 tiles and regenerate output
223    o pkg/seaice: fix bug in new version of S/R SEAICE_LSR.F:
224       update lab_sea.lsr & lab_sea.salt_plume results
225    o pkg/seaice:
226      Major overhaul of the LSOR solver: change discretization from original
227      finite differences to finite volume. The old routine is still in the
228      file seaice_lsr.F and can be turned on with the CPP-flag
229      SEAICE_OLD_AND_BAD_DISCRETIZATION
230      Along with this change, the discretization of the strain rates is changed.
231      The same flag turns on the old discretization in seaice_calc_strainrates.F
232      All verification experiments with seaice are affected (except for
233      seaice_exf_offline); global_ocean.cs32x15/input.icedyn recovers the old
234      results by setting the above flag.
235      Further changes involve changes in defaults:
236      - ALLOW_FLOODING is defined and turned on by default
237      - SEAICEadvSnow = .true. is now the default
238      - SEAICEadvSalt = .true. is now the default
239      - SEAICEadvAge  = .true. is now the default
240      - SEAICE_clipVelocities  = .false. is now the default
241                                         (as per J. Zhang's recommendation)
242      - B-grid, and thus not tested: SEAICE_TEST_ICE_STRESS_1/EXPLICIT_SSH_SLOPE
243        is defined, SEAICE_TEST_ICE_STRESS_1 is renamed into SEAICE_BICE_STRESS
244      - seaice_growth: replace computation of UG by a simple copy from wspeed
245    
246    checkpoint61k
247    o pkg/seaice: a few bug fixes mostly concerning the metric terms:
248      seaice_lsr: fixed a few metric terms,
249                  affects lab_sea.lsr, lab_sea.salt_plume, seaice_obcs
250                  add more debugging diagnostics
251      seaice_evp: removed superfluous many metric terms of FV discretization
252                  of stress divergence
253                  affects lab_sea, lab_sea.hb87
254      seaice_ocean_stress: removed superfluous many metric terms of FV discretization
255                  of stress divergence for non-evp solver
256                  (this part is not tested anywhere)
257      seaice_calc_strainrates: fixed a bug in no-slip boundary conditions,
258                  they were effectively only half slip
259                  affects lab_sea.hb87
260      updated all relevant verification experiments (except lab_sea//output_adm.txt
261        which was left over, but fails at level 7 & 6)
262    
263    checkpoint61j
264    o pkg/seaice: turn a few hard wired parameters into run time parameters
265    o pkg/seaice/exf/monitor: turn netcdf coordinate variable T into model
266      time in seconds as promised by the usints in seaice_monitor and exf_monitor
267    o pkg/mnc: fix the coordinate variables for multi-tile curvi-linear grids
268      that are not cubed-sphere grids
269    o pkg/seaice: fix a bug for the rare case of SEAICE_VECTORIZE_LSR
270    o pkg/ctrl:
271      fix bug in ctrl_depth_ini.F (EXCH_XYZ was applied to 2-D field xx_r_low)
272    o non-hydrostatic (calc_gw.F):
273      fix viscA4 on CS-grid + fix diagnostics: Wm_Diss & Wm_Advec
274      update output of exp. global_ocean.cs32x15.viscA4
275    o remove backward compatibility for thetaMin,phiMin
276      (replaced by xgOrigin,ygOrigin on Jan 27, 2009)
277    o pkg/exf: fix a serious bug in exf_set_uv.F: do a rotation of the
278      interpolated wind field not only for curvilineargrid=.true., but also
279      for a rotated grid (rotateGrid=.true.)
280    o pkg/exf: add capability to bypass interpolation for runoff
281      by setting yet another CPP-flag
282    o eesupp/src, B-grid exchange:
283      add simple (6 tiles, non MPI) CS exch + standard interface S/R.
284    o autodiff
285      new CPP option ALLOW_TAMC_SINGLEPREC_COMLEV
286      allows to switch to single-prec common blocks;
287      corresponding switch occurs in tamc.h, integer isbyte
288      TAF store directivs "byte" and "kind" seem equivalent;
289      Most tamc.h have been updated correspondingly.
290    
291    checkpoint61i
292    o autodiff:
293      fix autodiff_sore/_restore indices for obcs
294      (spotted by I. Hoteit, M. Mazloff)
295    o pkg/autodiff, pkg/ctrl:
296      Add flag to reduce all tapelev I/O to single-prec.
297      In data.ctrl set
298        doSinglePrecTapelev = .TRUE.
299    o model/src/ini_curvilinear_grid.F, new grid-file format without EXCH2:
300      check in hack for SX8 compiler (for dxV, dyU); remove it afterwards.
301    o pkg/shelfice: fix bug with SHELFICEconserve flag, prepare for
302      real fresh water flux (not yet)
303    o pkg/mdsio & rw:
304     - add argument & option to S/R MDS_WRITEVEC_LOC (avoid Open/Close & write-meta)
305     - add simple interface S/R to MDS_WRITEVEC_LOC in pkg/rw: WRITE_GLVEC_RX.
306    o pkg/flt:
307     - replace FLT_INIT by 3-steps standard initialization S/R: readparms,
308       init_fixed & ini_varia.
309     - clean-up & simplify linear interpolation S/R:
310       move flt_bilinear.F -> flt_interp_linear.F ; add arg. myThid.
311     - fix some indices (mainly vertical index) ; add IMPLICIT NONE.
312     - more efficient I/O, using new options of MDS_WRITEVEC_LOC S/R.
313     - add vert. coordinate & horiz. indices to output files.
314     - change vertical position in initial condition, from index to depth [m].
315     - store horizontal (decimal) indices i,j in common block (instead of coord. x,y)
316       (easier to get to a curvilinear grid) ; delX,delY no longer used.
317    o pkg/ggl90: a few modifications:
318     - fix bottom boundary condition, and interpolation of KappaE
319       near top and bottom
320     - use same set of viscosity/diffusivity for TKE/momnetum/tracer
321     - add smoothing option as in Blanke and Delecluse (1993) (--> ALLOW_GGL90_SMOOTH)
322     - add a few options for mixing-length upper limit (--> mxlMaxFlag)  
323    o rename thetaMin,phiMin -> xgOrigin,ygOrigin
324      (maintain temporary backward compatibility, until next checkpoint)
325    
326    checkpoint61h
327    o pkg/diagnostics with "periodic averaging diagnostics":
328      fix S/R DIAGNOSTICS_COUNT ; hack diagnostics_out.F to get NetCDF output.
329    o pkg/my82: add missing "_d 0" ; update output of exp. vermix.my82
330    o pkg/obcs: fix loop ranges for balancing code
331    o pkg/flt:
332     - improve criteria for finding which tile owns a float.
333     - fix FLT_EXCH for case without MPI ; apply cyclic BC + fix also the MPI case.
334    o Add well-mixed CO2 box in aim_v23 (which includes changes in pkg
335      aim_compon_interf, ocn_compon_interf, gchem, and dic to get coupling
336      with ocean biogeo)
337    o pkg/flt:
338     - do some cleaning and fix the restart ;
339     - add output.txt in flt_example/results to get this pkg tested.
340    o genmake2: add an EXTENDED_SRC option and flag.
341       for now, cannot be turned on on command line (only through genmake_local)
342    o pkg/mdsio: read/write vector:
343     - singleCpuIO: fix RL/RS in MDSWRITEVECTOR by calling new S/R MDS_BUFFERtoRL/RS
344     - singleCpuIO=F: fix arrType & filePrec & _BYTESWAPIO in both MDSREADVECTOR
345       & MDSWRITEVECTOR by calling new S/R MDS_RD/WR_VEC_RL/RS.
346     - affects gradient of test-exp "bottom_ctrl_5x5" (not due to _BYTESWAPIO)
347      (drop to 4 digits matching for grad). But this test does not look very
348      reliable (was already failling on eddy, fc5, at 6.digits). Update output.
349    
350    checkpoint61g
351    o pkg/flt:
352      - move FLT_CPPOPTIONS.h to FLT_OPTIONS.h (standard name)
353       + make OPTIONS file more standard.
354      - rename S/R MDSREADVECTOR_FLT to FLT_MDSREADVECTOR (to mach src file name)
355      - fix _BYTESWAPIO CPP-Flag
356      - standard reading of file-parameter "data.flt"
357      - fix argument list of S/R FLT_RESTART
358      still need to fix _BYTESWAPIO in pkg/mdsio/mdsio_writevector.F
359  o pkg/obcs: fix a CPP-flag, bug reported by Jeff Blundell  o pkg/obcs: fix a CPP-flag, bug reported by Jeff Blundell
360  o gather_2d.F/scatter_2d.F: suppress vectorisation on NEC SX in respond to  o gather_2d.F/scatter_2d.F: suppress vectorisation on NEC SX in respond to
361    a still unfixed compiler/optimization bug    a still unfixed compiler/optimization bug
# Line 2291  o add two new packages Line 2643  o add two new packages
2643    - my82 (Mellor and Yamada, 1982) level 2 turbulence closure scheme    - my82 (Mellor and Yamada, 1982) level 2 turbulence closure scheme
2644    - these scheme operate in a similar manner as KPP (without ghat, of    - these scheme operate in a similar manner as KPP (without ghat, of
2645      course), however there is only one diffusivity that is used for both      course), however there is only one diffusivity that is used for both
2646      temperature and salinity (and passive tracers)ç      temperature and salinity (and passive tracers)
2647  o add offline package  o add offline package
2648  o add "verification/dic_example" which contains a dic-gchem-ptracer setup  o add "verification/dic_example" which contains a dic-gchem-ptracer setup
2649    provided by Stephanie Dutkiewicz    provided by Stephanie Dutkiewicz

Legend:
Removed from v.1.1045  
changed lines
  Added in v.1.1132

  ViewVC Help
Powered by ViewVC 1.1.22