/[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.1090 by jmc, Wed Apr 22 19:25:06 2009 UTC revision 1.1183 by mlosch, Fri Oct 23 08:11:44 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 tools/gemake2: delete modifs from March 27-29 2009 (with option "-mpi":  o pkg/seaice/autodiff:
8     - add CPP-flag SEAICE_ALLOW_CLIPVELS (undefined by default)
9     - change seaice_calc_viscosities/strainraties for better
10       vectorization after AD by TAF
11     - rearrange common blocks in SEAICE.h once again for easier dumps of
12       adjoint output, adjust adcommon.h and g_common.h
13     - change flux form of default advection scheme for better vectorization
14       of adjoint code
15     - modify lab_sea/input.lsr to test the flux form of default advection
16       scheme
17    o autodiff: fix seaice fields in g_common.h (left over from 3 to 1 level
18      transition in checkpoint61r)
19    
20    checkpoint61x
21    o adjoint: adding two new verification experiments
22      * tutorial_global_biogeo: an online dic sensitivity experiment
23        (work with G. Forget)
24      * tutorial_dic_adjoffline: an offline dic sensitivity experiment
25        (an update of the Dutkiewicz et al., 2006 paper)
26      Both use the same CO2FLUX cost function, but different controls
27      Gradient checks look good, but some recomp. could still be traightened out
28      Should add xx_dic control variable to test at some point
29    o model/src:
30      add run-time parameters to check and/or mask initial Temp & Salt
31    o pkg/ptracers:
32      apply zonal filter (if used) to all passive tracers
33    
34    checkpoint61w
35    o adjoint:
36      * update and finish work on adjoint of dic package (G. Forget)
37      * start adjoint of offline package
38      * bug fix for cube-sphere TLM code (active_file_g.F)
39    o tutorial_global_oce_optim:
40      now tested with 4 tiles (change data.grdchk and update output_adm.txt)
41      and add SIZE.h_mpi for MPI testing
42    o pkg autodiff
43     - fix bug in autodiff_restore
44     - make genmake2 remove comments and empty lines from ad_input_code.f in
45       order to reduce file size (new script file remove_comments_sed)
46    o allows to specify vertical profile of vertical viscosity (viscArNr)
47    o pkg autodiff
48      autodiff_ini_model_io: comment out writing of grid information to netcdf
49      files, because at this time the grid fields are a bunch of zeros
50      (this is analogous to ini_mnc_vars.F)
51    o pkg seaice/autodiff:
52     left-over from removing the 3-time-level fields:
53     - small rearrangement of common blocks (SEAICE.h)
54     - fix field declaration in order to get meaningful output in
55       adseaice.*.nc (adcommon.h)
56    o pkg obcs:
57     - go through NH code only if nonHydrostatic=T
58     - fix restart for Orlanski (Hydrostatic; still to fix NH+Orlanski restart)
59     - move call to OBCS_CALC before SEAICE_MODEL which needs seaice-obc fields
60       => Update output of experiment seaice_obcs
61     - useOBCSbalance: code is wrong with multi-tile / proc: add a stop.
62     - fix NH+Orlanski restart (add wvel OBC value in pickup file)
63     - fix restart for useOBCSprescribe + useEXF
64    o pkg ctrl
65     - make obcs as control parameter work also with useSingleCPUio
66     - replace a few sny and snx by Ny and Nx to be consistent with
67       ctrl_set_globfld_x/yz.F
68    o pkg seaice
69     - re-introduce global field frWtrAtm for the case of
70       ALLOW_MEAN_SFLUX_COST_CONTRIBUTION defined (SEAICE.h, seaice_growth.F)
71     - catch the case of an empty filename in seaice_cost_init_fixed.F
72    o pkg ecco:
73     - small bug fixes in the_main_loop
74      . fix order of monitor arguments
75      . change SEAICE_MULTILEVEL to SEAICE_MULTICATEGORY
76      . move include DYNVARS.h up outside the ALLOW_AUTODIFF_TAMC block so
77        that ALLOW_ECCO_FORWARD_RUN works with NONLIN_FRSURF
78     - replace cos(yc*deg2rad) by _rA in cost_mean_heat/saltflux.F, fix the
79       imbalance diagnostics for wmean_h/sflux = 0
80     - add actual values to f_obcsn/s/w/e in ecco_cost_final.F
81     - catch the case of empty filenames in ecco_cost_init_fixed.F
82       and ecco_readparms.F
83     - fix the logic of wsalt=0.and.wsalt2=0 to wsalt=0.or.wsalt2=0 (and theta)
84       to make this statement meaningfull in ecco_cost_weights.F
85    o pkg zonal_filt:
86      - change ZONAL_FILTER S/R interface (allows to filter 2-D fields)
87      - when using rStar, filter uVel*hFacW instead of uVel
88       (tested without topography: fix a Pb of growing instability with NLFS)
89    o pkg/aim/thsice:
90     -Add capacity to read a Q-Flux to be applied in slab-mixed layer thsice_slab_ocean.F
91      (works as SST restoring)
92    o pkg/thsice:
93     -add specific parameter for salt restoring in slab mixed-layer
94       (by default equal to the temperature restoring)
95    
96    checkpoint61v
97    o pkg/kpp:
98     -add double diffusive contributions as a hack, turned off by default
99      for now, and the code can be excluded with a CPP-flag
100      EXCLUDE_KPP_DOUBLEDIFF just as EXCLUDE_KPP_SHEAR_MIX
101     -add a verification experiments to "vermix" to test this code
102    o pkg/mom_common:
103      Move computation of length scales to mom_init_fixed.F (avoid recomputation
104      at each time-steps*levels).
105    o pkg/layers:
106      Created a new package for computing volume fluxes in temperature or
107      [not yet implemented] density layers, e.g. for calculating residual
108      overturning circulations.
109    o pkg/diagnostics:
110      implement a RS type version for the set of DIAGNOSTICS-FILL subroutines
111    o fix few RS/RL type mismatch (check with ifort 11.1 "-warn all" option)
112    o pkg/mdsio (& pkg/rw): rework high level S/R interface:
113      To avoid mixing type (RS/RL) of input/output array argument,
114      replace single mixed-type array with a pair from each type (RS/RL).
115    o setting pickupSuff force to start from a pickup (even if nIter0=0): this
116      is now implemented for all pkgs (previously only done for the main pickup)
117    
118    checkpoint61u
119    o pkg/obcs:
120      - change OB?eta from _RS to _RL and put them into a separate commen block
121      - add code to read OB?eta from a file via obcs_prescribe_read: the code
122        compiles and does not destroy any other experiments, the data is read
123        properly, but the code is not tested in actual applications
124    o ini_mnc_vars.F: if available (ALLOW_CAL), use startTime_1 and 2 to
125      create a meaningful unit attribute for the time coordinate variable
126    o pkg/ebm: add plenty of missing "_d 0" ;
127       update results of experiment global_ocean_ebm (fwd + adjoint tests)
128      Notes: several pieces of code only works for (very) specific set up
129       (domain size, lat range, etc ...); should be documented & stop if one
130       of those assumptions is not meet.
131    o pkg/seaice: add missing "_d 0" in HSNOW initialisation
132       update results of experiments: lab_sea (3 fwd tests + 2 adjoint tests)
133          and offline_exf_seaice (seaicetd + adjoint).
134    o add missing "_d 0" in ShortWave flux penetration (S/R SWFRAC):
135       update results of experiments: lab_sea (4 fwd tests + 3 adjoint),
136        seaice_obcs, global_ocean.cs32x15 (icedyn & thsice) and natl_box (x2).
137    o few changes to pass when compiling with strick checking of arguments
138       across S/R (allow to find few bugs in S/R calls).
139    o pkg/seaice:
140     - seaice_evp: move if-statement out of loop (TEM) at the cost of 11 new
141       2d-fields, for slightly better performance
142    
143    checkpoint61t
144    o pkg/seaice:
145     - remove unused variables in seaice_evp
146     - reduce number of SQRT in seaice_evp and seaice_calc_viscosities
147    o eeset_parms.F, open_copy_data_file.F & ini_parms.F:
148       named file replacing SCRATCH files (for TARGET_BGL & TARGET_CRAYXT):
149       each processor opens a different file (a fix from Matt).
150    o cg2d:
151       use simple EXCH (overlap size = 1 and ignore corners), like in cg3d.F :
152       this reduces number of EXCH calls by 2 when using exch2.
153       For now, keep previous size array for gc2d_nsa.F to avoid Pb with store
154       directives and missing simple EXCH (exch_s3d_rx) adjoint code.
155    o pkg/thsice:
156       add code for Near Infra-Red albedo (from Jeff)
157       (turned on with run-time flag: thSIce_calc_albNIR, in data.ice).
158    o ini_model_io:
159       move (from packages_readparms.F & initialise_fixed.F) to ini_model_io.F
160       MNC & MONITOR initialisation calls; move MNC calls from {pkg}_readparms.F
161       to {pkg}_init_fixed.F for {pkg}= kpp, seaice & shelfice.
162    
163    checkpoint61s
164    o pkg/exch2 (exch2_rxN_cube*):
165       remove last BARRIER (no need to synchronise after getting data from shared
166       buffer (get) as long as any change to buffer (put,recv) is between BARRIER)
167    o change tiling of natl_box (from 2 to 4 tiles) to get it tested with MPI+MTH
168      and generate new output (for the 2 tests) using default rotationPeriod.
169    o pkg/exch2:
170      - always call exch2_*_cube, not exch-1 anymore, if useCubedSphereExchange=F
171        (was already the case with AUTODIFF_EXCH2 defined)
172      - add bj in exch2 arrays and S/R.
173      - can now use pkg/exch2 for regular (non CS) set-up; switch test experiment
174        MLAdjust to exch2 (identical results except with MPI -> truncation diff);
175    o add pkg/longstep for ptracer timestep a multiple of U/V/T/S timestep
176      - adds deltaT argument to the subroutines TIMESTEP_TRACER DWNSLP_APPLY
177        GAD_ADVECTION GAD_CALC_RHS GAD_IMPLICIT_R GAD_SOM_ADVECT
178      - verification experiment natl_box.longstep
179    o build_options: ifort+mpi_aces optfile (following Constantinos advise):
180      allow different versions of ifort (8 or 9) and different versions
181       of mpi (mpich or mpich2) to be used with the same optfile.
182      requires to load a netcdf module and set MPI_INC_DIR (see comments inside
183       optfile). Change testing scripts (example_scripts/ACESgrid) accordingly.
184    o pkg/seaice:
185      fix previous modif (3-level field removed) for old-pickup reading
186    
187    checkpoint61r
188    o pkg/seaice:
189     - clean up the 3-time levels of UICE,VICE,HEFF,AREA in three steps:
190      1. introduce fields for time step (N-1)
191      2. do not use levels 2 and 3 anymore (but keep them)
192      3. replace 3D versions of UICE,VICE,HEFF,AREA by 2D versions. This
193         includes changes to pkg/thsice/thsice_get_velocity.F,
194         model/src/external_forcing_surf.F, pkg/obcs/obcs_apply_seaice.F, and
195         pkg/obcs/obcs_apply_uvice.F.
196     - completely remove code for leap-frog time stepping for second order
197       advection
198     - adjusted adjoint accordingly;
199       gradient checks have improved by an order of magnitude!
200    
201    checkpoint61q
202    o pkg/seaice: fix bug in free slip boundary conditons for LSOR solver
203     - the boundary conditions along open boundaries are unclear, because the
204       associated masks are set to zero and are thus treated as closed when
205       applying the slip-boundary conditions.
206     - move initialisation of grid variables from seaice_init_varia to _fixed
207    o adding doc/diags_changes.txt
208    o model/src:
209      do_the_model_io: rename few pkg S/R which are called from do_the_model_io,
210       to avoid "_diags" naming if no relation to pkg/diagnostics.
211      ini_theta,ini_salt: fix (at least compile the code) reading from NetCDF
212       file and only call PLOT_FIELD_XYZ if debugMode=T.
213      remove few PLOT_FIELD_XY calls outside initialisation.
214      remove few unnecessary BARRIER.
215    o pkg/salt_plume: add SaltPlumeSouthernOcean to allow disabling
216      of salt_plume package in Southern Ocean
217    o pkg/thsice:
218     - fix hOceMxL in overlap (needed for sea-ice advection);
219       update output of experiment global_ocean.cs32x15.icedyn.
220    o pkg/seaice:
221     - make the (placeholder-)tracer IceAge a little more meaningful
222      by avoiding ridging of age due to convergent ice flow
223     - add diagnostics for surface temperature, atmospheric and oceanic
224       heat flux, and atmospheric and oceanic ice growth
225    
226    checkpoint61p
227    o eesupp (global_max & global_sum):
228     - add one more element to buffer (start at index 0) for GL_MAX/SUM output;
229     - remove starting & endding barrier (no longer needed)
230    o pkg/mdsio, writelocal:
231       works also in multi-threaded when LOCBIN_IO_THREAD_SAFE is undefined
232       (remove the STOP): uses shared buffer IO to store data from all threads
233       and then let master-thread write nThreads tiles.
234    o pkg/mdsio (read/write field):
235     - tiled IO done in 1 piece (all levels at a time);
236     - new header file "MDSIO_BUFF_3D.h" with 3-D buffers.
237     - multi-threaded: allow to read/write local (non-shared) array
238       (was already working with singleCpuIO ; now works also without);
239       => no longer needs barrier call after reading a file (added inside
240       mds_read/write field).
241     - move barrier calls outside gather/scatter_2d to mds_read/write_field
242    o pkg/obcs/exf
243      - exf_getffieldrec.F: add a few ". _d 0" and
244                            add code to catch the case fldperiod = 0.  
245      - obcs_precribe_read.F: change if statements so that they match
246                              the variable type (_RL)
247      - add a separate exf_iprec_obcs to exf_set_obcs (by default equal to
248        exf_iprec) as part of namelist EXF_NML_01
249    o pkg/mdsio:
250      - read/write tiled (local) files: read/write 1-level tile chunk at a time
251        (instead of segment of length sNx); expected to speed up tiled IO.
252    o pkg/exch2:
253      - take buffer copy from/to array out of S/R exch2_send/recv into new
254        S/R exch2_put/get ; Exch of local variable now works with MPI+MTH
255        (tested by removing commom block statement in SOLVE_FOR_PRESSURE.h).
256      - ad version of send/recv no longer needed (but needs exch2_ad_put/ad_get)
257      - simplify argument list of S/R exch2_get_uv_bounds and exch2_get_scal_bounds
258        (which replaces exch2_get_recv_bounds & exch2_get_send_bounds) by
259        including Topology header file (and common blocks).
260      - implement EXCH_IGNORE_CORNERS in scalar exchange (rx1);
261        remove 2nd exch call in exch2_s3d_rx (no longer needed).
262      - document S/R arg. list (+ update some comments)
263    o pkg/seaice:
264      - added SEAICE_CLIM_AIR code, which permits modulation of surface air
265        temperature and humidity over sea ice based on climatological values.
266      - clean up computation of Hibler+Bryan (1987) stress coupling for the case
267        of LSR (this change is expected to change the results slightly because now
268        it uses slightly difference moduli from the second last LSR solution, but
269        that is more consistent with the stress computations; this part of the code
270        is not tested in the verification experiments)
271      - fix diagnostics for flux of ice in case of B-grid
272      - turn strain rates eij and press into global fields within global
273        common block in SEAICE.h for more straighforward diagnostics
274      - fix diagnostics for principal stress components, press, zeta, eta
275        to give something meaning full
276    o change tiling of global_ocean.cs32x15 (12 tiles, 32x16) and generate
277      new output (all 4 tests) using default rotationPeriod.
278    o pkg KPP:
279      - add flag (LimitHblStable) to avoid limiting BL depth under stable conditions
280      - fix vertical indices in RI_IWMIX (match between KPP grid from 0 to Nr+1
281        and normal grid from 1 to Nr)
282    o Fix computations of L4rdt in mom_calc_visc.F
283    o eesupp:
284      new S/R ALL_PROC_DIE (like EEDIE, without irrelevant threads err msg):
285      should be called before a "stop" if we know that all process are going
286      to stop; if one Proc does call this S/R, will hang in MPI_FINALISE.
287    
288    checkpoint61o
289    o pkg/exch2: change check_depths to issue just a warning (instead of error+stop)
290    o model/src/ini_curvilinear_grid.F:
291       compute angleSin/CosC when they are not read from a file
292    o pkg/seaice:
293       add pseudo time stepping for LSOR
294    o pkg/exch2:
295       temporary fix for exch2_s3d_rx; update results of exp. hs94.cs-32x32x5.impIGW
296    o pkg/seaice:
297       Adding separate sea ice albedo, drag coefficients, and HO in Southern
298       Ocean from those in Arctic Ocean.  See *_south in SEAICE_PARAMS.h
299    o SingleCpuIO:
300       fix bug put before chkp61n (May 11), when more than 2 MPI cpu are used.
301    o pkg/seaice:
302       put exchanged variable in (local) common block for MPI+Multi-Threads
303       to work (necessary with exch1, and for now with exch2 too).
304    o verification: to get more experiments tested with mpi, mth and mpi+mth,
305      change tiling of:
306       aim.5l_Equatorial_Channel, aim.5l_LatLon, dome, exp4, global_ocean.90x40x15,
307       hs94.cs-32x32x5, MLAdjust, offline_exf_seaice & tutorial_plume_on_slope
308      update output of exp:
309       exp4, global_ocean.90x40x15(+.dwslp), MLAdjust(all 5),
310       tutorial_plume_on_slope and hs94.cs-32x32x5 (use default rotationPeriod)
311      + ideal_2D_ocean: use default rotation period + test Global_Sum_SINGLECPU
312                       (note: output.txt generated on faulks using only 1 tile)
313    o eesupp & mdsio (SingleCpuIO):
314     - remove MISSING_TILE_IO (from MDSIO_OPTIONS.h), replaced by run-time
315       parameter "W2_useE2ioLayOut" (read from file "data.exch2").
316       Note: was not effective for globalFile & useSingleCpuIO=F; now fixed.
317     - move MDSIO_SCPU.h (pkg/mdsio) to EEBUFF_SCPU.h (eesupp/inc).
318     - fix GLOBAL_SUM_SINGLECPU when using Exch2; re-use same buffers and same
319       gather/scatter S/R as with SingleCpuIO (=> 1 less 2D global RL array).
320    o pkg/exch2:
321     - new code to set-up W2-Exch2 topology (replace matlab-topology-generator):
322       read parameter file "data.exch2" if it exists; otherwise try default regular
323       cube without blank-tile; new header file "W2_EXCH2_SIZE.h";
324     - update matlab-topology-generator (exch2_setup_cs6_print.m) to be compatible
325       with new code (not to be used, just as a backup solution).
326     - add examples (e.g., for lat-lon-cap_120) of parameter file "data.exch2"
327       in utils/exch2/input.
328    
329    checkpoint61n
330    o eesupp & mdsio (SingleCpuIO):
331      - move mapping to global io-buffer inside gather_2d/scater_2d: this save
332        memory (1 less 2D global RL array); skip 1 buffer copy; and only send/
333        receive real*4 array (instead of real*8) to read/write 32.bit files.
334        (new template for gather_2d/scater_2d, new set of S/R MDS_PASS_R4,8toRL,S
335         and remove new mapping S/R)
336    o pkg/mdsio: in preparation for topology-generator within pkg/exch2:
337      - Change declaration of SingleCpuIO buffer, reduced to only 1 dim.:
338        xyBuffer_size, which is set to W2_ioBufferSize from W2_EXCH2_PARAMS.h
339        when using exch2. Add S/R to map global io-buffer to global model array.
340      - uses exch2_global_Nx,y instead of exch2_domain_nxt,nyt.
341    
342    checkpoint61m
343    o write_utils : change WRITE_1D_R8 to WRITE_1D_RL
344      and add 3 other S/R for other types (R4,R8,RS).
345      Change calls to WRITE_0D/1D_R[4,8,S,L] to match the type of argument.
346    o eesupp/inc/CPP_EEMACROS.h (EXCH & GLOBAL_SUM/MAX macros):
347     - ensure suffix (_R4,_R8,_RS,_RL) correspondance between macro and S/R ;
348     - reduce number of macro definitions (keep only those which are used);
349       can always call directly the corrresponding S/R if needed.
350     - everywhere in the code: change macros (EXCH & GLOBAL_SUM/MAX) suffix
351       _R4/_R8 to _RS/_RL when applied to _RS/_RL variable.
352     Fixed wrong type of arg in:
353       atm2d/put_ocnvars.F cheapaml/cheapaml.F ocn_compon_interf/ocn_apply_import.F
354       ctrl/ctrl_map_ini_ecco.F seaice/seaice_do_diags.F seaice/seaice_model.F
355     To fix: global_sum of var in common block is wrong (multi-threaded) in:
356       dic/dic_atmos.F ebm/ebm_area_t.F ebm/ebm_zonalmean.F sbo/sbo_calc.F
357     To fix: aim_v23/aim_do_co2.F (BARRIER & GLOBAL_SUM within bi,bj loops)
358    o adjustment.cs-32x32x1: test "blank tiles" when compiling with mpi; 2nd
359      "nlfs" mpi test is going to stop when checking for unconnected tile egdes.
360    o eesupp/src: new S/R to stop everyone when 1 (or more) Proc. find an error
361    o pkg/exch2: new S/R EXCH2_CHECK_DEPTHS to check that unconnected tile
362      edges are closed (zero depth); might happen with "blank tiles";
363      for now, disabled if using OBCS.
364    o tools/gemake2: delete modifs from March 27-29 2009 (with option "-mpi":
365    link *_mpi files from "MODS" dir) to return to a strait-forward behavior.    link *_mpi files from "MODS" dir) to return to a strait-forward behavior.
366  o move read_write.F to write_utils.F ; new parameter (maxLengthPrt1D)  o move read_write.F to write_utils.F ; new parameter (maxLengthPrt1D)
367    to reduce length of 1.D array print in STDOUT.    to reduce length of 1.D array print in STDOUT.
368  o adjustment.cs-32x32x1: changed to oceanic test with 1 continent  o adjustment.cs-32x32x1: changed to oceanic test with 1 continent
369    + use 48 tiles (16x8): this will allow to test "blank tiles".    + use 48 tiles (16x8): this will allow to test "blank tiles".

Legend:
Removed from v.1.1090  
changed lines
  Added in v.1.1183

  ViewVC Help
Powered by ViewVC 1.1.22