/[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.1793 by jmc, Thu Dec 27 17:20:52 2012 UTC revision 1.1853 by gforget, Thu Apr 11 14:52:12 2013 UTC
# Line 4  $Name$ Line 4  $Name$
4      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
5      ==============================      ==============================
6    
7    o verification/offline_exf_seaice : revise data.seaice and
8      update results accordingly in input.thermo and input_ad.
9        HO = 0.2, SEAICE_areaLossFormula = 2,
10        SEAICE_growMeltByConv = .TRUE.,
11        SEAICE_doOpenWaterMelt = .TRUE.,
12      Improve consistency with input.thsice and
13      fixes the spurious 500m ice corners.
14    o pkg/cd_code:
15      - swtich off Adams-Bashforth on the 1rst iteration. Update output of
16        experiment lab_sea.hb87
17    o tools/genmake2:
18      - change default makedepend from the system default to tools/xmakedepend:
19        When no makedepend is specified, select tools/xmakedepend it the test pass,
20        otherwise select the system default (if exist and pass the test) or
21        try to build cyrus-makedepend ; stop if none of the 3 works.
22    o pkg/exf & pkg/thsice:
23      - In case Energy-Reference-Level (ERL) is used (temp_EvPrRn=0), account
24        for energy content of Precip + RunOff & Evap. Assumes:
25       1) Rain has same temp as Air.
26       2) Snow has no heat capacity (consistent with seaice & thsice pkgs)
27       3) Run-Off comes at the temp of surface water (with same Cp)
28       4) Evap over ocean is released to the Atmos @ surf-temp (=SST); should be
29        put directly into Latent Heat flux & should use water-vapor heat capacity.
30       5) Evap over snow/seaice (pkg/thsice) assumes constant Lf+Lv (consistent
31        with snow having no Cp) => water vapor is released to Atmos @ 0^oC (=ERL).
32    o pkg/seaice: simplify the use of CPP flags in JFNK-related files
33        (compile when SEAICE_ALLOW_JFNK is defined)
34    o pkg/thsice:
35      - fix oceanic flux diagnostics & timeave (left from new sequence of calls
36        implementation): move THSICE_AVE call after the advection (called directly
37        from thsice_do_advect.F).
38    o verification/global_ocean.90x40x15 and OpenAD:
39      New experiment for OpenAD is identical to the one in input_ad except that
40        multiDimAdvection=.FALSE.,
41      There is still an issue with this.
42    o pkg/thsice:
43      - account for ice velocity (ocean relative velocity) in turbulent heat flux
44        at the base of the ice; update output of global_ocean.cs32x15 (icedyn)
45        and offline_exf_seaice (main test).
46    
47    checkpoint64f (2013/03/31)
48    o verification/OpenAD:
49      - new (smaller, without OAD_DEBUG) output_oadm.txt generated using gfortran.
50    o verification/isomip, pkg/shelfice and OpenAD:
51      - add OpenAD setup to this exp and make shelfice package work with OpenAD
52    o verification/tutorial_global_oce_biogeo:
53      - update AD and TLM output after turning off ALLOW_FE (in DIC_OPTIONS.h)
54    o verification/offline_exf_seaice:
55      - replace adjoint experiment (now similar to fwd input.thermo/ setup).
56    o verification/tutorial_tracer_adjsens, pkg/ptracers and OpenAD:
57      - add OpenAD setup to this exp and make ptracers package work with OpenAD
58    o pkg/seaice/seaice_evp.F:
59       mask the asymmetric contribution of the ice-ocean stress (when
60       the turningAngle is non-zero) over open water to avoid stripes
61       in the EVP-solutions after many EVP-substeps. This does not change
62       the verification experiment (because turningAngle=0 in that experiment),
63       but will affect other EVP solutions with turningAngle.ne.0
64    o verification, AD & TLM output:
65      - new TAF version 2.3.10 (previously 2.3.8) add more parenthesis in multiple
66        additions and multiplications which affect results at machine truncation
67        level. Update few AD output (1D_ocean_ice_column, bottom_ctrl_5x5,
68        lab_sea.noseaicedyn & offline_exf_seaice) and TLM output
69        (1D_ocean_ice_column & tutorial_tracer_adjsens.som81) which had only 12
70        digits of matching.
71    o pkg/seaice:
72       modify lsr-algorithm so that it vectorizes AND converges by using a
73       "zebra" method: split vector loop into two with step size 2; this
74       solves a general convergence problem for the JFNK solver on a vector
75       computer. It has the potential to even improve convergence for
76       non-vector machines, but it does change the results so it is turned
77       off by default. To turn on define CPP-flags SEAICE_VECTORIZE_LSR and
78       SEAICE_VECTORIZE_LSR_ZEBRA in SEAICE_OPTIONS.h
79    
80    checkpoint64e (2013/03/05)
81    o pkg/ecco and several verification experiments:
82       Remove special versions of the_model_main.F, the_main_loop.F, forward_step.F
83       from ecco package, and revert to unified routine in model/src/
84       Adjust verification experiments accordingly, especially gloal_oce_llc90
85       Should conclude long-standing work by Gael to reconsile the_main_loop with
86       ECCO and will enable OpenAD's revolve loop in ECCO setup.
87    o pkg/generic_advdiff:
88      - implement compressible flow method for multi-dim advection (similar to
89        gad_som_advect.F); add new option "GAD_MULTIDIM_COMPRESSIBLE" (since
90        TAF generates many recomputations) to use it; no yet coded with implicit
91        vertical advection.
92      - change experiment advect_cs (new output) and advect_xy (same output)
93        to test compressible flow multi-dim advection.
94    o pkg/generic_advdiff:
95      - fix SOM advection for unusual cases (linear Free-Surf in r*, other
96        cases where linear Free-Surf is not @ k=1, e.g., useShelfIce)
97    o pkg/seaice:
98      - fix bug in the no-slip implementation in seaice_calc_strainrates.F, this
99        bug affects two verification experiments (lab_sea and lab_sea.hb87),
100        recover old (wrong) results by setting noSlipFac = 2. in
101        seaice_calc_strainrates.F;
102      - add another averaging method for etaZ (etaZmethod=3), a mix of 2 and 0
103      - add runtime parameter SEAICE_JFNKepsilon for that finite difference
104        gradient in the Jacobian times vector operation (s/r seaice_jacvec)
105      - fix line search in seaice_jfnk.F (s/r seaice_jfnk_update)
106    o pkg/seaice: reorganize seaice_lsr.F to be able to reuse more code:
107      - new subroutines seaice_lsr_tridiagu/v
108      - introduce loop bounds i/jmin/max to be able to use a restricted
109        additive Schwarz method, but not used for now.
110      - rename SEAICE_CALC_LSR_COEFFS into SEAICE_LSR_CALC_COEFFS
111      - use new subroutines seaice_lsr_tridiagu/v in s/r seaice_preconditioner
112    o OpenAD support
113       Reorganize S/R call for overall time-stepping loop to facilitate
114       OpenAD revolve loop replacement
115       (new S/R main_do_loop.F)
116    o model/src:
117      - add 3-D array 3rd dimension to argument list of S/R ADAMS_BASHFORTH2/3;
118      - pkg/cheapaml: remove local S/R ADAMS2d and use instead standard S/R
119        ADAMS_BASHFORTH2.
120    o pkg/cheapaml:
121      - add new parameter (default=F, same as before) to disable full domain
122        (including land) periodicity in X & Y dir.
123      - new S/R CHEAPAML_COPY_EDGES to copy relaxation field into edges of
124        prognostics field ; merge the 2 S/R cycle_2d_tracer and timestep_2d_tracer
125        into single S/R CHEAPAML_TIMESTEP.
126    
127    checkpoint64d (2013/02/17)
128    o pkg/monitor (angular momentum):
129      - do not add Eta contribution to the total AM if using rigid-lid;
130      - correct for Adams-Bashforth in zonal-wind part (instead of in AM_eta).
131      - horizontal discretisation: compute Zonal wind contribution separately for
132        each component (instead of from cell centered averaged) and use the 4
133        grid-cell corner averaged value of Omega*(R*cos_Phi)^2 in AM_eta part.
134      - Check that AM is conserved with linear dynamics (momAdvection=F and linear
135        FreeSurf) if using vectorInvariant momentum.
136    o model/src (horizontal grid):
137      - move initialisation of all horizontal grid arrays to null value
138        from specific horizontal grid S/R to main S/R INI_GRID;
139      - add 2 new arrays to project model C-grid velocity on to zonal direction
140        (similar to AngleCosC & -AngleSinC but at U and V point location) and
141        compute them in S/R CALC_GRID_ANGLES.
142    o pkg/seaice: replace a hard coded parameter by a runtime parameter
143        SEAICE_JFNK_tolIter
144    o pkg/seaice: introduce a restricted additive Schwarz (RAS) method in
145        preconditioner for JFNK (run LSR into the overlaps as far as possible,
146        i.e OLx/y-2, but overwrite the solution in the overlaps)
147      - extend loop range in seaice_calc_viscosity
148      - add new parameters SEAICE_OLx/y (default = OLx/y-2)
149      - add flexible loop ranges to seaice_preconditioner, seaice_calc_lsr_coeffs
150      - slightly modify loops in tridiagonal solver lsr-preconditioner
151      - by default remove the extra exchange in seaice_preconditioner with a
152        replacement of an undocumented CPP flag by SEAICE_PRECOND_EXTRA_EXCHANGE
153      - modify offline_exf_seaice so that old results without RAS are reproduced
154        for output.dyn_jfnk.txt
155      in seaice_fgmres.F (unrelated to RAS)
156      - add a cpp-flag option for a re-ordered mapping (potentially more
157        efficient), but because this changes the results, turn it off by
158        default
159      - replace the save statement by local common blocks
160    o utils/matlab:
161      - fix tile-offset in rdmnc.m when output was produced using EXCH2.
162      - update cs_grid/bk_line scripts to use load_grid.m (in utils/matlab)
163        instead of local script "load_cs.m".
164    o pkg/seaice:
165      - fix a serious bug in seaice_fgmres.F, s/r seaice_scalprod, connected
166        to the scalarProductMetric and the unrolling. By chance this does not
167        change the verification experiment, because of the particular domain
168        decomposition
169    o pkg/diagnostics:
170      - change default missing_value from UNSET_RL to -999.
171    o pkg/layers: improve vectorization of layers_fluxcalc.F
172      - step 1: split loops, make TatU/V 2D arrays. This alone reduces the cpu
173        time spent in this routine by more than a factor of 3 on an NEC-SX8
174      - step 2: replace search algorithm by a bisection algorithm found in
175        Press et al. Numerical Recipes in Fortran. This algorithm can be
176        vectorized and finds the correct bin in approx. log2(Nlayers) steps.
177        For now, use this algorithm only for TARGET_NEC_SX defined.
178      - moved the search algorithm into a separate subroutine, that is kept
179        within the same file for automatic inlining.
180      - modified catching potential errors
181    o pkg/ecco:
182      - add compile switches ALLOW_SHALLOW_ALTIMETRY and ALLOW_HIGHLAT_ALTIMETRY.
183    o model/src, pkg/autodiff, pkg/ctrl, pkg/grdchk :
184      - modifications for ALLOW_GENTIM2D_CONTROL
185      -- forward_step.F, the_main_loop.F : pass options and arrays for adjoint tapes.
186      -- load_fields_driver.F : add call to CTRL_MAP_GENTIM2D( myTime, myIter, myThid )
187      -- add xx_gentim2d0/1 in ctrl_ad_check_lev1_dir.h checkpoint_lev1_directives.h etc
188      -- use ctrl_swapffields.F (new) instead of exf_SwapFFields.F in ctrl_get_gen.F
189    o model/src:
190      - account for moist thermodynamics in (dry) potential temperature equation;
191        update results of experiments: aim.5l_cs (x2), fizhi-cs (std + aqualev)
192        and cpl_aim+ocn.
193      - also switch to standard compiler (gfortran, was using ifort previously)
194        for fizhi-cs experiments reference output.
195    o pkg/ctrl:
196      - 1st step toward generic time-varying control variable (2D) via CPP option
197        ALLOW_GENTIM2D_CONTROL
198      - note new namelist ctrl_nml_genarr in data.ctrl
199    o verification/global_ocean.cs32x15:
200      - switch off OLD_THSICE_CALL_SEQUENCE (both forward and AD code) to use new
201        thsice calling sequence and update "output.icedyn.txt" output file.
202    o pkg/ptracers:
203      - from Ryan A.: code to allow to restart a tracer at regular time interval
204        specifying a reset frequency and a reset-phase (for each tracers).
205    o pkg/thsice:
206      - implement new sequence of calls for thsice+seaice:
207          previously:   ice-Dyn,ice-Advect,ice-Thermo(thsice)
208          new sequence: ice-Thermo(thsice),ice-Dyn,ice-Advect
209      - allows (with temporary CPP option "#define OLD_THSICE_CALL_SEQUENCE"
210        in CPP_OPTIONS.h) to recover old sequence;
211      - update diagnostics for Intermediate State thsice diagnostics: now filled
212        just before advection (instead of just after advection in old sequence);
213      - update output of experiment offline_exf_seaice (standard test);
214      - in global_ocean.cs32x15 exp., postpone output.icedyn.txt update
215        by setting #define OLD_THSICE_CALL_SEQUENCE in CPP_OPTIONS.h
216    o pkg/aim_v23:
217      - use tiled array for aim_CO2 and surface wind speed for export to coupler;
218      - use standard pkg/mdsio routines to read & write CO2 pickup file;
219      - in-line S/R AIM_AIM2DYN_EXCHANGES in aim_do_physics.F;
220      - Implement new sequence of calls for seaice advection (with thsice and
221        seaice diffusion): do ice advection after ice-Thermodynamic (previously
222        was done before); change aim_do_physics.F (split the bi,bj loops block);
223      - use same CPP option (OLD_THSICE_CALL_SEQUENCE) to recover old sequence;
224    o pkg/thsice:
225      - add a missing EXCH of sIceLoad in thsice_main.F (ocean set-up) for the case
226        where just seaice diffusion is used without seaice-dynamics.
227    
228    checkpoint64c (2013/01/18)
229    o pkg/seaice:
230      - add a parameter JFNKres_tFac that can replace JFNKres_t
231        (= JFNKres_tFac*JFNKresidual before JFNK solver)
232      - add a line search option to the JFNK solver
233      - move update of u/vIce and computation of initial residual into a separate
234        routine (same file); makes the code a little easier to read (?)
235      - the line search requires computation of the residual and for some reason
236        I cannot do the update of u/vIce at the beginning of the Newton loop
237        without changing the results, so we move computation of the new residual
238        also to the end of the Newton loop, even without the line search
239    o pkg/mdsio:
240      - write missing value (corresponding to undefined array value) to meta files;
241        skip this if missing value is one (UNSET_RL would be a better value for the
242        case of undefined missing value but it's current the default missing-value);
243        Also update mdsio_read_meta.F (reading timeInterval & missing-value).
244    o pkg/exf & thsice:
245      - for diagnostics purpose, fill-up and update Qsw array (+ fill the
246        diagnostics) even when SHORTWAVE_HEATING is #undef;
247    o verification/cfc_example:
248      - turn on pkg/layers with pkg/diagnostics (to test pkg/layers and to
249        illustrate how to use it).
250    o verification/offline_exf_seaice:
251      in input: add SEAICErestoreUnderIce=.TRUE., and reduce size of forcing
252        files to only one time-record.
253      in input.thermo: narrow the 2 bands of low initial ice-concentration near
254        N & S edges; switch on SEAICErestoreUnderIce=T; use pkg/diagnostics for
255        snapshot output; same mcPheePiston value as in thsice with
256        surf-current=0.2m/s.
257      in input.thsice: similar set-up as in input.thermo (channel with same seaice
258        initial conditions).
259    o verification/cfc_example:
260      - update file "data" with more standard choice (default rotationPeriod,
261        exactConserv=T, reduce ivdc_kappa to 10, remove multiDimAdvection=F,
262        set deltaTFreeSurf to deltaTtracer) and update the output.
263    o pkg/layers:
264      - fix bolus part so that it compiles with #undef GM_BOLUS_ADVEC;
265      - ensure that "layers_bolus" is switched to False when useGMRedi=False.
266    o verification/offline_exf_seaice:
267      - set all {field}period to zero for constant-in-time forcing (not tested
268        anywhere else).
269    o verification/offline_exf_seaice:
270      - dyn_lsr  : use pkg/seaice S/R (instead of pkg/thsice) to advect sea-ice;
271      - dyn_jfnk : use pkg/thsice S/R (instead of pkg/seaice) to advect sea-ice.
272    o model/src:
273      - implement Finite-Volume method for (hydrostatic) presure gradient
274        from S.-J. Lin (QJRMS 1997), for atmosphere using sigma-coordinate.
275    o pkg/seaice: implement EVP* by Lemieux et al. (2012).
276      - In order to be compatible with the default EVP code, this requires changing
277        the EVP code around a little; and because any small change immediately
278        breaks the lab_sea.hb87 verification experiment, the reference output is
279        updated two times (each time the remaining aggrement was only 2 digits):
280        - after replacing 1/dt by 1*recip_dt with a precomputed recip_dt
281        - after replacing 1/(1+0.5*dt/tau) by a/(a+1), with a=2*tau/dt
282        The adjoint test lab_sea.evp also needed to be updated.
283      - Turn on EVP* by setting SEAICE_deltaTevp as before, and setting
284        SEAICEuseEVPstar=.true.; for EVP* SEAICE_deltaTevp does not need to be
285        a factor of SEAICE_deltaTdyn, and you can set the number of EVP steps
286        explicitly with SEAICEnEVPstarSteps (thereby overriding the default of
287        INT(SEAICE_deltaTdyn/SEAICE_deltaTevp)). Controlling EVP* is deliberately
288        kept obscure, because I need to test it a little more.
289    o pkg/cheapaml:
290      - fix few incorrect type (_RL/_RS) in subroutine calls.
291    o pkg/seaice:
292      - new S/R SEAICE_MAP_RS2VEC (similar to SEAICE_MAP2VEC) to map _RS arrays
293        to _RL vector; this fixes JFNK code when compiling with "-use_r4" option.
294    o model/src:
295      - add a warning if using select_rStar=2 with nonlinFreeSurf <> 4;
296        change (in input_ad*/data) select_rStar from 2 to 1 for AD verification
297        experiments global_ocean.cs32x15 & tutorial_tracer_adjsens (+ update output)
298    o pkg/cheapaml:
299      - replace "273.1 _d 0" by celsius2K in pkg/cheapaml/cheapaml.F and
300        update output of experiment cheapAML_box.
301      - store in CHEAPAML.h coeffs ssq[0:2] for saturation specific humidity,
302        and set them only in one place (cheapaml_readparms.F).
303      - replace "273.16 _d 0" by celsius2K & replace "3.14159 _d 0" by constant "PI"
304      - remove re-setting of main-model parameter "celsius2K" (in cheapaml_readparms)
305        and re-setting of cheapAML params "lath" & "p0" (in cheapaml_init_varia).
306    o model/src/the_main_loop.F:
307      - fix revision 129 : need both ALLOW_ECCO and ALLOW_ECCO_EVOLUTION.
308    o pkg/seaice:
309      - remove legacy branch code.
310      - retire SEAICE_GROWTH_LEGACY and SEAICE_CAP_HEFF accordingly.
311    o verification/*/code*/SEAICE_OPTIONS.h:
312      - remove documentation for retired CPP flags (SEAICE_GROWTH_LEGACY,
313        SEAICE_CAP_HEFF).
314    o pkg/cheapaml:
315      - use standard routine (S/R GET_PERIODIC_INTERVAL) to decide when it is time
316        to load new forcing fields.
317      - add missing "_d 0" where it's needed; + some editing (improve indentation).
318      - remove switch to "-r8" compiler option (no longer needed) in verification
319        cheapAML_box, keep same output (small trucation error differences).
320  o verification/offline_exf_seaice:  o verification/offline_exf_seaice:
321    - add thermodynamics only test.    - add thermodynamics only test.
322  o verification/cheapAML_box:  o verification/cheapAML_box:

Legend:
Removed from v.1.1793  
changed lines
  Added in v.1.1853

  ViewVC Help
Powered by ViewVC 1.1.22