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