/[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.2172 by jmc, Thu Jul 16 23:26:38 2015 UTC revision 1.2225 by jmc, Sat Jan 16 22:05:05 2016 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/dic:
8      - use the updated seaice fraction (from seaice pkgs) for DIC forcing
9        instead of the previous time-step value ; done in new S/R DIC_FIELDS_UPDATE
10        called from the top of dic_biotic_forcing.F (instead of in DIC_FIELDS_LOAD);
11      - move also update with values from ATM-OCN Coupler in the same new S/R
12        (previously done in ocn_apply_import.F).
13    o tools/adjoint_options:
14      - specify explicitly "-f77" or "-f90" in default ad-optfile since TAF (default)
15        has been changed on Jan 14, 2016.
16    o pkg/ggl90:
17      - add alternative discretisation of vertical shear (off by default), to
18        compute the mean (@ grid-cell center) of vertical shear compon instead of
19        vertical shear of mean flow (@ grid-cell center).
20      - add correct parenthesis in ALLOW_GGL90_SMOOTH code to get the same truncated
21        results on adjacent faces of cubed-sphere grids; this fix the restart.
22    o pkg/atm_compon_interf:
23      - skip time-ave over cpl_atmSendFrq of Sea-Level Pressure (atmSLPr) and just
24        pass last value, stored directly before export to CPL ;
25        Note: this affects results of Coupled set-up where OCN uses useImportSLP=T
26      - remove atmSLPr from pickup files (no longer needed for restart)
27    
28    checkpoint65s (2016/01/13)
29    o eesupp & genmake2:
30       - first check for Function ETIME and, if it fails, then check for subroutine
31         ETIME and consistent output; use the appropriate version in timers.F
32    o model/src:
33      - to check restart test, add option to apply EXCH to U,V earlier, just before
34        integr_continuity; turned off by default.
35    o pkg/dic:
36      - remove CPP option USE_ATMOSCO2 (not needed)
37    o tools/tst_2+2:
38      - try to ignore difference between +0 and -0 in binary file, causing test
39        to fail for no good reason.
40    o eesupp & genmake2:
41      - put back modif from Jinxuan Zhu - change etime fct call to S/R call -
42        but also check (in genmake2) that test runs and returns meaningful results.
43    o coupler pkgs:
44      Finish implementation of optionally-exchanged variable switches (started on
45       Nov 12) with addition of few more exchanged fields, including 2-way exchange
46       of pkg/thsice vars to allow to use seaice dynamics (from pkg/seaice) & seaice
47       advection in OCN component; this requires sequential coupling.
48      coupler (pkg/atm_ocn_coupler):
49      - add argument msgUnit and iter number to CPL_SEND/RECV_ATM/OCN_FIELDS
50        routines (useful for debug).
51      - only export/import optionally exchanged fields (RunOff, seaice, Salt-Plume
52         flux, DIC fields) if corresponding switch is ON.
53      - add option to export RunOff to OCN (to use for DIC river input)
54      - add 2-way thSIce vars exchange; check for consistent cpl_sequential setting.
55      components (pkg/atm_ & ocn_compon_interf):
56      - to save memory space in ATM & OCN comp: add some #ifdef ALLOW_${PKG} around
57        optionally exchanged fields; add a check & stop in cpl_readparms.F.
58        ATMCPL.h: ALLOW_LAND for RunOff, ALLOW_THSICE for seaice, ALLOW_AIM for CO2;
59        OCNCPL.h: ALLOW_SALT_PLUME, ALLOW_DIC for DIC, ALLOW_THSICE for 2-W thsice.
60      - rename:   fracIce     -> sIceFrac_cpl (<- match coupler var name)
61        + in ATM: fracIceTime -> sIceFracTime ; seaIceTime -> sIceMassTime.
62      - add 2-way thSIce vars exchange ; add Salt-Plume flux ; add RunOff in OCN.
63      - hold in common block coupler optionally-exchanged variable switch
64        and set default "useImportXXX" accordingly; + check & stop if inconsistent.
65      - add diagnostics for all imported fields from coupler (in both OCN & ATM);
66        this is intended to replace the OCN cplFldtave.* output files (not updated
67        to contain new additional exchanged fields).
68    o pkg/dic:
69      - fix an other bug (also from 2014/08/30) in dic_atmos.F, in dic_int1=3 code.
70    o pkg/thsice:
71      - when coupled to ATM (useCoupler) with thSIce_skipThermo=T, skip reading
72        and writing of thsice pickup files.
73    o pkg/atm_compon_interf:
74      - move all field individual export-call out of cpl_export_import_data.F
75        into new S/R ATM_EXPORT_FIELDS
76    o pkg/atm_compon_interf:
77      - store state variables (for now, only thSIce fields) just before beeing
78        exported (in cpl_export_import_data.F, like on ocean side) instead of from
79        atm_store_my_data.F like all exported fluxes.
80      - new version of pickup_cpl: write field names into meta file and read meta
81        file to find field postion in pickup. Also only write/read fields meeded for
82        restart (same logic as main pickup file ; also using pickupStrictlyMatch)
83    o pkg/dic:
84      - fix bug (added on 2014/08/30) in dic_atmos.F for the case dic_int1= 2 or 3.
85    o pkg/seaice: add new parameter SEAICEpressReplFac to choose between
86        pressure replacement method (=1., default) or original Hibler (1979) (=0.)
87        value can be from domain [0,1] to average between the two methods
88    o pkg/atm_phys:
89      - allow to provide 2-d map of albedo and mixed-layer depth from files (but
90        hold fixed in time).
91    
92    checkpoint65r (2015/12/21)
93    o pkg/seaice:
94        fix bug in tridiagonal solvers for SEAICE_OLx/y>0, this affects the JFNK
95        solver (offline.dyn_jfnk) and global_ocean.cs32x15.seaice -> updated
96    o pkg/seaice:
97      - fix a bug in computing areaS in (J-J -> J-1) that affects runs
98        when SEAICEscaleSurfStress=.TRUE.,
99      - apply areaW and areaS in S/R seaice_calc_lhs
100    o pkg/atm_phys:
101      - add diagnostics for 3-D, Nr+1 levels diagnostics of radiative fluxes
102        (shortwave & longwave, net & downward) and diagnostics of temp. tendency
103        for individual contribution (convection, LSC cond, radiation)
104      - add Ruth Geen etal, GRL 2016, radiation optionis (see paper supp.
105        information); activated with: solar_exponent=0. and wv_exponent=-1.
106    o pkg/ggl90:
107      - within IDEMIX add upper limit for horizontal group velocity to
108        avoid instability problems in the explicit horizontal diffusion term for
109        long time steps or high resolution
110    o pkg/diagnostics:
111      - simplify setting of local-array size in diagstats_output.F: uses "numLevels"
112        (from DIAGNOSTICS_SIZE.h) to allow to write out more than Nr levels.
113    o testreport:
114      - remove "-v" (verbose option) in default mpirun command.
115    o pkg/seaice/seaice_growth:
116      - catch potential division by zero in ITD code, does not change
117        verification results
118    o pkg/ecco:
119      - retire old codes to the Attic. The calls to retired functions are
120        left in a commented form in case they ever need to be re-activated.
121      - ecco_check.F: issue warning if any old CPP option is
122        defined that corresponds to retired codes -- then stop.
123    o pkg/atm_compon_interf:
124      - change seaice fraction (use for OCN dic surf.flux) stored and sent to
125        coupler: use last time-step value instead of average over cplSendFrq_iter
126        time-step, similar to what is done for seaIceMass.
127        This will affect results (pkg/dic) of coupled ATM-OCN+dic set-up.
128    o eesupp & genmake2:
129      - from Jinxuan Zhu: change etime fct call to S/R call (fix gfortran ETIME);
130        <= remove until a fix is found for some ifort version on ia32 machine.
131    o pkg/seaice/seaice_evp:
132      - refine residual computation
133    
134    checkpoint65q (2015/11/18)
135    o pkg/seaice/seaice_evp:
136      - fix residual computations so that they no longer affect the actual
137        simulation code
138      - fix seaice_readparams to be able to actually set new parameters
139        SEAICEaEVPcStar and SEAICEaEVPalphaMin, also fix logic a little
140    o model/src:
141      - fix wrong truncation (single prec) from "FLOAT" function in expressions
142        like: deltaT*FLOAT(myIter) when myIter is large:
143         a) replaced by DFLOAT (if not seen by TAF);
144         b) remove "FLOAT" and use implicit type conversion (if seen by TAF).
145    o pkg/aim_v23:
146      - write/read the expected CO2 pickups (from S/R AIM_DO_CO2): fix reading
147        when pickupSuff is set ; do write rolling pickups + at the end of the run.
148    o pkg/atm_ocn_coupler + ocn_compon_interf:
149      - in prep for using pkg/seaice Dyn in coupled set-up, implement sequential
150        coupling option-1: ATM -(cpl)-> OCN and then OCN -(cpl)-> ATM ;
151        add corresponding call in DO_OCEANIC_PHYS.
152    o coupler pkgs:
153      Start a new way to set-up optionally-exchanged variables (previously, only
154       there for exchange_DIC): switch on/off from 1 place, at run-time, in coupler
155       data.cpl (new) instead of at compile time in ATMIDS.h & OCNIDS.h ;
156       send switch to ATM & OCN compon before reading their own "data.cpl" (so that
157       compononent param default can be set accordingly).
158      Note: this new switch (cpl_exchange_[xxx]) is not a per variable switch
159       but rather applies to a group a fields (going both ways like for DIC)
160       and corresponds to an option of coupling config.
161      coupler (pkg/atm_ocn_coupler):
162      - change parameter file from "data" to "data.cpl" (in cpl_read_params.F),
163        and read namelist directly (without eessup S/R calls) instead of
164        in the very un-flexible old way, 1 param per line, in fixed order;
165      - add error-counter to track errors (e.g., missing "data.cpl") and to stop
166        cleanly in CPL_CHECK_CPLCONFIG
167      - add new coupling-exchange selectors (some not yet active) including
168        cpl_exchange_DIC (to set to 3 to replace atmCpl_exchange_DIC=T in ATMIDS.h
169                                             and ocnCpl_exchange_DIC=T in OCNIDS.h);
170      - add 2 new S/R (CPL_SEND_ATM/OCN_CPLPARMS) to send coupler-params to each
171        component ; called directly from coupler.F, after comp. registrations;
172      - rename few variable-ID names (in ATMIDS.h & OCNIDS.h):
173         atm/ocnSeaIceName --> atm/ocnSIceMassName ; atm/ocnFCO2Name ->
174          -> atm/ocnCO2FluxName ; atm/ocnWSpdName --> atm/ocnWindSpdName
175      components (pkg/atm_ & ocn_compon_interf):
176      - add ATM/OCN_CPL_OPTIONS.h in component pkgs and include the corresponding
177        CPP header file in all pkg source files.
178      - move checking of comp. configs out of cpl_exch_configs.F into new S/R
179        ATM/OCN_CHECK_CPLCONFIG;
180      - add error-counter to track errors (e.g., in "data.cpl") and to stop
181        cleanly in ATM/OCN_CHECK_CPLCONFIG;
182      - add new S/R CPL_IMPORT_CPLPARMS to import coupling-exchange selectors
183        from coupler ; called from initialise_fixed.F, after PACKAGES_BOOT and
184        before PACKAGES_READPARMS;
185      verification/cpl_aim+ocn/:
186      - update ATMIDS.h & OCNIDS.h (in shared_code) and use new parameter file
187        "data.cpl" for coupler (in input_cpl).
188    o eesupp/inc:
189      - change the MPI-global-sum default (in CPP_EEOPTIONS.h) to:
190        #define GLOBAL_SUM_ORDER_TILES (only affects MPI built of global_sum_tile.F)
191        to always cumulate tile local-sum in the same order (still uses
192        MPI_Allreduce like standard code).
193      - in verification experiment code, remove local CPP_EEOPTIONS.h_mpi files
194        that are now identical to standard version (in eesupp/inc).
195    o pkg/shelfice:
196      - add SHELFICETransCoeffTFile for initializing spatially varying
197        heat and salt transfer coefficients
198    o pkg/atm_compon_interf & ocn_compon_interf:
199      - combine into 1 (new) S/R CPL_EXPORT_IMPORT_DATA the two previous
200        S/R CPL_EXPORT_MY_DATA & CPL_IMPORT_EXTERNAL_DATA
201      - move bi,bj loops inside S/R OCN_STORE_MY_DATA
202    o pkg/atm_ocn_coupler:
203      - initialise to zero ATMVARS.h & OCNVARS.h arrays.
204    o pkg/rbcs:
205      - select RBCS-mask for Temp & Salt the same way as for pTracer ; this allows
206        to compile with maskLEN=1 when all tracers use the same RBCS-mask.
207      - add a short report (to STDOUT) about RBCS-mask settings.
208    
209    checkpoint65p (2015/10/23)
210    o pkg/seaice/seaice_evp.F: fix store directives and key computations, since
211        there is no adjoint test for EVP, we don't know what it does to the adjoint
212        solution, but at least there are no more recomputation warnings left
213    o pkg/seaice/seaice_evp.F: slight modification of averaging so that EVP stands
214        a chance of giving the same results as a fully converged VP (JFNK) solver
215    o tools/adjoint_options/adjoint_diva:
216      - fix bug (typo) in Tang-Linear FLAGS setting for MPI + DIVA
217    o pkg/streamice:
218      - spotted bug in streamice_driving_stress using dyG where should have been dxG
219    o pkg/exf:
220      - remove STOP statement that prevented to use wind or (wind-)stress
221        interpolation for constant-in-time forcing fields (wind/stress period = 0).
222    o pkg/seaice:
223      - replace SEAICE_EVP_USE_ADAPTATION by run time parameters, requires
224        5 additionals 2D fields and 6 new store directives
225      - change logic in seaice_readparms.F: EVP code can now be turned on in
226        various ways, setting SEAICE_deltaTev is no longer required,
227        alternatively set SEAICE_evpAlpha/Beta, SEAICEuseEVPstar, SEAICEuseEVPrev,
228        or SEAICEaEVPcoeff
229    o pkg/shelfice/shelfice_thermodynamics.F
230      - redefining uLoc/vLoc to u/v, instead of 2u/v
231    o pkg/seaice/seaice_evp.F: introduce "adaptive EVP" by Madlen Kimmritz
232      - this is an intermediate state of a truely converging EVP that has
233        the potential of outperforming implicit VP methods, I do not recommend
234        using it yet.
235      - for now, "adaptive EVP" is turned on by setting SEAICE_EVP_USE_ADAPTATION
236        and all parameters are hardwired
237      - for convergence it still requires very smooth regularisations, that means,
238        if SEAICE_DELTA_SMOOTHREG is defined, deltaCreg=deltaC+deltaMin, which is
239        different from the VP case with deltaCreg=sqrt(deltaC**2+deltaMin**2); also
240        specific averaging appears to be required, still subject to tuning
241      - still needs cleaning and I need to turn parameter into runtime parameters
242      - add code to compute and print residuals of iteration, compile with
243        defining ALLOW_SEAICE_EVP_RESIDUAL to enable it
244    o pkg/seaice/seaice_evp.F: fix a bug (factor of four in front of e12Csq)
245      that entered with r1.54 on Sep04, 2015, -> update results
246    
247    checkpoint65o (2015/09/14)
248    o pkg/mom_vecinv & mom_common:
249      - add deep-atmosphere geometry factor in mom_calc_visc.F + in routines
250        called from there.
251      - add anelastic and deep-atmosphere geometry factor in pkg/mom_vecinv ; this
252        allows to use Vector-Invariant form in deep atmos and anelastic formulation
253    o pkg/seaice: introduce fast ice parameterisation following Itkin et al. (2015)
254      and Koenig-Beatty+Holland (2012)
255      - two new parameters SEAICE_tensilFac and SEAICE_tensilDepth
256      - global 2D field tensileStrength, computed in seaice_calc_ice_strength
257      - requires an extra input argument for S/R seaice_calc_viscosities
258      - tensileStrength's depth dependence is different from Itkin et al. (2015)
259        (to be changed and played with, for now exponential)
260    o pkg/seaice: introduce new parameter SEAICEscaleSurfStress (default = .FALSE.)
261      - if TRUE scale ice-ocean and ice-atmosphere stress acting on sea ice
262        by AREA according to Connelly et al. (2004), JPO.
263      - For EVP and the LSR solver, the implicit part of the drag term is
264        the only non-zero term in the denominator when there is no ice (which
265        prevented running the model with zero ice-ocean drag). If
266        SEAICEscaleSurfStress = .TRUE., this results in a division by zero
267        (or zero main diagonals BU/BV) which need to be caught. A practical
268        consequence is that for open water, the momentum equation reduce
269        to trivially 0 = 0 + 0 (for EVP). For LSR, BU/BV are reset to 1, if
270        they happen to be zero, often leading to a non-zero solution over
271        open water.
272    o pkg/seaice: revise EVP code (seaice_evp.F) in preparation for a more
273        efficient EVP method.
274      - introduce area-weighted averages for e12 (analogously to
275        seaice_calc_viscosities.F); this totally changes lab_sea.hb87 (2 digits
276        remain, but also changing the order of terms in the previous average
277        changed the results as much) -> update experiment
278      - add SEAICE_DELTA_SMOOTHREG code
279      - disentangle computation of zetaC/deltaC and zetaZ/deltaZ, also for
280        TEM-code
281      - adjust loop bounds so that only required fields are computed
282      - add new averaging code for zetaZ/deltaZ (again similar to
283        seaice_calc_viscosities.F), inactive by default for now (will change
284        results again), but will replace old code eventually
285    o eesupp/src:
286      - add an other alternative for multi-proc GLOBAL_SUM (in global_sum_tile.F)
287        to always cumulate tile local-sum in the same order
288        (option #define GLOBAL_SUM_ORDER_TILES in CPP_EEOPTIONS.h).
289        Still uses MPI_Allreduce (like standard code) but on array of tiles
290        and should be faster than GLOBAL_SUM_SEND_RECV code.
291      - switch almost all verification experiments that use GLOBAL_SUM_SEND_RECV
292        (defined in CPP_EEOPTIONS.h_mpi) to GLOBAL_SUM_ORDER_TILES instead.
293      - add CPP_EEOPTIONS.h_mpi with GLOBAL_SUM_ORDER_TILES defined in 14 more exp.:
294         aim.5l_cs, cheapAML_box, deep_anelastic, dome, fizhi-cs-32x32x40,
295         fizhi-cs-aqualev20, fizhi-gridalt-hs, global_ocean_ebm, matrix_example,
296         natl_box, solid-body.cs-32x32x1, tutorial_advection_in_gyre,
297         tutorial_baroclinic_gyre and tutorial_held_suarez_cs.
298    
299    checkpoint65n (2015/07/29)
300    o pkg/thsice:
301      - add a warning when used with pkg/aim but without updating SST
302        (stepFwd_oceMxL=F & useCoupler=F & tauRelax_MxL<>-1 ).
303    o pkg/offline:
304      - change where Offline-Fields are loaded, from load_fields_driver.F to
305        directly within FORWARD_STEP, just after updating iter number and time;
306      - move setting of GM and KPP diffusivity out of OFFLINE_FIELDS_LOAD into new
307        S/R OFFLINE_GET_DIFFUS that is called towards the end of DO_OCEANIC_PHYS;
308       These changes allow:
309        a) to recover exactly on-line results using deltaToffline=deltaTClock,
310        tested (using offlineIter0= nIter0 -1 and offlineTimeOffset= -deltaT/2)
311        both with Lin & Non-Lin Free-Surf, both with staggerTimeStep & synchronous
312        time-step, both with GM & KPP files and with recomputing these fields;
313        b) no needs to switch on staggerTimeStep anymore;
314        c) simplify rescaling of offline horiz. velocity with Non-Lin Free-Surf.
315      - update output of experiment tutorial_dic_adjoffline (only 3 digits left).
316  o pkg/offline:  o pkg/offline:
317    - change (rarely used) offset-iteration parameter to an offset time (in s):    - change (rarely used) offset-iteration parameter to an offset time (in s):
318      can easily recover previous offset-iter effect by setting offset-time      can easily recover previous offset-iter effect by setting offset-time

Legend:
Removed from v.1.2172  
changed lines
  Added in v.1.2225

  ViewVC Help
Powered by ViewVC 1.1.22