/[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.1836 by jmc, Mon Mar 4 19:08:25 2013 UTC revision 1.1934 by jmc, Fri Oct 25 18:44:25 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 eesupp/src & pkg/mdsio:
8      - remove gather/scatter_vector.F (from eesupp) and mdsio_read/writevector.F
9        (from pkg/mdsio) since they are no longer used for tape IO.
10      - move simple C subroutine "setdir.c" from pkg/compon_communic to eesupp/src
11        (+ keep a copy in pkg/atm_ocn_coupler).
12    o verification/isomip:
13      - split domain in 8 squared tiles (instead of just 2) for FWD, AD & OAD tests.
14    o tools/build_options:
15      - optfile linux_amd64_ifort11, fast (no ieee) case:
16        to get a working restart for all experiments, add "-fp-model source"
17        to FOPTIM (also remove "-opt-streaming-stores auto" <- this is the default).
18    
19    checkpoint64p (2013/10/20)
20    o pkg/autodiff:
21      - adread_adwrite.F: when type and length do not match ALLOW_AUTODIFF_WHTAPEIO
22        criteria, return to simple (default) tape IO format.
23    o pkg/mdsio (+ eesupp/src):
24      - rewrite/simplify tape-IO S/R (new S/R MDS_READ/WRITE_TAPE replace
25        previous MDSREAD/WRITEVECTOR) with 2 array argument from each type (R4/R8);
26      - fix globalFile and singleCpuIO options using simpler global mapping.
27      - new pairs of scatter/gather vector S/R (one for each type, generated from
28        template) corresponding to simpler global mapping.
29      - adread_write.F: switch to new S/R (but leave old calls commented) and write
30        to globalFile when parameter globalFiles is set to True.
31      - switch on useSingleCpuIO in AD experiments: global_ocean.90x40x15 & isomip.
32    o pkg/autodiff:
33      - adread_adwrite.F: add a hack for RS arrays (when RS=real*4);
34        with ALLOW_AUTODIFF_WHTAPEIO defined, since corresponding (real*4) is
35        missing, just add a STOP if array type do not match.
36    o pkg/exf:
37      - to allow to use seaice-fraction from pkg/exf in other pkgs than pkg/seaice:
38         rename CPP option EXF_ALLOW_SEAICE_RELAX to EXF_SEAICE_FRACTION
39         and rename array obsSIce to exf_iceFraction.
40    o tutorial_deep_convection:
41      - add 1rst version of isotropic 3-D Smagorinsky code (from Louis-Philippe),
42        for now all in this exp. code dir.
43      - add secondary test using Smag-3D viscosity (keep very low backgroud visc for
44        testing) and flux-limter advection scheme for Temp. (no explicit diffusion).
45    o eesupp/src:
46      - added SINGLE_DISK_IO option because reading of parameter files on
47        very large number of processors did not work with existing set-up
48    
49    checkpoint64o (2013/09/25)
50    o pkg/atm_phys & atm_compon_interf:
51      - allows for coupling of gray-atmosphere (atm_phys) with ocean.
52    o tools/build_options:
53      - switch from "-fPIC" to "-mcmodel=medium" in optfiles: linux_amd64_gfortran
54        & linux_amd64_ifort11 and in linux_amd64_g77 & linux_amd64_open64 also
55        remove CPP setting (was identical to genmake2 default).
56    o pkg/gmredi & model/src: Some updating and work done on the residual model.
57      - created new file model/src/calc_eddy_stress.F to calculate the eddy stress.
58      - fixed up the calculation of the eddy stress and how it is added to the
59        external forcing via taueddy_external_forcing.F
60      - at the moment, it is very experimental and only works with the variable eddy
61        diffusivity calculated in gmredi_k3d.F
62    o pkg/seaice/seaice_check.F: add a check to catch the case of specifying
63        lw/swfluxfile instead of lw/swdownfile in data.exf together with pkg/seaice
64    o testreport (& tools/do_tst_2+2):
65      - for MPI test, set default COMMAND to 'mpirun -v -np TR_NPROC ./{executable}'
66    o tools/build_options: linux_amd64_gfortran & linux_amd64_ifort11:
67      - remove CPP='cpp -traditional -P' (this is genmake2 default)
68      - add switch (for big setups) to use 'medium' memory-model (but leave it
69        commented out for now). Note: needed with FC 19 to compile big setups.
70    o pkg/atm_compon_interf:
71      - regroup per field origin (model & pkgs) the coupling storage routines
72         (which accumalate in time each coupling field)
73      - move bi,bj loops inside atm_store_my_data.F and store also wind-stress;
74        update accordingly + simplify aim_do_physics.F
75    
76    checkpoint64n (2013/08/26)
77    o pkg/bbl:
78      - update formulation to reduce dispersion of tracer properties in bbl
79      - replaced horizontal and vertical transport parameters with velocities
80      - changes results for verification/global_with_exf.yearly
81    o verification/halfpipe_streamice:
82      - change to matrix construction; update both Fwd & AD output.
83    o verification/tutorial_tracer_adjsens:
84      - OpenAD test: switch to vector-invariant momentum since fluxform with
85        rStar causes problems with OpenAD; update output.
86    o tools:
87      - genmake2: define a common prefix for all temp-file names ; move them
88        to local /tmp disk (but would be easy to move them back);
89      - xmakedepend: apply the same switch for temp-files.
90    o pkg/offline:
91      - changes for use with rStar (NLFS): assume that we continue to provide
92        UVELMASS/hFacW_ini & VVELMASS/hFacS_ini (hFacWi,S_ini as written in grid
93        output files) and apply r* scaling at each time-step (function of the time
94        varying Eta field); use normal pickup files (needed for restart with NLFS).
95    o adjoint with OpenAD:
96      - move OpenAD-specific files from verification/OpenAD/code_oad_all to a new
97        package openad.  This package has to be compiled for all runs with OpenAD
98        instead of specifying code_oad_all as a modification directory.
99      - link the fortran templates from tools/OAD_support to the build directory
100        for preprocessing to avoid conflicts between different setups compiled
101        simultaneously and simplify cleanup.
102    o pkg/exf:
103      - in exf_getforcing.F: comment out most of the EXCH calls (fields used by
104        main model or other pkgs are exchanged in EXF_MAPFIELDS).
105    o model/src:
106      - do_stagger_fields_exchanges.F: skip EXCH calls when using pkg/offline
107        (similar to do_fields_blocking_exchanges.F).
108    o tools/genmake2:
109      - remove TAF output src code just before calling staf (when not doing a fresh
110        built, will stop if things go wrong with staf).
111      - fix target "adobj" (left from TAF version-2 transition)
112    o pkg/diagnostics:
113      - add a parameter to track the status of pkg/diagnostics activation;
114        this replace/extend the use of logical param "settingDiags" (to check
115        when adding diag to the list) to also check any DIAGNOSTICS_[]FILL* call.
116        New S/R to print error message and stop when pkgStatus is not right.
117      - to try: add a warning when "data.diagnostics" exist but useDiagnostics=F.
118    
119    checkpoint64m (2013/08/13)
120    o model/inc:
121      - move FCORI_MACROS.h from PARAMS(_MACROS).h to GRID_MACROS.h (left from
122        8 yrs ago modif when Coriolis arrays were moved from PARAMS.h to GRID.h)
123    o pkg/kpp & aim_v23:
124      - only call DIAGNOSTICS_FILL when useDiagnostics=TRUE.
125    o pkg/seaice:
126      - restore (last time on 2011/12/24) ability to compile pkg/seaice with
127        undefined pkg/exf ALLOW_ATM_TEMP or ALLOW_RUNOFF CPP-options.
128      - only call DIAGNOSTICS_FILL when useDiagnostics=TRUE.
129    o eesupp/src:
130      - fix check for overlap-size argument (myOLn.GT.MAX_OLY_EXCH)
131        in exch1_rx[_ad].template (thanks to Francois Lefeuvre report).
132    o switch to TAF version-2:
133      - add flow directives to "freeze" naming convention for ad-routines that
134        are used in I/O routines of adjoint state and associated common blocks
135        that are used in adcommon.h (backward compatible with version-1)
136      - update all ad-optfiles (in tools/adjoint_options + few local to some
137        verification exp) to switch to latest TAF version (currently 2.5.5);
138        merge "*_mpi" adoptfiles into the corresponding non-mpi (but commented
139        out except for diva).
140      - the major change with TAF version-2 is that both the differentiated
141        routines and the forward routines (new) are returned as part of the TAF
142        output source code. Update genmake2 accordingly.
143      - new output for AD experiment tutorial_global_oce_biogeo (better AD gradient
144        with version-2 when compared to TLM or FWD gradient); other experiments
145        output show minor changes (machine truncation level); update the few AD &
146        TLM output with 12 matching digits only (some prior to version-2 switch).
147    
148    checkpoint64l (2013/08/05)
149    o model/src:
150      - in dynamics.F, skip the call to CALC_VISCOSITY if momViscosity=F
151      - main_do_loop.F: no longer includes additional (AUTODIFF) header files
152        when using TAF but only if using OpenAD.
153    o pkg/mom_common+fluxform+vecinv:
154      - always set horiz. viscosity arrays to background value before calling
155        MOM_CALC_VISC (in MOM_FLUXFORM & MOM_VECINV) and call S/R MOM_CALC_VISC
156        only when using variable horiz. viscosity (useVariableVisc=T);
157      - simplify mom_vecinv.F (only 1 block for momViscosity).
158    o pkg/mom_common: variable horiz. viscosity in NH code
159      - move horiz-viscosity 3-D arrays viscAh_W & viscA4_W (used in calc_gw.F)
160        from model/inc/NH_VARS.h to pkg/mom_common/MOM_VISC.h
161      - store in common block (in MOM_VISC.h) logical parameters: useHarmonicVisc,
162        useBiharmonicVisc & useVariableVisc.
163      - fix old bug (already in 1rst implementation of variable horizontal viscosity
164        in calc_gw.F, version 1.27, 2006/06/20):
165        horizontal bi-harmonic viscosity in gW was missing when wiscA4W=0
166        even if variable bi-harmonic viscosity (e.g., viscA4Grid) was non-zero.
167      - change background horiz viscosity in vertical momentum (gW) to always be
168        viscAhW,viscA4W (instead of switching to viscAhD, viscA4D when using
169        any type of variable viscosity, e.g., viscAh/4Grid)
170      - change default value of viscAhW,viscA4W to viscAhD,viscA4D (instead of
171        viscAh,viscA4) which default to viscAh,viscA4 :  this avoid missing
172        horiz. viscosity in W eq when just setting viscAhD/Z,viscA4D/Z (since
173        default viscAh,viscA4 is zero).
174    o model/src:
175      - extend length of drC & recip_drC from Nr to Nr+1.
176    o Divided Adjoint (DIVA) with MPI:
177      - genmake2: do not create local ./mpi_headers dir when no file to copy into.
178      - remove option ALLOW_DIVIDED_ADJOINT_MPI (replaced by ALLOW_DIVIDED_ADJOINT
179        and ALLOW_USE_MPI) from model/src code and from ECCO_CPPOPTIONS.h.
180    
181    checkpoint64k (2013/07/23)
182    o tools/OAD_support:
183      - avoid the need for file "oadtempflile" in insertTemplateDir.bash
184        and remove it from genmake2.
185    o tools/build_options:
186      - remove MPI_HEADER_FILES_INC setting (no longer used) and comment out
187        MPI_HEADER_FILES list (let genmake2 find the missing headers)
188    o verification/halfpipe_streamice:
189      - change to driving stress formulation; update both Fwd & AD output.
190    o tools/genmake2 and Divided Adjoint with MPI:
191      - remove default setting of MPIINCLUDEDIR from $MPIHOME (always been broken);
192      - only create local ./mpi_headers dir if using MPI and MPIINCLUDEDIR is set;
193      - add new variable (LOCAL_MPI_HEADERS, replace MPI_HEADER_FILES_INC) in
194        Makefile, which is set by genmake2, and check which mpi-headers are needed.
195        Note: need to update adjoint_options/adjoint_diva_mpi accordingly.
196      - update target "CLEAN" (to also remove ./mpi_headers/ dir).
197    o testreport:
198      - changes for Divided Adjoint Run: when file "run_ADM_DIVA" is present,
199        perform additional runs (number of additional runs is "add_DIVA_runs" taken
200        from file run_ADM_DIVA ; should match "nchklev_3" in file "code_ad/tamc.h").
201    o verification/lab_sea (AD test):
202      - change to test Divided Adjoint Run in this set-up (both standard AD test
203        and 2 secondary AD test: noseaice + noseaicedyn); drop the evp AD test.
204    o pkg atm_ & ocn_compon_interf:
205      - remove cpl_earlyExpImpCall=F option (was mainly a hack added for the old
206        staggerTimeStep implementation) and disable (retired) corresponding
207        parameter.
208    o pkg/offline:
209      - comment out unused field-array: ICEM, ConvectCount, hflx[01], icem[01];
210      - add S/R OFFLINE_CHECK to print offline parameter summary
211        and stop if loading a file that is not used;
212      - enable to calculate GMRedi tensor, KPP mixing and convective index
213        when the corresponding file(s) is/are not specified (but pkg is used).
214      - comment out resetting of monitorFreq (to endTime) since it still useful
215        for monitoring CFL numbers and checking offline fields.
216    o OpenAD:
217      - do not define OAD_DEBUG (in tools/adjoint_options/adjoint_oad) to get
218        smaller standard output) and generate new (smaller) output (on baudelaire
219        using gfortran with -ieee) for exp. tutorial_tracer_adjsens
220    o model/src:
221      - call EXTERNAL_FORCING_SURF outside bi,bj loops (in do_oceanic_phys.F),
222        (where external_forcing_adjust was called) and change list of arguments;
223      - merge external_forcing_adjust.F into external_forcing_surf.F
224      - put surface relaxation of SST & SSS in new S/R (forcing_surf_relax.F)
225        which also contains balancing surface relaxation (ALLOW_BALANCE_RELAX,
226        previously in file "balance_relax.F", now removed); remove corresponding
227        arrays in FFIELDS.h.
228      - EXTERNAL_FORCING_SURF is now called when using pkg/offline; implement
229        a hack for getting back surfaceForcingS (loaded from file) when using
230        ALLOW_OLD_VIRTUALFLUX (pkg/dic, in DIC_OPTIONS.h) in offline set-up.
231      - update output of AD test experiment global_ocean.cs32x15.seaice_dynmix
232        (was getting only 11 digits for AD-grad, due to different truncation
233         error in AD code).
234    
235    checkpoint64j (2013/07/04)
236    o pkg/offline:
237      - to enable the use of other features with pkg offline (e.g., obcs),
238        also call S/R DO_OCEANIC_PHYS in offline run - more modifs to come.
239      - comment out "allowFreezing" (now active) in 2 offline testing exp.
240    o pkg/exf ALLOW_RUNOFTEMP
241      - fix sign and make compatible with temp_EvPrRn .NE. UNSET_RL
242      - ahdded RUNOFTEMP test to verification/global_ocean.cs32x15/input.seaice
243      - changes verification/global_ocean.cs32x15/results/output.seaice.txt
244    o pkg/atm_phys:
245      - allows to add some damping on uVel,vVel (strength only function of k).
246    o verification, ad-tests:
247      - reduce length of the simulation for AD test: global_ocean.cs32x15.seaice,
248        global_ocean_ebm and offline_exf_seaice.thsice.
249    o verification/halfpipe_streamice:
250      - found a bug in s/r pkg/streamice/streamice_advect_thickness which caused
251        flux input boundary conditions to be ignored; a fix for this bug changed
252        verification results for forward test.
253      - add adjoint test.
254    o pkg/ptracers:
255      - allow to restart without SOM pickup (resetting SOM to zero).
256    o pkg/gmredi:
257      - New parameterisation for the calculation of K and (GM_PsiX,GM_PsiY) added.
258      - It uses mixing length theory and suppression by steering level effects to
259        calculate a diffusivity used in the Redi tensor (see Bates et al., submitted
260        to JPO, for more details).
261      - It also uses the diffusivity to calculate an eddy induced streamfunction
262        based on a PV closure (rather than a buoyancy closure as in GM).
263      - New files added to pkg/gmredi:
264        gmredi_k3d.F calculates the diffusivity for the Redi tensor and the eddy
265          induced overturning streamfunction (called from gmredi_calc_tensor)
266        gmredi_calc_eigs.F calculates the eigenvalues and eigenvectors required
267          for use in gmredi_k3d
268        gmredi_calc_urms.F calculates the the rms eddy velocity used to calculate
269          the eddy diffusivity
270    
271    checkpoint64i (2013/06/18)
272    o pkg/cheapaml:
273      - take src files from MITgcm_contrib/verification_other/offline_cheapaml/code
274        that enable the use of seaice (pkg/thsice thermo & pkg/seaice dynamics)
275        with pkg/cheapaml.
276    o utils/python: add subpackage for cubed sphere function, first example is
277        simple cs.pcol to plot 2D fields on a map
278    o pkg/streamice:
279      - This is a new package in MITgcm, which solves dynamic equations for ice
280         sheets and ice streams
281      - At the moment it is independent of atmosphere/ocean dynamics, or any
282        other physical packages in MITgcm, even ShelfIce (although it is anticipat-
283        ed that this will change)
284      - There is a new verification exercise that tests it, halfpipe_streamice,
285        evolves a 40x80 km ice stream/ice shelf (at coarse resolution) in a
286        rectangular embayment
287      - halfpipe_streamice is forward-only for now, will have TAF test soon
288      - for questions, email dgoldber@mit.edu
289    o pkg/thsice:
290      - fix filling of snow-precip diagnostic (SIsnwPrc) (left from the removing
291        of snowPrc from common bloc, on May 2nd, thsice_step_fwd.F revision 1.40);
292      - remove diagnostics 'SI_FrcFx' ( now identical to 'SI_Fract' with new
293         call sequence: ice-Thermo(thsice),ice-Dyn,ice-Advect );
294      - allow precip from pkg/cheapaml to make snow over seaice.
295    o pkg/cfc:
296      - move in separated header file (CFC_ATMOS.h) time-serie arrays of Atmos CFC
297      - remove atmos CFC 2-D arrays from common bloc (now local variables)
298      - change some hard-coded constant to run time parameter for reading and
299        and interpolating atmospheric CFC ; keep same default as previously.
300      - more flexible reading of Atmos CFC (number of records no longer hard-coded)
301      - put arguments in standard order: myTime, myIter, myThid
302      - use standard S/R GET_PERIODIC_INTERVAL for time interpolation of
303        atmospheric CFC ; fix for the case myIter < PTRACERS_iter0.
304      - fix calculation of surf tendencys for partial-cell and/or Non-Lin Free-Surf
305      - new routine to print summary of CFC params and check for consistency.
306    o pkg/seaice: add uice/vice to the parameter list of s/r seaice_advdiff,
307      similar to thsice_advdiff.
308      Note that for the B-grid case the routine parameters uc/vc are overwritten.
309    o pkg/cfc:
310      - change time interpolation for atmospheric CFC-11 (now same as for CFC-12);
311        update output (CFC-11) of experiment cfc_example and tutorial_cfc_offline.
312      - fix interpolation of atmphspheric CFC between Sorthern and Northern
313        hemisphere value.
314    o pkg/seaice:
315      - allow to use seaice dynamics without pkg/exf (e.g., using pkg/cheapaml)
316    o pkg/seaice: add Adams-Bashforth2 time discretization for ice dynamics;
317                  so far only for JFNK-solver, turn on with SEAICEuseAB2 = .TRUE.
318      - requires an additional time level that is stored in the pickup
319      - you can start with AB2 from a pickup without this time level with
320        pickupStrictlyMatch = .false.
321    
322    checkpoint64h (2013/05/26)
323    o pkg/cheapaml:
324      - replace one ()**.5 with SQRT + use temp var "cdu = tau/Du" (coare3_flux);
325        affect machine truncation (but exp cheapAML_box still pass with 13 digits).
326      - consistent treatment of uWind,vWind location: assume everywhere in
327        pkg/cheapaml that they are on C-grid (@ uVel,vVel location)
328        (Note: was already the case for Tair,Qair advection);
329        add option (INCONSISTENT_WIND_LOCATION) to recover old results.
330      - compute diabatic forcing only in the interior (excluding overlap regions)
331      - fill in array "Qsw" (for short-wave heating);
332      - add diagnostic for precip ; fill-in diagnostics for cheapAML state vars
333        before any update (consistent with MITgcm state var diagnostics);
334      - cheapaml_fields_load.F & cheapaml_init_varia.F : only fill outer edges if
335        non periodic ; + adjusted for uWind,vWind (C-grid).
336    o pkg/seaice:
337      - remove dependence of pkg/seaice versus pkg/exf (file: pkg/pkg_depend) and
338        add explicitly pkg/exf in experiments that currently use pkg/seaice.
339      - allows to compile pkg/seaice without pkg/exf.
340    o tools/OAD_support:
341      - use a local copy (ourRWCPinBd) to ensure index ourRWCP stays within bounds
342        of array "ourStepOf" in function rvNextAction (revolve.F90).
343    o pkg/aim_compon_interf renamed to pkg/atm_compon_interf
344      - change pkg name (to be used with other atmospheric physics pkg)
345        but by-pass CVS when making the new pkg so that we keep the CVS history.
346    o pkg/offline:
347      - call OFFLINE_FIELDS_LOAD from pkg/offline/offline_init_varia.F
348        (previously called from pkg/dic_surfforcing_init.F)
349    o model/src:
350      - re-order sequence of calls within packages_init_variables.F, using
351        standard order, except for:
352         DIAGNOSTICS_INIT_VARIA and OFFLINE_INIT_VARIA (called @ the top),
353         CTRL_MAP_INI and OBCS_INIT_VARIABLES (called @ the bottom)
354    o new pkg: atm_phys
355      - add entry point and source code (previously in contrib
356          verification_other/atm_gray/code_gray/) for new pkg "atm_phys"
357        (atmospheric physics pkg from P. O'Gorman and T. Schneider, JCl, 2008).
358      - re-order sequence of calls within packages_readparms.F,
359         packages_init_fixed.F and packages_write_pickup.F
360        (but skip packages_init_variables.F because of TAF issues).
361      - update contrib/*/atm_gray/ experiment.
362    o pkg/seaice (ITD code):
363      - simplify and fix the writing of pickup-files for ITD case
364    o pkg/seaice (ITD code):
365      - introduce floe size dependent lateral melt for ITD case;
366      - cleaning (avoid division by zero, change to f77 syntax, remove tabs)
367      - update seaice_itd experiment main output.
368    o pkg/thsice:
369      - remove snowPrc from common block THSICE_FLUX, in THSICE_VARS.h and
370        declared it locally (without bi,bj) in thsice_main.F & aim_do_physics.F;
371        changed accordingly pkg/atm2d.
372    
373    checkpoint64g (2013/05/01)
374    o pkg/exf ALLOW_RUNOFTEMP
375      - if defined,  river and glacier runoff temperature can be read-in from
376        runoftempfile; runoftempfile needs to have same format as runofffile
377    o pkg/thsice (& pkg/bulk_force):
378      - In case Energy-Reference-Level (ERL) is used (temp_EvPrRn=0), account
379        for energy content of Precip + RunOff & Evap (like with pkg/exf).
380      - change computation domain to always be the interior (without overlap).
381      - fix bug (wrong overlap) when using EXF and stressReduction > 0.
382    o verification/tutorial_global_oce_biogeo:
383      - add preliminary OpenAD setup (still has problems)
384    o pkg/seaice:
385      - turn some parameters for choosing the convergence criterion of
386        the inexact Newton method (JFNK) into runtime parameters for convenience
387      - clean up seaice_calc_viscosities.F: remove/retire two options for
388        SEAICetaZmethod (1,2), introduce area weighting for method 3; this changes
389        the results of lab_sea/input because a division has been replaced by
390        a multiplication with the reciprocal (sumNorm)
391      - modify offline_exf_seaice/input.dyn_jfnk/data.seaice to improve
392        convergence of solver (changes results, of course), and to test some
393        new features
394    o model/src:
395      - remove unused (since Aug 2006) S/R SET_WRITE_GLOBAL_PICKUP and its call.
396      - remove unused forcing field "dQdT" and associated parameter "dQdTFile".
397      - move forcing adjustment (balancing surface forcing) in specific S/R
398        where might also be applied exchanges (if needed) to forcing arrays.
399    o pkg/shelfice:
400      - allow to specify the Ice-Shelf mass per unit area (instead of
401        shelficeLoadAnomaly) and do the conversion from one to the other
402      - move loading from file of shelficeLoadAnomaly from SHELFICE_INIT_VARIA
403        to SHELFICE_INIT_FIXED: (a) avoid some TAF recompution issues and
404        (b) will be loaded in time (before calling INI_PRESSURE) in case it is
405        needed in phi0surf for initial pressure calculation when useDynP_inEos_Zc
406      - update isomip experiment forward output (get machine truncation error
407        from conversion shelficeLoadAnomaly -> shelficeMass -> shelficeLoadAnomaly)
408    o verification/offline_exf_seaice:
409      revise data.seaice in input.thermo and input_ad and update results accordingly
410        HO = 0.2, SEAICE_areaLossFormula = 2,
411        SEAICE_growMeltByConv = .TRUE.,
412        SEAICE_doOpenWaterMelt = .TRUE.,
413      Improve consistency with input.thsice and fix the spurious 500m ice corners
414    o verification/global_ocean.cs32x15:
415      - set Energy Reference Level to liq. water @ 0 degC: temp_EvPrRn=0.,
416        in icedyn (thsice+seaice-dyn) test and update output.
417    o verification/offline_exf_seaice:
418      - adjust list of tested variables (with testreport) and update
419        standard test output with better slection of monitored fields.
420    o pkg/cd_code:
421      - swtich off Adams-Bashforth on the 1rst iteration. Update output of
422        experiment lab_sea.hb87
423    o tools/genmake2:
424      - change default makedepend from the system default to tools/xmakedepend:
425        When no makedepend is specified, select tools/xmakedepend it the test pass,
426        otherwise select the system default (if exist and pass the test) or
427        try to build cyrus-makedepend ; stop if none of the 3 works.
428    o pkg/exf & pkg/thsice:
429      - In case Energy-Reference-Level (ERL) is used (temp_EvPrRn=0), account
430        for energy content of Precip + RunOff & Evap. Assumes:
431       1) Rain has same temp as Air.
432       2) Snow has no heat capacity (consistent with seaice & thsice pkgs)
433       3) Run-Off comes at the temp of surface water (with same Cp)
434       4) Evap over ocean is released to the Atmos @ surf-temp (=SST); should be
435        put directly into Latent Heat flux & should use water-vapor heat capacity.
436       5) Evap over snow/seaice (pkg/thsice) assumes constant Lf+Lv (consistent
437        with snow having no Cp) => water vapor is released to Atmos @ 0^oC (=ERL).
438    o pkg/seaice: simplify the use of CPP flags in JFNK-related files
439        (compile when SEAICE_ALLOW_JFNK is defined)
440    o pkg/thsice:
441      - fix oceanic flux diagnostics & timeave (left from new sequence of calls
442        implementation): move THSICE_AVE call after the advection (called directly
443        from thsice_do_advect.F).
444    o verification/global_ocean.90x40x15 and OpenAD:
445      New experiment for OpenAD is identical to the one in input_ad except that
446        multiDimAdvection=.FALSE.,
447      There is still an issue with this.
448    o pkg/thsice:
449      - account for ice velocity (ocean relative velocity) in turbulent heat flux
450        at the base of the ice; update output of global_ocean.cs32x15 (icedyn)
451        and offline_exf_seaice (main test).
452    
453    checkpoint64f (2013/03/31)
454    o verification/OpenAD:
455      - new (smaller, without OAD_DEBUG) output_oadm.txt generated using gfortran.
456    o verification/isomip, pkg/shelfice and OpenAD:
457      - add OpenAD setup to this exp and make shelfice package work with OpenAD
458    o verification/tutorial_global_oce_biogeo:
459      - update AD and TLM output after turning off ALLOW_FE (in DIC_OPTIONS.h)
460    o verification/offline_exf_seaice:
461      - replace adjoint experiment (now similar to fwd input.thermo/ setup).
462    o verification/tutorial_tracer_adjsens, pkg/ptracers and OpenAD:
463      - add OpenAD setup to this exp and make ptracers package work with OpenAD
464    o pkg/seaice/seaice_evp.F:
465       mask the asymmetric contribution of the ice-ocean stress (when
466       the turningAngle is non-zero) over open water to avoid stripes
467       in the EVP-solutions after many EVP-substeps. This does not change
468       the verification experiment (because turningAngle=0 in that experiment),
469       but will affect other EVP solutions with turningAngle.ne.0
470    o verification, AD & TLM output:
471      - new TAF version 2.3.10 (previously 2.3.8) add more parenthesis in multiple
472        additions and multiplications which affect results at machine truncation
473        level. Update few AD output (1D_ocean_ice_column, bottom_ctrl_5x5,
474        lab_sea.noseaicedyn & offline_exf_seaice) and TLM output
475        (1D_ocean_ice_column & tutorial_tracer_adjsens.som81) which had only 12
476        digits of matching.
477    o pkg/seaice:
478       modify lsr-algorithm so that it vectorizes AND converges by using a
479       "zebra" method: split vector loop into two with step size 2; this
480       solves a general convergence problem for the JFNK solver on a vector
481       computer. It has the potential to even improve convergence for
482       non-vector machines, but it does change the results so it is turned
483       off by default. To turn on define CPP-flags SEAICE_VECTORIZE_LSR and
484       SEAICE_VECTORIZE_LSR_ZEBRA in SEAICE_OPTIONS.h
485    
486    checkpoint64e (2013/03/05)
487    o pkg/ecco and several verification experiments:
488       Remove special versions of the_model_main.F, the_main_loop.F, forward_step.F
489       from ecco package, and revert to unified routine in model/src/
490       Adjust verification experiments accordingly, especially gloal_oce_llc90
491       Should conclude long-standing work by Gael to reconsile the_main_loop with
492       ECCO and will enable OpenAD's revolve loop in ECCO setup.
493  o pkg/generic_advdiff:  o pkg/generic_advdiff:
494    - implement compressible flow method for multi-dim advection (similar to    - implement compressible flow method for multi-dim advection (similar to
495      gad_som_advect.F); add new option "GAD_MULTIDIM_COMPRESSIBLE" (since      gad_som_advect.F); add new option "GAD_MULTIDIM_COMPRESSIBLE" (since
# Line 30  o pkg/seaice: reorganize seaice_lsr.F to Line 516  o pkg/seaice: reorganize seaice_lsr.F to
516    - rename SEAICE_CALC_LSR_COEFFS into SEAICE_LSR_CALC_COEFFS    - rename SEAICE_CALC_LSR_COEFFS into SEAICE_LSR_CALC_COEFFS
517    - use new subroutines seaice_lsr_tridiagu/v in s/r seaice_preconditioner    - use new subroutines seaice_lsr_tridiagu/v in s/r seaice_preconditioner
518  o OpenAD support  o OpenAD support
519    Reorganize S/R call for overall time-stepping loop to facilitate     Reorganize S/R call for overall time-stepping loop to facilitate
520    OpenAD revolve loop replacement     OpenAD revolve loop replacement
521    (new S/R main_do_loop.F)     (new S/R main_do_loop.F)
522  o model/src:  o model/src:
523    - add 3-D array 3rd dimension to argument list of S/R ADAMS_BASHFORTH2/3;    - add 3-D array 3rd dimension to argument list of S/R ADAMS_BASHFORTH2/3;
524    - pkg/cheapaml: remove local S/R ADAMS2d and use instead standard S/R    - pkg/cheapaml: remove local S/R ADAMS2d and use instead standard S/R

Legend:
Removed from v.1.1836  
changed lines
  Added in v.1.1934

  ViewVC Help
Powered by ViewVC 1.1.22