/[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.392 by mlosch, Tue Nov 23 20:22:12 2004 UTC revision 1.614 by jmc, Fri Feb 24 00:57:55 2006 UTC
# Line 1  Line 1 
1  $Header$  $Header$
2    $Name$
3    
4      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
5      ==============================      ==============================
6    
7    o 1rst implementation of Implicit IGW using the 3-D solver (use3Dsolver=T)
8      and based on a reference stratification. Tested within hs94.cs-32x32x5.
9    o seaice advection: add code for choosing between the default (2nd order
10      central differences) and multiDim advection schemes provided by
11      pkg/generic_advdiff
12    o seaice:
13      - removed all references to SEAICE_GRID.h fields from
14        dynsolver, advect, diffus, lsr, without changing the results of lab_sea
15        (cubed-sphere configuration will change a little)
16      - grid computations in seaice_init are obsolete now.
17      - cleaned up code a little so that is easier to read (hopefully).
18      - added new (better?) discretizations of default advection and diffusion
19        (turn on with SEAICEuseFluxForm)
20      - replace 1 in 3d global fields by kSurface in growth.F
21      
22    checkpoint58a_post
23    o seaice:
24      - bug fix in advect.F
25    o shelfice:
26      - separate shelfice load anomaly from pload (breaks with time
27        dependent forcing), => introduce constant field shelficeLoadAnomaly.
28        Its default is 0., but it may be computed more cleverly from (unknown)
29        t- and s-profiles (tRef, sRef) and actual EOS. For now this has to be
30        done offline. A good approximation of the pressure load anomaly is
31        necessary to avoid large initial adjustment processes underneath
32        deep-reaching shelfice.
33      - small fix in ISOMIP shelfice_thermodynamics
34      - add more sophisticated thermodynamics of Hellmer (1989), this is
35        now the default
36      - give a few constants more intuitive names
37    o add code to balance EmPmP and Qnet at the end of do_ocean_physics. Useful
38      if bulk formulae are used in long integration (especially EmPmR). Turn
39      on with balanceEmPmR = .true. or balanceQnet = .true. in data, PARM01
40      if balancePrintMean, the imbalance that is substracted is printed to
41      STDOUT.
42    o fix bug (typo) in obcs_apply_eta.F
43    o add a new package shelfice which enables us to model the circulation
44      underneath shelf-ice with draughts of hundreds of meters.
45     - works
46     - to do: + include improved thermodynamics
47              + check whether kSurfC (from model/inc/SURFACE.h) can be used
48                instead of new field kTopC
49              + can we use topoFile, Ro_surf and existing code in ini_depth.F
50                instead of new field R_shelfIce and new code in ini_depth.F
51                and ini_mask_etc.F?
52              + adjust vertical mixing schemes?
53    o seaice pkg:
54     - add a cpp-flag to include a simple formula for a salinity dependent
55       freezing point of sea water at the surface
56     - change the contents of SEAICE_TEST_ICE_STRESS_1 so that zonal and
57       meridional stresses are now different (and hopefully correct)
58    o AIM pkg: change threshold on large-scale condensation to conserve Heat.
59       update results of aim.5l_LatLon test ; fix title & units of AIM diagnostics
60       and  change 3 diagnostics names (OSR,OSRCLR,KM -> TSR,TSRCLR,DRAG)
61    o pkg/diagnostics: add options and parameters for regional-statistics
62       diagnostics ; define regions by reading a region-mask from a file
63       (might add other options later)
64    o pkg/bulk_force:
65     - add option and parameters to use AIM surface-flux formulae.
66     - Change loading part: S/R BULKF_FIELDS_LOAD only take care of bulkf_fields ;
67       others forcing fields (defined in PARM05, parameter file "data") are loaded
68       from S/R EXTERNAL_FIELDS_LOAD, as usually (with useBulkforce = T or F).
69     - initialise all bulkf_fields in bulkf_init.F ; do in-lining of exf_bulkcdn.F
70     - use the right EXCH call for uwind,vwind (to work on CS-grid).
71     - re-arrange header files (move parameters from BULKF.h to BULKF_PARAMS.h)
72       and parameters (note: calcWindStress replaces .NOT.readwindstress).
73     - cleaner interface with pkg/thSIce (new S/R THSICE_GET_PRECIP).
74    
75    checkpoint58
76    o add "--embed-source" option to genmake2 which, when the supporting
77      tools can be compiled (as determined by genmake2) will turn on the
78      embed_files package which then embeds the entire MITgcm source code
79      (*.[fFhc] + Makefile) used for the build within the executable
80      - requested by CNH but off by default
81      - adds a paltry <9MB to the mitgcmuv executable in most cases
82      - only writes the output when useEMBEDSRC in eedata is true
83    o clean up 3-dim. relaxation codes:
84      * Fix I/O inconsistency in pkg/rbcs: replace precFloat32 by readBinaryPrec
85      * Remove 3-dim. relaxation code from pkg/exf (now use only pkg/rbcs)
86      * Thanks to Tom Haine for testing!
87    
88    checkpoint57z_post
89    o make CG3D solver compatible with Free-surface at k > 1 (p-coordinate):
90      affects output (@ truncation level) of exp. plume_on_slope ; updated.
91    o fix the pickup behavior for obcs in combination with ptracers
92    o fix a bug in ptracers_external_fields_load
93    o no AB-2 at the first iteration for Gw (consistent with Gu,Gv,Gt,Gs ...)
94      affects results of NH experiments (but all our NH tests have Gw[it=1]==0 !!)
95    o transfered surface NH pressure to eta field (if exactConserv).
96    o gmredi: fix bugs in slope tapering ldd97 (depthZ had the wrong sign)
97         and update lab_sea output.txt file.
98    o completing 3-dim. relaxation code: additional term
99         1/tau*( theta(i,j,k) - thetaStar(i,j,k) )
100      in S/R EXTERNAL_FORCING_T (file external_forcing.F)
101      * At compile time, in "EXF_OPTIONS.h" set
102          #define  ALLOW_CLIMTEMP_RELAXATION
103          #define  ALLOW_CLIMSALT_RELAXATION
104      * At runtime in "data", under &PARM03 set (e.g., units in seconds)
105          tauThetaClimRelax3Dim = 15552000.0,
106          tauSaltClimRelax3Dim  = 15552000.0,
107      * 3 cases:
108          1) Use a constant-in-time 3-dim. field
109          2) Provide 12 3-dim. fields corresponding to Jan., ... , Dec.
110          3) Real-time relax. to time-evolving prescribed field
111    o Adding Laplacian-type smoothness constraint on
112      initial and boundary (forcing) controls (ic, bc)
113      * enable at compile-time via
114        ALLOW_SMOOTH_IC_COST_CONTRIBUTION, ALLOW_SMOOTH_BC_COST_CONTRIBUTION
115      * enable at run-time via mult_smooth_ic = 1, mult_smooth_bc = 1
116      * Absolute scale somewhat unclear so far
117      * no area weighting was included since no other cost term does so far
118    
119    checkpoint57y_post
120    First step for a NLFS adjoint
121    o initially suppress rStar (new flag DISABLE_RSTAR_CODE)
122    o new init. routines for calc_r_star, calc_surf_dr
123    o still need to deal with ini_masks_etc
124    o testreport seemed happy
125    checkpoint57y_pre
126    
127    o obcs: always reset wVel & eta (to zero or prescribe it) at the OBs.
128    o add new package rbcs (for relaxing ptracers at boundaries/depth/surface etc)
129    o add new offset feature to offline pkg (plus minor cleanup)
130    o cfc fixes for final year of integration
131    o clean up GCHEM_OPTIONS.h/DIC_OPTIONS.h (all options should
132      be in GCHEM_OPTIONS).
133    
134    checkpoint57x_post
135    o fix testreport Pb in multiple test experiment.
136    o cleanup our namelist syntax in verification:
137      - make it more F95 standards compliant (some work still needed)
138      - allows the vast majority of our tests to work with gfortran
139    o Non-Hydrostatic:
140     - put all NH variables (formerly in DYNVARS.h & GW.h) in NH_VARS.h
141     - fix exchange calls for CS-grid using the new EXCH_S3D_RL.
142    o Allow to apply AB on T,S rather than on AB(gT,gS):
143     - implemented within #ifdef ALLOW_ADAMSBASHFORTH_3
144     - use the same arrays (gtNm,gsNm) to hold tracer field at previous
145       time-steps (if AB(T,S)) and tendencies (if AB(gT,gS)).
146     - need a "clever" read_pickup to restart from the other AB formulation
147    o move local commom bloc /GlobalLo/ (used for singleCpuIO) in EESUPPORT.h
148    o remove some unused variables (reduces number of compiler warnings)
149    o adjoint
150      * Disable useVariableViscosity for adjoint for time being
151        (there are a number of extensive recomputations)
152      * Initialize certain variables needed for adjoint
153    
154    checkpoint57w_post
155    o pkg/mdsio: use MAX_LEN_FNAM (instead of hard coded 80) for file-names
156                 fix the singleCpuIO RS version of MDSWRITEFIELD
157       and add a WARNING if using globalFiles=T with multi-processors.
158    o write grid-file(s) in only one format (MDS or MNC): default=MNC (if
159      using MNC); can force MDS output by setting writegrid_mnc=F in data.mnc
160    o tracer advection scheme:
161      add DST-2 (=Lax-Wendroff) & 1rst order upwind schemes (the way to get
162       those 2 schemes before, changing GAD_FLUX_LIMITER.h + to recompile,
163       was not very convenient).
164      add implicit version for DST-3 & DST3_Flux_Limit vertical advection;  
165       (further tests needed for the implicit DST & Flux-Limiter schemes).
166    o PTRACERS: added namelist parameter PTRACERS_ref in analogy to tRef
167      and sRef for convenience, handy, when you want to initialize
168      ptracers in a large domain with a constant non-zero value. Also,
169      untangle ptracers_init.F a little.
170    
171    checkpoint57v_post
172    o change dome exp. to use variable viscosity in side-drag terms.
173    o remove all the files from the cfc_offline test since they are now
174      incorporated into tutorial_cfc_offline which has been producing
175      identical results for a few days of testing
176    o add new parameter: sideDragFactor to enable half-slip-side BC
177      and change side-drag S/R to use viscosity from mom_calc_visc
178      (can switch back by setting sideDragFactor=0.); affects results at
179      truncation level; regenerate output for exp.: global_ocean.cs32x10(x3),
180       global_ocean_ebm, global_ocean_pressure, matrix_example & MLAdjust(x5).
181    o vorticity & strain are now computed without masking, to account for
182      no-slip-side BC in variable viscosity calculation ; mask is applied
183      anyway later on, before computing advection/coriolis terms.
184    o modify exp4 to test new obcs support for reading from external files
185      and for passive tracers
186    o fix behavior of ptracers_monitor when monitor output is written to a
187      netcdf-file; still not independent of monitor.F (writes in the same
188      file).
189    o change default value of tRef (put a uniform, meaningful value for ocean
190      & atmos) and initialize sRef (was missing, as Martin mentioned).
191    o Redefine viscC?Leith*.  Change the coefficients of Leith viscosity
192      to be proportional to enstrophy dissipation length.  Thus,
193       viscC?leith?_new=pi*(viscC?leith?_old)^(1/3)
194      This is now consistent with Smagorinsky and with the manual.
195    
196    checkpoint57u_post
197    o OBCS and PTRACERS: add open boundary support for passive tracers
198      - either use homogenous (pseudo) v.Neumann conditions or prescribe
199        OB-values from file; this is not different from the way theta and salinity
200        are treated
201      - however, Orlanski-radiation conditions are not supported, and the model
202        will stop if you use pTracers and Orlanski at the same time.
203      - beefed up the rountine obcs_external_fields_load: now only those open
204        boundary values are overwritten with values from files for which there
205        are really files, otherwise the OB-fields remain untouched. This makes
206        it possible to use different OBs at different ends of the domain (as
207        with EXF)
208    o reduce number of exch calls in shap_filt (s2).
209    o CS-grid: fix divergence contribution to Leith viscosity. Affects (@
210      truncation level) 4/5 MLAdjust test results (11&12 digits match).
211    
212    checkpoint57t_post
213    o test biharmonic viscosity + no-slip with exp: global_ocean.cs32x15.viscA4
214    o finish adding brackets (coriolis S/R, vecinv) in multi-terms sum to get
215      same truncation, independent of the face number: affects affects several
216      test. exp.: Update global_ocean.cs32x15 (+.thsice) & MLAdjust (all 5).
217    o fix bug in mom_v_sidedrag (biharmonic viscosity was counted 2 times);
218      note that biharmonic viscosity + no-slip was not tested in any experiment.
219    o fix bug in non-linear free-surf side-drag ; generate new output for
220      global_ocean.cs32x15(+.thsice) & global_ocean_pressure.
221    o make mnc honor the writeBinaryPrec flag for most kinds of output
222    o eesupp: split EEPARAMS into EESIZE/EEPARAMS and add a vector global sum
223    o minor change in averaging W (mom_vi_?_vertshear.F) that affects truncation
224      error; update output from global_ocean.cs32x15(+.thsice) & MLAdjust (all 5).
225    o add brackets in multi-terms sum to get same truncation, independent of
226      the face number (manage to get exactly same value of gU,gV at the cs-edge)
227      affects several test. exp. Update global_ocean.cs32x15 (+.thsice).
228    o Clean, new version of mom_calc_visc.  Adds a new feature, viscAhRemax
229      and viscA4Remax, which allow one to limit the minimum viscosity by specifying
230      a maximum grid-scale Reynolds number.  Also, changes all length calculations
231      in the case of grid anisotropy to geometric mean of dx**2 and dy**2.
232      Suggested values of viscosity parameters found in mom_calc_visc.
233      Many new viscosity diagnostics available.
234    
235    checkpoint57s_post
236    o Add clearing-house for viscosity calculations, mom_calc_visc and
237      correct numerous bugs in LeithD and Smagorinsky discretizations.  
238      mom_calc_visc now allows for Smagorinsky, Leith and LeithD to be
239      used simultaneously, and to be used in mom_vecinv and mom_fluxform.
240      Mom_vecinv can be used with standard viscous calculation (mom_vi_hdissip.F)
241      and with Strain-Tension viscous term calculation (mom_hdissip.F).  
242      The Strain-Tension form is now called with useStrainTensonVisc=.true.
243      ViscAhStrain and viscAhTension are now retired.
244    o diagnostics "VISCAHD" & "VISCA4D" replace "VISCAH" & "VISCA4".
245    o In response to JMC's bug report, genmake2 now ignores any soft-links
246      that it finds in the *current* directory when searching for files.
247      Soft-links found in any other source directories will still work as
248      they did before.
249    o add HAVE_SETRLSTK define and useSETRLSTK flag (eedata) which calls
250      a C routine to unlimit the stack size
251    o mnc-ify pkg/aim snapshot vars
252    o mnc-ify pkg/land as requested by Daniel
253    o various changes to mnc including:
254      - all files use the new "BASENAME[[.ITER].{t|f}NUM].nc" format
255      - output can now be grouped so that all files within a group
256          change the ITER portion of their names in lock-step together
257      - can now read ("global") PER-FACE (in addition to PER-TILE) files
258          and works with both EXCH1 and EXCH2 (but needs more testing)
259      - writing works for all verification test cases w/ g77 on Linux
260    
261    checkpoint57r_post
262    o make MNC variable initialization happen earlier
263    o add diagnostics for (almost) each term of the momentum Eq.
264    o remove SST,SSS,UVEL_k2,VVEL_k2 diagnostics.
265    o move pkg timeave out of pkg-groups "oceanic" & "atmospheric" (pkg/pkg_groups)
266    o remove MINIMAL_TAVE_OUTPUT option ;
267      add "dumpInitAndLast" to replace "nodumps" (and MINIMAL_TAVE_OUTPUT).
268    o move some basic S/R (used elsewhere in the code) from pkg/mdsio to eesupp/src
269      allows to compile without pkg/mdsio (& test successfully hs94.1x64x5)
270    o add MNC output for adjoint of EXF (adexf), SEAICE (adseaice)
271    o add time-dependent SST, SSS control
272    o cost terms in ecco/ modularized (new cost_generic.F)
273    o cost_averages routines are modularized (new cost_averagesgeneric.F)
274    o new infrastructure for seaice model vs. obs misfits in seaice/
275      * driver from the_main_loop.F is seaice_cost_driver.F
276      * first term is AREA vs. SMR observation misfit (ALLOW_SEAICE_COST_SMR_AREA)
277      * gradient check in Exp. 6a of natl_box_adjoint/
278    o natl_box_adjoint input_seaice for Lab. Sea now available online via wget
279      (see README file)
280    
281    checkpoint57q_post
282    o fix initialisation Pb in calc_gw (get NANs in exp5 with g77)
283    o fix annoying months-old bug in mnc/mdsio output for seaice pkg
284      - bug was reported by Jean-Michel and Ian Fenty
285    o fix initialization Pb in AIM-Ocean coupled model (aim_sice2aim.F)
286      update exp. cpl_aim+ocn output.
287    o Implement 4-level checkpointing: AUTODIFF_4_LEVEL_CHECKPOINT
288      Needed for high-res. runs with large global files
289      to be able to limit tapelev size
290    o Changed some tests for adjoint:
291      * carbon/ now tests AUTODIFF_2_LEVEL_CHECKPOINT
292      * natl_box_adjoint/code_bulk/ tests AUTODIFF_4_LEVEL_CHECKPOINT
293      * global_ocean.90x40x15/ back to default 3-level checkpointing
294    o Modifications to enable scripted COMMON to MODULE conversion
295      * verification/OpenAD/code/ has necessary scripts and description;
296      * converts, compiles and runs successfully;
297    
298    checkpoint57p_post
299    o reading of delXfile & delYfile moved from ini_parms.F to ini_grid.F
300    o move ${PKG}_MNC_INIT from ${PKG}_READ_PARAMS to ${PKG}_INIT_FIXED
301      (already the case for some pkgs, including recent MNC_init, e.g. thsice)
302    o call CALC_GW from DYNAMICS (instead of from FORWARD_STEP)
303    o Dissipation & phiHyd gradient are always added to gU,gV in timestep.F
304      (was already the case for dissipation with mom_vecinv,
305       and also the case for grad.PhiHyd when staggered-timeStep)
306      This will allow to put dissipation out-off the AB time-stepping.
307     Unfortunately, affects truncation error. Update output of exp. that "fails"
308      (cg2d matching number) and leave the others for later: front_relax(10),
309      global_ocean_pressure(10),ideal_2D_oce(10),internal_wave(9),lab_sea(12),
310      matrix_example(7!),natl_box(12),plume_on_slope(11),vermix(12x3,11x2).
311    
312    checkpoint57o_post
313    o fix MNC/MDSIO flags in PTRACERS so that (1) it works, and
314      (2) its ready to be documented during the 20050801 DocFest
315    o MNC-ify the gmredi package -- tested and works w/ and w/o mnc
316      for ideal_2D_oce and lab_sea on linux_ia32_g77
317    o set to #undef CTRL_SET_OLD_MAXCVARS_30
318      in global1x1_tot/code_seaice_adj to access indices 32, 33
319      (worth another tag)
320    
321    checkpoint57n_post
322    o adding precip control xx_precip 32
323    o adding swflux control xx_swflux 33
324    o Replacing 10 forcing penalty cost routines by 1 generic.
325    o Implementing partitioning in time-mean and variable adjustment
326      for atmos. state.
327    o Time-mean constant error is now read in data.ecco (instead of hard-coded)
328    o in natl_box_adjoint:
329      * updated exp 1 to 10
330        (but does not have time-mean vs. variable split results yet)
331    o updated global1x1/*_seaice_adj/
332      to incorporate latest changes
333    o enable useSingleCpuIO for mdsio vector routines
334      (applies in particular to tapelev I/O)
335    o removed inefficient K-loop in all active_file*control*.F
336    
337    checkpoint57m_post
338    o re-arrange exp. adjustment.cs-32x32x1 (adjust_nlfs.cs-32x32x1 merge into)
339      and advect_cs (use EXCH-2 and standard cs-32 grid files).
340    o pkg/exch2: filling of face-corner halo regions is optional
341      (ifdef W2_FILL_NULL_REGIONS) and is turned off by default.
342    
343    checkpoint57l_post
344     tag before updating tutorials (pre-release-2)
345    o add cos & sin of grid-direction orientation angle (curvilinear grid only,
346       for output and interpolation of input fields onto model grid)
347    o fix radius in mom_v_metric_cylinder.F (wrong by dy/2).
348       Affects results of exp. rotating_tank ; update output.txt.
349    o pkg/diagnostics: 2 new S/R to fill a diagnostics using a scaling factor
350          and taking the square quantity and/or using a fraction-weight field.
351    o fix aim_diagnostics_init.F mate numbers.
352    o add mnc output capability to diagnostics/diagstat and update
353      our cvsignore files for the ACSII output generated
354    
355    checkpoint57k_post
356    o lab_sea now uses diagnostics package and provides a suite of
357      new EXF diagnostics output.
358    
359    checkpoint57j_post
360    o exf
361      * bug fix in exf_bulkformulae for psimh (spotted by Eliyahu Biton):
362        Bracket missing in eqn. for psimh (eqn. 7 of Large & Pond, 1982);
363        terms -2*ATAN(x) + pi/2 were added to both stable and unstable case.
364      * extended diagnostics output
365        - now includes all EXF variables
366        - hs, hl moved to common block
367        - separate S/R for diagnostics, called after ctrl vector is added
368        - name changed for EXFsw -> EXFswnet, EXFlw -> EXFlwnet
369        - EXFqnet contains true Qnet (i.e. includes qsw)
370      * updated verif. lab_sea, global_with_exf
371      * Rearranging hflux (exf_mapfields -> exf_getforcing)
372        to ensure that diagnostics sees a true EXFqnet
373    o add matlab scripts for CS-grid vectors (in utils/cs_grid)
374    o diagnostics pkg: change pointers so that 1 diagnostic can be used
375      several times with different frequencies.
376    o mnc-ify the thsice package as requested by Daniel Enderton
377      - the monitor--MNC output needs work
378      - many attributes need to be added (most are currently blank)
379    o mnc-ify theaim_v23 package as requested by Daniel Enderton
380      - many attributes need to be added (most are currently blank)
381    o remove "groundAtK1" & "rkFac" (replaced by usingPCoords & -rkSign)
382    o add Sadourny like discretization to the 4th.O vorticity advection scheme
383    o fix bug in rStar code (calc_r_star.F) when used with topo/bathy.
384      update output.txt from: aim.5l_cs (& .thSI), cpl_aim+ocn,
385         fizhi-gridalt-hs and global_ocean.cs32x15 (& .thsice).
386    o fix few little things to avoid floating error with -ftrapuv (ifort v9)
387     (ini_spherical_polar_grid.F, calc_r_star.F, monitor/mon_surfcor.F)
388    
389    checkpoint57i_post
390    o Diagnostics package was changed: diagnostic allocation is now deactivated
391      by setting an empty file name in the namelist (used to be by setting 0 frequency)
392    o momentum advection (VI form) options: activate 4th.O vorticity Coriolis ;
393      add 1rst Order upwind vertical Shear.
394    o add a quick fix for A/B checkpointing in OFFLINE mode
395    o add IGNORE_TIME define and -ignore_time option to genmake2 as an
396      end-run around the C/Fortran linking horror in eesup/src/timers.F
397    o added pkg/runclock: allows user to terminate run based on wall clock
398    o added non-hydrostatic scaling parameter: nh_am2
399    o pkg/my82: futile attempt to make results compiler-independent by adding
400      some "_d 0", etc. no effect.
401    o for mnc output, fill the 'T' coordinate var with myTime and create a
402      separate 'iter' variable for iteration count
403    o pkg/diagnostics:
404      add new capability: output of Global/Regional & per-level statistics ;
405      tested against monitor output (including land & thSIce) on various exp.
406    o make seaice-pkg I/O type flags (MNC/MDS) more standard.
407    
408    checkpoint57h_done
409      tag after changing user customized files (DIAGNOSTICS_SIZE.h,
410      data.diagnostics, data.gchem, data.ptracers) is done.
411    
412    o pkg/diagnostics: prepare for Global/Regional statistics.
413      add 2nd namelist in data.diagnostics and update DIAGNOSTICS_SIZE.h
414    o remove "baseTime" (no used) from arg. list of DIFF_BASE_MULTIPLE
415      and rename it: DIFFERENT_MULTIPLE
416    
417    checkpoint57h_post
418    o move tIter0 from data.gchem to PTRACERS_Iter0 in data.ptracers
419    o pkg/diagnostics: change output frequency, now in seconds (instead
420      of Nb of iterations). This requires to change data.diagnostics.
421    checkpoint57h_pre
422    
423    o in mnc, let the user create the output dir if both mnc_outdir_date
424      and mnc_outdir_num are false
425    o add fflags() to pkg/diagnostics which allows one to specify, on a
426      per-file basis, the file precision [default behavior is unchanged]
427    o added diagnostics to pkg/kpp, including computation of mixed layer
428      depth based on a temperature/density criterion
429    o updated verification/natl_box to test the new pkg/kpp diagnostics
430    o bugs fixed (with useAbsVorticity=T & with useJamartMomAdv=T)
431    o test AB-3 code with 4th.Order advection scheme (advect_xy/input.ab3_c4)
432    o fix mnc flags in pkg/seaice so that they are "owned" by seaice,
433      not mnc -- this mistake was found by JMC and is now corrected
434    
435    checkpoint57g_post
436    o add Adams-Bashforth 3rd Order time stepping:
437      - default remains AB-2 (#undef ALLOW_ADAMSBASHFORTH_3)
438      - not yet implemented with NH, CD_Scheme, pTracers, MNC_pickup ...
439    checkpoint57g_pre
440    
441    o make lambda[Theta,Salt]ClimRelax spatially varying
442      (default remains spatially uniform constant based on tau)
443    o added new control variables
444      * init. uVel, vVel, etanN
445      * lambda[Theta,Salt]ClimRelax
446    o separate masks used for ctrl_pack/unpack 'from write_grid' output
447      (suggested by G. Forget)
448    
449    checkpoint57f_post
450    o add parameter "baseTime" = model base time = time at iteration zero
451      new function DIFF_BASE_MULTIPLE (replace DIFFERENT_MULTIPLE) to
452      account for baseTime shift
453    o move Coriolis arrays from PARAMS.h to GRID.h
454    checkpoint57f_pre
455    
456    o mnc-ifying the seaice package
457    o introduced an alternative length scale calculation for use in the limiting
458      maximum viscosity (off by default, see mom_vi_hdissip.F for details)
459    o S/R ptracers_calc_diff: argument list bug for CALL KPP_CALC_DIFF_S
460    o added modified Leith scheme with viscosity proportional to grad(Hdiv)
461      as well as grad(vort3).  Made it compatible with both approximate form
462      for the cubed sphere and full form of these gradients.  Added flag
463      useFullLeith to select between approximate and full form.  
464    o add "d0" in getcon.F to reduce differences between compilers/platforms
465      update output.txt of exp. fizhi-gridalt-hs
466    o change port_rand funtion to allow to set the seed.
467    o fix obscure MNC bug involving conversion between reals and integers
468      - it caused no problems with certain (eg. gcc, ifc) compilers but the
469        new version is now more standards-compliant
470    o adjoint-related modifications
471      * gad_dst3_adv_?: to be able to use advscheme=30 in adjoint
472      * gad_calc_rhs: fool adjoint, revert to advscheme=30 in reverse when using 33
473      * mom_vi_hdissip: to be able to use Leith scheme in adjoint
474    
475    checkpoint57e_post
476    o Adding eddy stress controls a la Ferreira et al.
477    checkpoint57e_pre (actually was tagged: eckpoint57e_pre)
478    
479    o fizhi pkg: fix AtoC and CtoA when called with same input & output array
480      update output.txt of exp. fizhi-gridalt-hs
481    o add more "Z" dimensions and associated coordinate vars to MNC
482    o pkg/exf: added twoDigitYear capability for option useExfYearlyFields
483    o parameter calendarDumps: when set, approximate months (30-31 days) and years
484      (360-372 days) for parameters chkPtFreq, pChkPtFreq, taveFreq, SEAICE_taveFreq,
485      KPP_taveFreq, and freq in pkg/diagnostics are converted to exact calendar
486      months and years.  Requires pkg/cal.
487    o diagnostics: add pickup functionality
488      - off by default for backwards compatibility
489      - current version only uses MDSIO -- MNC will be added
490    o allow to use pkg diagnostics for LAND pkg output & Shapiro Tendencies.
491    o pkg/shap_filt: fix for multi-tiles per face (CS-grid with EXCH2)
492    
493    checkpoint57d_post
494    o pkg/thsice used in the ocean: fix tile Pb related to sIceLoad forcing
495    o modify verification/ideal_2D_oce so that it tests the creation of mnc
496      output from the diagnostics package
497    o temporary fix in EXCH2_UV: add a 3rd call to fix the 24.tiles CS-grid set-up
498    o make genmake2 more portable by using ".$FS" instead of ".f" for more
499      of the test compilations (eg. IBM xlf on Mac OSX w/HFS filesys)
500    o add coordinate variables to the MNC output from the diagnostics package
501    o change MNC-flags default value: true when useMNC=TRUE
502    o add a "readgrid_mnc" flag to MNC and give model/src/ini_curvilinear_grid.F
503      the ability (off by default) to read the grid information using MNC
504    o lab_sea: use LSR_ERROR=1.E-12 for testing ; update output.txt
505    o diagnostics: fix a bug (sub-set of levels); keep full precision when averaging
506    o do snap-shot diagnostics at the middle of the time-interval (instead of
507      at the end).
508    o update exp. aim.5l_cs (use default LAND_OPTIONS.h file) and add
509      new test for AIM with thSIce + slab-ocean
510    o cpl_aim+ocn test: update output (optfile & default GMREDI_OPTIONS.h)
511    o move state variable diagnostics to the beginning of the time step ;
512      (note: still needs some work specially for snap-shot diagnostics)
513      duplicate & replace UVELVVEL diagnostics with UV_VEL_C & UV_VEL_Z
514    
515    checkpoint57c_post
516    o Change testing option (-ieee) to get closer agreement between g77 & ifort
517    o generate new output.txt files on faulks with g77 & the new -ieee option
518    checkpoint57c_pre
519    
520    o add "_d 0" in freeze ; change some "FLOAT" to "DFLOAT" and went back (TAF);
521      update output.txt of exp. global_ocean.90x40x15
522    o fix in ini_eos.F (one "_d 0" was missing): now g77 & ifc give much closer
523      results; update global_ocean_pressure & global_ocean.cs32x15 output.txt
524    o add saltFlux(+File) forcing capability with periodicExternalForcing
525    o allow to use pkg diagnostics for GM-Redi & ThSIce output
526    o diagnostics_fill: fill in up to sNx+1 or sNy+1 (if U,V or Z pt) for mnc.
527    o improvements for the MNC output from DIAGNOSTICS:
528      - correct handling of the unlimitied dimension
529      - all dims now have associated coordinate variables
530      - variables located at U, V, and Corner ("Z") points are now written
531        with the correct horzontal dimensions (eg. 'Xp1', 'Yp1') so that
532        *all* data points are output
533    o added seaice_summary.F and removed obsolete ALLOW_SEAICE's from pkg/seaice
534    o added pkg/seaice/seaice_monitor.F
535    o mdsio_readfield.F can now read files visible only to master MPI process.
536    o pkg/seaice: some test code for ice-modified ocean stress
537    o exf_getffields interpolation defaults to bilinear for all scalar forcing
538      fields (remains bicubic for wind velocity and stress).  This avoids, e.g.,
539      spurious negative numbers for precipitation and humidity.  Will cause
540      some small numerical differences for integrations using pkg/exf/exf_interp.F.
541    o pkg/ctrl/
542      * Add ctrlvec diagnostics in pack/unpack for nondimensional I/O
543        May be enabled via doPackDiag
544      *  introduce z-scaling of
545         + gradient (doZscalePack) and
546         + control (doZscaleUnpack)
547    o pkg/mdsio/
548      * enable useSingleCpuIO for state estimation (ctrl_pack/unpack) mdsio_gl.F
549      * N.B.: leave tiled I/O for tapelev (i.e. mds...vector)
550    
551    checkpoint57b_post
552    o simplify options in diagnostics_fill arguments (remove weird combinations)
553    o fix MNC coordinate variables so that they now contain per-face index
554      values (instead of xC,xG,yC,yG values) when run with EXCH2 since, in
555      general, EXCH2 grids are only aligned to a "local" coordinate system
556    o test diagnostics pkg in dic_example (no output) and ideal_2D_oce (with ouput)
557    o Implement instantaneous output within pkg diagnostics:
558      setting a negative frequency in data.diagnostics produces
559     "snap-shot" output (as opposed to time-average if freq > 0)
560    o move out off EQUIVALENCE all diagnostics that fizhi does not use ;
561      NOTE: change ptracer diagnostics name: use 2 characters to work with
562            more than 9 tracers. Add also diagnostics of advect & diffus flux.
563    o more CF compliance: convert all MNC 'iter' variables to 'T' so that it
564      becomes the coordinate variable for the 'T' dimension
565    o add CF-style coordinate variables to MNC
566      - just a first cut:  numbers are meaningless on the cubesphere and
567        missing along the T axis but otherwise it works!
568    o add surface flux diagnostics: TAUX, TAUY, TFLUX, TRELAX, TICE, SFLUX, SRELAX
569    o pkg/thSIce: read initial state from 6 files ; allow monitor to use MNC.
570    o add diagnostics of advective & diffusive flux for T & S (pTracers coming soon)
571      this is based directly on what the model compute, so that the content of
572      each diagnostics might be affected by switching options (e.g., implicitDiffusion)
573    o small changes for diagnostics (gad_biharm_?.F, arguments of S/R impldiff.F)
574    o clean up build options files for SunOS, in particular
575      - remove -traditional from xmakedepend
576      - change -xarch=v9 to -xarch=native for f77
577      - restore use of make (not gmake)
578      - add CC=gcc if using g77
579    o changes per discussion of AM, JMC, & EH3 today:
580      - move redundant "grid" variables out of "state" -- in preparation of
581        creating proper CF-style coordinate variables and in recognition
582        that any/all grid variables can be easily copied (eg. w/NCO) from
583        the grid files into any of the other NetCDF files
584      - *always* output grid variables with MDSIO (this will eventually be
585        changed)
586      - output MNC grid file if useMNC is true
587    o every instance of _END_MASTER() has been found and replaced with
588      _END_MASTER( myThid ) in order to satisfy certain picky Sun
589      preprocessors
590    o re-arrange diagnostics pkg initialisation:
591     - allow each package to extend the available diagnostics list
592     - add some checking and fix small problems (multi-threaded, ...)
593     - new S/R diagnostics_fill (replace fill_diagnostics):
594       look through the short list of active diag. (instead of the long list)
595       create function DIAGNOSTICS_IS_ON to tell if a diagnostics is active
596    
597    checkpoint57a_post
598    o this set of changes restores TAMC compatibility
599      tested for global_ocean.90x40x15 adjoint on trough
600      (what a pain)
601    
602    checkpoint57a_pre
603    o some modifs in tools/build_options for sunos
604    o call ctrl_pack for ctrl at end for optimcycle=0
605    
606    checkpoint57
607    o add horizontal bi-harmonic viscosity for vertical velocity
608      - new parameter viscA4W (defaults to viscA4)
609      - affects only non-hydrostatic code
610      - to be done: no-slip lateral boundary conditions
611    o switch back to "sflux = sflux - runoff" in exf_getforcing.F
612    o implement Implicit Vertical advection for pTracers
613    o depth convergence accelerator: replace deltaTtracer by dTtracerLev(k)
614      (implemented also for implicit vertical diffusion & advection)
615    o OBCS: extend application of tracer open boundary conditions
616      to a strip beyond the boundary of width Olx/Oly
617    o add call monitor to TAF set of required calls (flow directives)
618      and include in natl_box_adjoint
619    o untangle a few ifdefs related to cost, ctrl, autodiff
620      untangel a few ALLOW_ECCO_OPTIMIZATION
621    o Apply Cooper & Haines representative estimates also at sub-surface (k=Nr)
622      (theta2, salt2 in ecco_cost_weights)
623    
624    checkpoint56c_post
625    o OBCS as control variables
626      - update ctrl_ad_diff.list, obcs_ad_diff.list
627      - remove balance of obcs controls from default
628      - fix index bug nobcs in ctrl_init
629      - fix dummy fields filen in ctrl_pack
630      - add dummy weights for obcs
631    o lsopt: comment out fort.94 output
632    o optim:
633      - remove unnecessary header files
634      - adjusted namelists
635      - add xerbla.F to Makefile
636    o exf:
637      - multiple modifications of exf_check_range
638      - unit change in exf_init_runoff from [m/year] to [m/s]
639    o mdsio:
640      - i/o-vector change debugLevel from debLevA to debLevB
641    
642    checkpoint56b_post
643    o collect all ptracers output in an new routine (PTRACERS_OUTPUT) that
644      is called from DO_THE_MODEL_IO
645      - rename PTRACERS_STATV_WRITE to PTRACERS_WRITE_TIMEAVE
646      - introduce new parameter PTRACERS_dumpFreq (defaults to dumpFreq)
647      - clean up PTRACERS_WRITE_TIMEAVE a little
648    o add another verification experiment that tests both the cfc-pkg and
649      the gchem-pkg with GCHEM_SEPARATE_FORCING undefined, that is with
650      reactive tendencies treated inline with regular timestepping
651    o testreport previously did a lousy job of comparing small numbers, so
652      we (Ed and Martin) changed the  way the cut-off of 1.e-12 is
653      implemented in the little c-progam tmp_cmpnum.c which is embedded in
654      testreport.
655    o ptracers_monitor: bug fix for MPI
656    
657    checkpoint56a_post
658    o GCHEM: finish reorganizating the package
659      - forward_step calls GCHEM_CALC_TENDENDY, which computes gchemTendency
660        (introduces another 3D-array for each passive tracer, but only if
661        GCHEM_SEPARATE_FORCING is undefined. For GCHEM_SEPARATE_FORCING
662        gchemTendency is not needed because the timestep is done separately)
663      - gchemTendency is added to gPtr in GCHEM_ADD_TENDENCY
664      - GCHEM is now---more or less cleanly---separate from PTRACERS
665      - to bo done: handling of surfaceForcingPtr that is needed for KPP not
666        so clear to me how that should be done in a general way.
667    o CFC: improve formatting of output in CFC_ATMOS
668    o PTRACERS:
669      - rename GCHEM_MONITOR to PTRACERS_MONITOR and call it from MONITOR, so
670        that ALL experiments with ptracers enable can be checked. This makes
671        GCHEM_MONITOR obsolete.
672      - include a runtime parameter PTRACERS_monitorFreq that defaults to
673        monitorFreq
674      - set default PTRACERS_write_mdsio to false if PTRACERS_write_mnc is true
675    o add PTRACERS (PTR0[1-5]) to testreport output
676  o exf_check_range:  o exf_check_range:
677    * default changed to .TRUE.    * default changed to .TRUE.
678    * will check values for niter0 if debLevA, every timestep if debLevB    * will check values for niter0 if debLevA, every timestep if debLevB
# Line 24  o GCHEM: Line 694  o GCHEM:
694      be replaced by a 3D array that is computed before thermodynamics; for      be replaced by a 3D array that is computed before thermodynamics; for
695      this, the tracer ID iTracer had to be included into the parameter list      this, the tracer ID iTracer had to be included into the parameter list
696      of ptracers_forcing      of ptracers_forcing
697    - set up forward_step for a call to gchem_forcing, that will replace    - set up forward_step to call gchem_forcing, that will replace
698      gchem_forcing_int in ptracers_forcing      gchem_forcing_int in ptracers_forcing
699    - cleaning up    - cleaning up
700  o small cleaning of exch2: remove files that are also in eesupp/src  o small cleaning of exch2: remove files that are also in eesupp/src

Legend:
Removed from v.1.392  
changed lines
  Added in v.1.614

  ViewVC Help
Powered by ViewVC 1.1.22