/[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.1738 by jmc, Mon Sep 17 22:32:32 2012 UTC revision 1.1758 by mlosch, Fri Nov 2 07:53:56 2012 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: simplify seaice_evp by avoiding use of extra fields u/vIceC
8    o pkg/autodiff
9      - Add inAdMode switch for SEAICEuseDYNAMICS:
10        SEAICEuseDYNAMICSinAdMode
11        (extend logic in ADAUTODIFF_INADMODE_SET)
12    o pkg/seaice
13      - Step 1 of merging ice-thickness distribution (ITD) code from
14        MITgcm_contrib/torge/itd/ (author: Torge Martin)
15        The code implements concept laid out in Thorndike et al. (1975);
16        code is enabled via CPP flag SEAICE_ITD
17        nITD (in SEAICE_SIZE.h) defines number of ice thickness categories;
18        code is not complete yet!!
19      - SEAICE_SIZE.h now included in almost all seaice routines
20    o verification (using seaice):
21      - switch back seaice-dynamics (+ lower SOLV_MAX_ITERS) in ADM exp.
22        global_ocean.cs32x15.seaice ; update results.
23      - compile B-grid dynamics in FW test exp 1D_ocean_ice_column ;
24      - compile C-grid dynamics with #define SEAICE_ALLOW_JFNK in FW test exp
25        offline_exf_seaice (previously compiling B-grid).
26    o pkg/exf:
27      - fix some STOP conditions in exf_check.F (from Aug 28, modifs did not allow
28        to specify neither just 1 comp of wind or wind-stress, nor no surface momentum
29        forcing): just replace CPP option ALLOW_ATM_WIND with run-time param useAtmWind.
30    o pkg/layers (adjust changes from Sept 19):
31      - do snap-shot output and fill diagnostics from inside iLa loop (layers_calc.F)
32        this allows to remove "layers_maxNum" dimension from all layers arrays.
33    o pkg/layers (adjust changes from Sept 19):
34      - remove ALLOW_LAYERS_OUTPUT ; fix snap-shot and time-ave output ;
35      - snap-shot extended to all type of layers (not just the 1rst type);
36      - remove layers variables when corresponding CPP option is undef.
37    o pkg/layers:
38      - rename pkg/layers variable and add new output (layer exist & un-weighted
39        velocity in layer)
40    o pkg/seaice:
41      - simplify (?) ocean-drag coefficent computations by using
42        the same routine in all three solvers (saves some code lines)
43      - remove global fields DRAGA and DRAGS from C-grid code and
44        declare them locally in lsr.F (the only place where they are
45        still used)
46      - update AD output of lab_sea.evp (very sensitive anyway).
47    o pkg/seaice: add Jacobian-free Newton Krylov solver for ice dynamics
48      (Lemieux et al. 2010, 2012). Turn on with the usual SEAICE_ALLOW_JFNK
49      and SEAICEuseJFNK=.true. A few things to note:
50      - involves a FGMRES, that is parallel but very inefficient,
51        because of many scalar products (global sums) in the Gram-Schmidt
52        orthogonalization method. Needs to be replaced. Original code
53        from J.-F. Lemieux
54      - preconditioner for FGMRES is a stripped down version of LSOR
55      - SOLV_MAX_ITERS is now used for LSOR as a solver and as a
56        preconditioner. For the preconditioner, the default is set to 10.
57      - Jacobian times vector is computed with finite differences
58      - marginally related to JFNK-solver: add a smooth regularization for
59        zeta with tanh that will improve convergence (Lemieux and Tremblay,
60        2009) for any implicit solver (not implemented for EVP because EVP
61        does not need this regularization in the first place). Define
62        SEAICE_ZETA_SMOOTHREG in SEAICE_OPTIONS.h to use.
63    o eesupp/src:
64      - check in genmake2 that intrinsic S/R flush is available.
65      - add S/R to flush IO unit (if intrinsic S/R flush is available)
66      - flush standard & error msg unit before MPI termination
67      - if debugMode, flush IO unit after printing msg
68    o pkg/ctrl:
69      - ctrl_unpack.F, ALLOW_GENARR3D_CONTROL code: fix missing argument in call
70        to ctrl_set_unpack_xyz.
71    
72    checkpoint64 (2012/10/09)
73    o pkg/shelfice
74      - in mom_fluxform.F, mom_vecinv.F:
75        CALL SHELFICE_V_DRAG(bi,bj,k,vFld,KE,KappaRU,vF,myThid)
76        Bug fix KappaRU -> KappaRV
77        does not change verification/isomip but may change other setups
78        (presumably constant KappaRU=KappaRV)
79      - fix diagnostics SHIU/Vdrag (thanks to Tom Millgate)
80    o pkg/obcs:
81      - fixed some indexing bugs in obcs_seaice_sponge.F
82      - added obcs relaxation subroutines for HSALT and HSNOW
83      - changes verification/seaice_obcs/results/obcs_seaice_sponge.F
84    o pkg/seaice:
85      - add "multicategory" parameterization for snow in analogy to ice,
86        default is off for now.
87    o pkg/obcs:
88      - added obcs_seaice_sponge.F for obcs relaxation of seaice variables;
89        see verification/seaice_obcs/input.seaiceSponge for example application.
90      - remove local commom bloc "OBCS_GLOBAL" in obcs_readparms.F and do all the OB
91        index setting by master-thread. This fix multi-threads Pb in previous modifs.
92    o pkg/layers :
93      - add diagnostics to pkg/layers
94      - move the meat of pkg/layers from layers_calc.F to layers_fluxcalc.F
95        (from R.Tulloch). Layers_calc.F now calls layers_fluxcalc.F accordingly.
96    o pkg/generic_advdiff/gad_implicit_r.F :
97      - fix ADVr diagnostics in case of both implicit vertical advection and diffusion.
98    o pkg/seaice/seaice_growth.F:
99      - fix SItflux diagnostics in virtual salt flux + balanceEmpmr case.
100    o pkg/obcs:
101      - change default setting of OBCS_UVICE_OLD from defined to undef;
102      - use new parameter OB_indexNone for null index value (instead of hard-coded 0)
103        in all pkg/obcs code and in other pkgs; remove k loop in obcs_check_depths.F;
104      - pass OB_indexNone as (new) argument to S/R CTRL_MASK_SET_[X,Y]Z;
105      - change OB null index value (OB_indexNone) to -99 ; this fix the problem
106        with Southern & Western OB at a tile edge (at least in a simple set-up).
107  o pkg/obcs:  o pkg/obcs:
108    - change settings of OB index in case OB coincides with tile edge (shared by    - change settings of OB index in case OB coincides with tile edge (shared by
109      2 tiles), so that common OB location is registered by both tiles (instead of      2 tiles), so that common OB location is registered by both tiles (instead of
110      only 1). However, for Southern and Western OB, index value = 0 is not      only 1). However, for Southern and Western OB, index value = 0 is not
111      recognized as valid OB index (=default null value) and needs to be fixed.      recognized as valid OB index (=default null value) and needs to be fixed.
112    - modify accordingly OBCS_BALANCE_FLOW and OBCS_MONITOR (avoid double counting)    - modify accordingly OBCS_BALANCE_FLOW and OBCS_MONITOR (avoid double counting)
113    - new parameter "OB_indexNone" for null index value (for now, still zero)    - new parameter "OB_indexNone" for null index value (for now, still zero)
114      to replace hard-coded null index value of zero.      to replace hard-coded null index value of zero.
115  o pkg/exch2:  o pkg/exch2:
116    - improve W2_USE_E2_SAFEMODE (faster) and W2_E2_DEBUG_ON:    - improve W2_USE_E2_SAFEMODE (faster) and W2_E2_DEBUG_ON:
Legend:
Removed from v.1.1738  
changed lines
  Added in v.1.1758
  ViewVC Help
Powered by ViewVC 1.1.22