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

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

  ViewVC Help
Powered by ViewVC 1.1.22