/[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.2161 by mlosch, Mon Jun 15 15:48:41 2015 UTC revision 1.2184 by dimitri, Fri Oct 2 01:02:41 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/shelfice/shelfice_thermodynamics.F
8      - redefining uLoc/vLoc to u/v, instead of 2u/v
9    o pkg/seaice/seaice_evp.F: introduce "adaptive EVP" by Madlen Kimmritz
10      - this is an intermediate state of a truely converging EVP that has
11        the potential of outperforming implicit VP methods, I do not recommend
12        using it yet.
13      - for now, "adaptive EVP" is turned on by setting SEAICE_EVP_USE_ADAPTATION
14        and all parameters are hardwired
15      - for convergence it still requires very smooth regularisations, that means,
16        if SEAICE_DELTA_SMOOTHREG is defined, deltaCreg=deltaC+deltaMin, which is
17        different from the VP case with deltaCreg=sqrt(deltaC**2+deltaMin**2); also
18        specific averaging appears to be required, still subject to tuning
19      - still needs cleaning and I need to turn parameter into runtime parameters
20      - add code to compute and print residuals of iteration, compile with
21        defining ALLOW_SEAICE_EVP_RESIDUAL to enable it
22    o pkg/seaice/seaice_evp.F: fix a bug (factor of four in front of e12Csq)
23      that entered with r1.54 on Sep04, 2015, -> update results
24    
25    checkpoint65o (2015/09/14)
26    o pkg/mom_vecinv & mom_common:
27      - add deep-atmosphere geometry factor in mom_calc_visc.F + in routines
28        called from there.
29      - add anelastic and deep-atmosphere geometry factor in pkg/mom_vecinv ; this
30        allows to use Vector-Invariant form in deep atmos and anelastic formulation
31    o pkg/seaice: introduce fast ice parameterisation following Itkin et al. (2015)
32      and Koenig-Beatty+Holland (2012)
33      - two new parameters SEAICE_tensilFac and SEAICE_tensilDepth
34      - global 2D field tensileStrength, computed in seaice_calc_ice_strength
35      - requires an extra input argument for S/R seaice_calc_viscosities
36      - tensileStrength's depth dependence is different from Itkin et al. (2015)
37        (to be changed and played with, for now exponential)
38    o pkg/seaice: introduce new parameter SEAICEscaleSurfStress (default = .FALSE.)
39      - if TRUE scale ice-ocean and ice-atmosphere stress acting on sea ice
40        by AREA according to Connelly et al. (2004), JPO.
41      - For EVP and the LSR solver, the implicit part of the drag term is
42        the only non-zero term in the denominator when there is no ice (which
43        prevented running the model with zero ice-ocean drag). If
44        SEAICEscaleSurfStress = .TRUE., this results in a division by zero
45        (or zero main diagonals BU/BV) which need to be caught. A practical
46        consequence is that for open water, the momentum equation reduce
47        to trivially 0 = 0 + 0 (for EVP). For LSR, BU/BV are reset to 1, if
48        they happen to be zero, often leading to a non-zero solution over
49        open water.
50    o pkg/seaice: revise EVP code (seaice_evp.F) in preparation for a more
51        efficient EVP method.
52      - introduce area-weighted averages for e12 (analogously to
53        seaice_calc_viscosities.F); this totally changes lab_sea.hb87 (2 digits
54        remain, but also changing the order of terms in the previous average
55        changed the results as much) -> update experiment
56      - add SEAICE_DELTA_SMOOTHREG code
57      - disentangle computation of zetaC/deltaC and zetaZ/deltaZ, also for
58        TEM-code
59      - adjust loop bounds so that only required fields are computed
60      - add new averaging code for zetaZ/deltaZ (again similar to
61        seaice_calc_viscosities.F), inactive by default for now (will change
62        results again), but will replace old code eventually
63    o eesupp/src:
64      - add an other alternative for multi-proc GLOBAL_SUM (in global_sum_tile.F)
65        to always cumulate tile local-sum in the same order
66        (option #define GLOBAL_SUM_ORDER_TILES in CPP_EEOPTIONS.h).
67        Still uses MPI_Allreduce (like standard code) but on array of tiles
68        and should be faster than GLOBAL_SUM_SEND_RECV code.
69      - switch almost all verification experiments that use GLOBAL_SUM_SEND_RECV
70        (defined in CPP_EEOPTIONS.h_mpi) to GLOBAL_SUM_ORDER_TILES instead.
71      - add CPP_EEOPTIONS.h_mpi with GLOBAL_SUM_ORDER_TILES defined in 14 more exp.:
72         aim.5l_cs, cheapAML_box, deep_anelastic, dome, fizhi-cs-32x32x40,
73         fizhi-cs-aqualev20, fizhi-gridalt-hs, global_ocean_ebm, matrix_example,
74         natl_box, solid-body.cs-32x32x1, tutorial_advection_in_gyre,
75         tutorial_baroclinic_gyre and tutorial_held_suarez_cs.
76    
77    checkpoint65n (2015/07/29)
78    o pkg/thsice:
79      - add a warning when used with pkg/aim but without updating SST
80        (stepFwd_oceMxL=F & useCoupler=F & tauRelax_MxL<>-1 ).
81    o pkg/offline:
82      - change where Offline-Fields are loaded, from load_fields_driver.F to
83        directly within FORWARD_STEP, just after updating iter number and time;
84      - move setting of GM and KPP diffusivity out of OFFLINE_FIELDS_LOAD into new
85        S/R OFFLINE_GET_DIFFUS that is called towards the end of DO_OCEANIC_PHYS;
86       These changes allow:
87        a) to recover exactly on-line results using deltaToffline=deltaTClock,
88        tested (using offlineIter0= nIter0 -1 and offlineTimeOffset= -deltaT/2)
89        both with Lin & Non-Lin Free-Surf, both with staggerTimeStep & synchronous
90        time-step, both with GM & KPP files and with recomputing these fields;
91        b) no needs to switch on staggerTimeStep anymore;
92        c) simplify rescaling of offline horiz. velocity with Non-Lin Free-Surf.
93      - update output of experiment tutorial_dic_adjoffline (only 3 digits left).
94    o pkg/offline:
95      - change (rarely used) offset-iteration parameter to an offset time (in s):
96        can easily recover previous offset-iter effect by setting offset-time
97        (offlineTimeOffset) to offset-iter x deltaToffline;
98      - do not turn-off Temp & Salinity surface relaxation when using KPP (used
99        in KPP as part of surface buoyancy forcing) in offline run.
100      - fix ptracers KPP non-local term in the case where KPP is re-computed.
101    o OpenAD
102      - Adding preliminary ggl90 test to OpenAD verification.
103        Config. with ggl90 currently diverges after 64 time steps
104        so will need more work to become more stable
105      - Adding preliminary kpp test to OpenAD verification.
106        required small tweakes in kpp_routines.F
107      - Merge and update from Krishna Narayanan's contrib area:
108        * Add one variable for preliminary OpenAD DIVA
109        * Add replacement time-stepping for OpenAD DIVA
110        * genmake2 flag -diva (but only for OpenAD; preliminary)
111        * required modifs for OAD_support
112        successful preliminary test for tutorial_global_oce_biogeo/
113    o pkg/seaice:
114      - initialise deltaC in seaice_lsr, cosmetic changes in seaice_lsr
115      - add cpp-flag SEAICE_DELTA_SMOOTHREG for option of regularising
116        delta with a smooth function in s/r seaice_calc_viscosity
117        (no effect on EVP)
118      - rename local suffix "sqr" into "sq" for more consistent variable names
119    o pkg/seaice:
120        add new parameter SEAICE_deltaMin that is used ***only***
121        for regularizing Delta (and nothing else, like the parameter
122        which is also used for all sorts of things SEAICE_EPS).
123        Defaults to SEAICE_EPS for backward compatibility only.
124    o pkg/seaice: replace erroneously storing e12Csqr on the wrong tape with
125        the wrong key by inialising it before use -> fixes a recomputation
126        problem but does not fix the gradient
127    o pkg/layers:
128      - move resetting of flux arrays to outside the iLa loop
129        (this caused layers_diapycnal to not work with more than one layer coord.)
130      - fix bug in calculation of layers_recip_delta in layers_init_fixed.F
131    o pkg/layers:
132      - remove from common block 3D layers fields (used to diagnose transport
133        in tracer coords) and declare them as local arrays in layers_calc.F
134  o pkg/seaice:  o pkg/seaice:
135    - Rename SEAICE_VECTORIZE_LSR_ZEBRA to SEAICE_LSR_ZEBRA and define it    - Rename SEAICE_VECTORIZE_LSR_ZEBRA to SEAICE_LSR_ZEBRA and define it
136      in lab_sea forward experiment (changes results -> updated)      in lab_sea forward experiment (changes results -> updated)
137  o pkg/seaice: modify seaice_lsr.F in order to improve the gradient    - update adjoint experiment lab_sea after cleaning up seaice_lsr_tridiagu/v
138    computations; for SEAICE_VECTORIZE_LSR  o pkg/seaice: modify seaice_lsr.F in order to improve the gradient
139       computations; for SEAICE_VECTORIZE_LSR
140    - move the loops over which the tridiagonal solvers (seaice_lsr_tridiagu/v)    - move the loops over which the tridiagonal solvers (seaice_lsr_tridiagu/v)
141      sweep to the outside,      sweep to the outside,
142    - remove store directives and add "CADJ loop sequential" directives    - remove store directives and add "CADJ loop sequential" directives
143      in analogy to model/src/solve_tridiagonal.F      in analogy to model/src/solve_tridiagonal.F
144    - replace many "#ifdef SEAICE_VECTORIZE_LSR_ZEBRA" by variable loop    - replace many "#ifdef SEAICE_VECTORIZE_LSR_ZEBRA" by variable loop
145      boundaries to yield more readable code. This has also the charming side      boundaries to yield more readable code. This has also the charming side
146      effect that your can use SEAICE_VECTORIZE_LSR_ZEBRA without      effect that your can use SEAICE_VECTORIZE_LSR_ZEBRA without
147      SEAICE_VECTORIZE_LSR (i.e. adjoint without recomputations in      SEAICE_VECTORIZE_LSR (i.e. adjoint without recomputations in
148      seaice_lsr_tridiagu/v now requires either of these flags, vectorization      seaice_lsr_tridiagu/v now requires either of these flags, vectorization
149      still requires SEAICE_VECTORIZE_LSR with SEAICE_VECTORIZE_LSR_ZEBRA as      still requires SEAICE_VECTORIZE_LSR with SEAICE_VECTORIZE_LSR_ZEBRA as
150      an additional option)      an additional option)
151    The tridiagonal solvers are now completely analoguous to solve_tridiagonal.F     The tridiagonal solvers are now completely analoguous to solve_tridiagonal.F
152    but the adjoint gradients (which are affected by this change) still explode.     but the adjoint gradients (which are affected by this change) still explode.
153  o coupler (atm & ocn _compon_interf):  o coupler (atm & ocn _compon_interf):
154    - avoid using Exch2 IO mapping if not suitable.    - avoid using Exch2 IO mapping if not suitable.
155    
 >>>>>>> 1.2160  
156  checkpoint65m (2015/06/15)  checkpoint65m (2015/06/15)
157  o pkg/layers:  o pkg/layers:
158    - clean-up old parameters setting (from Sep 2012):    - clean-up old parameters setting (from Sep 2012):
Legend:
Removed from v.1.2161  
changed lines
  Added in v.1.2184
  ViewVC Help
Powered by ViewVC 1.1.22