/[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.1999 by jmc, Sat May 24 13:26:28 2014 UTC revision 1.2032 by jmc, Sun Jul 27 18:04:14 2014 UTC
# Line 4  $Name$ Line 4  $Name$
4      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
5      ==============================      ==============================
6    
7    checkpoint65a (2014/07/27)
8    o model/src:
9      - Non-Lin Free-Surf in p/z coords: allows hFacW,S to be larger than
10        surrounding hFacC=1 @ edge of a step with different kSurfC on either side
11        (topo in p-coords, ice-shelf in z-coords); update output of experiments
12        tutorial_global_oce_in_p and isomip (both AD and OpenAD output).
13    o model/src:
14      - collect momentum forcing tendency into local 2-D array (passed as argument
15        to apply_forcing S/R); this simplifies mom. Forcing diagnostics and CD_CODE
16        parts in timestep.F (without CD_CODE_NO_AB_MOMENTUM defined = the default);
17        affects results at machine truncation level (e.g., if using OBCS_SPONGE).
18      - collect T & S forcing tendency into local 2-D array (passed as argument to
19        apply_forcing S/R); add 3-D diagnostics for temp & salt forcing tendency.
20        affects results at machine truncation level (e.g., ocean: if using
21        SHORTWAVE_HEATING and atmos:, if using FrictionHeating or (atm_Rq<>0)
22        moist correction due to vertical advection).
23      - update output of experiments: fizhi-cs-32x32x40 (13 digits),
24        global_ocean.cs32x15 (fwd: icedyn, seaice, thsice, all 3 @ 11d; adm & tml:
25        seaice_dynmix & thsice, fwd grd @ 3 & 8), hs94.1x64x5 (ad grd @ 13),
26        obcs_ctrl (ad & fwd grd @ 12 & 7).
27    o pkg/ptracers & diagnostics:
28      - use ptracers (short) name (instead of long one) to make ptracers diagnostics
29        title (to avoid some unexpected truncation); sort out ptracers diags units.
30      - collect ptracer forcing tendency into local 2-D array (passed as argument
31        to apply_forcing S/R); add 3-D diagnostics for ptracers forcing tendency.
32    
33    checkpoint65 (2014/07/18)
34    o pkg/fizhi:
35      - add parenthesis around fizhi tendencies in fizhi_tendency_apply.F: this
36        affects results at machine truncation level. Update output of experiments
37        fizhi-cs-32x32x40 (12 digits for cg2d) & fizhi-cs-aqualev20 (13 digits).
38    o pkg/seaice: more clean up
39      - rename TAF common block tapes to be more consistent and remove
40        unnecessary duplicates
41      - add store directives to remove more recomputations, but use smaller
42        tapes (and adjust keys) when SEAICE_LSR_ADJOINT_ITER is undefined.
43        The idea behind this flag is to reduce the memory footprint when the
44        adjoint of the LSR is not required. As of now, defining the flag
45        SEAICE_LSR_ADJOINT_ITER increases 4 common block tapes for 2D field by
46        a factor of SOLV_MAX_FIXED (=500 by default in SEAICE_SIZE.h).
47        For experiment offline_exf_seaice, this means on a linux box with gfortran
48        nearly a factor of 10 (but keep in mind there is mostly seaice code in
49        this experiment):
50        > size mitgcmuv_ad
51        text    data        bss     dec        hex  
52        5100659 1216  220458952  225560827  d71c8fb
53        5339425 1200 1868783272 1874123897 6fb4dc79 (SEAICE_LSR_ADJOINT_ITER)
54        but only a moderate increase of 36% for global_ocean.cs32x15:
55        text    data        bss     dec        hex
56        6695027 1648 1014578664 1021275339 3cdf6ccb
57        6695203 1648 1382505320 1389202171 52cd8afb (SEAICE_LSR_ADJOINT_ITER)
58      - note: there are still recomputations in global_ocean.cs32x15 and
59        offline_exf_seaice and the gradient agreement changes (gets a
60        little worse) for global_ocean.cs32x15 (-> update output_adm.seaice.txt)
61      - in general, the AD gradients still are wrong
62    o OpenAD
63      - clean up code_oad/ directories somewhat by moving files
64        dontCompile, dontTransform, keepOriginal to OAD_support/
65    o model/src:
66      - using MOM fluxform and Non-Lin Free-Surf with momDissip_In_AB=F
67        (not currently tested): add dissipation terms before NLFS rescaling of
68        tendencies (previously done after).
69    o pkg/diagnostics:
70      - write stats-diags in ascii file with more digits (same as in monitor).
71    o pkg/down_slope:
72      - drop the bi,bj dimension of argument recip_hFac in S/R DWNSLP_APPLY; with
73        NLFS (not tested), use new recip_hFac value in order to conserve tracer.
74    o model/src:
75      - new file "apply_forcing.F" containing all the code previously in
76        external_forcing.F, but with new argument list: pass, as new argument,
77        the current level tendency array to update (instead of a direct update
78        of the common bloc array). Change the corresponding calls.
79      - add CPP option USE_OLD_EXTERNAL_FORCING (default is undef) to use
80        EXTERNAL_FORCING_U,V,T,S subroutines (for backward compatibility).
81    o pkg/seaice/seaice_lsr, a little bit of clean up
82      - rearrange some code blocks, so that now all calculations for
83        setting up the right-hand-sides done within one single block
84        (I find that easier to read and understand)
85      - rearrange freedrift-code to reduce the number of ifdefs
86      - adjust some indentations and comments
87    
88    checkpoint64z (2014/07/10)
89    o model/src and several pkgs:
90      - change argument list of all S/R called from external_forcing.F and
91        called from ptracers_forcing.F (renamed to ptracers_apply_forcing.F)
92        to pass, as argument, the current level tendency array to update
93        (instead of a direct update of the common bloc array).
94      - change also (in the same way) argument list of PTRACERS_APPLY_FORCING;
95      - rename taueddy_external_forcing.f to taueddy_tendency_apply.F
96    o pkg/seaice/seaice_lsr.F
97      - always initialize local arrays in seaice_lsr_tridiagu/v
98      - add more local tapes to seaice_lsr_tridiagu/v
99      - move storing u/vice(:,:,bi,bj) out of seaice_lsr_tridiagu/v, but
100        keep it in bi/bj loop; this seems to be essential
101      - gradients are still wrong but only by 4-7 orders of magnitude
102        instead of 15 orders of magnitude, lab_sea improves a little
103    o pkg/seaice
104      - (pkg/autodiff) remove storing/restoring of seaiceMaskU/V that
105        generates unnecessary code in seaice_lsr_ad
106      - add more store directives in seaice_lsr to avoid unnecessary code
107        in seaice_lsr_ad
108      - exclude some residual computations in seaice_lsr when
109        ALLOW_AUTODIFF_TAMC is defined in order to avoid unnecessary code
110        in seaice_lsr_ad
111      - add more storing, a new tape (model/src/the_main_loop.F) and
112        another key for SEAICE_LSR_VECTORIZE_ZEBRA (the local keys ilockey
113        and klockey should be checked by a specialist)
114      - modify ( define SEAICE_VECTORIZE_LSR and SEAICE_LSR_ADJOINT_ITER )
115        and update lab_sea after previous changes in pkg/seaice below
116    o pkg/autodiff:
117      - introduce useSmoothCorrel2DinAdMode, SEAICEapproxLevInAd, and
118        viscFacInAd : run time parameters allowing for adjoint switches
119        note : logic for useSmoothCorrel2DinAdMode is in ctrl_readparms.F
120    o pkg/ecco:
121      - introduce cost_gencost_boxmean.F that computes a box mean
122        temperature (for use in adjoint sensitivity studies).
123    o model/src/set_defaults.F : initialize viscFacAdj to 1.
124    o pkg/mom_common/mom_calc_visc.F:
125      - ifdef ALLOW_AUTODIFF then include viscFacAdj factor.
126    o pkg/seaice/seaice_summary.F:
127      - statement for SEAICEadjMODE now in autodiff_readparms.F
128    o pkg/seaice:
129      - remove some store directives and add new ones in an effort to get
130        finally rid off the recomputation warnings. This is successful
131        when SEAICE_VECTORIZE_LSR ***and*** SEAICE_LSR_ADJOINT_ITER are
132        defined (i.e. no extensive recomputation warnings left), but there
133        are still too many recomputations on the solver iteration level
134      - define SEAICE_VECTORISZE_LSR in offline_exf_seaice/code_ad/SEAICE_OPTIONS.h
135      - move definition of MPSEUDOTIMESTEPS into TAMC-block in
136        SEAICE_SIZE.h
137      - add a store directive in S/R seaice_calc_viscosities
138    o pkg/seaice (ITD code):
139      - fix definition of Hlimit for multi-threaded case
140      - use double precision arithmetic to compute Hlimit: this affects results;
141        update output of experiment seaice_itd.
142    o replace content of verification/seaice_itd
143      - instead of a 1D column use the configuration of offline_exf_seaice
144        to test the ice thickness distribution code with ridging
145      - this is a first shot at this experiment
146      - change of some defaults of ridging parameters SEAICEhStar (from
147        100 to 25) and SEAICEmuRidging (from 5 to 3) in seaice_readparms.F
148        and lot of new comments.
149    
150    checkpoint64y (2014/06/20)
151    o pkg/seaice/seaice_reg_ridge.F:
152      - modify the regularization of very small volumes of ice; only check
153        total volume HEFF, and do not apply correction to individual ice
154        classes because that leads to funny patterns in the heat and fw-fluxes
155        into the ocean. Adjust seaice_itd/results/output.txt
156      - replace nITD by SEAICE_multDim in loop-ranges
157    o pkg/exf:
158      - add code to process the case field period = -12. with UV interpolation.
159      - move debug print (formerly active when exf_verbose=T) from EXF_getFFieldRec
160        into S/R that calls it (i.e., exf_set_gen.F, exf_set_uv.F & obcs_exf_load.F)
161    o verification/global_with_exf:
162      - switch all fields period to -12. (was previously untested) in input/data.exf
163        (but keep secondary test ".yearly" unchanged); update the reference output.
164    o pkg/cal:
165      - new calendar type "noLeapYear" for 365 days calendar without any leap year.
166    o pkg/exf:
167      - new integer parameter "exf_debugLev" (default set to main model debugLevel)
168        to replace logical switch "exf_verbose"
169      - report in STDOUT when loading a new record (if exf_debugLev >= 3)
170    o pkg/sbo:
171      - update from Katherine Quinn: evaluate mass directly from density and SSH
172        (budget using Finite Volume method) instead of from bottom pressure.
173      - allow testreport to collect few selected SBO output ; tested in exp.
174        global_ocean.90x40x15 (after updating the results, to include SBO output).
175      - fix some multi-threading issues (in output part).
176    o model/src & most packages:
177      - always call {PKG}_READPARMS even if corresponding use{PKG}=F;
178      - in {PKG}_READPARMS, in case use{PKG}=F: add a call to new S/R
179        PACKAGES_UNUSED_MSG to print a weak warning when parameter file
180        "data.{pkg}" exist ; and leave the S/R (return) just after.
181      - more standard LAND pkg interface (land_init_fixed.F, land_check.F).
182        + call LAND_READPARMS from packages_readparms.F
183    o pkg/seaice:
184      - remove code that resets useThSice in S/R seaice_readparms and
185        replace it by a warning in S/R seaice_check
186      - remove parameter MULTDIM, replace with nITD and runtime
187        parameter SEAICE_multDim (still hard-wired to nITD if
188        SEAICE_ITD is defined)
189    
190  checkpoint64x (2014/05/24)  checkpoint64x (2014/05/24)
191  o verification/lab_sea:  o verification/lab_sea:
192    - update adjoint and TLM output following changes to pkg/salt_plume    - update adjoint and TLM output following changes to pkg/salt_plume
Legend:
Removed from v.1.1999  
changed lines
  Added in v.1.2032
  ViewVC Help
Powered by ViewVC 1.1.22