/[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.1774 by mlosch, Fri Nov 23 15:54:18 2012 UTC revision 1.1938 by jahn, Fri Nov 8 20:04:16 2013 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/cost:
8      - define a "final cost function" for each tile, and use GLOBAL_SUM_TILE_RL
9        (instead of _GLOBAL_SUM_RL) to get the final cost from all procs & threads:
10        o this fixes global-sum value if multi-threads (since GLOBAL_SUM_TILE_RL
11          cannot be applied to var in common block).
12        o allows to get results independent of tile-per-proc setting (if using
13          slower GLOBAL_SUM_SEND_RECV option).
14      - update accordingly: seaice_cost_final.F, shelfice_cost_final.F,
15        thsice_cost_final.F and ecco_cost_final.F (this one affects lab_sea
16        forward gradients: update output).
17      - update */code_oad/externalDummies.F by adding dummy GLOBAL_SUM_TILE_RL
18        S/R and, for most exp., also comment out GLOBAL_SUM_R8.
19    o pkg/mom_common:
20      - move to pkg/mom_common and model/src (previously in tutorial_deep_convection
21        code) 2nd version of isotropic 3-D Smagorinsky code interface: strain and
22        viscosity are locally declared in dynmics.F and pass as argument to CALC_GW;
23        ensure that all field value that are used are set.
24    o tutorial_deep_convection (smag3d):
25      - add missing scaling in mom_calc_smag_3d.F;
26      - set smag3D_coeff in input.smag3d such that results do not change.
27    o pkg/dic:
28      - change dic_cost.F to compute cost on each tile first and then apply
29        GLOBAL_SUM_TILE_RL: this affects forward gradient of AD experiments
30        tutorial_dic_adjoffline & tutorial_global_oce_biogeo; update output.
31      - remove kLev argument and use GLOBAL_SUM_TILE_RL in tracer_meanarea.
32    o eesupp/src & pkg/mdsio:
33      - remove gather/scatter_vector.F (from eesupp) and mdsio_read/writevector.F
34        (from pkg/mdsio) since they are no longer used for tape IO.
35      - move simple C subroutine "setdir.c" from pkg/compon_communic to eesupp/src
36        (+ keep a copy in pkg/atm_ocn_coupler).
37    o verification/isomip:
38      - split domain in 8 squared tiles (instead of just 2) for FWD, AD & OAD tests.
39    o tools/build_options:
40      - optfile linux_amd64_ifort11, fast (no ieee) case:
41        to get a working restart for all experiments, add "-fp-model source"
42        to FOPTIM (also remove "-opt-streaming-stores auto" <- this is the default).
43    
44    checkpoint64p (2013/10/20)
45    o pkg/autodiff:
46      - adread_adwrite.F: when type and length do not match ALLOW_AUTODIFF_WHTAPEIO
47        criteria, return to simple (default) tape IO format.
48    o pkg/mdsio (+ eesupp/src):
49      - rewrite/simplify tape-IO S/R (new S/R MDS_READ/WRITE_TAPE replace
50        previous MDSREAD/WRITEVECTOR) with 2 array argument from each type (R4/R8);
51      - fix globalFile and singleCpuIO options using simpler global mapping.
52      - new pairs of scatter/gather vector S/R (one for each type, generated from
53        template) corresponding to simpler global mapping.
54      - adread_write.F: switch to new S/R (but leave old calls commented) and write
55        to globalFile when parameter globalFiles is set to True.
56      - switch on useSingleCpuIO in AD experiments: global_ocean.90x40x15 & isomip.
57    o pkg/autodiff:
58      - adread_adwrite.F: add a hack for RS arrays (when RS=real*4);
59        with ALLOW_AUTODIFF_WHTAPEIO defined, since corresponding (real*4) is
60        missing, just add a STOP if array type do not match.
61    o pkg/exf:
62      - to allow to use seaice-fraction from pkg/exf in other pkgs than pkg/seaice:
63         rename CPP option EXF_ALLOW_SEAICE_RELAX to EXF_SEAICE_FRACTION
64         and rename array obsSIce to exf_iceFraction.
65    o tutorial_deep_convection:
66      - add 1rst version of isotropic 3-D Smagorinsky code (from Louis-Philippe),
67        for now all in this exp. code dir.
68      - add secondary test using Smag-3D viscosity (keep very low backgroud visc for
69        testing) and flux-limter advection scheme for Temp. (no explicit diffusion).
70    o eesupp/src:
71      - added SINGLE_DISK_IO option because reading of parameter files on
72        very large number of processors did not work with existing set-up
73    
74    checkpoint64o (2013/09/25)
75    o pkg/atm_phys & atm_compon_interf:
76      - allows for coupling of gray-atmosphere (atm_phys) with ocean.
77    o tools/build_options:
78      - switch from "-fPIC" to "-mcmodel=medium" in optfiles: linux_amd64_gfortran
79        & linux_amd64_ifort11 and in linux_amd64_g77 & linux_amd64_open64 also
80        remove CPP setting (was identical to genmake2 default).
81    o pkg/gmredi & model/src: Some updating and work done on the residual model.
82      - created new file model/src/calc_eddy_stress.F to calculate the eddy stress.
83      - fixed up the calculation of the eddy stress and how it is added to the
84        external forcing via taueddy_external_forcing.F
85      - at the moment, it is very experimental and only works with the variable eddy
86        diffusivity calculated in gmredi_k3d.F
87    o pkg/seaice/seaice_check.F: add a check to catch the case of specifying
88        lw/swfluxfile instead of lw/swdownfile in data.exf together with pkg/seaice
89    o testreport (& tools/do_tst_2+2):
90      - for MPI test, set default COMMAND to 'mpirun -v -np TR_NPROC ./{executable}'
91    o tools/build_options: linux_amd64_gfortran & linux_amd64_ifort11:
92      - remove CPP='cpp -traditional -P' (this is genmake2 default)
93      - add switch (for big setups) to use 'medium' memory-model (but leave it
94        commented out for now). Note: needed with FC 19 to compile big setups.
95    o pkg/atm_compon_interf:
96      - regroup per field origin (model & pkgs) the coupling storage routines
97         (which accumalate in time each coupling field)
98      - move bi,bj loops inside atm_store_my_data.F and store also wind-stress;
99        update accordingly + simplify aim_do_physics.F
100    
101    checkpoint64n (2013/08/26)
102    o pkg/bbl:
103      - update formulation to reduce dispersion of tracer properties in bbl
104      - replaced horizontal and vertical transport parameters with velocities
105      - changes results for verification/global_with_exf.yearly
106    o verification/halfpipe_streamice:
107      - change to matrix construction; update both Fwd & AD output.
108    o verification/tutorial_tracer_adjsens:
109      - OpenAD test: switch to vector-invariant momentum since fluxform with
110        rStar causes problems with OpenAD; update output.
111    o tools:
112      - genmake2: define a common prefix for all temp-file names ; move them
113        to local /tmp disk (but would be easy to move them back);
114      - xmakedepend: apply the same switch for temp-files.
115    o pkg/offline:
116      - changes for use with rStar (NLFS): assume that we continue to provide
117        UVELMASS/hFacW_ini & VVELMASS/hFacS_ini (hFacWi,S_ini as written in grid
118        output files) and apply r* scaling at each time-step (function of the time
119        varying Eta field); use normal pickup files (needed for restart with NLFS).
120    o adjoint with OpenAD:
121      - move OpenAD-specific files from verification/OpenAD/code_oad_all to a new
122        package openad.  This package has to be compiled for all runs with OpenAD
123        instead of specifying code_oad_all as a modification directory.
124      - link the fortran templates from tools/OAD_support to the build directory
125        for preprocessing to avoid conflicts between different setups compiled
126        simultaneously and simplify cleanup.
127    o pkg/exf:
128      - in exf_getforcing.F: comment out most of the EXCH calls (fields used by
129        main model or other pkgs are exchanged in EXF_MAPFIELDS).
130    o model/src:
131      - do_stagger_fields_exchanges.F: skip EXCH calls when using pkg/offline
132        (similar to do_fields_blocking_exchanges.F).
133    o tools/genmake2:
134      - remove TAF output src code just before calling staf (when not doing a fresh
135        built, will stop if things go wrong with staf).
136      - fix target "adobj" (left from TAF version-2 transition)
137    o pkg/diagnostics:
138      - add a parameter to track the status of pkg/diagnostics activation;
139        this replace/extend the use of logical param "settingDiags" (to check
140        when adding diag to the list) to also check any DIAGNOSTICS_[]FILL* call.
141        New S/R to print error message and stop when pkgStatus is not right.
142      - to try: add a warning when "data.diagnostics" exist but useDiagnostics=F.
143    
144    checkpoint64m (2013/08/13)
145    o model/inc:
146      - move FCORI_MACROS.h from PARAMS(_MACROS).h to GRID_MACROS.h (left from
147        8 yrs ago modif when Coriolis arrays were moved from PARAMS.h to GRID.h)
148    o pkg/kpp & aim_v23:
149      - only call DIAGNOSTICS_FILL when useDiagnostics=TRUE.
150    o pkg/seaice:
151      - restore (last time on 2011/12/24) ability to compile pkg/seaice with
152        undefined pkg/exf ALLOW_ATM_TEMP or ALLOW_RUNOFF CPP-options.
153      - only call DIAGNOSTICS_FILL when useDiagnostics=TRUE.
154    o eesupp/src:
155      - fix check for overlap-size argument (myOLn.GT.MAX_OLY_EXCH)
156        in exch1_rx[_ad].template (thanks to Francois Lefeuvre report).
157    o switch to TAF version-2:
158      - add flow directives to "freeze" naming convention for ad-routines that
159        are used in I/O routines of adjoint state and associated common blocks
160        that are used in adcommon.h (backward compatible with version-1)
161      - update all ad-optfiles (in tools/adjoint_options + few local to some
162        verification exp) to switch to latest TAF version (currently 2.5.5);
163        merge "*_mpi" adoptfiles into the corresponding non-mpi (but commented
164        out except for diva).
165      - the major change with TAF version-2 is that both the differentiated
166        routines and the forward routines (new) are returned as part of the TAF
167        output source code. Update genmake2 accordingly.
168      - new output for AD experiment tutorial_global_oce_biogeo (better AD gradient
169        with version-2 when compared to TLM or FWD gradient); other experiments
170        output show minor changes (machine truncation level); update the few AD &
171        TLM output with 12 matching digits only (some prior to version-2 switch).
172    
173    checkpoint64l (2013/08/05)
174    o model/src:
175      - in dynamics.F, skip the call to CALC_VISCOSITY if momViscosity=F
176      - main_do_loop.F: no longer includes additional (AUTODIFF) header files
177        when using TAF but only if using OpenAD.
178    o pkg/mom_common+fluxform+vecinv:
179      - always set horiz. viscosity arrays to background value before calling
180        MOM_CALC_VISC (in MOM_FLUXFORM & MOM_VECINV) and call S/R MOM_CALC_VISC
181        only when using variable horiz. viscosity (useVariableVisc=T);
182      - simplify mom_vecinv.F (only 1 block for momViscosity).
183    o pkg/mom_common: variable horiz. viscosity in NH code
184      - move horiz-viscosity 3-D arrays viscAh_W & viscA4_W (used in calc_gw.F)
185        from model/inc/NH_VARS.h to pkg/mom_common/MOM_VISC.h
186      - store in common block (in MOM_VISC.h) logical parameters: useHarmonicVisc,
187        useBiharmonicVisc & useVariableVisc.
188      - fix old bug (already in 1rst implementation of variable horizontal viscosity
189        in calc_gw.F, version 1.27, 2006/06/20):
190        horizontal bi-harmonic viscosity in gW was missing when wiscA4W=0
191        even if variable bi-harmonic viscosity (e.g., viscA4Grid) was non-zero.
192      - change background horiz viscosity in vertical momentum (gW) to always be
193        viscAhW,viscA4W (instead of switching to viscAhD, viscA4D when using
194        any type of variable viscosity, e.g., viscAh/4Grid)
195      - change default value of viscAhW,viscA4W to viscAhD,viscA4D (instead of
196        viscAh,viscA4) which default to viscAh,viscA4 :  this avoid missing
197        horiz. viscosity in W eq when just setting viscAhD/Z,viscA4D/Z (since
198        default viscAh,viscA4 is zero).
199    o model/src:
200      - extend length of drC & recip_drC from Nr to Nr+1.
201    o Divided Adjoint (DIVA) with MPI:
202      - genmake2: do not create local ./mpi_headers dir when no file to copy into.
203      - remove option ALLOW_DIVIDED_ADJOINT_MPI (replaced by ALLOW_DIVIDED_ADJOINT
204        and ALLOW_USE_MPI) from model/src code and from ECCO_CPPOPTIONS.h.
205    
206    checkpoint64k (2013/07/23)
207    o tools/OAD_support:
208      - avoid the need for file "oadtempflile" in insertTemplateDir.bash
209        and remove it from genmake2.
210    o tools/build_options:
211      - remove MPI_HEADER_FILES_INC setting (no longer used) and comment out
212        MPI_HEADER_FILES list (let genmake2 find the missing headers)
213    o verification/halfpipe_streamice:
214      - change to driving stress formulation; update both Fwd & AD output.
215    o tools/genmake2 and Divided Adjoint with MPI:
216      - remove default setting of MPIINCLUDEDIR from $MPIHOME (always been broken);
217      - only create local ./mpi_headers dir if using MPI and MPIINCLUDEDIR is set;
218      - add new variable (LOCAL_MPI_HEADERS, replace MPI_HEADER_FILES_INC) in
219        Makefile, which is set by genmake2, and check which mpi-headers are needed.
220        Note: need to update adjoint_options/adjoint_diva_mpi accordingly.
221      - update target "CLEAN" (to also remove ./mpi_headers/ dir).
222    o testreport:
223      - changes for Divided Adjoint Run: when file "run_ADM_DIVA" is present,
224        perform additional runs (number of additional runs is "add_DIVA_runs" taken
225        from file run_ADM_DIVA ; should match "nchklev_3" in file "code_ad/tamc.h").
226    o verification/lab_sea (AD test):
227      - change to test Divided Adjoint Run in this set-up (both standard AD test
228        and 2 secondary AD test: noseaice + noseaicedyn); drop the evp AD test.
229    o pkg atm_ & ocn_compon_interf:
230      - remove cpl_earlyExpImpCall=F option (was mainly a hack added for the old
231        staggerTimeStep implementation) and disable (retired) corresponding
232        parameter.
233    o pkg/offline:
234      - comment out unused field-array: ICEM, ConvectCount, hflx[01], icem[01];
235      - add S/R OFFLINE_CHECK to print offline parameter summary
236        and stop if loading a file that is not used;
237      - enable to calculate GMRedi tensor, KPP mixing and convective index
238        when the corresponding file(s) is/are not specified (but pkg is used).
239      - comment out resetting of monitorFreq (to endTime) since it still useful
240        for monitoring CFL numbers and checking offline fields.
241    o OpenAD:
242      - do not define OAD_DEBUG (in tools/adjoint_options/adjoint_oad) to get
243        smaller standard output) and generate new (smaller) output (on baudelaire
244        using gfortran with -ieee) for exp. tutorial_tracer_adjsens
245    o model/src:
246      - call EXTERNAL_FORCING_SURF outside bi,bj loops (in do_oceanic_phys.F),
247        (where external_forcing_adjust was called) and change list of arguments;
248      - merge external_forcing_adjust.F into external_forcing_surf.F
249      - put surface relaxation of SST & SSS in new S/R (forcing_surf_relax.F)
250        which also contains balancing surface relaxation (ALLOW_BALANCE_RELAX,
251        previously in file "balance_relax.F", now removed); remove corresponding
252        arrays in FFIELDS.h.
253      - EXTERNAL_FORCING_SURF is now called when using pkg/offline; implement
254        a hack for getting back surfaceForcingS (loaded from file) when using
255        ALLOW_OLD_VIRTUALFLUX (pkg/dic, in DIC_OPTIONS.h) in offline set-up.
256      - update output of AD test experiment global_ocean.cs32x15.seaice_dynmix
257        (was getting only 11 digits for AD-grad, due to different truncation
258         error in AD code).
259    
260    checkpoint64j (2013/07/04)
261    o pkg/offline:
262      - to enable the use of other features with pkg offline (e.g., obcs),
263        also call S/R DO_OCEANIC_PHYS in offline run - more modifs to come.
264      - comment out "allowFreezing" (now active) in 2 offline testing exp.
265    o pkg/exf ALLOW_RUNOFTEMP
266      - fix sign and make compatible with temp_EvPrRn .NE. UNSET_RL
267      - ahdded RUNOFTEMP test to verification/global_ocean.cs32x15/input.seaice
268      - changes verification/global_ocean.cs32x15/results/output.seaice.txt
269    o pkg/atm_phys:
270      - allows to add some damping on uVel,vVel (strength only function of k).
271    o verification, ad-tests:
272      - reduce length of the simulation for AD test: global_ocean.cs32x15.seaice,
273        global_ocean_ebm and offline_exf_seaice.thsice.
274    o verification/halfpipe_streamice:
275      - found a bug in s/r pkg/streamice/streamice_advect_thickness which caused
276        flux input boundary conditions to be ignored; a fix for this bug changed
277        verification results for forward test.
278      - add adjoint test.
279    o pkg/ptracers:
280      - allow to restart without SOM pickup (resetting SOM to zero).
281    o pkg/gmredi:
282      - New parameterisation for the calculation of K and (GM_PsiX,GM_PsiY) added.
283      - It uses mixing length theory and suppression by steering level effects to
284        calculate a diffusivity used in the Redi tensor (see Bates et al., submitted
285        to JPO, for more details).
286      - It also uses the diffusivity to calculate an eddy induced streamfunction
287        based on a PV closure (rather than a buoyancy closure as in GM).
288      - New files added to pkg/gmredi:
289        gmredi_k3d.F calculates the diffusivity for the Redi tensor and the eddy
290          induced overturning streamfunction (called from gmredi_calc_tensor)
291        gmredi_calc_eigs.F calculates the eigenvalues and eigenvectors required
292          for use in gmredi_k3d
293        gmredi_calc_urms.F calculates the the rms eddy velocity used to calculate
294          the eddy diffusivity
295    
296    checkpoint64i (2013/06/18)
297    o pkg/cheapaml:
298      - take src files from MITgcm_contrib/verification_other/offline_cheapaml/code
299        that enable the use of seaice (pkg/thsice thermo & pkg/seaice dynamics)
300        with pkg/cheapaml.
301    o utils/python: add subpackage for cubed sphere function, first example is
302        simple cs.pcol to plot 2D fields on a map
303    o pkg/streamice:
304      - This is a new package in MITgcm, which solves dynamic equations for ice
305         sheets and ice streams
306      - At the moment it is independent of atmosphere/ocean dynamics, or any
307        other physical packages in MITgcm, even ShelfIce (although it is anticipat-
308        ed that this will change)
309      - There is a new verification exercise that tests it, halfpipe_streamice,
310        evolves a 40x80 km ice stream/ice shelf (at coarse resolution) in a
311        rectangular embayment
312      - halfpipe_streamice is forward-only for now, will have TAF test soon
313      - for questions, email dgoldber@mit.edu
314    o pkg/thsice:
315      - fix filling of snow-precip diagnostic (SIsnwPrc) (left from the removing
316        of snowPrc from common bloc, on May 2nd, thsice_step_fwd.F revision 1.40);
317      - remove diagnostics 'SI_FrcFx' ( now identical to 'SI_Fract' with new
318         call sequence: ice-Thermo(thsice),ice-Dyn,ice-Advect );
319      - allow precip from pkg/cheapaml to make snow over seaice.
320    o pkg/cfc:
321      - move in separated header file (CFC_ATMOS.h) time-serie arrays of Atmos CFC
322      - remove atmos CFC 2-D arrays from common bloc (now local variables)
323      - change some hard-coded constant to run time parameter for reading and
324        and interpolating atmospheric CFC ; keep same default as previously.
325      - more flexible reading of Atmos CFC (number of records no longer hard-coded)
326      - put arguments in standard order: myTime, myIter, myThid
327      - use standard S/R GET_PERIODIC_INTERVAL for time interpolation of
328        atmospheric CFC ; fix for the case myIter < PTRACERS_iter0.
329      - fix calculation of surf tendencys for partial-cell and/or Non-Lin Free-Surf
330      - new routine to print summary of CFC params and check for consistency.
331    o pkg/seaice: add uice/vice to the parameter list of s/r seaice_advdiff,
332      similar to thsice_advdiff.
333      Note that for the B-grid case the routine parameters uc/vc are overwritten.
334    o pkg/cfc:
335      - change time interpolation for atmospheric CFC-11 (now same as for CFC-12);
336        update output (CFC-11) of experiment cfc_example and tutorial_cfc_offline.
337      - fix interpolation of atmphspheric CFC between Sorthern and Northern
338        hemisphere value.
339    o pkg/seaice:
340      - allow to use seaice dynamics without pkg/exf (e.g., using pkg/cheapaml)
341    o pkg/seaice: add Adams-Bashforth2 time discretization for ice dynamics;
342                  so far only for JFNK-solver, turn on with SEAICEuseAB2 = .TRUE.
343      - requires an additional time level that is stored in the pickup
344      - you can start with AB2 from a pickup without this time level with
345        pickupStrictlyMatch = .false.
346    
347    checkpoint64h (2013/05/26)
348    o pkg/cheapaml:
349      - replace one ()**.5 with SQRT + use temp var "cdu = tau/Du" (coare3_flux);
350        affect machine truncation (but exp cheapAML_box still pass with 13 digits).
351      - consistent treatment of uWind,vWind location: assume everywhere in
352        pkg/cheapaml that they are on C-grid (@ uVel,vVel location)
353        (Note: was already the case for Tair,Qair advection);
354        add option (INCONSISTENT_WIND_LOCATION) to recover old results.
355      - compute diabatic forcing only in the interior (excluding overlap regions)
356      - fill in array "Qsw" (for short-wave heating);
357      - add diagnostic for precip ; fill-in diagnostics for cheapAML state vars
358        before any update (consistent with MITgcm state var diagnostics);
359      - cheapaml_fields_load.F & cheapaml_init_varia.F : only fill outer edges if
360        non periodic ; + adjusted for uWind,vWind (C-grid).
361    o pkg/seaice:
362      - remove dependence of pkg/seaice versus pkg/exf (file: pkg/pkg_depend) and
363        add explicitly pkg/exf in experiments that currently use pkg/seaice.
364      - allows to compile pkg/seaice without pkg/exf.
365    o tools/OAD_support:
366      - use a local copy (ourRWCPinBd) to ensure index ourRWCP stays within bounds
367        of array "ourStepOf" in function rvNextAction (revolve.F90).
368    o pkg/aim_compon_interf renamed to pkg/atm_compon_interf
369      - change pkg name (to be used with other atmospheric physics pkg)
370        but by-pass CVS when making the new pkg so that we keep the CVS history.
371    o pkg/offline:
372      - call OFFLINE_FIELDS_LOAD from pkg/offline/offline_init_varia.F
373        (previously called from pkg/dic_surfforcing_init.F)
374    o model/src:
375      - re-order sequence of calls within packages_init_variables.F, using
376        standard order, except for:
377         DIAGNOSTICS_INIT_VARIA and OFFLINE_INIT_VARIA (called @ the top),
378         CTRL_MAP_INI and OBCS_INIT_VARIABLES (called @ the bottom)
379    o new pkg: atm_phys
380      - add entry point and source code (previously in contrib
381          verification_other/atm_gray/code_gray/) for new pkg "atm_phys"
382        (atmospheric physics pkg from P. O'Gorman and T. Schneider, JCl, 2008).
383      - re-order sequence of calls within packages_readparms.F,
384         packages_init_fixed.F and packages_write_pickup.F
385        (but skip packages_init_variables.F because of TAF issues).
386      - update contrib/*/atm_gray/ experiment.
387    o pkg/seaice (ITD code):
388      - simplify and fix the writing of pickup-files for ITD case
389    o pkg/seaice (ITD code):
390      - introduce floe size dependent lateral melt for ITD case;
391      - cleaning (avoid division by zero, change to f77 syntax, remove tabs)
392      - update seaice_itd experiment main output.
393    o pkg/thsice:
394      - remove snowPrc from common block THSICE_FLUX, in THSICE_VARS.h and
395        declared it locally (without bi,bj) in thsice_main.F & aim_do_physics.F;
396        changed accordingly pkg/atm2d.
397    
398    checkpoint64g (2013/05/01)
399    o pkg/exf ALLOW_RUNOFTEMP
400      - if defined,  river and glacier runoff temperature can be read-in from
401        runoftempfile; runoftempfile needs to have same format as runofffile
402    o pkg/thsice (& pkg/bulk_force):
403      - In case Energy-Reference-Level (ERL) is used (temp_EvPrRn=0), account
404        for energy content of Precip + RunOff & Evap (like with pkg/exf).
405      - change computation domain to always be the interior (without overlap).
406      - fix bug (wrong overlap) when using EXF and stressReduction > 0.
407    o verification/tutorial_global_oce_biogeo:
408      - add preliminary OpenAD setup (still has problems)
409    o pkg/seaice:
410      - turn some parameters for choosing the convergence criterion of
411        the inexact Newton method (JFNK) into runtime parameters for convenience
412      - clean up seaice_calc_viscosities.F: remove/retire two options for
413        SEAICetaZmethod (1,2), introduce area weighting for method 3; this changes
414        the results of lab_sea/input because a division has been replaced by
415        a multiplication with the reciprocal (sumNorm)
416      - modify offline_exf_seaice/input.dyn_jfnk/data.seaice to improve
417        convergence of solver (changes results, of course), and to test some
418        new features
419    o model/src:
420      - remove unused (since Aug 2006) S/R SET_WRITE_GLOBAL_PICKUP and its call.
421      - remove unused forcing field "dQdT" and associated parameter "dQdTFile".
422      - move forcing adjustment (balancing surface forcing) in specific S/R
423        where might also be applied exchanges (if needed) to forcing arrays.
424    o pkg/shelfice:
425      - allow to specify the Ice-Shelf mass per unit area (instead of
426        shelficeLoadAnomaly) and do the conversion from one to the other
427      - move loading from file of shelficeLoadAnomaly from SHELFICE_INIT_VARIA
428        to SHELFICE_INIT_FIXED: (a) avoid some TAF recompution issues and
429        (b) will be loaded in time (before calling INI_PRESSURE) in case it is
430        needed in phi0surf for initial pressure calculation when useDynP_inEos_Zc
431      - update isomip experiment forward output (get machine truncation error
432        from conversion shelficeLoadAnomaly -> shelficeMass -> shelficeLoadAnomaly)
433    o verification/offline_exf_seaice:
434      revise data.seaice in input.thermo and input_ad and update results accordingly
435        HO = 0.2, SEAICE_areaLossFormula = 2,
436        SEAICE_growMeltByConv = .TRUE.,
437        SEAICE_doOpenWaterMelt = .TRUE.,
438      Improve consistency with input.thsice and fix the spurious 500m ice corners
439    o verification/global_ocean.cs32x15:
440      - set Energy Reference Level to liq. water @ 0 degC: temp_EvPrRn=0.,
441        in icedyn (thsice+seaice-dyn) test and update output.
442    o verification/offline_exf_seaice:
443      - adjust list of tested variables (with testreport) and update
444        standard test output with better slection of monitored fields.
445    o pkg/cd_code:
446      - swtich off Adams-Bashforth on the 1rst iteration. Update output of
447        experiment lab_sea.hb87
448    o tools/genmake2:
449      - change default makedepend from the system default to tools/xmakedepend:
450        When no makedepend is specified, select tools/xmakedepend it the test pass,
451        otherwise select the system default (if exist and pass the test) or
452        try to build cyrus-makedepend ; stop if none of the 3 works.
453    o pkg/exf & pkg/thsice:
454      - In case Energy-Reference-Level (ERL) is used (temp_EvPrRn=0), account
455        for energy content of Precip + RunOff & Evap. Assumes:
456       1) Rain has same temp as Air.
457       2) Snow has no heat capacity (consistent with seaice & thsice pkgs)
458       3) Run-Off comes at the temp of surface water (with same Cp)
459       4) Evap over ocean is released to the Atmos @ surf-temp (=SST); should be
460        put directly into Latent Heat flux & should use water-vapor heat capacity.
461       5) Evap over snow/seaice (pkg/thsice) assumes constant Lf+Lv (consistent
462        with snow having no Cp) => water vapor is released to Atmos @ 0^oC (=ERL).
463    o pkg/seaice: simplify the use of CPP flags in JFNK-related files
464        (compile when SEAICE_ALLOW_JFNK is defined)
465    o pkg/thsice:
466      - fix oceanic flux diagnostics & timeave (left from new sequence of calls
467        implementation): move THSICE_AVE call after the advection (called directly
468        from thsice_do_advect.F).
469    o verification/global_ocean.90x40x15 and OpenAD:
470      New experiment for OpenAD is identical to the one in input_ad except that
471        multiDimAdvection=.FALSE.,
472      There is still an issue with this.
473    o pkg/thsice:
474      - account for ice velocity (ocean relative velocity) in turbulent heat flux
475        at the base of the ice; update output of global_ocean.cs32x15 (icedyn)
476        and offline_exf_seaice (main test).
477    
478    checkpoint64f (2013/03/31)
479    o verification/OpenAD:
480      - new (smaller, without OAD_DEBUG) output_oadm.txt generated using gfortran.
481    o verification/isomip, pkg/shelfice and OpenAD:
482      - add OpenAD setup to this exp and make shelfice package work with OpenAD
483    o verification/tutorial_global_oce_biogeo:
484      - update AD and TLM output after turning off ALLOW_FE (in DIC_OPTIONS.h)
485    o verification/offline_exf_seaice:
486      - replace adjoint experiment (now similar to fwd input.thermo/ setup).
487    o verification/tutorial_tracer_adjsens, pkg/ptracers and OpenAD:
488      - add OpenAD setup to this exp and make ptracers package work with OpenAD
489    o pkg/seaice/seaice_evp.F:
490       mask the asymmetric contribution of the ice-ocean stress (when
491       the turningAngle is non-zero) over open water to avoid stripes
492       in the EVP-solutions after many EVP-substeps. This does not change
493       the verification experiment (because turningAngle=0 in that experiment),
494       but will affect other EVP solutions with turningAngle.ne.0
495    o verification, AD & TLM output:
496      - new TAF version 2.3.10 (previously 2.3.8) add more parenthesis in multiple
497        additions and multiplications which affect results at machine truncation
498        level. Update few AD output (1D_ocean_ice_column, bottom_ctrl_5x5,
499        lab_sea.noseaicedyn & offline_exf_seaice) and TLM output
500        (1D_ocean_ice_column & tutorial_tracer_adjsens.som81) which had only 12
501        digits of matching.
502    o pkg/seaice:
503       modify lsr-algorithm so that it vectorizes AND converges by using a
504       "zebra" method: split vector loop into two with step size 2; this
505       solves a general convergence problem for the JFNK solver on a vector
506       computer. It has the potential to even improve convergence for
507       non-vector machines, but it does change the results so it is turned
508       off by default. To turn on define CPP-flags SEAICE_VECTORIZE_LSR and
509       SEAICE_VECTORIZE_LSR_ZEBRA in SEAICE_OPTIONS.h
510    
511    checkpoint64e (2013/03/05)
512    o pkg/ecco and several verification experiments:
513       Remove special versions of the_model_main.F, the_main_loop.F, forward_step.F
514       from ecco package, and revert to unified routine in model/src/
515       Adjust verification experiments accordingly, especially gloal_oce_llc90
516       Should conclude long-standing work by Gael to reconsile the_main_loop with
517       ECCO and will enable OpenAD's revolve loop in ECCO setup.
518    o pkg/generic_advdiff:
519      - implement compressible flow method for multi-dim advection (similar to
520        gad_som_advect.F); add new option "GAD_MULTIDIM_COMPRESSIBLE" (since
521        TAF generates many recomputations) to use it; no yet coded with implicit
522        vertical advection.
523      - change experiment advect_cs (new output) and advect_xy (same output)
524        to test compressible flow multi-dim advection.
525    o pkg/generic_advdiff:
526      - fix SOM advection for unusual cases (linear Free-Surf in r*, other
527        cases where linear Free-Surf is not @ k=1, e.g., useShelfIce)
528    o pkg/seaice:
529      - fix bug in the no-slip implementation in seaice_calc_strainrates.F, this
530        bug affects two verification experiments (lab_sea and lab_sea.hb87),
531        recover old (wrong) results by setting noSlipFac = 2. in
532        seaice_calc_strainrates.F;
533      - add another averaging method for etaZ (etaZmethod=3), a mix of 2 and 0
534      - add runtime parameter SEAICE_JFNKepsilon for that finite difference
535        gradient in the Jacobian times vector operation (s/r seaice_jacvec)
536      - fix line search in seaice_jfnk.F (s/r seaice_jfnk_update)
537    o pkg/seaice: reorganize seaice_lsr.F to be able to reuse more code:
538      - new subroutines seaice_lsr_tridiagu/v
539      - introduce loop bounds i/jmin/max to be able to use a restricted
540        additive Schwarz method, but not used for now.
541      - rename SEAICE_CALC_LSR_COEFFS into SEAICE_LSR_CALC_COEFFS
542      - use new subroutines seaice_lsr_tridiagu/v in s/r seaice_preconditioner
543    o OpenAD support
544       Reorganize S/R call for overall time-stepping loop to facilitate
545       OpenAD revolve loop replacement
546       (new S/R main_do_loop.F)
547    o model/src:
548      - add 3-D array 3rd dimension to argument list of S/R ADAMS_BASHFORTH2/3;
549      - pkg/cheapaml: remove local S/R ADAMS2d and use instead standard S/R
550        ADAMS_BASHFORTH2.
551    o pkg/cheapaml:
552      - add new parameter (default=F, same as before) to disable full domain
553        (including land) periodicity in X & Y dir.
554      - new S/R CHEAPAML_COPY_EDGES to copy relaxation field into edges of
555        prognostics field ; merge the 2 S/R cycle_2d_tracer and timestep_2d_tracer
556        into single S/R CHEAPAML_TIMESTEP.
557    
558    checkpoint64d (2013/02/17)
559    o pkg/monitor (angular momentum):
560      - do not add Eta contribution to the total AM if using rigid-lid;
561      - correct for Adams-Bashforth in zonal-wind part (instead of in AM_eta).
562      - horizontal discretisation: compute Zonal wind contribution separately for
563        each component (instead of from cell centered averaged) and use the 4
564        grid-cell corner averaged value of Omega*(R*cos_Phi)^2 in AM_eta part.
565      - Check that AM is conserved with linear dynamics (momAdvection=F and linear
566        FreeSurf) if using vectorInvariant momentum.
567    o model/src (horizontal grid):
568      - move initialisation of all horizontal grid arrays to null value
569        from specific horizontal grid S/R to main S/R INI_GRID;
570      - add 2 new arrays to project model C-grid velocity on to zonal direction
571        (similar to AngleCosC & -AngleSinC but at U and V point location) and
572        compute them in S/R CALC_GRID_ANGLES.
573    o pkg/seaice: replace a hard coded parameter by a runtime parameter
574        SEAICE_JFNK_tolIter
575    o pkg/seaice: introduce a restricted additive Schwarz (RAS) method in
576        preconditioner for JFNK (run LSR into the overlaps as far as possible,
577        i.e OLx/y-2, but overwrite the solution in the overlaps)
578      - extend loop range in seaice_calc_viscosity
579      - add new parameters SEAICE_OLx/y (default = OLx/y-2)
580      - add flexible loop ranges to seaice_preconditioner, seaice_calc_lsr_coeffs
581      - slightly modify loops in tridiagonal solver lsr-preconditioner
582      - by default remove the extra exchange in seaice_preconditioner with a
583        replacement of an undocumented CPP flag by SEAICE_PRECOND_EXTRA_EXCHANGE
584      - modify offline_exf_seaice so that old results without RAS are reproduced
585        for output.dyn_jfnk.txt
586      in seaice_fgmres.F (unrelated to RAS)
587      - add a cpp-flag option for a re-ordered mapping (potentially more
588        efficient), but because this changes the results, turn it off by
589        default
590      - replace the save statement by local common blocks
591    o utils/matlab:
592      - fix tile-offset in rdmnc.m when output was produced using EXCH2.
593      - update cs_grid/bk_line scripts to use load_grid.m (in utils/matlab)
594        instead of local script "load_cs.m".
595    o pkg/seaice:
596      - fix a serious bug in seaice_fgmres.F, s/r seaice_scalprod, connected
597        to the scalarProductMetric and the unrolling. By chance this does not
598        change the verification experiment, because of the particular domain
599        decomposition
600    o pkg/diagnostics:
601      - change default missing_value from UNSET_RL to -999.
602    o pkg/layers: improve vectorization of layers_fluxcalc.F
603      - step 1: split loops, make TatU/V 2D arrays. This alone reduces the cpu
604        time spent in this routine by more than a factor of 3 on an NEC-SX8
605      - step 2: replace search algorithm by a bisection algorithm found in
606        Press et al. Numerical Recipes in Fortran. This algorithm can be
607        vectorized and finds the correct bin in approx. log2(Nlayers) steps.
608        For now, use this algorithm only for TARGET_NEC_SX defined.
609      - moved the search algorithm into a separate subroutine, that is kept
610        within the same file for automatic inlining.
611      - modified catching potential errors
612    o pkg/ecco:
613      - add compile switches ALLOW_SHALLOW_ALTIMETRY and ALLOW_HIGHLAT_ALTIMETRY.
614    o model/src, pkg/autodiff, pkg/ctrl, pkg/grdchk :
615      - modifications for ALLOW_GENTIM2D_CONTROL
616      -- forward_step.F, the_main_loop.F : pass options and arrays for adjoint tapes.
617      -- load_fields_driver.F : add call to CTRL_MAP_GENTIM2D( myTime, myIter, myThid )
618      -- add xx_gentim2d0/1 in ctrl_ad_check_lev1_dir.h checkpoint_lev1_directives.h etc
619      -- use ctrl_swapffields.F (new) instead of exf_SwapFFields.F in ctrl_get_gen.F
620    o model/src:
621      - account for moist thermodynamics in (dry) potential temperature equation;
622        update results of experiments: aim.5l_cs (x2), fizhi-cs (std + aqualev)
623        and cpl_aim+ocn.
624      - also switch to standard compiler (gfortran, was using ifort previously)
625        for fizhi-cs experiments reference output.
626    o pkg/ctrl:
627      - 1st step toward generic time-varying control variable (2D) via CPP option
628        ALLOW_GENTIM2D_CONTROL
629      - note new namelist ctrl_nml_genarr in data.ctrl
630    o verification/global_ocean.cs32x15:
631      - switch off OLD_THSICE_CALL_SEQUENCE (both forward and AD code) to use new
632        thsice calling sequence and update "output.icedyn.txt" output file.
633    o pkg/ptracers:
634      - from Ryan A.: code to allow to restart a tracer at regular time interval
635        specifying a reset frequency and a reset-phase (for each tracers).
636    o pkg/thsice:
637      - implement new sequence of calls for thsice+seaice:
638          previously:   ice-Dyn,ice-Advect,ice-Thermo(thsice)
639          new sequence: ice-Thermo(thsice),ice-Dyn,ice-Advect
640      - allows (with temporary CPP option "#define OLD_THSICE_CALL_SEQUENCE"
641        in CPP_OPTIONS.h) to recover old sequence;
642      - update diagnostics for Intermediate State thsice diagnostics: now filled
643        just before advection (instead of just after advection in old sequence);
644      - update output of experiment offline_exf_seaice (standard test);
645      - in global_ocean.cs32x15 exp., postpone output.icedyn.txt update
646        by setting #define OLD_THSICE_CALL_SEQUENCE in CPP_OPTIONS.h
647    o pkg/aim_v23:
648      - use tiled array for aim_CO2 and surface wind speed for export to coupler;
649      - use standard pkg/mdsio routines to read & write CO2 pickup file;
650      - in-line S/R AIM_AIM2DYN_EXCHANGES in aim_do_physics.F;
651      - Implement new sequence of calls for seaice advection (with thsice and
652        seaice diffusion): do ice advection after ice-Thermodynamic (previously
653        was done before); change aim_do_physics.F (split the bi,bj loops block);
654      - use same CPP option (OLD_THSICE_CALL_SEQUENCE) to recover old sequence;
655    o pkg/thsice:
656      - add a missing EXCH of sIceLoad in thsice_main.F (ocean set-up) for the case
657        where just seaice diffusion is used without seaice-dynamics.
658    
659    checkpoint64c (2013/01/18)
660    o pkg/seaice:
661      - add a parameter JFNKres_tFac that can replace JFNKres_t
662        (= JFNKres_tFac*JFNKresidual before JFNK solver)
663      - add a line search option to the JFNK solver
664      - move update of u/vIce and computation of initial residual into a separate
665        routine (same file); makes the code a little easier to read (?)
666      - the line search requires computation of the residual and for some reason
667        I cannot do the update of u/vIce at the beginning of the Newton loop
668        without changing the results, so we move computation of the new residual
669        also to the end of the Newton loop, even without the line search
670    o pkg/mdsio:
671      - write missing value (corresponding to undefined array value) to meta files;
672        skip this if missing value is one (UNSET_RL would be a better value for the
673        case of undefined missing value but it's current the default missing-value);
674        Also update mdsio_read_meta.F (reading timeInterval & missing-value).
675    o pkg/exf & thsice:
676      - for diagnostics purpose, fill-up and update Qsw array (+ fill the
677        diagnostics) even when SHORTWAVE_HEATING is #undef;
678    o verification/cfc_example:
679      - turn on pkg/layers with pkg/diagnostics (to test pkg/layers and to
680        illustrate how to use it).
681    o verification/offline_exf_seaice:
682      in input: add SEAICErestoreUnderIce=.TRUE., and reduce size of forcing
683        files to only one time-record.
684      in input.thermo: narrow the 2 bands of low initial ice-concentration near
685        N & S edges; switch on SEAICErestoreUnderIce=T; use pkg/diagnostics for
686        snapshot output; same mcPheePiston value as in thsice with
687        surf-current=0.2m/s.
688      in input.thsice: similar set-up as in input.thermo (channel with same seaice
689        initial conditions).
690    o verification/cfc_example:
691      - update file "data" with more standard choice (default rotationPeriod,
692        exactConserv=T, reduce ivdc_kappa to 10, remove multiDimAdvection=F,
693        set deltaTFreeSurf to deltaTtracer) and update the output.
694    o pkg/layers:
695      - fix bolus part so that it compiles with #undef GM_BOLUS_ADVEC;
696      - ensure that "layers_bolus" is switched to False when useGMRedi=False.
697    o verification/offline_exf_seaice:
698      - set all {field}period to zero for constant-in-time forcing (not tested
699        anywhere else).
700    o verification/offline_exf_seaice:
701      - dyn_lsr  : use pkg/seaice S/R (instead of pkg/thsice) to advect sea-ice;
702      - dyn_jfnk : use pkg/thsice S/R (instead of pkg/seaice) to advect sea-ice.
703    o model/src:
704      - implement Finite-Volume method for (hydrostatic) presure gradient
705        from S.-J. Lin (QJRMS 1997), for atmosphere using sigma-coordinate.
706    o pkg/seaice: implement EVP* by Lemieux et al. (2012).
707      - In order to be compatible with the default EVP code, this requires changing
708        the EVP code around a little; and because any small change immediately
709        breaks the lab_sea.hb87 verification experiment, the reference output is
710        updated two times (each time the remaining aggrement was only 2 digits):
711        - after replacing 1/dt by 1*recip_dt with a precomputed recip_dt
712        - after replacing 1/(1+0.5*dt/tau) by a/(a+1), with a=2*tau/dt
713        The adjoint test lab_sea.evp also needed to be updated.
714      - Turn on EVP* by setting SEAICE_deltaTevp as before, and setting
715        SEAICEuseEVPstar=.true.; for EVP* SEAICE_deltaTevp does not need to be
716        a factor of SEAICE_deltaTdyn, and you can set the number of EVP steps
717        explicitly with SEAICEnEVPstarSteps (thereby overriding the default of
718        INT(SEAICE_deltaTdyn/SEAICE_deltaTevp)). Controlling EVP* is deliberately
719        kept obscure, because I need to test it a little more.
720    o pkg/cheapaml:
721      - fix few incorrect type (_RL/_RS) in subroutine calls.
722    o pkg/seaice:
723      - new S/R SEAICE_MAP_RS2VEC (similar to SEAICE_MAP2VEC) to map _RS arrays
724        to _RL vector; this fixes JFNK code when compiling with "-use_r4" option.
725    o model/src:
726      - add a warning if using select_rStar=2 with nonlinFreeSurf <> 4;
727        change (in input_ad*/data) select_rStar from 2 to 1 for AD verification
728        experiments global_ocean.cs32x15 & tutorial_tracer_adjsens (+ update output)
729    o pkg/cheapaml:
730      - replace "273.1 _d 0" by celsius2K in pkg/cheapaml/cheapaml.F and
731        update output of experiment cheapAML_box.
732      - store in CHEAPAML.h coeffs ssq[0:2] for saturation specific humidity,
733        and set them only in one place (cheapaml_readparms.F).
734      - replace "273.16 _d 0" by celsius2K & replace "3.14159 _d 0" by constant "PI"
735      - remove re-setting of main-model parameter "celsius2K" (in cheapaml_readparms)
736        and re-setting of cheapAML params "lath" & "p0" (in cheapaml_init_varia).
737    o model/src/the_main_loop.F:
738      - fix revision 129 : need both ALLOW_ECCO and ALLOW_ECCO_EVOLUTION.
739    o pkg/seaice:
740      - remove legacy branch code.
741      - retire SEAICE_GROWTH_LEGACY and SEAICE_CAP_HEFF accordingly.
742    o verification/*/code*/SEAICE_OPTIONS.h:
743      - remove documentation for retired CPP flags (SEAICE_GROWTH_LEGACY,
744        SEAICE_CAP_HEFF).
745    o pkg/cheapaml:
746      - use standard routine (S/R GET_PERIODIC_INTERVAL) to decide when it is time
747        to load new forcing fields.
748      - add missing "_d 0" where it's needed; + some editing (improve indentation).
749      - remove switch to "-r8" compiler option (no longer needed) in verification
750        cheapAML_box, keep same output (small trucation error differences).
751    o verification/offline_exf_seaice:
752      - add thermodynamics only test.
753    o verification/cheapAML_box:
754      - split the domain in 8 (2x4) tiles; add MPI and multi-threaded test;
755        keep same output (same 14 digits for cg2d).
756      - switch on "-r8" compiler option (some numerical values are missing
757        the "_d 0") and update output (only 6 digits for cg2d were preserved).
758    
759    checkpoint64b (2012/12/24)
760    o pkg/cheapaml:
761      - remove un-used variables (Qnet,EmPmR) from pickup files.
762      - store tendency of advected variables into cheapaml pickup file
763        and fix the AB-2 restart.
764    o verification/offline_exf_seaice and lab_sea adjoint:
765      - switch to new code (#undef SEAICE_GROWTH_LEGACY).
766      - switch to 4 time steps in offline, consistent with lab_sea.
767      - match parameters from SEAICE_GROWTH_LEGACY branch in data.seaice.
768      - update results after change of experiment seaice.
769    o pkg/seaice:
770      - remove option SEAICE_MULTICATEGORY; has been replaced with run-time
771        parameter SEAICE_multDim. Issue a warning if SEAICE_MULTICATEGORY is
772        still defined since SEAICE_multDim needs to be set to MULTDIM value (from
773        SEAICE_SIZE.h) to recover previous #define SEAICE_MULTICATEGORY default.
774    o model/src:
775      - add run-time parameter to select CG2D solver minimum-residual solution
776        (in case of a poor convergence).
777    o pkg/seaice/seaice_lsr.F
778      - move calculation of matrix coefficients into new
779        S/R SEAICE_CALC_LSR_COEFFS; for now this routine remains in the same
780        file; the main purpose of this change is to be able to reuse the code
781        in S/R SEAICE_PRECONDITIONER
782      - modify the computation of the residual so that the mass of sea-ice
783        is not used twice; this simplifies the calculations, but since
784        the residual is never used except when LSR_mixIniGuess.ge.4---and this is
785        not tested in testreport---this does not affect the testreport results.
786    o verification/seaice_obcs:
787      - disable SEAICE_GROWTH_LEGACY; reset input files to match those of
788        verification/lab_sea/input.salt_plume; changes results/output.txt*
789    o model/src:
790      - adapt quasi-hydrostatic formulation to work with p or p* coordinate.
791    o pkg/seaice:
792      - improve description of parameters in seaice_fgmres.F and start
793        cleaning up the code
794      - change name of SCALPROD to SEAICE_SCALPROD and modify meaning of
795        second argument
796      - add a metric based on grid cell area to SEAICE_SCALPROD; update
797        results of offline_exf_seaice.dyn_jfnk
798      - compute initial residual norm in seaice_jfnk based on SEAICE_SCALPROD
799    o pkg/aim_v23:
800      - fix an other bug in dissipation heating (grid-cell area was missing).
801    o model/src + pkg/seaice:
802      - add global (multi-tile) linear solver for pair of tri-diagonal system
803        along X and Y lines, respectively for U and V component.
804        Note: 1) MPI and cube-exchange not yet coded.
805            2) probably not accurate for poorly conditioned / large size problem.
806      - add option for SEAICE_LSR to use global (multi-tile) tri-diagonal solver;
807        default: not compiled and not used.
808    o pkg/seaice:
809      - first attempt at a multi-threaded fgmres solver for the JFNK solver:
810        introduce bi/bj loops and use global_sum_tile_rl instead of
811        mpi_allreduce in seaice_fgmres.F
812      - improve diagnostics output for fgmres solver (seaice_fgmres.F)
813      - replace a few unnecessary goto statements by if-blocks  (seaice_fgmres.F)
814      - add a verification experiment for JFNK in offline_exf_seaice
815    o verification:
816      - change (for now, only fwd test) offline_exf_seaice experiment:
817        use a re-entrant channel with idealized coastline (trapezoidal shape)
818        at southern boundary to explore seaice dynamics and thermodynamics
819        with simple forcing. Main set-up has dynamics+thermo and SST relaxation;
820        secondary test (input.thsice) use thermodynamics only with SST relaxation.
821      - add seaice-dynamics only secondary test using LSR and pkg/thsice advection.
822    o pkg/exf:
823      - disable "sstExtrapol" when Nr=1
824    o verification:
825      - add new experiment "cheapAML_box" (simple 1 layer set-up from Bruno Deremble)
826        to test pkg/cheapaml.
827    o pkg/aim_v23 & pkg/shap_filt:
828      - fix dissipation heating (erroneous gravity multiplication factor).
829    o pkg/seaice:
830      - remove unused field pressPre from jfnk solver
831      - rename two subroutines in seaice_preconditioner.F
832    o eesupp/src:
833       - S/R MDSFINDUNIT avoids reserved unit numbers (scrUnit1, eeDataUnit ...)
834  o pkg/seaice/seaice_preconditioner.F: fix preconditioner once again  o pkg/seaice/seaice_preconditioner.F: fix preconditioner once again
835     - remove wrong parts of the rhs and the lhs (were in there     - remove wrong parts of the rhs and the lhs (were in there
836       because I did not understand the concept of this preconditioner)       because I did not understand the concept of this preconditioner)
Legend:
Removed from v.1.1774  
changed lines
  Added in v.1.1938
  ViewVC Help
Powered by ViewVC 1.1.22