/[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.2173 by jmc, Sat Jul 18 22:15:07 2015 UTC revision 1.2191 by mlosch, Thu Oct 22 15:51:59 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: fix store directives and key computations, since
8      there is no adjoint test for EVP, we don't know what it does to the adjoint
9      solution, but at least there are no more recomputation warnings left
10    o pkg/seaice/seaice_evp.F: slight modification of averaging so that EVP stands
11      a chance of giving the same results as a fully converged VP (JFNK) solver
12    o tools/adjoint_options/adjoint_diva:
13      - fix bug (typo) in Tang-Linear FLAGS setting for MPI + DIVA
14    o pkg/streamice:
15      - spotted bug in streamice_driving_stress using dyG where should have been dxG
16    o pkg/exf:
17      - remove STOP statement that prevented to use wind or (wind-)stress
18        interpolation for constant-in-time forcing fields (wind/stress period = 0).
19    o pkg/seaice:
20      - replace SEAICE_EVP_USE_ADAPTATION by run time parameters, requires
21        5 additionals 2D fields and 6 new store directives
22      - change logic in seaice_readparms.F: EVP code can now be turned on in
23        various ways, setting SEAICE_deltaTev is no longer required,
24        alternatively set SEAICE_evpAlpha/Beta, SEAICEuseEVPstar, SEAICEuseEVPrev,
25        or SEAICEaEVPcoeff
26    o pkg/shelfice/shelfice_thermodynamics.F
27      - redefining uLoc/vLoc to u/v, instead of 2u/v
28    o pkg/seaice/seaice_evp.F: introduce "adaptive EVP" by Madlen Kimmritz
29      - this is an intermediate state of a truely converging EVP that has
30        the potential of outperforming implicit VP methods, I do not recommend
31        using it yet.
32      - for now, "adaptive EVP" is turned on by setting SEAICE_EVP_USE_ADAPTATION
33        and all parameters are hardwired
34      - for convergence it still requires very smooth regularisations, that means,
35        if SEAICE_DELTA_SMOOTHREG is defined, deltaCreg=deltaC+deltaMin, which is
36        different from the VP case with deltaCreg=sqrt(deltaC**2+deltaMin**2); also
37        specific averaging appears to be required, still subject to tuning
38      - still needs cleaning and I need to turn parameter into runtime parameters
39      - add code to compute and print residuals of iteration, compile with
40        defining ALLOW_SEAICE_EVP_RESIDUAL to enable it
41    o pkg/seaice/seaice_evp.F: fix a bug (factor of four in front of e12Csq)
42      that entered with r1.54 on Sep04, 2015, -> update results
43    
44    checkpoint65o (2015/09/14)
45    o pkg/mom_vecinv & mom_common:
46      - add deep-atmosphere geometry factor in mom_calc_visc.F + in routines
47        called from there.
48      - add anelastic and deep-atmosphere geometry factor in pkg/mom_vecinv ; this
49        allows to use Vector-Invariant form in deep atmos and anelastic formulation
50    o pkg/seaice: introduce fast ice parameterisation following Itkin et al. (2015)
51      and Koenig-Beatty+Holland (2012)
52      - two new parameters SEAICE_tensilFac and SEAICE_tensilDepth
53      - global 2D field tensileStrength, computed in seaice_calc_ice_strength
54      - requires an extra input argument for S/R seaice_calc_viscosities
55      - tensileStrength's depth dependence is different from Itkin et al. (2015)
56        (to be changed and played with, for now exponential)
57    o pkg/seaice: introduce new parameter SEAICEscaleSurfStress (default = .FALSE.)
58      - if TRUE scale ice-ocean and ice-atmosphere stress acting on sea ice
59        by AREA according to Connelly et al. (2004), JPO.
60      - For EVP and the LSR solver, the implicit part of the drag term is
61        the only non-zero term in the denominator when there is no ice (which
62        prevented running the model with zero ice-ocean drag). If
63        SEAICEscaleSurfStress = .TRUE., this results in a division by zero
64        (or zero main diagonals BU/BV) which need to be caught. A practical
65        consequence is that for open water, the momentum equation reduce
66        to trivially 0 = 0 + 0 (for EVP). For LSR, BU/BV are reset to 1, if
67        they happen to be zero, often leading to a non-zero solution over
68        open water.
69    o pkg/seaice: revise EVP code (seaice_evp.F) in preparation for a more
70        efficient EVP method.
71      - introduce area-weighted averages for e12 (analogously to
72        seaice_calc_viscosities.F); this totally changes lab_sea.hb87 (2 digits
73        remain, but also changing the order of terms in the previous average
74        changed the results as much) -> update experiment
75      - add SEAICE_DELTA_SMOOTHREG code
76      - disentangle computation of zetaC/deltaC and zetaZ/deltaZ, also for
77        TEM-code
78      - adjust loop bounds so that only required fields are computed
79      - add new averaging code for zetaZ/deltaZ (again similar to
80        seaice_calc_viscosities.F), inactive by default for now (will change
81        results again), but will replace old code eventually
82    o eesupp/src:
83      - add an other alternative for multi-proc GLOBAL_SUM (in global_sum_tile.F)
84        to always cumulate tile local-sum in the same order
85        (option #define GLOBAL_SUM_ORDER_TILES in CPP_EEOPTIONS.h).
86        Still uses MPI_Allreduce (like standard code) but on array of tiles
87        and should be faster than GLOBAL_SUM_SEND_RECV code.
88      - switch almost all verification experiments that use GLOBAL_SUM_SEND_RECV
89        (defined in CPP_EEOPTIONS.h_mpi) to GLOBAL_SUM_ORDER_TILES instead.
90      - add CPP_EEOPTIONS.h_mpi with GLOBAL_SUM_ORDER_TILES defined in 14 more exp.:
91         aim.5l_cs, cheapAML_box, deep_anelastic, dome, fizhi-cs-32x32x40,
92         fizhi-cs-aqualev20, fizhi-gridalt-hs, global_ocean_ebm, matrix_example,
93         natl_box, solid-body.cs-32x32x1, tutorial_advection_in_gyre,
94         tutorial_baroclinic_gyre and tutorial_held_suarez_cs.
95    
96    checkpoint65n (2015/07/29)
97    o pkg/thsice:
98      - add a warning when used with pkg/aim but without updating SST
99        (stepFwd_oceMxL=F & useCoupler=F & tauRelax_MxL<>-1 ).
100  o pkg/offline:  o pkg/offline:
101    - change where Offline-Fields are loaded, from load_fields_driver.F to    - change where Offline-Fields are loaded, from load_fields_driver.F to
102      directly within FORWARD_STEP, just after updating iter number and time;      directly within FORWARD_STEP, just after updating iter number and time;
# Line 16  o pkg/offline: Line 109  o pkg/offline:
109      time-step, both with GM & KPP files and with recomputing these fields;      time-step, both with GM & KPP files and with recomputing these fields;
110      b) no needs to switch on staggerTimeStep anymore;      b) no needs to switch on staggerTimeStep anymore;
111      c) simplify rescaling of offline horiz. velocity with Non-Lin Free-Surf.      c) simplify rescaling of offline horiz. velocity with Non-Lin Free-Surf.
112    - update output of experiments tutorial_cfc_offline (11 to 13 digits left)    - update output of experiment tutorial_dic_adjoffline (only 3 digits left).
     and tutorial_dic_adjoffline (only 3 digits left).  
113  o pkg/offline:  o pkg/offline:
114    - change (rarely used) offset-iteration parameter to an offset time (in s):    - change (rarely used) offset-iteration parameter to an offset time (in s):
115      can easily recover previous offset-iter effect by setting offset-time      can easily recover previous offset-iter effect by setting offset-time
Legend:
Removed from v.1.2173  
changed lines
  Added in v.1.2191
  ViewVC Help
Powered by ViewVC 1.1.22