/[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.1296 by dimitri, Thu Sep 16 20:43:00 2010 UTC revision 1.1325 by dimitri, Tue Nov 2 21:23:50 2010 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/generic_advdiff:
8      Change OBCS in multiDim advection (gad_advection.F):
9       Use maskInC to prevent intermediate update of tracer-field in OB regions,
10       thus keeping current-time-step OB values; and remove the call to OBCS_APPLY
11       which were applying future time-step OB values.
12      This way, all explicit calculations use current-time-step OB values
13       (simpler and more consistent).
14      Affects results when set-up uses time-varying OBC and multi-dim advection;
15       update output of experiment seaice_obcs.
16    o pkg/seaice:
17      move OPTIONS (formerly locally defined) to SEAICE_OPTIONS.h:
18       this includes SEAICE_GROWTH_LEGACY (+ SEAICE_DO_OPEN_WATER_GROWTH &
19         SEAICE_OCN_MELT_ACT_ON_AREA ) and SEAICE_SOLVE4TEMP_LEGACY.
20      Default (pkg/seaice/SEAICE_OPTIONS.h) is to define SEAICE_SOLVE4TEMP_LEGACY
21       but leave SEAICE_GROWTH_LEGACY undef ; For now, keep verification exp
22       unchanged (SEAICE_GROWTH_LEGACY defined in */code*/SEAICE_OPTIONS.h)
23      For OBCS, close southern and western boundaries for sea ice solver
24       based on myYGlobalLo and myXGlobalLo instead of OB_Js and OB_Iw.
25    o pkg aim (2010/10/26):
26      - Introduce ALLOW_INSOLATION, moving toward a more general incoming solar radiation
27        formulation with obliquity, eccentricity, ... as run-time parameters,
28      - Introduce diagnostics for top-of-the-atmaophere incoming solar radiation and
29        upward shortwave radition at the surface (all-sky and clear-sky).
30    o pkg/obcs:
31      - writing sections with globalFiles=T does not work: switch to tiled-files
32        when writing sections for OBCS (change pkg/rw/write_rec.F accordingly).
33      - fix index bug (OB_Ie -> OB_Iw_) in obcs_apply_uvice.F
34      - apply OBCS_SEAICE_COMPUTE_UVICE to both uice and vice at all obcs
35    
36    checkpoint62m (2010/10/21)
37    o pkg/seaice: vectorize seaice_solve4temp.F
38      - move iterations loops outside of i/j-loops
39      - requires 2D arrays for TB and effConduct
40      - add 2D iceOrNot to replace HICE_ACTUAL.GT.0, as HICE_ACTUAL is
41        modified in this routine
42      - a little bit of clean up, adjusting indentation, etc.
43      - replaced a few TEN**a by exp(a*log(10.)) in new code for better
44        performance and '**TWO' by '**2' (because TWO is a float),
45        replaced ONE, ZERO
46    o eesupp/src:
47       fix flow dir for exch1_z_RX_cube, exch1_uv_RX_cube & exch1_bg_RX_cube
48    o pkg/seaice/seaice_growth.F
49      Revision 1.89 of seaice_growth.F completes the overall october 2010 revision and
50      merging process. The merged routine has two branches, and one activates the 'evolution'
51      branch by setting SEAICE_GROWTH_LEGACY to undef (at the top of the routine for now).
52      In summary, here is what the '10/10 revision' did consist of
53        part 1: seaice_growth.F revisions 1.71 to 1.75
54          - Preliminary seaice_growth.F clean-up, keeping test results strictly unchanged.
55            The emphasis here was on improving readability (e.g. by trying to use meaningful
56            variable names) and modularity (e.g. by adding arrays to store heff etc. updates).
57        part 1.5:
58          - revising the lab_sea experiment to make it a more useful (i.e. stable
59            and representative of an actual model run). This was done by generating a
60            (winter) pickup to start from, switching back the main experiment dynamic
61            solver from evp to lsr, having one of the experiments use the free drift solution.
62        part 2: seaice_growth.F revisions 1.76 to 1.85
63          - Sorting out hard-coded constants that prevented heat and fresh water conservation.
64          - Sorting out the various fields units and sign (set early on and reset at the end)
65            to avoid numerous scalings/re-scalings. By v1.85 the basic unit for ice, heat and
66            fresh water stock is 'effective ice meters', and >0 implies more ice/snow.
67          - Using tendency terms arrays (d_HEFF* etc.) to alleviate remaining code locks,
68            and simplifying logical expressions to further improve readibility.
69        part 3: seaice_growth.F revisions 1.86 to 1.89
70          - Re-organize the blocks of code more logically.
71          - Merging the '_if.F' pieces that add something new.
72    o pkg/thsice:
73      - vectorize thsice_solve4temp and thsice_get_exf.
74        Problem: thsice_solve4temp has an unvectorizable statement that
75        might also be problematic for TAF
76      - thsice_get_exf: add the ALLOW_BULK_LARGEYEAGER04 CPP-flag. This
77        does not change the verifcation experiment, but will change all
78        experiments that have this flag undefined (not tested in verification);
79        a little reorganization to save a few 2D fields and recomputations
80      - 1st step of revisiting thsice adjoint (major recomp. sorted out)
81      - 2nd step:  fix some store dirs
82        o seems to work somewhat, i.e. using ifort but not g77
83        o likely cause are incomplete init. or recomputation
84        o start verification/global_ocean.cs32x15/input_ad.thsice
85    o pkg/obcs: add support for exch2
86      - the position of the boundary and prescribed values are specified using
87        a global domain with exch2 facets stacked
88        - in x for N,S boundaries (like W2_mapIO=-1)
89        - in y for E,W boundaries (so E,W boundaries do not overlap)
90    o pkg/mdsio: add support for writing global sections with exch2
91        file layout is as for obcs described above:
92        - for x-z sections, facets are stacked in x (like W2_mapIO=-1)
93        - for y-z sections, facets are stacked in y
94    o pkg/thsice:
95       - change a few routines for better vectorization: thsice_albedo.F,
96         thsice_step_temp.F, and in thsice_calc_thickn.F. Results are
97         not affected, but in particular the latter routine has changed
98         dramatically (not necessarily easier to read now).
99       - re-arrange calculation in THSICE_SOLVE4TEMP: move outside of i,j loops
100       iteration loop and call to S/R THSICE_GET_BULKF & THSICE_GET_EXF.
101    o verification/lab_sea:
102      - overhaul of the forward model experiments. The main motivation
103        for this revision is that lab_sea has proven quite fragile (when
104        changing the sea-ice code, and switching platforms) to a point
105        when failures were very hard to attribute to specific code problems.
106      - the suspected reason for this fragility is noise amplification by evp so,
107        to start with, we decided to replace evp with lsr in the main experiment.
108      - the input.lsr experiment thus became obsolete and was removed.
109      - it was replaced with the input.fd experiment that uses an
110        analytical "free drift" solution as the ice velocity field.
111      - Taken together, input (lsr), input.fd (free drift) and input.hb87 (evp)
112        will hopefully allow us to discriminate problems due to dynamic solvers
113        from problems due to the thermodynamic code. input.salt_plume (lsr) allows
114        us to test additional options, along with the salt plume package.
115      - Finally a pickup was generated (at the end of a one year run) to
116        make the tests more representative of an actual model run. Previously
117        the tests starting point was a 1m thickness of ice everywhere.
118    o pkg/seaice:
119      - introducing S\R seaice_freedrift.F that compute the
120        analytical solution to an approximate ice momentum
121        equation (i.e. the simple stress/coriolis balance).
122      - define SEAICE_ALLOW_FREEDRIFT to trigger the computation
123        of the free drift velocity field, and enable its use.
124        To actually use the free drift solution as the model ice
125        velocity field set SEAICEuseFREEDRIFT to .TRUE. in data.seaice.
126    o move "useOASIS" from PARAMS.h to EEPARAMS.h and read it from "eedata"
127    
128    checkpoint62l (2010/10/04)
129    o pkg/seaice:
130      -fix the USE_ORIGINAL_SBI version (for variable freezing point);
131      -use "celsius2K" in SEAICE_VARIABLE_FREEZING_POINT and in new
132       version (with USE_ORIGINAL_SBI undef);
133    o add hooks for new packages OASIS, the package will follow
134      (hooks exclude the seaice pkg for now)
135    o pkg/autodiff: ALLOW_AUTODIFF_WHTAPEIO
136      Option to use mdsio_read_whalos.F/mdsio_write_whalos.F,
137      rather than mdsio_readvector.F/mdsio_writevector.F, within
138      adread_adwrite.F to read/write adjoint-related tapes to disk.
139    o eesupp and pkg/mdsio: ALLOW_WHIO
140      Option to write/read 2D fields to files including tiles halos.
141      The main purpose is for adjoint related "tape I/O".
142      The secondary purpose is for debugging phases.
143      A word on those rather specifically intended purposes.
144      The code was meant to be stripped down to the bare minimum.
145      Typically, there is a unique and basic mapping to the global
146      buffer (tile by tile, then proc by proc, then level by level), and
147      byteswaps and meta-files are omitted. Indeed the (*_WH*) code is not
148      meant to generate user friendly outputs of various shapes and forms.
149      Rather it is intended, for the knowledgeabale developers,
150      to be largely independent of the many mdsio options, to allow
151      outputs of tiles+halos in debugging phases and in adjoint runs,
152      and to allow some flexibility in tuning (adoint) I/O depending on each
153      file system behavior. With regard to tuning, most of the potential knobs
154      (to leave files open, to allocate the I/O indep. of usesinglecpuio,
155      to specify directories, to specify the typical size of I/O calls, etc.)
156      are not yet included, and will be added if they prove useful.
157    o pkg/seaice:
158       fix uninitialised arrays in USE_ORIGINAL_SBI version of seaice_solve4temp.F
159    o pkg/seaice:
160       new S/R seaice_solve4temp.F: this is the merged version of Ian Fenty
161       code (seaice_budget_ice_if.F) which enable to recover old standard
162       code (seaice_budget_ice.F) when option USE_ORIGINAL_SBI is defined.
163       temporary: add "#define USE_ORIGINAL_SBI" at the top of seaice_solve4temp.F.
164    o model/src (rStar):
165       use simple average (not area weighted) of rStarFac at U,V point when using
166       vectorInvariant and selectKEscheme=1,3 (consistent with mom. vertical advection)
167       update output of exp. fizhi-cs-32x32x40 and fizhi-cs-aqualev20.
168    
169    checkpoint62k (2010/09/17)
170    o adjoint
171       * changed verification/global_ocean.cs32x15 ; adjoint now tests
172         Adams-Bashforth 3 + NLFS + RSTAR
173       * some more modifs required for ECCO-GODAE v4
174  o pkg/seaice  o pkg/seaice
175     add SEAICE_AGE_VOL option to use volume instead of area for     add SEAICE_AGE_VOL option to use volume instead of area for
176     integrating SEAICE_AGE     integrating SEAICE_AGE
# Line 276  o tools/build_options: Line 443  o tools/build_options:
443      but keep the former optfile for few weeks for smooth transition      but keep the former optfile for few weeks for smooth transition
444  o pkg/ecco  o pkg/ecco
445    - allow user defined cost function contributions    - allow user defined cost function contributions
446          CPP: ALLOW_USERCOST_CONTRIBUTION       CPP: ALLOW_USERCOST_CONTRIBUTION
447          USAGE: compute contributions in cost_usercost_all.F       USAGE: compute contributions in cost_usercost_all.F
448    - note: this is identical to what ALLOW_GENCOST_CONTRIBUTION once was.    - note: this is identical to what ALLOW_GENCOST_CONTRIBUTION once was.
449        the two options are complementary.        the two options are complementary.
450  o pkg/smooth  o pkg/smooth
Legend:
Removed from v.1.1296  
changed lines
  Added in v.1.1325
  ViewVC Help
Powered by ViewVC 1.1.22