/[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.1108 by dfer, Wed May 27 03:12:04 2009 UTC revision 1.1199 by mlosch, Tue Nov 24 14:54:25 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 mnc: add horizontal grid coordinates (XC,YC) and (XG,YG) to netcdf output
8      for complex grids, for now only to "state", "adstate", "adexf", "adseaice"
9    o model: add Christopher Wolfe's implemetation of the single reduction cg2d
10     (d'Azevedo, Eijkhout, and Romine, 1999) with CPP-flag ALLOW_SRCG and
11     runtime flag useSRCGSolver. Modify verification experiment global_with_exf
12     to test this code
13    o pkg/cd_code:
14      Change defaut D-grid time-stepping to get the same for D-grid momentum as for
15       C-grid momentum. D-grid velocity used to be stepped forward in time with
16       Adams-Bashforth only on surface pressure term. Tests show that using AB on
17       D-grid coriolis term improves stability (as expected from CD-scheme paper).
18      To reproduce old results, add 2 options (in CD_CODE_OPTIONS.h).
19      Update results of exp: global_with_exf (x2), ideal_2D_oce, isomip (x2),
20       lab_sea (3 Fwd +3 AD), natl_box (x2), tutorial_global_oce_biogeo (Fwd+AD),
21       tutorial_global_oce_latlon, tutorial_global_oce_optim & tutorial_tracer_adjsens.
22      Keep old results by adding CD_CODE_OPTIONS.h in exp: cfc_example, exp2,
23           lobal_ocean.90x40x15(+AD), global_ocean_ebm(+AD) and OpenAD.
24    o pkg/obcs:
25     - fix bug in initialisation (broken since Oct 1rst modification)
26       for useOBCSprescribe without useEXF.
27    
28    checkpoint61y
29    o pkg/seaice
30     - change (once again) the discretization for the evp solver: rather than
31       squaring the averages of strain rates, average the squares to give
32       a better approximation to VP behavior (in simple test); this breaks
33       verification experiments that use EVP (lab_sea and lab_sea.hb87)
34        --> update outputs
35    o add CPP-flag (in mom_calc_visc.F) and TAF directive (in gmredi_slope_limit)
36      to improve vectorization of adjoint code on TARGET_NEC_SX
37    o verification/vermix:
38      - change vermix to a more "classical" 1d column with high vertical
39        resolution, seasonal heat flux forcing and constant wind forcing
40        --> update outputs
41    o pkg/seaice
42     - fix diagnostics of principal stress components for EVP
43     - add optional diagnostics of non-linear convergence for LSOR solver
44     - make maximum zeta a runtime parameter
45    o model/src (3-d solver):
46      Undo modif from Dec 12, 2005 (solve_for_pressure.F cvs version v.49),
47       which consisted in removing the Surf. Non-Hydrostatic pressure from the
48       whole column and adding it to Eta (inspired from Casualli, 1999):
49       Only affects 3-D solver initial guess and only if using ExactConserv,
50       but tests did show a degradation of solver convergence. For this reason,
51       decide to disable this piece of code.
52       to recover previous version, just set zeroPsNH = exactConserv (line 78).
53      Update output of exp. global_ocean.cs32x15.viscA4 & hs94.cs-32x32x5.impIGW
54    o pkg/ggl90:
55     - fix previous modif (unfortunatly has not been documented here).
56    o pkg/seaice/autodiff:
57     - add CPP flag to include clipping of zeta in EVP code;
58       although clipping of zeta is not necessary for EVP and not recommended
59       it was done by default, now clipping is disabled by default
60     - modify computation and store zetaC/Z in seaice_evp for TAF
61     - add CPP-flag SEAICE_ALLOW_CLIPVELS (undefined by default)
62     - change seaice_calc_viscosities/strainrates for better
63       vectorization after AD by TAF
64     - rearrange common blocks in SEAICE.h once again for easier dumps of
65       adjoint output, adjust adcommon.h and g_common.h
66     - change flux form of default advection scheme for better vectorization
67       of adjoint code
68     - modify lab_sea/input.lsr to test the flux form of default advection
69       scheme
70    o autodiff: fix seaice fields in g_common.h (left over from 3 to 1 level
71      transition in checkpoint61r)
72    
73    checkpoint61x
74    o adjoint: adding two new verification experiments
75      * tutorial_global_biogeo: an online dic sensitivity experiment
76        (work with G. Forget)
77      * tutorial_dic_adjoffline: an offline dic sensitivity experiment
78        (an update of the Dutkiewicz et al., 2006 paper)
79      Both use the same CO2FLUX cost function, but different controls
80      Gradient checks look good, but some recomp. could still be traightened out
81      Should add xx_dic control variable to test at some point
82    o model/src:
83      add run-time parameters to check and/or mask initial Temp & Salt
84    o pkg/ptracers:
85      apply zonal filter (if used) to all passive tracers
86    
87    checkpoint61w
88    o adjoint:
89      * update and finish work on adjoint of dic package (G. Forget)
90      * start adjoint of offline package
91      * bug fix for cube-sphere TLM code (active_file_g.F)
92    o tutorial_global_oce_optim:
93      now tested with 4 tiles (change data.grdchk and update output_adm.txt)
94      and add SIZE.h_mpi for MPI testing
95    o pkg autodiff
96     - fix bug in autodiff_restore
97     - make genmake2 remove comments and empty lines from ad_input_code.f in
98       order to reduce file size (new script file remove_comments_sed)
99    o allows to specify vertical profile of vertical viscosity (viscArNr)
100    o pkg autodiff
101      autodiff_ini_model_io: comment out writing of grid information to netcdf
102      files, because at this time the grid fields are a bunch of zeros
103      (this is analogous to ini_mnc_vars.F)
104    o pkg seaice/autodiff:
105     left-over from removing the 3-time-level fields:
106     - small rearrangement of common blocks (SEAICE.h)
107     - fix field declaration in order to get meaningful output in
108       adseaice.*.nc (adcommon.h)
109    o pkg obcs:
110     - go through NH code only if nonHydrostatic=T
111     - fix restart for Orlanski (Hydrostatic; still to fix NH+Orlanski restart)
112     - move call to OBCS_CALC before SEAICE_MODEL which needs seaice-obc fields
113       => Update output of experiment seaice_obcs
114     - useOBCSbalance: code is wrong with multi-tile / proc: add a stop.
115     - fix NH+Orlanski restart (add wvel OBC value in pickup file)
116     - fix restart for useOBCSprescribe + useEXF
117    o pkg ctrl
118     - make obcs as control parameter work also with useSingleCPUio
119     - replace a few sny and snx by Ny and Nx to be consistent with
120       ctrl_set_globfld_x/yz.F
121    o pkg seaice
122     - re-introduce global field frWtrAtm for the case of
123       ALLOW_MEAN_SFLUX_COST_CONTRIBUTION defined (SEAICE.h, seaice_growth.F)
124     - catch the case of an empty filename in seaice_cost_init_fixed.F
125    o pkg ecco:
126     - small bug fixes in the_main_loop
127      . fix order of monitor arguments
128      . change SEAICE_MULTILEVEL to SEAICE_MULTICATEGORY
129      . move include DYNVARS.h up outside the ALLOW_AUTODIFF_TAMC block so
130        that ALLOW_ECCO_FORWARD_RUN works with NONLIN_FRSURF
131     - replace cos(yc*deg2rad) by _rA in cost_mean_heat/saltflux.F, fix the
132       imbalance diagnostics for wmean_h/sflux = 0
133     - add actual values to f_obcsn/s/w/e in ecco_cost_final.F
134     - catch the case of empty filenames in ecco_cost_init_fixed.F
135       and ecco_readparms.F
136     - fix the logic of wsalt=0.and.wsalt2=0 to wsalt=0.or.wsalt2=0 (and theta)
137       to make this statement meaningfull in ecco_cost_weights.F
138    o pkg zonal_filt:
139      - change ZONAL_FILTER S/R interface (allows to filter 2-D fields)
140      - when using rStar, filter uVel*hFacW instead of uVel
141       (tested without topography: fix a Pb of growing instability with NLFS)
142    o pkg/aim/thsice:
143     -Add capacity to read a Q-Flux to be applied in slab-mixed layer thsice_slab_ocean.F
144      (works as SST restoring)
145    o pkg/thsice:
146     -add specific parameter for salt restoring in slab mixed-layer
147       (by default equal to the temperature restoring)
148    
149    checkpoint61v
150    o pkg/kpp:
151     -add double diffusive contributions as a hack, turned off by default
152      for now, and the code can be excluded with a CPP-flag
153      EXCLUDE_KPP_DOUBLEDIFF just as EXCLUDE_KPP_SHEAR_MIX
154     -add a verification experiments to "vermix" to test this code
155    o pkg/mom_common:
156      Move computation of length scales to mom_init_fixed.F (avoid recomputation
157      at each time-steps*levels).
158    o pkg/layers:
159      Created a new package for computing volume fluxes in temperature or
160      [not yet implemented] density layers, e.g. for calculating residual
161      overturning circulations.
162    o pkg/diagnostics:
163      implement a RS type version for the set of DIAGNOSTICS-FILL subroutines
164    o fix few RS/RL type mismatch (check with ifort 11.1 "-warn all" option)
165    o pkg/mdsio (& pkg/rw): rework high level S/R interface:
166      To avoid mixing type (RS/RL) of input/output array argument,
167      replace single mixed-type array with a pair from each type (RS/RL).
168    o setting pickupSuff force to start from a pickup (even if nIter0=0): this
169      is now implemented for all pkgs (previously only done for the main pickup)
170    
171    checkpoint61u
172    o pkg/obcs:
173      - change OB?eta from _RS to _RL and put them into a separate commen block
174      - add code to read OB?eta from a file via obcs_prescribe_read: the code
175        compiles and does not destroy any other experiments, the data is read
176        properly, but the code is not tested in actual applications
177    o ini_mnc_vars.F: if available (ALLOW_CAL), use startTime_1 and 2 to
178      create a meaningful unit attribute for the time coordinate variable
179    o pkg/ebm: add plenty of missing "_d 0" ;
180       update results of experiment global_ocean_ebm (fwd + adjoint tests)
181      Notes: several pieces of code only works for (very) specific set up
182       (domain size, lat range, etc ...); should be documented & stop if one
183       of those assumptions is not meet.
184    o pkg/seaice: add missing "_d 0" in HSNOW initialisation
185       update results of experiments: lab_sea (3 fwd tests + 2 adjoint tests)
186          and offline_exf_seaice (seaicetd + adjoint).
187    o add missing "_d 0" in ShortWave flux penetration (S/R SWFRAC):
188       update results of experiments: lab_sea (4 fwd tests + 3 adjoint),
189        seaice_obcs, global_ocean.cs32x15 (icedyn & thsice) and natl_box (x2).
190    o few changes to pass when compiling with strick checking of arguments
191       across S/R (allow to find few bugs in S/R calls).
192    o pkg/seaice:
193     - seaice_evp: move if-statement out of loop (TEM) at the cost of 11 new
194       2d-fields, for slightly better performance
195    
196    checkpoint61t
197    o pkg/seaice:
198     - remove unused variables in seaice_evp
199     - reduce number of SQRT in seaice_evp and seaice_calc_viscosities
200    o eeset_parms.F, open_copy_data_file.F & ini_parms.F:
201       named file replacing SCRATCH files (for TARGET_BGL & TARGET_CRAYXT):
202       each processor opens a different file (a fix from Matt).
203    o cg2d:
204       use simple EXCH (overlap size = 1 and ignore corners), like in cg3d.F :
205       this reduces number of EXCH calls by 2 when using exch2.
206       For now, keep previous size array for gc2d_nsa.F to avoid Pb with store
207       directives and missing simple EXCH (exch_s3d_rx) adjoint code.
208    o pkg/thsice:
209       add code for Near Infra-Red albedo (from Jeff)
210       (turned on with run-time flag: thSIce_calc_albNIR, in data.ice).
211    o ini_model_io:
212       move (from packages_readparms.F & initialise_fixed.F) to ini_model_io.F
213       MNC & MONITOR initialisation calls; move MNC calls from {pkg}_readparms.F
214       to {pkg}_init_fixed.F for {pkg}= kpp, seaice & shelfice.
215    
216    checkpoint61s
217    o pkg/exch2 (exch2_rxN_cube*):
218       remove last BARRIER (no need to synchronise after getting data from shared
219       buffer (get) as long as any change to buffer (put,recv) is between BARRIER)
220    o change tiling of natl_box (from 2 to 4 tiles) to get it tested with MPI+MTH
221      and generate new output (for the 2 tests) using default rotationPeriod.
222    o pkg/exch2:
223      - always call exch2_*_cube, not exch-1 anymore, if useCubedSphereExchange=F
224        (was already the case with AUTODIFF_EXCH2 defined)
225      - add bj in exch2 arrays and S/R.
226      - can now use pkg/exch2 for regular (non CS) set-up; switch test experiment
227        MLAdjust to exch2 (identical results except with MPI -> truncation diff);
228    o add pkg/longstep for ptracer timestep a multiple of U/V/T/S timestep
229      - adds deltaT argument to the subroutines TIMESTEP_TRACER DWNSLP_APPLY
230        GAD_ADVECTION GAD_CALC_RHS GAD_IMPLICIT_R GAD_SOM_ADVECT
231      - verification experiment natl_box.longstep
232    o build_options: ifort+mpi_aces optfile (following Constantinos advise):
233      allow different versions of ifort (8 or 9) and different versions
234       of mpi (mpich or mpich2) to be used with the same optfile.
235      requires to load a netcdf module and set MPI_INC_DIR (see comments inside
236       optfile). Change testing scripts (example_scripts/ACESgrid) accordingly.
237    o pkg/seaice:
238      fix previous modif (3-level field removed) for old-pickup reading
239    
240    checkpoint61r
241    o pkg/seaice:
242     - clean up the 3-time levels of UICE,VICE,HEFF,AREA in three steps:
243      1. introduce fields for time step (N-1)
244      2. do not use levels 2 and 3 anymore (but keep them)
245      3. replace 3D versions of UICE,VICE,HEFF,AREA by 2D versions. This
246         includes changes to pkg/thsice/thsice_get_velocity.F,
247         model/src/external_forcing_surf.F, pkg/obcs/obcs_apply_seaice.F, and
248         pkg/obcs/obcs_apply_uvice.F.
249     - completely remove code for leap-frog time stepping for second order
250       advection
251     - adjusted adjoint accordingly;
252       gradient checks have improved by an order of magnitude!
253    
254    checkpoint61q
255    o pkg/seaice: fix bug in free slip boundary conditons for LSOR solver
256     - the boundary conditions along open boundaries are unclear, because the
257       associated masks are set to zero and are thus treated as closed when
258       applying the slip-boundary conditions.
259     - move initialisation of grid variables from seaice_init_varia to _fixed
260    o adding doc/diags_changes.txt
261    o model/src:
262      do_the_model_io: rename few pkg S/R which are called from do_the_model_io,
263       to avoid "_diags" naming if no relation to pkg/diagnostics.
264      ini_theta,ini_salt: fix (at least compile the code) reading from NetCDF
265       file and only call PLOT_FIELD_XYZ if debugMode=T.
266      remove few PLOT_FIELD_XY calls outside initialisation.
267      remove few unnecessary BARRIER.
268    o pkg/salt_plume: add SaltPlumeSouthernOcean to allow disabling
269      of salt_plume package in Southern Ocean
270    o pkg/thsice:
271     - fix hOceMxL in overlap (needed for sea-ice advection);
272       update output of experiment global_ocean.cs32x15.icedyn.
273    o pkg/seaice:
274     - make the (placeholder-)tracer IceAge a little more meaningful
275      by avoiding ridging of age due to convergent ice flow
276     - add diagnostics for surface temperature, atmospheric and oceanic
277       heat flux, and atmospheric and oceanic ice growth
278    
279    checkpoint61p
280    o eesupp (global_max & global_sum):
281     - add one more element to buffer (start at index 0) for GL_MAX/SUM output;
282     - remove starting & endding barrier (no longer needed)
283    o pkg/mdsio, writelocal:
284       works also in multi-threaded when LOCBIN_IO_THREAD_SAFE is undefined
285       (remove the STOP): uses shared buffer IO to store data from all threads
286       and then let master-thread write nThreads tiles.
287    o pkg/mdsio (read/write field):
288     - tiled IO done in 1 piece (all levels at a time);
289     - new header file "MDSIO_BUFF_3D.h" with 3-D buffers.
290     - multi-threaded: allow to read/write local (non-shared) array
291       (was already working with singleCpuIO ; now works also without);
292       => no longer needs barrier call after reading a file (added inside
293       mds_read/write field).
294     - move barrier calls outside gather/scatter_2d to mds_read/write_field
295    o pkg/obcs/exf
296      - exf_getffieldrec.F: add a few ". _d 0" and
297                            add code to catch the case fldperiod = 0.  
298      - obcs_precribe_read.F: change if statements so that they match
299                              the variable type (_RL)
300      - add a separate exf_iprec_obcs to exf_set_obcs (by default equal to
301        exf_iprec) as part of namelist EXF_NML_01
302    o pkg/mdsio:
303      - read/write tiled (local) files: read/write 1-level tile chunk at a time
304        (instead of segment of length sNx); expected to speed up tiled IO.
305    o pkg/exch2:
306      - take buffer copy from/to array out of S/R exch2_send/recv into new
307        S/R exch2_put/get ; Exch of local variable now works with MPI+MTH
308        (tested by removing commom block statement in SOLVE_FOR_PRESSURE.h).
309      - ad version of send/recv no longer needed (but needs exch2_ad_put/ad_get)
310      - simplify argument list of S/R exch2_get_uv_bounds and exch2_get_scal_bounds
311        (which replaces exch2_get_recv_bounds & exch2_get_send_bounds) by
312        including Topology header file (and common blocks).
313      - implement EXCH_IGNORE_CORNERS in scalar exchange (rx1);
314        remove 2nd exch call in exch2_s3d_rx (no longer needed).
315      - document S/R arg. list (+ update some comments)
316    o pkg/seaice:
317      - added SEAICE_CLIM_AIR code, which permits modulation of surface air
318        temperature and humidity over sea ice based on climatological values.
319      - clean up computation of Hibler+Bryan (1987) stress coupling for the case
320        of LSR (this change is expected to change the results slightly because now
321        it uses slightly difference moduli from the second last LSR solution, but
322        that is more consistent with the stress computations; this part of the code
323        is not tested in the verification experiments)
324      - fix diagnostics for flux of ice in case of B-grid
325      - turn strain rates eij and press into global fields within global
326        common block in SEAICE.h for more straighforward diagnostics
327      - fix diagnostics for principal stress components, press, zeta, eta
328        to give something meaning full
329    o change tiling of global_ocean.cs32x15 (12 tiles, 32x16) and generate
330      new output (all 4 tests) using default rotationPeriod.
331    o pkg KPP:
332      - add flag (LimitHblStable) to avoid limiting BL depth under stable conditions
333      - fix vertical indices in RI_IWMIX (match between KPP grid from 0 to Nr+1
334        and normal grid from 1 to Nr)
335  o Fix computations of L4rdt in mom_calc_visc.F  o Fix computations of L4rdt in mom_calc_visc.F
336  o eesupp:  o eesupp:
337    new S/R ALL_PROC_DIE (like EEDIE, without irrelevant threads err msg):    new S/R ALL_PROC_DIE (like EEDIE, without irrelevant threads err msg):

Legend:
Removed from v.1.1108  
changed lines
  Added in v.1.1199

  ViewVC Help
Powered by ViewVC 1.1.22