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

Legend:
Removed from v.1.390  
changed lines
  Added in v.1.605

  ViewVC Help
Powered by ViewVC 1.1.22