/[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.1724 by jmc, Sat Aug 25 00:07:09 2012 UTC revision 1.1791 by jmc, Mon Dec 24 15:55:10 2012 UTC
# Line 4  $Name$ Line 4  $Name$
4      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
5      ==============================      ==============================
6    
7    checkpoint64b (2012/12/24)
8    o pkg/cheapaml:
9      - remove un-used variables (Qnet,EmPmR) from pickup files.
10      - store tendency of advected variables into cheapaml pickup file
11        and fix the AB-2 restart.
12    o verification/offline_exf_seaice and lab_sea adjoint:
13      - switch to new code (#undef SEAICE_GROWTH_LEGACY).
14      - switch to 4 time steps in offline, consistent with lab_sea.
15      - match parameters from SEAICE_GROWTH_LEGACY branch in data.seaice.
16      - update results after change of experiment seaice.
17    o pkg/seaice:
18      - remove option SEAICE_MULTICATEGORY; has been replaced with run-time
19        parameter SEAICE_multDim. Issue a warning if SEAICE_MULTICATEGORY is
20        still defined since SEAICE_multDim needs to be set to MULTDIM value (from
21        SEAICE_SIZE.h) to recover previous #define SEAICE_MULTICATEGORY default.
22    o model/src:
23      - add run-time parameter to select CG2D solver minimum-residual solution
24        (in case of a poor convergence).
25    o pkg/seaice/seaice_lsr.F
26      - move calculation of matrix coefficients into new
27        S/R SEAICE_CALC_LSR_COEFFS; for now this routine remains in the same
28        file; the main purpose of this change is to be able to reuse the code
29        in S/R SEAICE_PRECONDITIONER
30      - modify the computation of the residual so that the mass of sea-ice
31        is not used twice; this simplifies the calculations, but since
32        the residual is never used except when LSR_mixIniGuess.ge.4---and this is
33        not tested in testreport---this does not affect the testreport results.
34    o verification/seaice_obcs:
35      - disable SEAICE_GROWTH_LEGACY; reset input files to match those of
36        verification/lab_sea/input.salt_plume; changes results/output.txt*
37    o model/src:
38      - adapt quasi-hydrostatic formulation to work with p or p* coordinate.
39    o pkg/seaice:
40      - improve description of parameters in seaice_fgmres.F and start
41        cleaning up the code
42      - change name of SCALPROD to SEAICE_SCALPROD and modify meaning of
43        second argument
44      - add a metric based on grid cell area to SEAICE_SCALPROD; update
45        results of offline_exf_seaice.dyn_jfnk
46      - compute initial residual norm in seaice_jfnk based on SEAICE_SCALPROD
47    o pkg/aim_v23:
48      - fix an other bug in dissipation heating (grid-cell area was missing).
49    o model/src + pkg/seaice:
50      - add global (multi-tile) linear solver for pair of tri-diagonal system
51        along X and Y lines, respectively for U and V component.
52        Note: 1) MPI and cube-exchange not yet coded.
53            2) probably not accurate for poorly conditioned / large size problem.
54      - add option for SEAICE_LSR to use global (multi-tile) tri-diagonal solver;
55        default: not compiled and not used.
56    o pkg/seaice:
57      - first attempt at a multi-threaded fgmres solver for the JFNK solver:
58        introduce bi/bj loops and use global_sum_tile_rl instead of
59        mpi_allreduce in seaice_fgmres.F
60      - improve diagnostics output for fgmres solver (seaice_fgmres.F)
61      - replace a few unnecessary goto statements by if-blocks  (seaice_fgmres.F)
62      - add a verification experiment for JFNK in offline_exf_seaice
63    o verification:
64      - change (for now, only fwd test) offline_exf_seaice experiment:
65        use a re-entrant channel with idealized coastline (trapezoidal shape)
66        at southern boundary to explore seaice dynamics and thermodynamics
67        with simple forcing. Main set-up has dynamics+thermo and SST relaxation;
68        secondary test (input.thsice) use thermodynamics only with SST relaxation.
69      - add seaice-dynamics only secondary test using LSR and pkg/thsice advection.
70    o pkg/exf:
71      - disable "sstExtrapol" when Nr=1
72    o verification:
73      - add new experiment "cheapAML_box" (simple 1 layer set-up from Bruno Deremble)
74        to test pkg/cheapaml.
75    o pkg/aim_v23 & pkg/shap_filt:
76      - fix dissipation heating (erroneous gravity multiplication factor).
77    o pkg/seaice:
78      - remove unused field pressPre from jfnk solver
79      - rename two subroutines in seaice_preconditioner.F
80    o eesupp/src:
81       - S/R MDSFINDUNIT avoids reserved unit numbers (scrUnit1, eeDataUnit ...)
82    o pkg/seaice/seaice_preconditioner.F: fix preconditioner once again
83       - remove wrong parts of the rhs and the lhs (were in there
84         because I did not understand the concept of this preconditioner)
85       - add an exchange between u- and v-sweep. This increases communication
86         by a factor of two but improves the convergence of jfnk. Ideally,
87         one would like the u- and v-sweep to be independent. So far within
88         undocumented ifdefs
89       - reorganize the extra computations for printing the residuals
90    o tools/genmake2:
91      - comment out the "mpif.h" include in NetCDF test programm (for MPI built)
92        since it's not necessary (and was preventing to use MNC in some cases).
93    
94    checkpoint64a (2012/11/16)
95    o pkg/autodiff:
96      - improve pkg-internal logical parameter switch (inAd/inFwd mode) so that
97        autodiff_readparms.F no longer depends on calling order (within
98        packages_readparms.F) and no longer requires to include ${PKG}_PARMS.h):
99        rename param turnFreeDriftInAdMode to SEAICEuseFREEDRIFTswitchInAd
100        and change param SEAICEuseDYNAMICSinAdMode to SEAICEuseDYNAMICSswitchInAd.
101    o pkg/obcs:
102      - adding tidal velocity forcing capability to obcs
103        for usage example, see verification/seaice_obcs/input.tides
104    o pkg/seaice:
105      - turn the LSR-relaxation parameters WFAU/V into runtime parameters
106        SEAICE_LSRrelaxU/V
107      - introduce runtime parameter SEAICEetaZmethod (=0 by default)
108    o pkg/aim_v23:
109      - implement heating from friction (bottom friction and stratospheric drag).
110    o pkg/shap_filt:
111      - fix diagnostic "SHAP_dKE" (was 2 time smaller); implement heating from
112        dissipation but only for Shap_funct=2, active when addFrictionHeating=T.
113    o model/src:
114      - move AddMass common block from DYNVARS.h to FFIELDS.h; allow to read-in
115        time-constant addMass field from file (addMassFile);
116      - support for heating due to friction (and Kinetic Energy dissipation):
117        add option, run-time switch and 3-D field array + support code ; each pkg
118        or S/R can add heating contribution (if there is a sink of KE); the sum is
119        accounted for as a Temperature tendency term (EXTERNAL_FORCING_T).
120    o pkg/seaice:
121      - complete merge of ITD code into seaice_growth.F
122        from MITgcm_contrib/torge/itd/code/
123        no verification so far, but Torge Martin tested
124        "1-D test case and Beaufort set up;
125        deliver okay results after 1 model year."
126      - complete merge of SEAICE_SIZE.h into all relevant routines
127      - move computation of eta at Z-points into calc_viscosities
128        + introduce global field etaZ
129        + add different methods for computing etaZ, for now the old
130          method (simple averages of eta to Z-points) is hardwired
131        + replace etaMeanZ by etaZ in seaice_lsr, seaice_calc_residual,
132          seaice_preconditioner
133      - simplify computations for SEAICE_ZETA_SMOOTHREG
134      - seaice_preconditioner:
135        + replace one remaining press by pressLoc (mistake)
136        + adjust printing of residuals
137      - improve I/O of diagnostics of JFNK in seaice_jfnk.F:
138        + add (SEAICE.h) and rename some counters (initiliased in
139          seaice_init_fixed.F)
140        + add _BEGIN/END_MASTER (myThid)  around print statements
141        + print total number of iterations and failures over SEAICE_monFreq
142          interval to STDOUT
143      - move control over output-messages by fgmres to seaice_jfnk to save
144        some operations
145      - move resetting initial guess for fgmres outside of seaice_fgmres_driver
146        in order to make restarts with sol .ne. zero (fgmres with restarts)
147        (seaice_jfnk.F, seaice_fgmres.F)
148      - seaice_calc_l/rhs.F: replace maskW/S by maskInW/S; this fixes
149        a serious convergence problem when open boundaries are involved
150    o pkg/ggl90: fix computation of control volumes for the case of
151      ALLOW_GGL90_HORIZDIFF (probably not used very often): move control volumes
152      up half a grid cell from C-centered to W-centered (pointed out by C. Eden)
153    o pkg/seaice: simplify seaice_evp by avoiding use of extra fields u/vIceC
154    o pkg/autodiff
155      - Add inAdMode switch for SEAICEuseDYNAMICS:
156        SEAICEuseDYNAMICSinAdMode
157        (extend logic in ADAUTODIFF_INADMODE_SET)
158    o pkg/seaice
159      - Step 1 of merging ice-thickness distribution (ITD) code from
160        MITgcm_contrib/torge/itd/ (author: Torge Martin)
161        The code implements concept laid out in Thorndike et al. (1975);
162        code is enabled via CPP flag SEAICE_ITD
163        nITD (in SEAICE_SIZE.h) defines number of ice thickness categories;
164        code is not complete yet!!
165      - SEAICE_SIZE.h now included in almost all seaice routines
166    o verification (using seaice):
167      - switch back seaice-dynamics (+ lower SOLV_MAX_ITERS) in ADM exp.
168        global_ocean.cs32x15.seaice ; update results.
169      - compile B-grid dynamics in FW test exp 1D_ocean_ice_column ;
170      - compile C-grid dynamics with #define SEAICE_ALLOW_JFNK in FW test exp
171        offline_exf_seaice (previously compiling B-grid).
172    o pkg/exf:
173      - fix some STOP conditions in exf_check.F (from Aug 28, modifs did not allow
174        to specify neither just 1 comp of wind or wind-stress, nor no surface momentum
175        forcing): just replace CPP option ALLOW_ATM_WIND with run-time param useAtmWind.
176    o pkg/layers (adjust changes from Sept 19):
177      - do snap-shot output and fill diagnostics from inside iLa loop (layers_calc.F)
178        this allows to remove "layers_maxNum" dimension from all layers arrays.
179    o pkg/layers (adjust changes from Sept 19):
180      - remove ALLOW_LAYERS_OUTPUT ; fix snap-shot and time-ave output ;
181      - snap-shot extended to all type of layers (not just the 1rst type);
182      - remove layers variables when corresponding CPP option is undef.
183    o pkg/layers:
184      - rename pkg/layers variable and add new output (layer exist & un-weighted
185        velocity in layer)
186    o pkg/seaice:
187      - simplify (?) ocean-drag coefficent computations by using
188        the same routine in all three solvers (saves some code lines)
189      - remove global fields DRAGA and DRAGS from C-grid code and
190        declare them locally in lsr.F (the only place where they are
191        still used)
192      - update AD output of lab_sea.evp (very sensitive anyway).
193    o pkg/seaice: add Jacobian-free Newton Krylov solver for ice dynamics
194      (Lemieux et al. 2010, 2012). Turn on with the usual SEAICE_ALLOW_JFNK
195      and SEAICEuseJFNK=.true. A few things to note:
196      - involves a FGMRES, that is parallel but very inefficient,
197        because of many scalar products (global sums) in the Gram-Schmidt
198        orthogonalization method. Needs to be replaced. Original code
199        from J.-F. Lemieux
200      - preconditioner for FGMRES is a stripped down version of LSOR
201      - SOLV_MAX_ITERS is now used for LSOR as a solver and as a
202        preconditioner. For the preconditioner, the default is set to 10.
203      - Jacobian times vector is computed with finite differences
204      - marginally related to JFNK-solver: add a smooth regularization for
205        zeta with tanh that will improve convergence (Lemieux and Tremblay,
206        2009) for any implicit solver (not implemented for EVP because EVP
207        does not need this regularization in the first place). Define
208        SEAICE_ZETA_SMOOTHREG in SEAICE_OPTIONS.h to use.
209    o eesupp/src:
210      - check in genmake2 that intrinsic S/R flush is available.
211      - add S/R to flush IO unit (if intrinsic S/R flush is available)
212      - flush standard & error msg unit before MPI termination
213      - if debugMode, flush IO unit after printing msg
214    o pkg/ctrl:
215      - ctrl_unpack.F, ALLOW_GENARR3D_CONTROL code: fix missing argument in call
216        to ctrl_set_unpack_xyz.
217    
218    checkpoint64 (2012/10/09)
219    o pkg/shelfice
220      - in mom_fluxform.F, mom_vecinv.F:
221        CALL SHELFICE_V_DRAG(bi,bj,k,vFld,KE,KappaRU,vF,myThid)
222        Bug fix KappaRU -> KappaRV
223        does not change verification/isomip but may change other setups
224        (presumably constant KappaRU=KappaRV)
225      - fix diagnostics SHIU/Vdrag (thanks to Tom Millgate)
226    o pkg/obcs:
227      - fixed some indexing bugs in obcs_seaice_sponge.F
228      - added obcs relaxation subroutines for HSALT and HSNOW
229      - changes verification/seaice_obcs/results/obcs_seaice_sponge.F
230    o pkg/seaice:
231      - add "multicategory" parameterization for snow in analogy to ice,
232        default is off for now.
233    o pkg/obcs:
234      - added obcs_seaice_sponge.F for obcs relaxation of seaice variables;
235        see verification/seaice_obcs/input.seaiceSponge for example application.
236      - remove local commom bloc "OBCS_GLOBAL" in obcs_readparms.F and do all the OB
237        index setting by master-thread. This fix multi-threads Pb in previous modifs.
238    o pkg/layers :
239      - add diagnostics to pkg/layers
240      - move the meat of pkg/layers from layers_calc.F to layers_fluxcalc.F
241        (from R.Tulloch). Layers_calc.F now calls layers_fluxcalc.F accordingly.
242    o pkg/generic_advdiff/gad_implicit_r.F :
243      - fix ADVr diagnostics in case of both implicit vertical advection and diffusion.
244    o pkg/seaice/seaice_growth.F:
245      - fix SItflux diagnostics in virtual salt flux + balanceEmpmr case.
246    o pkg/obcs:
247      - change default setting of OBCS_UVICE_OLD from defined to undef;
248      - use new parameter OB_indexNone for null index value (instead of hard-coded 0)
249        in all pkg/obcs code and in other pkgs; remove k loop in obcs_check_depths.F;
250      - pass OB_indexNone as (new) argument to S/R CTRL_MASK_SET_[X,Y]Z;
251      - change OB null index value (OB_indexNone) to -99 ; this fix the problem
252        with Southern & Western OB at a tile edge (at least in a simple set-up).
253    o pkg/obcs:
254      - change settings of OB index in case OB coincides with tile edge (shared by
255        2 tiles), so that common OB location is registered by both tiles (instead of
256        only 1). However, for Southern and Western OB, index value = 0 is not
257        recognized as valid OB index (=default null value) and needs to be fixed.
258      - modify accordingly OBCS_BALANCE_FLOW and OBCS_MONITOR (avoid double counting)
259      - new parameter "OB_indexNone" for null index value (for now, still zero)
260        to replace hard-coded null index value of zero.
261    o pkg/exch2:
262      - improve W2_USE_E2_SAFEMODE (faster) and W2_E2_DEBUG_ON:
263        always check buffer-size (but outside the loop); move checking for valid
264        index from W2_USE_E2_SAFEMODE to W2_E2_DEBUG_ON; in W2_E2_DEBUG_ON,
265        print each tile and point connection only if |W2_printMsg|>= 2 and 3.
266    o pkg/ctrl (pack & unpack):
267      - in ctrl_pack.F & ctrl_unpack.F, use the same local file-name variable
268        (locally defined) for different variables (+ remove the variable-specific
269        file name in CTRL_GENARR.h). For now, only implemented for GENARR2/3D_CONTROL
270        (to be extended to other local file names in ctrl.h).
271    o pkg/flt:
272      - add a FLT_Iter0 option (similar to PTRACERS_Iter0).
273    o pkg/mdsio:
274      - fix missing STOP in mdsio_readvec_loc.F (case ioUnit=0) when file is missing
275    o tools/genmake2:
276      - add similar dependences for target EXE_AD & EXE_FTL as we have for EXECUTABLE
277        (added by Alistair in version 1.51 for building without make depend);
278      - remove Header files list (H_SRC_FILES) from the dependence of targets
279         ad_input_code.f & ftl_input_code.f (not clear why it was added in the first
280        place in version 1.23): this save going through staf when unrelated header
281        file has been modified.
282       Note: for an intermediate target (e.g., ad_input_code.f or adtaf) without
283        "make depend" step, still need to do "make links" before "make adtaf".
284    o eesupp/src (+ few pkg src files) & genmake2:
285      - finish to remove ALWAYS_USE_MPI in source code that TAF does not see;
286      - no longer set ALWAYS_USE_MPI in Makefile (was done by genmake2).
287      - improve genmake2 checking for ALLOW_${PKG} in CPP_OPTIONS.h files (distinguish
288         upper/lower case ; account for space between "#" and undef/define)
289    o pkg/flt (with exch2):
290      - cleaning: fix so that it compiles; move remaining MPI calls out of flt_exch2.F;
291      - simplify threads synchronisation; fix (many issues) multi-threaded with MPI.
292    
293    checkpoint63s (2012/09/05)
294    o pkg/ecco:
295      - remove #ifdef ALLOW_SMOOTH_CORREL3D brackets.
296      - add more relevant #ifdef ALLOW_SMOOTH ones.
297      - sort out useAtmWind, useSMOOTH, ctrlSmoothCorrel2D.
298    o pkg/ctrl:
299      - CTRL_OPTIONS.h : document CPPs related to smoothing; ALLOW_SMOOTH_CORREL2D/3D,
300        ALLOW_CTRL_SMOOTH, ALLOW_SMOOTH_CTRL2D, ALLOW_SMOOTH_CTRL3D.
301      - ctrl.h : document run time params doInitXX, doMainPack, doMainUnpack, doPackDiag,
302        doSinglePrecTapelev (audodiff param!), ctrlSmoothCorrel2D (new) and
303        ctrlSmoothCorrel3D (new).
304        The added ctrlSmoothCorrel2D/3D replace the old ALLOW_SMOOTH_CORREL2D/3D CPP.
305      - ctrl_readparms.F : initialise ctrlSmoothCorrel2D/3D and read from namelist.
306      - ctrl_init.F ctrl_pack.F ctrl_unpack.F : switch off ustress/vstress
307        controls when useAtmWind, and vice versa.
308      - ctrl_map_ini_ecco.F : UVEL0, VVEL0, edtaux, edtauy cannot be smoothed as tracers.
309      - ctrl_get_gen.F : replace #ifdef ALLOW_SMOOTH_CORREL2D with
310          #ifdef ALLOW_SMOOTH, IF ctrlSmoothCorrel2D, IF useSMOOTH.
311      - ctrl_map_ini_gen.F : replace #ifdef ALLOW_SMOOTH_CORREL3D with
312          #ifdef ALLOW_SMOOTH, IF ctrlSmoothCorrel3D, IF useSMOOTH.
313      - ctrl_set_pack_xy[z].F, ctrl_set_unpack_xy[z].F:
314        replace ALLOW_SMOOTH_CORREL2D/3D brackets with ctrlSmoothCorrel2D/3D.
315    o pkg/smooth:
316      - remove un-necessary includes.
317      - remove ALLOW_SMOOTH*D and ALLOW_SMOOTH_CORREL*D brakets.
318        Those CPP options were never defined, and not necessary.
319    o model/src:
320      - fix CPP brackets in ini_mixing.F (GM_background_K belongs to gmredi).
321      - standard treatment of useSMOOTH.
322    o pkg/autodiff:
323      - add LOGICAL useGGL90inAdMode, useGGL90inFwdMode, useSALT_PLUMEinAdMode,
324        useSALT_PLUMEInFwdMode allowing to switch off ggl90 and salt_plume adjoint.
325    o pkg/exch2:
326      - add "if usingMPI" where it was missing (in exch2_rx?_cube*.template)
327        + start to remove ALWAYS_USE_MPI in source code that TAF does not see.
328      - make new S/R W2_MAP_PROCS out of w2_eeboot.F to setup exch2 CommonFlag
329        and tiles to process mapping.
330      - re-arrange 3 header files (W2_EXCH2_PARAMS.h, _BUFFER.h & _TOPOLOGY.h)
331        moving /EXCH2_FILLVAL_RX/ & /W2_EXCH2_COMMFLAG/ to W2_EXCH2_TOPOLOGY.h;
332        rename "exch2_tProc" to "W2_tileProc" and move seting of e2FillValue_RX
333        from w2_eeboot.F to w2_readparms.F.
334    o eesupp/src:
335      - set exchNeedsMemsync & exchUsesBarrier according to "usingMPI" (instead of
336        ALWAYS_USE_MPI) + moved from ini_communication_patterns.F to eeset_parms.F;
337        also move setting of exchCollectStatistics from exch_init.F to eeset_parms.F
338      - in eedie.F and all_proc_die.F, always call MPI_FINALIZE when ALLOW_USE_MPI
339        is defined (remove condition on usingMPI) to match cases where MPI_INIT
340        has been called.
341      - fix bug I put in previous modif of exch_sm_3d_rx.template (not much used);
342      - change default (non-MPI) value of myPid (from 1 to 0) to match MPI settings
343        (so that on 1 proc, get same myPid independently of usingMPI)
344      - add a stop if multi-procs but usingMPI=F ; stop also if usingMPI=T and
345        useCubedSphereExchange with EXCH-1 (safer).
346      - add "if usingMPI" where it was missing (+ start to remove ALWAYS_USE_MPI).
347    o pkg/ggl90:
348      - fix bug in filling corner mask array "mskCor" (FILL_CS_CORNER_TR_RL call)
349        used in ALLOW_GGL90_SMOOTH code.
350    o pkg/flt:
351      - fix typo in FLT_WITHOUT_Y_PERIODICITY (the effect was to ignore stopping
352        floats at the Southern end of the domain).
353    o eesupp/inc:
354      - add a list of CPP-options which can be set by Makefile (DEF_IN_MAKEFILE.h);
355        but not used (and not to be used) anywhere, just for documentation.
356    o pkg/autodiff, exf, ptracers & seaice (AD & TL output):
357      - comment out the very nasty "writeBinaryPrec = writeStatePrec" line, since:
358        a) can damage previously written records in any other multi-record files.
359        b) if needed, can easily do it cleanly, by calling S/R WRITE_REC_3D_RL/S.
360    o tools/genmake2:
361      - define variable MAKE within Makefile and use it in place of 'make';
362      - update AD_CONFIG.h (if needed) before building any of the executable
363        (forward, adjoint or tangent-lin); this involves a 2nd invocation of make;
364        and add the usual #ifndef/#define AD_CONFIG_H in AD_CONFIG.h;
365      - always print usual warning comment 'this file is automatically generated'
366        from script convert_cpp_cmd2defines (+ remove this warning from genmake2).
367    o pkg/cheapaml:
368      - updated to correspond to the code used in accepted paper (Nico Wienders).
369    o pkg/exf:
370        added run time switch useAtmWind to replace ALLOW_ATM_WIND cpp switch.
371        ALLOW_ATM_WIND now just sets the useAtmWind default (see exf_readparms.F)
372        and force defines ALLOW_BULKFORMULAE (EXF_OPTIONS.h).
373    o pkg/exf, autodiff, ctrl, ecco and seaice:
374        remove ALLOW_ATM_WIND brackets, or replace them with useAtmWind ones.
375    o pkg/ctrl, ecco:
376        allow to compile both ALLOW_U/VSTRESS_CONTROL and ALLOW_U/VWIND_CONTROL.
377        Depending on useAtmWind, one is inactive, and the other is active (see
378        exf_getffields.F/exf_getsurfacefluxes.F).
379  o tools/genmake2:  o tools/genmake2:
380    - improve error msg when failing in 2nd test of check_for_broken_Ff    - improve error msg when failing in 2nd test of check_for_broken_Ff
381    - fix option '-makefile' short-name (-mf according to genmake2 -help)    - fix option '-makefile' short-name (-mf according to genmake2 -help)

Legend:
Removed from v.1.1724  
changed lines
  Added in v.1.1791

  ViewVC Help
Powered by ViewVC 1.1.22