/[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.2121 by jmc, Mon Feb 2 22:14:35 2015 UTC revision 1.2183 by mlosch, Thu Sep 24 08:17:32 2015 UTC
# Line 4  $Name$ Line 4  $Name$
4      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
5      ==============================      ==============================
6    
7    o pkg/seaice/seaice_evp.F: introduce "adaptive EVP" by Madlen Kimmritz
8      - this is an intermediate state of a truely converging EVP that has
9        the potential of outperforming implicit VP methods, I do not recommend
10        using it yet.
11      - for now, "adaptive EVP" is turned on by setting SEAICE_EVP_USE_ADAPTATION
12        and all parameters are hardwired
13      - for convergence it still requires very smooth regularisations, that means,
14        if SEAICE_DELTA_SMOOTHREG is defined, deltaCreg=deltaC+deltaMin, which is
15        different from the VP case with deltaCreg=sqrt(deltaC**2+deltaMin**2); also
16        specific averaging appears to be required, still subject to tuning
17      - still needs cleaning and I need to turn parameter into runtime parameters
18      - add code to compute and print residuals of iteration, compile with
19        defining ALLOW_SEAICE_EVP_RESIDUAL to enable it
20    o pkg/seaice/seaice_evp.F: fix a bug (factor of four in front of e12Csq)
21      that entered with r1.54 on Sep04, 2015, -> update results
22    
23    checkpoint65o (2015/09/14)
24    o pkg/mom_vecinv & mom_common:
25      - add deep-atmosphere geometry factor in mom_calc_visc.F + in routines
26        called from there.
27      - add anelastic and deep-atmosphere geometry factor in pkg/mom_vecinv ; this
28        allows to use Vector-Invariant form in deep atmos and anelastic formulation
29    o pkg/seaice: introduce fast ice parameterisation following Itkin et al. (2015)
30      and Koenig-Beatty+Holland (2012)
31      - two new parameters SEAICE_tensilFac and SEAICE_tensilDepth
32      - global 2D field tensileStrength, computed in seaice_calc_ice_strength
33      - requires an extra input argument for S/R seaice_calc_viscosities
34      - tensileStrength's depth dependence is different from Itkin et al. (2015)
35        (to be changed and played with, for now exponential)
36    o pkg/seaice: introduce new parameter SEAICEscaleSurfStress (default = .FALSE.)
37      - if TRUE scale ice-ocean and ice-atmosphere stress acting on sea ice
38        by AREA according to Connelly et al. (2004), JPO.
39      - For EVP and the LSR solver, the implicit part of the drag term is
40        the only non-zero term in the denominator when there is no ice (which
41        prevented running the model with zero ice-ocean drag). If
42        SEAICEscaleSurfStress = .TRUE., this results in a division by zero
43        (or zero main diagonals BU/BV) which need to be caught. A practical
44        consequence is that for open water, the momentum equation reduce
45        to trivially 0 = 0 + 0 (for EVP). For LSR, BU/BV are reset to 1, if
46        they happen to be zero, often leading to a non-zero solution over
47        open water.
48    o pkg/seaice: revise EVP code (seaice_evp.F) in preparation for a more
49        efficient EVP method.
50      - introduce area-weighted averages for e12 (analogously to
51        seaice_calc_viscosities.F); this totally changes lab_sea.hb87 (2 digits
52        remain, but also changing the order of terms in the previous average
53        changed the results as much) -> update experiment
54      - add SEAICE_DELTA_SMOOTHREG code
55      - disentangle computation of zetaC/deltaC and zetaZ/deltaZ, also for
56        TEM-code
57      - adjust loop bounds so that only required fields are computed
58      - add new averaging code for zetaZ/deltaZ (again similar to
59        seaice_calc_viscosities.F), inactive by default for now (will change
60        results again), but will replace old code eventually
61    o eesupp/src:
62      - add an other alternative for multi-proc GLOBAL_SUM (in global_sum_tile.F)
63        to always cumulate tile local-sum in the same order
64        (option #define GLOBAL_SUM_ORDER_TILES in CPP_EEOPTIONS.h).
65        Still uses MPI_Allreduce (like standard code) but on array of tiles
66        and should be faster than GLOBAL_SUM_SEND_RECV code.
67      - switch almost all verification experiments that use GLOBAL_SUM_SEND_RECV
68        (defined in CPP_EEOPTIONS.h_mpi) to GLOBAL_SUM_ORDER_TILES instead.
69      - add CPP_EEOPTIONS.h_mpi with GLOBAL_SUM_ORDER_TILES defined in 14 more exp.:
70         aim.5l_cs, cheapAML_box, deep_anelastic, dome, fizhi-cs-32x32x40,
71         fizhi-cs-aqualev20, fizhi-gridalt-hs, global_ocean_ebm, matrix_example,
72         natl_box, solid-body.cs-32x32x1, tutorial_advection_in_gyre,
73         tutorial_baroclinic_gyre and tutorial_held_suarez_cs.
74    
75    checkpoint65n (2015/07/29)
76    o pkg/thsice:
77      - add a warning when used with pkg/aim but without updating SST
78        (stepFwd_oceMxL=F & useCoupler=F & tauRelax_MxL<>-1 ).
79    o pkg/offline:
80      - change where Offline-Fields are loaded, from load_fields_driver.F to
81        directly within FORWARD_STEP, just after updating iter number and time;
82      - move setting of GM and KPP diffusivity out of OFFLINE_FIELDS_LOAD into new
83        S/R OFFLINE_GET_DIFFUS that is called towards the end of DO_OCEANIC_PHYS;
84       These changes allow:
85        a) to recover exactly on-line results using deltaToffline=deltaTClock,
86        tested (using offlineIter0= nIter0 -1 and offlineTimeOffset= -deltaT/2)
87        both with Lin & Non-Lin Free-Surf, both with staggerTimeStep & synchronous
88        time-step, both with GM & KPP files and with recomputing these fields;
89        b) no needs to switch on staggerTimeStep anymore;
90        c) simplify rescaling of offline horiz. velocity with Non-Lin Free-Surf.
91      - update output of experiment tutorial_dic_adjoffline (only 3 digits left).
92    o pkg/offline:
93      - change (rarely used) offset-iteration parameter to an offset time (in s):
94        can easily recover previous offset-iter effect by setting offset-time
95        (offlineTimeOffset) to offset-iter x deltaToffline;
96      - do not turn-off Temp & Salinity surface relaxation when using KPP (used
97        in KPP as part of surface buoyancy forcing) in offline run.
98      - fix ptracers KPP non-local term in the case where KPP is re-computed.
99    o OpenAD
100      - Adding preliminary ggl90 test to OpenAD verification.
101        Config. with ggl90 currently diverges after 64 time steps
102        so will need more work to become more stable
103      - Adding preliminary kpp test to OpenAD verification.
104        required small tweakes in kpp_routines.F
105      - Merge and update from Krishna Narayanan's contrib area:
106        * Add one variable for preliminary OpenAD DIVA
107        * Add replacement time-stepping for OpenAD DIVA
108        * genmake2 flag -diva (but only for OpenAD; preliminary)
109        * required modifs for OAD_support
110        successful preliminary test for tutorial_global_oce_biogeo/
111    o pkg/seaice:
112      - initialise deltaC in seaice_lsr, cosmetic changes in seaice_lsr
113      - add cpp-flag SEAICE_DELTA_SMOOTHREG for option of regularising
114        delta with a smooth function in s/r seaice_calc_viscosity
115        (no effect on EVP)
116      - rename local suffix "sqr" into "sq" for more consistent variable names
117    o pkg/seaice:
118        add new parameter SEAICE_deltaMin that is used ***only***
119        for regularizing Delta (and nothing else, like the parameter
120        which is also used for all sorts of things SEAICE_EPS).
121        Defaults to SEAICE_EPS for backward compatibility only.
122    o pkg/seaice: replace erroneously storing e12Csqr on the wrong tape with
123        the wrong key by inialising it before use -> fixes a recomputation
124        problem but does not fix the gradient
125    o pkg/layers:
126      - move resetting of flux arrays to outside the iLa loop
127        (this caused layers_diapycnal to not work with more than one layer coord.)
128      - fix bug in calculation of layers_recip_delta in layers_init_fixed.F
129    o pkg/layers:
130      - remove from common block 3D layers fields (used to diagnose transport
131        in tracer coords) and declare them as local arrays in layers_calc.F
132    o pkg/seaice:
133      - Rename SEAICE_VECTORIZE_LSR_ZEBRA to SEAICE_LSR_ZEBRA and define it
134        in lab_sea forward experiment (changes results -> updated)
135      - update adjoint experiment lab_sea after cleaning up seaice_lsr_tridiagu/v
136    o pkg/seaice: modify seaice_lsr.F in order to improve the gradient
137       computations; for SEAICE_VECTORIZE_LSR
138      - move the loops over which the tridiagonal solvers (seaice_lsr_tridiagu/v)
139        sweep to the outside,
140      - remove store directives and add "CADJ loop sequential" directives
141        in analogy to model/src/solve_tridiagonal.F
142      - replace many "#ifdef SEAICE_VECTORIZE_LSR_ZEBRA" by variable loop
143        boundaries to yield more readable code. This has also the charming side
144        effect that your can use SEAICE_VECTORIZE_LSR_ZEBRA without
145        SEAICE_VECTORIZE_LSR (i.e. adjoint without recomputations in
146        seaice_lsr_tridiagu/v now requires either of these flags, vectorization
147        still requires SEAICE_VECTORIZE_LSR with SEAICE_VECTORIZE_LSR_ZEBRA as
148        an additional option)
149       The tridiagonal solvers are now completely analoguous to solve_tridiagonal.F
150       but the adjoint gradients (which are affected by this change) still explode.
151    o coupler (atm & ocn _compon_interf):
152      - avoid using Exch2 IO mapping if not suitable.
153    
154    checkpoint65m (2015/06/15)
155    o pkg/layers:
156      - clean-up old parameters setting (from Sep 2012):
157        + remove old parameters setting from common block (declared locally);
158        + stop if any old parameter is mixed with new parameters;
159        + check for valid or consistent parameter setting (e.g., layers_name)
160      - check that layers_bounds for Pot.Density are in the right range;
161        this is supposed to catch un-updated setting after recently shifting
162        potential density field "prho" by -1000 (now in "sigma" scale).
163    o pkg/layers:
164      - Added missing call to layers_fill for keeping track of tottend in
165        pkg/diagnostics/diagnostics_fill_state.F
166    o pkg/layers:
167      - added CPP flux LAYERS_FINEGRID_DIAPYCNAL to disable fine-grid binning
168        of diapycnal fluxes
169    o pkg/layers:
170      - reverted changes to interpolation to previous treatment; there
171        was never actually a problem
172    o pkg/layers:
173      - fix interpolation of layers tracers in the presence of partial cells
174      - add lots of new functionality to compute water mass transformation
175        (enabled by CPP flag LAYERS_THERMODYNAMICS)
176      - add new calls from main code (diags_oceanic_surf_flux.F,
177        thermodynamics.F, impldiff.F) and gad to save the necessary fields
178        for layers_thermodynamics calculations
179    o pkg/ptracers:
180      - add capability of stepping a passive tracer forward without advection
181        by setting PTRACERS_advScheme to 0.
182    o pkg/diagnostics:
183      - add few (threads) _BARRIER around any update of "diag_pkgStatus".
184    o pkg/layers:
185      - shift potential density field "prho" by -1000 for density layer averaging so
186        that it corresponds to usual "sigma" definition (i.e., rho(p=const)-1000);
187        Note that density bin boundaries (layers_bounds in data.layers) have
188        to be set accordingly (i.e., shifted, in sigma scale).
189    o pkg/layers:
190      - fix snap-shot output of layer thickness @ V.pt and output of pot.density
191      - fix calculation of layer non-weighted velocity and layer probability
192      - add diagnostics for pot.density used for layer averaging
193    o pkg/seaice:
194      - seaice_calc_viscosity.F: change computation of deltaC**2 to ensure
195        positiveness, modify a few comments and improve variable names,
196      - the changed computation of deltaC affects most seaice-related
197        experiments at the truncation level.
198        seaice_itd.thermo is affected by far the most, with
199        only 7 digits of agreement remaining; update experiments.
200    
201    checkpoint65l (2015/05/04)
202    o pkg/shelfice:
203      - in shelfice_thermodynamics.F, if SHELFICEBoundaryLayer=.true., then set
204        factor drKp1 to zero where it is negative (ie. do not use
205        SHELFICEBoundaryLayer if hFacC>1)
206      - the same for factor drLoc in shelfice_forcing.F
207    o tools/genmake2: move FEXTRAFLAGS from FFLAGS to FOPTIM to avoid a problem
208        with modules in ptracers, when using testreport option "-small_f"
209    o pkg/seaice/seaice_growth.F (in case of SEAICE_ITD)
210      - replace tmpscal1**1.36 by faster exp(1.36*log(tmpscal1))
211    o verification/testreport
212      - add new option -small_f
213        with this option testreport first makes the target small_f (*.f files)
214        before running make (all). This is useful for explicit inlining when
215        all sources need to be full available.
216    o pkg/exf/exf_interp.F:
217        remove previous changes related to reverse (decreasing) input grids latitude
218    
219    checkpoint65k (2015/04/01)
220    o pkg/ecco:
221        fix missing cosphi initialization; use eccoiter in place of optimcycle;
222        if autodiff is not compiled then use READ_REC_XY_RL/READ_REC_XYZ_RL
223        instead of active read/write
224    o pkg/ctrl:
225        if autodiff is not compiled then use READ_REC_XY_RL/READ_REC_XYZ_RL
226        instead of active read/write.
227    o pkg/autodiff: make autodiff_readparms.F and autodiff_whtapeio_sync.F
228        return if .NOT.useAUTODIFF
229    o pkg/seaice: add CPP brackets (cost function related).
230    o model/src:
231        add useAUTODIFF run time switch.
232    o pkg/exf/exf_interp.F:
233        add code to allow input grids with latitude starting in the north (i.e. when
234        j=1 corresponds to northern edge of field); modify global_with_exf to test
235        this feature.
236    o verification/testreport:
237      - replace "grep -A3" with POSIX conformal sed command, because there are
238        implementations of grep that cannot do contextual searches
239      - replace "cat -n" by "nl" (two occurrences), because not all implementations
240        of "cat" have the option "-n"
241    o model/src:
242        add new file gsw_teos10.F that contains 3 functions, more or less a copy
243        of the gsw_toolbox.f90 of the TEOS-10 fortran interface, that are used to
244        convert from conservative to potential temperature. For now they are not
245        called anywhere.
246    o pkg/ggl90: make "useIDEMIX" useful
247      - re-organize code in s/r ggl90_calc a little to move computations that
248        are different for idemix into separate loops (changes verification
249        experiment idemix at truncation level, 12 digits remain -> updated)
250      - requires 2 new 2D fields KappaM, verticalShear
251      - make hFacI and recip_hFacI local 3D variables that are passed to
252        s/r ggl90_idemix
253      - fix and improve some comments
254    
255    checkpoint65j (2015/02/25)
256    o pkg/openad
257      - synchronize seawater.F with version in model/src: replace routine
258        body of sw_temp by calling sw_ptmp and fix/add a few comments.
259    o pkg/mnc:
260      - change default value of pickup_read/write_mnc params to FALSE
261        (since few features are missing in reading MNC pickup files).
262    o model/src:
263      - change default (background) vertical diffusivity for salt to be the
264        temperature (background) vertical diffusivity.
265      - change background vertical diffusivity in vertical mixing pkgs ggl90,
266        kl10, my82 and pp81 from temperature diffusivity to salinity diffusivity.
267        This makes ptracers default diffusivity (that uses salt diffKr) more
268        consistent with vertical mixing schemes.
269    o verification/testreport:
270      - fix typo in "-pcls" option ;
271      - put temp files in /tmp (instead of in local dir)
272      - add option "-ncad" to use genmake2 new option "-nocat4ad"
273      - allow to use "-j" option (for multi-procs make) for AD, OAD and TLM.
274    o pkg/ggl90 (idemix code):
275      - apply few fixes ; add some missing "_d 0" ; update output of exp.
276        global_ocean.90x40x15.idemix
277    o OpenAD
278      - enable active file I/O and begin removal of work-around code
279        for OpenAD for handling control variables (memory vs. active files)
280      - new CPP options file OPENAD_OPTIONS.h
281      - revert to genarr2d gradient check for hs94.1x64x5
282      - (genarr3d tested by tutorial_global_oce_biogeo, tutorial_tracer_adjsens)
283      - gentim2d works for halfpipe_streamice but not yet hs94.1x64x5
284    o tools/genmake2:
285      - add option for an alternative Makefile (from Martin) where multiple
286        source files are sent to TAF (instead of all in one file).
287    o pkg/ggl90: add IDEMIX (Olbers and Eden, 2013)
288      - code provided by Carsten Eden as an extension of ggl90
289      - so far the code is turned on within ggl90 by setting a CPP-flag at
290        compile time; a runtime flag implementation is not yet complete
291      - tested in new verification experiment global_ocean.90x40x15.idemix
292    o verification:
293      - reduce length of the AD test run by half for experiments
294        global_ocean.cs32x15 (std + thsice) and isomip
295      - enable active file I/O for OpenAD, adjust verifications accordingly
296      - halfpipe_streamice: change gradient check from genarr2d to gentim2d
297    o verification/testreport:
298      - fix two non-POSIX sed statements so that there are no error messages
299        on Mac OS X
300      - replace hostname -s with hostname | sed 's/\..*$//'
301    o model/src/seawater.F:
302      - simplify SW_TEMP to just calling SW_PTMP with P and PR exchanged
303        (following the matlab code), which would have avoided the bug caught
304        by Matt (see below) in the first place
305    o pkg/shelfice:
306      - implement alternative discretisation of quadratic drag (selectBotDragQuadr)
307        and true distance in no-slip BC (bottomVisc_pCell) similar to bottom-drag.
308      - fix sign of friction tendency (no-slip BC) and fix u,v drag for thin
309        ice-shelf case (i.e., within surface level).
310      - add config summary in shelfice_check.F
311      - remove global variable recip_SHELFICElatentHeat and replace by a local
312        variable recip_latentHeat
313      - introduce diagnostics for friction velocity uStar (in case of
314        SHELFICEuseGammaFrict)
315      - change verification experiment isomip.htd:
316        + start from pickup (as input_ad at nIter0=8640)
317        + start from SHELFICEmassFile instead of SHELFICEloadAnomaly file
318        + test SHELFICEuseGammaFrict = .TRUE.
319        + test SHELFICEadvDiffHeatFlux = .TRUE.
320    o model/src:
321      - fix bug (Thanks to Matt) in function SW_TEMP (to convert Pot.Temp to
322        in-situ Temp) and in routine SW_TEMP (in pkg/openad);
323        update output of experiments global_with_exf.yearly and isomip (all 6).
324  o pkg/shelfice:  o pkg/shelfice:
325    - change setting of kTopC to also be zero where there is no ice-shelf.    - change setting of kTopC to also be zero where there is no ice-shelf.
326    - zero out surface forcing arrays (for now, for T & S + EmPmR) where    - zero out surface forcing arrays (for now, for T & S + EmPmR) where
327      ice-shelf is.      ice-shelf is.
328    - first attempt to get RealFreshWaterFlux effects: for now, only without    - first attempt to get RealFreshWaterFlux effects: for now, only without
329      SHELFICEboundaryLayer ; not consistent if using synchronous time-stepping.      SHELFICEboundaryLayer ; not consistent if using synchronous time-stepping.

Legend:
Removed from v.1.2121  
changed lines
  Added in v.1.2183

  ViewVC Help
Powered by ViewVC 1.1.22