/[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.1303 by gforget, Fri Sep 24 23:25:45 2010 UTC revision 1.1319 by jmc, Sun Oct 24 21:35:13 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/obcs:
8       writing sections with globalFiles=T does not work: add a STOP for this case.
9    
10    checkpoint62m (2010/10/21)
11    o pkg/seaice: vectorize seaice_solve4temp.F
12      - move iterations loops outside of i/j-loops
13      - requires 2D arrays for TB and effConduct
14      - add 2D iceOrNot to replace HICE_ACTUAL.GT.0, as HICE_ACTUAL is
15        modified in this routine
16      - a little bit of clean up, adjusting indentation, etc.
17      - replaced a few TEN**a by exp(a*log(10.)) in new code for better
18        performance and '**TWO' by '**2' (because TWO is a float),
19        replaced ONE, ZERO
20    o eesupp/src:
21       fix flow dir for exch1_z_RX_cube, exch1_uv_RX_cube & exch1_bg_RX_cube
22    o pkg/seaice/seaice_growth.F
23      Revision 1.89 of seaice_growth.F completes the overall october 2010 revision and
24      merging process. The merged routine has two branches, and one activates the 'evolution'
25      branch by setting SEAICE_GROWTH_LEGACY to undef (at the top of the routine for now).
26      In summary, here is what the '10/10 revision' did consist of
27        part 1: seaice_growth.F revisions 1.71 to 1.75
28          - Preliminary seaice_growth.F clean-up, keeping test results strictly unchanged.
29            The emphasis here was on improving readability (e.g. by trying to use meaningful
30            variable names) and modularity (e.g. by adding arrays to store heff etc. updates).
31        part 1.5:
32          - revising the lab_sea experiment to make it a more useful (i.e. stable
33            and representative of an actual model run). This was done by generating a
34            (winter) pickup to start from, switching back the main experiment dynamic
35            solver from evp to lsr, having one of the experiments use the free drift solution.
36        part 2: seaice_growth.F revisions 1.76 to 1.85
37          - Sorting out hard-coded constants that prevented heat and fresh water conservation.
38          - Sorting out the various fields units and sign (set early on and reset at the end)
39            to avoid numerous scalings/re-scalings. By v1.85 the basic unit for ice, heat and
40            fresh water stock is 'effective ice meters', and >0 implies more ice/snow.
41          - Using tendency terms arrays (d_HEFF* etc.) to alleviate remaining code locks,
42            and simplifying logical expressions to further improve readibility.
43        part 3: seaice_growth.F revisions 1.86 to 1.89
44          - Re-organize the blocks of code more logically.
45          - Merging the '_if.F' pieces that add something new.
46    o pkg/thsice:
47      - vectorize thsice_solve4temp and thsice_get_exf.
48        Problem: thsice_solve4temp has an unvectorizable statement that
49        might also be problematic for TAF
50      - thsice_get_exf: add the ALLOW_BULK_LARGEYEAGER04 CPP-flag. This
51        does not change the verifcation experiment, but will change all
52        experiments that have this flag undefined (not tested in verification);
53        a little reorganization to save a few 2D fields and recomputations
54      - 1st step of revisiting thsice adjoint (major recomp. sorted out)
55      - 2nd step:  fix some store dirs
56        o seems to work somewhat, i.e. using ifort but not g77
57        o likely cause are incomplete init. or recomputation
58        o start verification/global_ocean.cs32x15/input_ad.thsice
59    o pkg/obcs: add support for exch2
60      - the position of the boundary and prescribed values are specified using
61        a global domain with exch2 facets stacked
62        - in x for N,S boundaries (like W2_mapIO=-1)
63        - in y for E,W boundaries (so E,W boundaries do not overlap)
64    o pkg/mdsio: add support for writing global sections with exch2
65        file layout is as for obcs described above:
66        - for x-z sections, facets are stacked in x (like W2_mapIO=-1)
67        - for y-z sections, facets are stacked in y
68    o pkg/thsice:
69       - change a few routines for better vectorization: thsice_albedo.F,
70         thsice_step_temp.F, and in thsice_calc_thickn.F. Results are
71         not affected, but in particular the latter routine has changed
72         dramatically (not necessarily easier to read now).
73       - re-arrange calculation in THSICE_SOLVE4TEMP: move outside of i,j loops
74       iteration loop and call to S/R THSICE_GET_BULKF & THSICE_GET_EXF.
75    o verification/lab_sea:
76      - overhaul of the forward model experiments. The main motivation
77        for this revision is that lab_sea has proven quite fragile (when
78        changing the sea-ice code, and switching platforms) to a point
79        when failures were very hard to attribute to specific code problems.
80      - the suspected reason for this fragility is noise amplification by evp so,
81        to start with, we decided to replace evp with lsr in the main experiment.
82      - the input.lsr experiment thus became obsolete and was removed.
83      - it was replaced with the input.fd experiment that uses an
84        analytical "free drift" solution as the ice velocity field.
85      - Taken together, input (lsr), input.fd (free drift) and input.hb87 (evp)
86        will hopefully allow us to discriminate problems due to dynamic solvers
87        from problems due to the thermodynamic code. input.salt_plume (lsr) allows
88        us to test additional options, along with the salt plume package.
89      - Finally a pickup was generated (at the end of a one year run) to
90        make the tests more representative of an actual model run. Previously
91        the tests starting point was a 1m thickness of ice everywhere.
92    o pkg/seaice:
93      - introducing S\R seaice_freedrift.F that compute the
94        analytical solution to an approximate ice momentum
95        equation (i.e. the simple stress/coriolis balance).
96      - define SEAICE_ALLOW_FREEDRIFT to trigger the computation
97        of the free drift velocity field, and enable its use.
98        To actually use the free drift solution as the model ice
99        velocity field set SEAICEuseFREEDRIFT to .TRUE. in data.seaice.
100    o move "useOASIS" from PARAMS.h to EEPARAMS.h and read it from "eedata"
101    
102    checkpoint62l (2010/10/04)
103    o pkg/seaice:
104      -fix the USE_ORIGINAL_SBI version (for variable freezing point);
105      -use "celsius2K" in SEAICE_VARIABLE_FREEZING_POINT and in new
106       version (with USE_ORIGINAL_SBI undef);
107    o add hooks for new packages OASIS, the package will follow
108      (hooks exclude the seaice pkg for now)
109  o pkg/autodiff: ALLOW_AUTODIFF_WHTAPEIO  o pkg/autodiff: ALLOW_AUTODIFF_WHTAPEIO
110    Option to use mdsio_read_whalos.F/mdsio_write_whalos.F,    Option to use mdsio_read_whalos.F/mdsio_write_whalos.F,
111    rather than mdsio_readvector.F/mdsio_writevector.F, within    rather than mdsio_readvector.F/mdsio_writevector.F, within
112    adread_adwrite.F to read/write adjoint-related tapes to disk.    adread_adwrite.F to read/write adjoint-related tapes to disk.
113  o eesup and pkg/mdsio: ALLOW_WHIO  o eesupp and pkg/mdsio: ALLOW_WHIO
114      Option to write/read 2D fields to files including tiles halos.
115    Option to write/read 2D fields to files including tiles halos.    The main purpose is for adjoint related "tape I/O".
   The main purpose is for adjoint related "tape I/O".  
116    The secondary purpose is for debugging phases.    The secondary purpose is for debugging phases.
   
117    A word on those rather specifically intended purposes.    A word on those rather specifically intended purposes.
118    The code was meant to be stripped down to the bare minimum.    The code was meant to be stripped down to the bare minimum.
119    Typically, there is a unique and basic mapping to the global    Typically, there is a unique and basic mapping to the global
120    buffer (tile by tile, then proc by proc, then level by level), and    buffer (tile by tile, then proc by proc, then level by level), and
121    byteswaps and meta-files are omitted. Indeed the (*_WH*) code is not    byteswaps and meta-files are omitted. Indeed the (*_WH*) code is not
122    meant to generate user friendly outputs of various shapes and forms.    meant to generate user friendly outputs of various shapes and forms.
123    Rather it is intended, for the knowledgeabale developers,    Rather it is intended, for the knowledgeabale developers,
124    to be largely independent of the many mdsio options, to allow    to be largely independent of the many mdsio options, to allow
125    outputs of tiles+halos in debugging phases and in adjoint runs,    outputs of tiles+halos in debugging phases and in adjoint runs,
126    and to allow some flexibility in tuning (adoint) I/O depending on each    and to allow some flexibility in tuning (adoint) I/O depending on each
127    file system behavior. With regard to tuning, most of the potential knobs    file system behavior. With regard to tuning, most of the potential knobs
128    (to leave files open, to allocate the I/O indep. of usesinglecpuio,    (to leave files open, to allocate the I/O indep. of usesinglecpuio,
129    to specify directories, to specify the typical size of I/O calls, etc.)    to specify directories, to specify the typical size of I/O calls, etc.)
130    are not yet included, and will be added if they prove useful.    are not yet included, and will be added if they prove useful.
Legend:
Removed from v.1.1303  
changed lines
  Added in v.1.1319
  ViewVC Help
Powered by ViewVC 1.1.22