/[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.1807 by mlosch, Wed Jan 16 21:20:56 2013 UTC revision 1.1827 by jmc, Sun Feb 17 23:10:07 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 pkg/seaice: add a line search option to the JFNK solver  checkpoint64d (2013/02/17)
8    o pkg/monitor (angular momentum):
9      - do not add Eta contribution to the total AM if using rigid-lid;
10      - correct for Adams-Bashforth in zonal-wind part (instead of in AM_eta).
11      - horizontal discretisation: compute Zonal wind contribution separately for
12        each component (instead of from cell centered averaged) and use the 4
13        grid-cell corner averaged value of Omega*(R*cos_Phi)^2 in AM_eta part.
14      - Check that AM is conserved with linear dynamics (momAdvection=F and linear
15        FreeSurf) if using vectorInvariant momentum.
16    o model/src (horizontal grid):
17      - move initialisation of all horizontal grid arrays to null value
18        from specific horizontal grid S/R to main S/R INI_GRID;
19      - add 2 new arrays to project model C-grid velocity on to zonal direction
20        (similar to AngleCosC & -AngleSinC but at U and V point location) and
21        compute them in S/R CALC_GRID_ANGLES.
22    o pkg/seaice: replace a hard coded parameter by a runtime parameter
23        SEAICE_JFNK_tolIter
24    o pkg/seaice: introduce a restricted additive Schwarz (RAS) method in
25        preconditioner for JFNK (run LSR into the overlaps as far as possible,
26        i.e OLx/y-2, but overwrite the solution in the overlaps)
27      - extend loop range in seaice_calc_viscosity
28      - add new parameters SEAICE_OLx/y (default = OLx/y-2)
29      - add flexible loop ranges to seaice_preconditioner, seaice_calc_lsr_coeffs
30      - slightly modify loops in tridiagonal solver lsr-preconditioner
31      - by default remove the extra exchange in seaice_preconditioner with a
32        replacement of an undocumented CPP flag by SEAICE_PRECOND_EXTRA_EXCHANGE
33      - modify offline_exf_seaice so that old results without RAS are reproduced
34        for output.dyn_jfnk.txt
35      in seaice_fgmres.F (unrelated to RAS)
36      - add a cpp-flag option for a re-ordered mapping (potentially more
37        efficient), but because this changes the results, turn it off by
38        default
39      - replace the save statement by local common blocks
40    o utils/matlab:
41      - fix tile-offset in rdmnc.m when output was produced using EXCH2.
42      - update cs_grid/bk_line scripts to use load_grid.m (in utils/matlab)
43        instead of local script "load_cs.m".
44    o pkg/seaice:
45      - fix a serious bug in seaice_fgmres.F, s/r seaice_scalprod, connected
46        to the scalarProductMetric and the unrolling. By chance this does not
47        change the verification experiment, because of the particular domain
48        decomposition
49    o pkg/diagnostics:
50      - change default missing_value from UNSET_RL to -999.
51    o pkg/layers: improve vectorization of layers_fluxcalc.F
52      - step 1: split loops, make TatU/V 2D arrays. This alone reduces the cpu
53        time spent in this routine by more than a factor of 3 on an NEC-SX8
54      - step 2: replace search algorithm by a bisection algorithm found in
55        Press et al. Numerical Recipes in Fortran. This algorithm can be
56        vectorized and finds the correct bin in approx. log2(Nlayers) steps.
57        For now, use this algorithm only for TARGET_NEC_SX defined.
58      - moved the search algorithm into a separate subroutine, that is kept
59        within the same file for automatic inlining.
60      - modified catching potential errors
61    o pkg/ecco:
62      - add compile switches ALLOW_SHALLOW_ALTIMETRY and ALLOW_HIGHLAT_ALTIMETRY.
63    o model/src, pkg/autodiff, pkg/ctrl, pkg/grdchk :
64      - modifications for ALLOW_GENTIM2D_CONTROL
65      -- forward_step.F, the_main_loop.F : pass options and arrays for adjoint tapes.
66      -- load_fields_driver.F : add call to CTRL_MAP_GENTIM2D( myTime, myIter, myThid )
67      -- add xx_gentim2d0/1 in ctrl_ad_check_lev1_dir.h checkpoint_lev1_directives.h etc
68      -- use ctrl_swapffields.F (new) instead of exf_SwapFFields.F in ctrl_get_gen.F
69    o model/src:
70      - account for moist thermodynamics in (dry) potential temperature equation;
71        update results of experiments: aim.5l_cs (x2), fizhi-cs (std + aqualev)
72        and cpl_aim+ocn.
73      - also switch to standard compiler (gfortran, was using ifort previously)
74        for fizhi-cs experiments reference output.
75    o pkg/ctrl:
76      - 1st step toward generic time-varying control variable (2D) via CPP option
77        ALLOW_GENTIM2D_CONTROL
78      - note new namelist ctrl_nml_genarr in data.ctrl
79    o verification/global_ocean.cs32x15:
80      - switch off OLD_THSICE_CALL_SEQUENCE (both forward and AD code) to use new
81        thsice calling sequence and update "output.icedyn.txt" output file.
82    o pkg/ptracers:
83      - from Ryan A.: code to allow to restart a tracer at regular time interval
84        specifying a reset frequency and a reset-phase (for each tracers).
85    o pkg/thsice:
86      - implement new sequence of calls for thsice+seaice:
87          previously:   ice-Dyn,ice-Advect,ice-Thermo(thsice)
88          new sequence: ice-Thermo(thsice),ice-Dyn,ice-Advect
89      - allows (with temporary CPP option "#define OLD_THSICE_CALL_SEQUENCE"
90        in CPP_OPTIONS.h) to recover old sequence;
91      - update diagnostics for Intermediate State thsice diagnostics: now filled
92        just before advection (instead of just after advection in old sequence);
93      - update output of experiment offline_exf_seaice (standard test);
94      - in global_ocean.cs32x15 exp., postpone output.icedyn.txt update
95        by setting #define OLD_THSICE_CALL_SEQUENCE in CPP_OPTIONS.h
96    o pkg/aim_v23:
97      - use tiled array for aim_CO2 and surface wind speed for export to coupler;
98      - use standard pkg/mdsio routines to read & write CO2 pickup file;
99      - in-line S/R AIM_AIM2DYN_EXCHANGES in aim_do_physics.F;
100      - Implement new sequence of calls for seaice advection (with thsice and
101        seaice diffusion): do ice advection after ice-Thermodynamic (previously
102        was done before); change aim_do_physics.F (split the bi,bj loops block);
103      - use same CPP option (OLD_THSICE_CALL_SEQUENCE) to recover old sequence;
104    o pkg/thsice:
105      - add a missing EXCH of sIceLoad in thsice_main.F (ocean set-up) for the case
106        where just seaice diffusion is used without seaice-dynamics.
107    
108    checkpoint64c (2013/01/18)
109    o pkg/seaice:
110      - add a parameter JFNKres_tFac that can replace JFNKres_t
111        (= JFNKres_tFac*JFNKresidual before JFNK solver)
112      - add a line search option to the JFNK solver
113      - move update of u/vIce and computation of initial residual into a separate
114        routine (same file); makes the code a little easier to read (?)
115      - the line search requires computation of the residual and for some reason
116        I cannot do the update of u/vIce at the beginning of the Newton loop
117        without changing the results, so we move computation of the new residual
118        also to the end of the Newton loop, even without the line search
119  o pkg/mdsio:  o pkg/mdsio:
120    - write missing value (corresponding to undefined array value) to meta files;    - write missing value (corresponding to undefined array value) to meta files;
121      skip this if missing value is one (UNSET_RL would be a better value for the      skip this if missing value is one (UNSET_RL would be a better value for the
122      case of undefined missing value but it's current the default missing-value);      case of undefined missing value but it's current the default missing-value);
123      Also update mdsio_read_meta.F (reading timeInterval & missing-value).      Also update mdsio_read_meta.F (reading timeInterval & missing-value).
124  o pkg/exf & thsice:  o pkg/exf & thsice:
125    - for diagnostics purpose, fill-up and update Qsw array (+ fill the diagnostics)    - for diagnostics purpose, fill-up and update Qsw array (+ fill the
126      even when SHORTWAVE_HEATING is #undef;      diagnostics) even when SHORTWAVE_HEATING is #undef;
127  o verification/cfc_example:  o verification/cfc_example:
128    - turn on pkg/layers with pkg/diagnostics (to test pkg/layers and to illustrate    - turn on pkg/layers with pkg/diagnostics (to test pkg/layers and to
129      how to use it).      illustrate how to use it).
130  o verification/offline_exf_seaice:  o verification/offline_exf_seaice:
131    in input: add SEAICErestoreUnderIce=.TRUE., and reduce size of forcing    in input: add SEAICErestoreUnderIce=.TRUE., and reduce size of forcing
132      files to only one time-record.      files to only one time-record.
133    in input.thermo: narrow the 2 bands of low initial ice-concentration near    in input.thermo: narrow the 2 bands of low initial ice-concentration near
134      N & S edges; switch on SEAICErestoreUnderIce=T; use pkg/diagnostics for      N & S edges; switch on SEAICErestoreUnderIce=T; use pkg/diagnostics for
135      snapshot output; same mcPheePiston value as in thsice with surf-current=0.2m/s.      snapshot output; same mcPheePiston value as in thsice with
136        surf-current=0.2m/s.
137    in input.thsice: similar set-up as in input.thermo (channel with same seaice    in input.thsice: similar set-up as in input.thermo (channel with same seaice
138      initial conditions).      initial conditions).
139  o verification/cfc_example:  o verification/cfc_example:
Legend:
Removed from v.1.1807  
changed lines
  Added in v.1.1827
  ViewVC Help
Powered by ViewVC 1.1.22