/[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.232 by jmc, Wed Dec 31 20:29:54 2003 UTC revision 1.456 by dimitri, Sun Feb 20 11:46:24 2005 UTC
# Line 2  $Header$ Line 2  $Header$
2    
3      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
4      ==============================      ==============================
5    o parameter calendarDumps: when set, approximate months (30-31 days) and years
6      (360-372 days) for parameters chkPtFreq, pChkPtFreq, taveFreq, SEAICE_taveFreq,
7      KPP_taveFreq, and freq in pkg/diagnostics are converted to exact calendar
8      months and years.  Requires pkg/cal.
9    o diagnostics: add pickup functionality
10      - off by default for backwards compatibility
11      - current version only uses MDSIO -- MNC will be added
12    o allow to use pkg diagnostics for LAND pkg output & Shapiro Tendencies.
13    o pkg/shap_filt: fix for multi-tiles per face (CS-grid with EXCH2)
14    
15    checkpoint57d_post
16    o pkg/thsice used in the ocean: fix tile Pb related to sIceLoad forcing
17    o modify verification/ideal_2D_oce so that it tests the creation of mnc
18      output from the diagnostics package
19    o temporary fix in EXCH2_UV: add a 3rd call to fix the 24.tiles CS-grid set-up
20    o make genmake2 more portable by using ".$FS" instead of ".f" for more
21      of the test compilations (eg. IBM xlf on Mac OSX w/HFS filesys)
22    o add coordinate variables to the MNC output from the diagnostics package
23    o change MNC-flags default value: true when useMNC=TRUE
24    o add a "readgrid_mnc" flag to MNC and give model/src/ini_curvilinear_grid.F
25      the ability (off by default) to read the grid information using MNC
26    o lab_sea: use LSR_ERROR=1.E-12 for testing ; update output.txt
27    o diagnostics: fix a bug (sub-set of levels); keep full precision when averaging
28    o do snap-shot diagnostics at the middle of the time-interval (instead of
29      at the end).
30    o update exp. aim.5l_cs (use default LAND_OPTIONS.h file) and add
31      new test for AIM with thSIce + slab-ocean
32    o cpl_aim+ocn test: update output (optfile & default GMREDI_OPTIONS.h)
33    o move state variable diagnostics to the beginning of the time step ;
34      (note: still needs some work specially for snap-shot diagnostics)
35      duplicate & replace UVELVVEL diagnostics with UV_VEL_C & UV_VEL_Z
36    
37    checkpoint57c_post
38    o Change testing option (-ieee) to get closer agreement between g77 & ifort
39    o generate new output.txt files on faulks with g77 & the new -ieee option
40    checkpoint57c_pre
41    
42    o add "_d 0" in freeze ; change some "FLOAT" to "DFLOAT" and went back (TAF);
43      update output.txt of exp. global_ocean.90x40x15
44    o fix in ini_eos.F (one "_d 0" was missing): now g77 & ifc give much closer
45      results; update global_ocean_pressure & global_ocean.cs32x15 output.txt
46    o add saltFlux(+File) forcing capability with periodicExternalForcing
47    o allow to use pkg diagnostics for GM-Redi & ThSIce output
48    o diagnostics_fill: fill in up to sNx+1 or sNy+1 (if U,V or Z pt) for mnc.
49    o improvements for the MNC output from DIAGNOSTICS:
50      - correct handling of the unlimitied dimension
51      - all dims now have associated coordinate variables
52      - variables located at U, V, and Corner ("Z") points are now written
53        with the correct horzontal dimensions (eg. 'Xp1', 'Yp1') so that
54        *all* data points are output
55    o added seaice_summary.F and removed obsolete ALLOW_SEAICE's from pkg/seaice
56    o added pkg/seaice/seaice_monitor.F
57    o mdsio_readfield.F can now read files visible only to master MPI process.
58    o pkg/seaice: some test code for ice-modified ocean stress
59    o exf_getffields interpolation defaults to bilinear for all scalar forcing
60      fields (remains bicubic for wind velocity and stress).  This avoids, e.g.,
61      spurious negative numbers for precipitation and humidity.  Will cause
62      some small numerical differences for integrations using pkg/exf/exf_interp.F.
63    o pkg/ctrl/
64      * Add ctrlvec diagnostics in pack/unpack for nondimensional I/O
65        May be enabled via doPackDiag
66      *  introduce z-scaling of
67         + gradient (doZscalePack) and
68         + control (doZscaleUnpack)
69    o pkg/mdsio/
70      * enable useSingleCpuIO for state estimation (ctrl_pack/unpack) mdsio_gl.F
71      * N.B.: leave tiled I/O for tapelev (i.e. mds...vector)
72    
73    checkpoint57b_post
74    o simplify options in diagnostics_fill arguments (remove weird combinations)
75    o fix MNC coordinate variables so that they now contain per-face index
76      values (instead of xC,xG,yC,yG values) when run with EXCH2 since, in
77      general, EXCH2 grids are only aligned to a "local" coordinate system
78    o test diagnostics pkg in dic_example (no output) and ideal_2D_oce (with ouput)
79    o Implement instantaneous output within pkg diagnostics:
80      setting a negative frequency in data.diagnostics produces
81     "snap-shot" output (as opposed to time-average if freq > 0)
82    o move out off EQUIVALENCE all diagnostics that fizhi does not use ;
83      NOTE: change ptracer diagnostics name: use 2 characters to work with
84            more than 9 tracers. Add also diagnostics of advect & diffus flux.
85    o more CF compliance: convert all MNC 'iter' variables to 'T' so that it
86      becomes the coordinate variable for the 'T' dimension
87    o add CF-style coordinate variables to MNC
88      - just a first cut:  numbers are meaningless on the cubesphere and
89        missing along the T axis but otherwise it works!
90    o add surface flux diagnostics: TAUX, TAUY, TFLUX, TRELAX, TICE, SFLUX, SRELAX
91    o pkg/thSIce: read initial state from 6 files ; allow monitor to use MNC.
92    o add diagnostics of advective & diffusive flux for T & S (pTracers coming soon)
93      this is based directly on what the model compute, so that the content of
94      each diagnostics might be affected by switching options (e.g., implicitDiffusion)
95    o small changes for diagnostics (gad_biharm_?.F, arguments of S/R impldiff.F)
96    o clean up build options files for SunOS, in particular
97      - remove -traditional from xmakedepend
98      - change -xarch=v9 to -xarch=native for f77
99      - restore use of make (not gmake)
100      - add CC=gcc if using g77
101    o changes per discussion of AM, JMC, & EH3 today:
102      - move redundant "grid" variables out of "state" -- in preparation of
103        creating proper CF-style coordinate variables and in recognition
104        that any/all grid variables can be easily copied (eg. w/NCO) from
105        the grid files into any of the other NetCDF files
106      - *always* output grid variables with MDSIO (this will eventually be
107        changed)
108      - output MNC grid file if useMNC is true
109    o every instance of _END_MASTER() has been found and replaced with
110      _END_MASTER( myThid ) in order to satisfy certain picky Sun
111      preprocessors
112    o re-arrange diagnostics pkg initialisation:
113     - allow each package to extend the available diagnostics list
114     - add some checking and fix small problems (multi-threaded, ...)
115     - new S/R diagnostics_fill (replace fill_diagnostics):
116       look through the short list of active diag. (instead of the long list)
117       create function DIAGNOSTICS_IS_ON to tell if a diagnostics is active
118    
119    checkpoint57a_post
120    o this set of changes restores TAMC compatibility
121      tested for global_ocean.90x40x15 adjoint on trough
122      (what a pain)
123    
124    checkpoint57a_pre
125    o some modifs in tools/build_options for sunos
126    o call ctrl_pack for ctrl at end for optimcycle=0
127    
128    checkpoint57
129    o add horizontal bi-harmonic viscosity for vertical velocity
130      - new parameter viscA4W (defaults to viscA4)
131      - affects only non-hydrostatic code
132      - to be done: no-slip lateral boundary conditions
133    o switch back to "sflux = sflux - runoff" in exf_getforcing.F
134    o implement Implicit Vertical advection for pTracers
135    o depth convergence accelerator: replace deltaTtracer by dTtracerLev(k)
136      (implemented also for implicit vertical diffusion & advection)
137    o OBCS: extend application of tracer open boundary conditions
138      to a strip beyond the boundary of width Olx/Oly
139    o add call monitor to TAF set of required calls (flow directives)
140      and include in natl_box_adjoint
141    o untangle a few ifdefs related to cost, ctrl, autodiff
142      untangel a few ALLOW_ECCO_OPTIMIZATION
143    o Apply Cooper & Haines representative estimates also at sub-surface (k=Nr)
144      (theta2, salt2 in ecco_cost_weights)
145    
146    checkpoint56c_post
147    o OBCS as control variables
148      - update ctrl_ad_diff.list, obcs_ad_diff.list
149      - remove balance of obcs controls from default
150      - fix index bug nobcs in ctrl_init
151      - fix dummy fields filen in ctrl_pack
152      - add dummy weights for obcs
153    o lsopt: comment out fort.94 output
154    o optim:
155      - remove unnecessary header files
156      - adjusted namelists
157      - add xerbla.F to Makefile
158    o exf:
159      - multiple modifications of exf_check_range
160      - unit change in exf_init_runoff from [m/year] to [m/s]
161    o mdsio:
162      - i/o-vector change debugLevel from debLevA to debLevB
163    
164    checkpoint56b_post
165    o collect all ptracers output in an new routine (PTRACERS_OUTPUT) that
166      is called from DO_THE_MODEL_IO
167      - rename PTRACERS_STATV_WRITE to PTRACERS_WRITE_TIMEAVE
168      - introduce new parameter PTRACERS_dumpFreq (defaults to dumpFreq)
169      - clean up PTRACERS_WRITE_TIMEAVE a little
170    o add another verification experiment that tests both the cfc-pkg and
171      the gchem-pkg with GCHEM_SEPARATE_FORCING undefined, that is with
172      reactive tendencies treated inline with regular timestepping
173    o testreport previously did a lousy job of comparing small numbers, so
174      we (Ed and Martin) changed the  way the cut-off of 1.e-12 is
175      implemented in the little c-progam tmp_cmpnum.c which is embedded in
176      testreport.
177    o ptracers_monitor: bug fix for MPI
178    
179    checkpoint56a_post
180    o GCHEM: finish reorganizating the package
181      - forward_step calls GCHEM_CALC_TENDENDY, which computes gchemTendency
182        (introduces another 3D-array for each passive tracer, but only if
183        GCHEM_SEPARATE_FORCING is undefined. For GCHEM_SEPARATE_FORCING
184        gchemTendency is not needed because the timestep is done separately)
185      - gchemTendency is added to gPtr in GCHEM_ADD_TENDENCY
186      - GCHEM is now---more or less cleanly---separate from PTRACERS
187      - to bo done: handling of surfaceForcingPtr that is needed for KPP not
188        so clear to me how that should be done in a general way.
189    o CFC: improve formatting of output in CFC_ATMOS
190    o PTRACERS:
191      - rename GCHEM_MONITOR to PTRACERS_MONITOR and call it from MONITOR, so
192        that ALL experiments with ptracers enable can be checked. This makes
193        GCHEM_MONITOR obsolete.
194      - include a runtime parameter PTRACERS_monitorFreq that defaults to
195        monitorFreq
196      - set default PTRACERS_write_mdsio to false if PTRACERS_write_mnc is true
197    o add PTRACERS (PTR0[1-5]) to testreport output
198    o exf_check_range:
199      * default changed to .TRUE.
200      * will check values for niter0 if debLevA, every timestep if debLevB
201      * STOP in exf_check_range if 'out of range'
202      * range chosen generous to allow for potential extema
203      TBD: a version with stricter range
204           permitted range would be runtime parameters
205    o GMredi:
206      - fix ldd97 slope limit ; extend valid domain of Psi-Bolus ;
207      - change S/R GRAD_SIGMA to be valid in the overlap on CS-grid.
208      - change the default GMREDI_OPTIONS.h: compile everything except Visbeck
209    o GCHEM:
210      - move call to GCHEM_FORCING_SEP, so that it is now called before
211        the blocking exchanges to avoid computing halos in GCHEM_FORCING_SEP
212      - introduced a useGCHEM runtime flag and applied it everywhere.
213      - renamed PTRACERS_SEPARATE_FORCING to GCHEM_SEPARATE_FORCING;
214        this flag is now only visible within GCHEM-pkg routines.
215      - moved gchem_forcing_int temporarily into ptracers_forcing, this will
216        be replaced by a 3D array that is computed before thermodynamics; for
217        this, the tracer ID iTracer had to be included into the parameter list
218        of ptracers_forcing
219      - set up forward_step to call gchem_forcing, that will replace
220        gchem_forcing_int in ptracers_forcing
221      - cleaning up
222    o small cleaning of exch2: remove files that are also in eesupp/src
223    
224    checkpoint56
225    o AIM: read land-Fraction earlier (from ini-fixed) and monitor land at nIter0
226    o mnc: cleanup "full" files--fixes out-of-memory problem
227    o ECCO-GODAE
228      * model/, pkg/, code/, input/ now complete w.r.t. ECCO-GODAE production setup
229      * fully merged and updated from SIO code
230        windstress sign change affects
231        - exf_inscal_[u/v]stress
232        - sign of xx_tau[u/v] adjustments
233        - signs in cost_scat
234      * ctrl_get_gen: avoid update of tau[u/v] for first few days
235      * natl_box_adjoint updated and compressed
236      * 2-level checkpoint switch added
237        (tested in EXP. 1 of natl_box_adjoint)
238      * new switch doInitXX to avoid init. of xx_... in optimcycle=0
239    o SVD
240      * first complete, non-hacked version for SVD calculations
241      * all interfaces clean w.r.t to handling of xx_..., adxx_..., g_xx_...
242        performed via direct access files admtlm_vector.it????
243        (via ctrl_[un]pack)
244      * optimcycle is used for ARPACK iteration counter
245        (note: adxx_... is iterated, rather than xx_...; latter is constant!)
246      * EV problem ok w.r.t. L2-norm in tangent space; need to check other norms
247        and add linear transforms where necessary
248      * controls are initial/final SST, SSS (can be extended; need to check
249        ctrl_set_[un]pack_xyz.F according to ctrl_set_[un]pack_xy.F
250    o mdsio: filename info to both stderr, stdout
251    o exf: exf_inscal_runoff: separated default conversion [m/year] -> [m/sec]
252      from sign change
253    
254    checkpoint55j_post
255    o isolate dissipation tendency (allow to keep it out off AB) in V.I form.
256      affects (truncation level) output.txt of global_ocean.cs32x15(&.thsice)
257    o add consistent build directories and .cvsignore files to verification
258    o added viscA4GridMin and viscA4GridMax for more control on biharmonic viscosity
259    o make all subroutines 30 or less characters in length for certain older
260      compilers to work (a thorough search only turns up two routines in
261      pkg/ptracers)
262    o allow (vecinv only) a different horiz. viscosity for Divergence and Vorticity;
263      affects (truncation level) output.txt of global_ocean.cs32x15(&.thsice)
264    
265    checkpoint55i_post
266    o change arrays for total vertical diffusivity (2.D if explicit,
267      3.D if implicit) ; each ptracer uses its own background value.
268    o added mnc/NetCDF output capability to ocn_compon_interf
269      - added bits to verification/cpl_aim+ocn so that it can run on
270        ITRDA and other clusters without the mpirun "-wd" option
271    o fixed wayward endif in mom_fluxform.F that messed up cyclindrical
272      configuration
273    o give MNC the ability to create additional files rather than
274      exceeding a specified file size limit "MNC_MAX_FILESIZE"
275      - the default is just less than 2GB which is a practical limit
276        on many systems
277    o Compute the non-local term using either KPP vertical diffusivity (=default)
278       or the total vertical diffusivity (when KPP_ghatUseTotalDiffus=T)
279    
280    checkpoint55h_post
281    o add verification test "vermix"
282    o remove all the *_ioinc flags and replace with the single global
283      outputTypesInclusive flag
284    o allow to set a vertical profile of vertical diffusivity for T & S
285    o call PACKAGE_BOOT & PACKAGES_READPARMS just after INI_PARMS
286    o useMNC moved from data.mnc to data.pkg
287    
288    checkpoint55g_post
289    o Renaming carbon AD experiments to comply with testreport conventions
290    o Updating adjoint of Held-Suarez hs94.1x64x5
291    o Fixing directives for non-ECCO forward_step, after ECCO directives
292      have been fixed to avoid recomp. of do_oceanic_phys
293      Works well for global_ocean.90x40x15,
294      but recomp. issue remains for carbon
295    o Updating adjoint output for natl_box_adjoint after
296      pkg/ecco has been modified (see taga c55d_post)
297    
298    checkpoint55f_post
299    o rename to timeave_init_fixed
300    
301    checkpoint55e_post
302    o some delicate re-shuffle of store directives to avoid one
303      extra call of do_oceanic_physics
304      affects ecco_the_main_loop.F, do_oceanic_physics.F
305      checkpoint_lev1_directives.h
306    o NB: this may break global_ocean adjoint temporarily,
307      but it is clear how to fix it. Will do later, need this now.
308    
309    checkpoint55d_post
310    o ECCO specific cost function terms (up-to-date with 1x1 runs)
311    o ecco_cost_weights is modified to 1x1 runs
312    o modifs to allow observations to be read in as
313      single file or yearly files
314    o enable to read exf forcing fields as either
315      single file or yearly files (flag useExfYearlyFields)
316    checkpoint55d_pre
317    
318    o switch between 2-level and 3-level checkpointing
319    o add ability of MNC to write local and "partial" (eg. 2D slices where
320      the full 3D field is never actually stored) arrays to NetCDF files
321      with the correct (that is, the complete multi-dimensional) set of
322      array indicies
323      - used in mom_vecinv() to write the diagFreq output
324      - tested (demonstrated) in verification/aim.5l_cs
325    o further mnc-ification
326      - add verification/lab_sea to the list with mnc setups
327      - fix undefined mnc v-names in cd_code
328      - add kpp to mnc output (both snapshot and timeave)
329    o non-hydrostatic code
330      - added horizontal viscosity parameter viscAhW for w, that replaces viscAh
331        in calc_gw.F
332      - defaults to viscAh
333      - this is useful, when I want to use viscAh=0 (e.g., only biharmonic
334        friction or Leith) in non-hydrostatic mode, as long as there is
335        no biharmonic friction and/or Leigth or equivalent for vertical velocity.
336    o pkg/cost
337      extend simple cost functions to treat time mean quantities
338      (T,U,V) via cost_accumulate_state from cost_tile
339    o multidimAdvection on CS-grid: implement the full 3 pass ; overlap needs to
340        be 2 times larger than the minimum size the advection scheme requires.
341        now get a perfect conservation of salt in advect_cs test-exp (updated)
342    o GGL90 (Gaspar's TKE mixed layer)
343      - add horizontal diffusion of TKE, requires exchanges in
344        do_fields_blocking_exchanges, horizontal diffusivity is zero by
345        default. In OPA there is no horizontal diffusion of TKE but the
346        mixing coefficients are computed from a horizontal average of TKE of
347        6 points or so. I think that diffusion has a little more physical
348        justification.
349      - clean up ggl90_calc in the hope of reducing memory usage (this hope
350        was in vain)
351      - mask tke-variable in ggl90_init
352    o fix mask Pb in multidimAdvection on CS-grid (still a small Pb with ipass)
353    o fix biharmonic diffusion on CS-grid.
354    o more obcs stuff
355      - added obcs_external_fields_load as an alternative to the EXF-pkg option
356        for reading boundary values. This routine is called from
357        obcs_prescribe_read
358      - updated obcs_check
359      - changed sponge layer relaxation parameter units for days to seconds
360      - turned off sponge layer by default
361      - changed argument list of OBCS_CALC to myTime+deltaTclock
362      - fixed about 35 bugs in obcs_external_fields_load
363    
364    checkpoint55c_post
365    o obcs stuff
366      - new runtime flags useOBCSprescribe, useOBCSbalance
367      - re-arranged obcs_calc
368      - moved exf_getobcs to obcs_prescibe_read
369        obcs_prescibe_read is called after ORLANSKI stuff
370        and independent of it (executable can contain all).
371      - updated obcs_check
372      - added exchanges for U,V,gU,gV for useOBCSprescribe=.TRUE.
373        These may be necessary only if tangential velocities are
374        prescribed at OB, not sure, needs testing.
375    o make the MNC_CW_ADD_VATTR_* subroutines all take scalar arguments
376      which simplifies the calls and passes the ifc check-bounds test
377    
378    checkpoint55b_post
379    o add a 2D scalar, Z location, exchange_2 routine.
380    o EXCH_UV for CS-grid: add one u,v in the corner-halo region, allows
381      to compute vorticity on a wider stencil [e.g., vort3(0,1)&(1,0)]
382      affects advect_cs results (for the wrong reason): update output.txt
383    o use a more standard discretization for biharmonic viscosity ;
384     (original version still available with #define MOM_VI_ORIGINAL_VISCA4 )
385      affects results of dome set-up: update output.txt
386    o fix gad_advection.F (CubedSphere part), back to version 1.23
387    o fix bug in mom_vi_del2uv.F (CubedSphere part).
388    
389    checkpoint55a_post
390    o merged code to
391      * prescribe/read time-dependent open boundaries
392        (works in conjunction with exf, cal)
393      * sponge layer code for open boundaries
394      * each boundary N/S/E/W now has its own CPP option
395        (healthy for the adjoint)
396    
397    checkpoint55
398    o up-to-date adjoint of previous breaks
399      - natl_box_adjoint: problem for exp 12, 13: cost = 0. needs fix
400      - global_ocean_ebm: adjoint still needs fix (TLM ok)
401      - global_ocean.90x40x15: extended GMREDI_OPTIONS.h
402      - carbon: * remove tr1-related adjoint verification
403                * add ptracer as control variable
404    o remove all tr1-related code (ALLOW_PASSIVE_TRACER)
405      (adjoint stuff still has some tr1 'names', but all use ptracer arrays)
406    
407    checkpoint54f_post
408    o in verification/carbin/ remove experiment related to
409      ALLOW_PASSIVE_TRACER (array tr1 etc)
410      code_ad/ input_ad/ results_ad
411    
412    o add two new packages
413      - ggl90, TKE-model of Gaspar et al. (1990), Blanke+Delecluse (1993)
414        + subject to optimization (memory)
415      - opps, OPPS (ocean penetrative plume scheme) of Paluszkiewicz+Romea (1997)
416        + subject to reordering (many GOTOs in the original code)
417        + so far opps is called for every (i,j) point, and then operates in
418          the vertical; one may want to reorder the scheme in order operate in
419          horizontal slabs at the cost of (a lot of) memory, at least ten new 3D
420          fields)
421    o overhaul of IO so that we now have flags for MDSIO and/or MNC
422      - all verification tests compile and run with linux_ia32_g77
423      - defaults are compatible with current input files--nothing
424        should change if you were not previously using MNC
425      - MNC output has been added in numerous places (eg. timeave)
426        but there are still a few writes not yet do-able with MNC
427        (this is in progress)
428      - flags now allow for either/or/both use of MDSIO and MNC and
429        documentation will soon follow
430      - numerous small formatting cleanups for ProTeX
431    o add version info (via BUILD_INFO.h) to MNC and stdout
432    
433    checkpoint54e_post
434    o add two new packages
435      - pp81 (Packanowski and Philander, 1981), Richardson number and
436        stratification dependent mixing
437      - my82 (Mellor and Yamada, 1982) level 2 turbulence closure scheme
438      - these scheme operate in a similar manner as KPP (without ghat, of
439        course), however there is only one diffusivity that is used for both
440        temperature and salinity (and passive tracers)ç
441    o add offline package
442    o add "verification/dic_example" which contains a dic-gchem-ptracer setup
443      provided by Stephanie Dutkiewicz
444    o genmake2: add ability to search for C compilers and cleanup test output
445      - fixes recent SX-6 problems and should help with others (eg. HP-UX)
446    
447    checkpoint54d_post
448    o data.diagnostics: if levels() is unspecified then ALL levels are diagnosed
449    o bug fix in pkg/cal/cal_addtime.F for TheCalendar='model' day count is wrong,
450      leading to additional shift by 5 days each year of integration
451       (i.e. leading to considerable offset over the years)
452    o update global_ocean.cs32x15/thsice: use z* + RealFW + staggerTimeStep
453    o start with no AB and no initial filtering (staggerTimeStep):
454      change the results of test-exp: hs94.128x64x5, hs94.1x64x5 & hs94.cs-32x32x5
455    o move ocean specific S/R calls from forward_step.F to do_oceanic_phys.F
456      (external loading & cpl import/export stay in forward_step.F)
457    o update coupling interface to work with new staggerTimeStep
458    o change MDSIO ptracers_pickup.
459    o fix a bug in Flux deriv. vs Ts (pkg/aim_v23/phy_suflux_sice.F);
460    o genmake2 fixes so netcdf now works on halem ("osf1_alpha_*")
461    o bugs fixed in kpp_calc.F (ustarY) and external_forcing.F (Qsw & Partial-cell)
462    
463    checkpoint54c_post
464    o replace surfaceTendency U,V,S,T,Tice,Ptr by surfaceForcing U,V,S,T,Tice,Ptr
465      this affects (truncation error) output.txt of test-exp:
466       global_ocean_ebm,_pressure,.cs32x15(&.thsice), lab_sea & natl_box
467    o pkg/kpp: replace pMask by maskC to reduce memory size (one 3.D array)
468    o change Qnet to always be the net heat flux, (+upward).
469      affects (truncation error) output of natl_box exp.: updated
470    o define the max number of passive tracers in PTRACERS_SIZE.h
471    o start to test mnc routinely in exp0
472    o add sea-ice loading (sIceLoad) as a new forcing field.
473    
474    checkpoint54b_post
475    o add testreport "-j JOBS" option for parallel builds (tested and pretty
476       impressive on Altix and Opteron systems)
477    o diagnostics can now write NetCDF files through the mnc package
478    o staggerTimeStep & NLFS: update pkg aim_v23 & thsice; update aim.5l_cs output
479    o staggerTimeStep & multi-Dim-Advec: Exch(U,V) added
480    o staggerTimeStep & NLFS: rescale T & S surface tendencies.
481      affects (truncation error) output of lab_sea & natl_box exp.: updated
482    o change pickup _cd & _ptracers (old pickup restart: usePickupBeforeC54=.T.)
483    
484    checkpoint54a_post
485    o re-write staggerTimeStep: step fwd U,V and then T,S (reverse the order).
486    o change pickup file (to restart from old pickup, set usePickupBeforeC54=.T.)
487    o get a perfect restart when using RealFreshWaterFlux + NLFS (not finished)
488    o update output of all exp. that use staggerTimeStep (except hs94.cs & 2D)
489      + internal_wave (shift in writing cg2d_ini_res)
490    checkpoint54a_pre
491    
492    o require forcing_In_AB=.FALSE. and remove surfaceTendencyTice from pkg/seaice
493      - modified verification/lab_sea and natl_box_adjoint accordingly
494    
495    checkpoint54
496    checkpoint53g_post
497    o prepare splitting of thermodynamics: store convect.counter in common block
498    o store d.etaH/dt (instead of Div.hV) in common ; affects (truncation error)
499      output of: global_ocean.cs32x15, global_ocean_pressure & internal_wave
500    o cnh's modifs to gad_advection.F for cube-sphere multi-dim advection
501    o new parameter "nShapS": allow to treat differently S from T (nShapT)
502    o separate Vert.Advec.Scheme from horizontal Advec.Scheme (T & S)
503    o restore adjointability
504    
505    checkpoint53f_post
506    o modifications that affect AIM experiment results (update output.txt):
507      - no advection of S at top interface Nr (pkg/generic_advdiff)
508      - air surf. temp: use lower value between T0 & T1 (phy_suflux_prep.F)
509    o AIM pkg:
510      - include stability function into surf.Flux derivative relative to Tsurf
511      - calculate clear-sky radiation & surface temp. change (for diagnostics)
512    
513    checkpoint53e_post
514    o merged cylindrical coord configuration
515    o added rotating_tank experiment
516    checkpoint53e_pre
517    
518    checkpoint53d_post
519    o faulks upgraded: re-generate output.txt that have changed:
520       global_ocean.cs32x15(.thsice), global_ocean_ebm, lab_sea, natl_box
521    o Added Bryan and Lewis, 1979, vertical diffusivity profile
522    o pkg/land: fix problem in land_impl_temp.F when Ts arround freezing.
523    o added Sadourny discretization of Coriolis
524      - moved some PARAMETERS from mom_*_coriolis.F to PARAMS.h
525      - re-enabled use of omega3 (useAbsVorticity)
526    o added Leith (Phys. Fluids, 1968, 10, 1409-1416) variable viscosity for
527      vector-invariant equations
528      - note: does *NOT* yet use same variable visc. for no-slip BC's.
529    o recoded the lateral boundary conditions for non-hydrostatic w, this time
530      hopefully without new bugs
531      - still to do: biharmonic friction for non-hydrostatic w
532    o bug fix for mom_*_sidedrag: rA -> rAs,rAw
533    o regenerated output for global_ocean_ebm (using -ieee)
534    o fix mom_common/mom_*_sidedrag.F to recover numeric output prior to creation
535      of mom_common/
536      - simple re-arrangement of instructions within computing loop
537    checkpoint53d_pre
538    
539    checkpoint53c_post
540    o new version of coupler for AIM+thSIce and Ocean models
541     - new fields (SSS,h_MixLayer ...) are exchanged through the coupler
542     - land: change Evaporation of snow
543     - aim with thSIce: fix bug in snow-precip.
544    o changes to calc_gw.F
545     - bug fix for w^2 term near the bottom boundary
546     - (hopefully) improve the lateral slip boundary condtions for use with
547        partial/looped cells
548      - because the bug fix changes two verification experiments anyway (exp5
549        and plume on slope), change the lateral boundary condition from half slip
550        to the value of no_slip_sides
551    o ebm package wasnt quite up to the parallel task. Fixed that.
552    o Adjoint of EBM package + verification
553      * TLM OK
554      * ADM not sufficient accuracy yet (needs further testing)
555    
556    checkpoint53b_post
557    o Commiting new energy balance model to repository
558      * package is pkg/ebm
559      * verif. is verification/global_ocean_ebm
560      * references are in ebm_driver.F
561    C      * X. Wang, P. Stone and J. Marotzke, 1999:
562    C        Global thermohaline circulation. Part I:
563    C        Sensitivity to atmospheric moisture transport.
564    C        J. Climate 12(1), 71-82
565    C      * X. Wang, P. Stone and J. Marotzke, 1999:
566    C        Global thermohaline circulation. Part II:
567    C        Sensitivity with interactive transport.
568    C        J. Climate 12(1), 83-91
569    C      * M. Nakamura, P. Stone and J. Marotzke, 1994:
570    C        Destabilization of the thermohaline circulation
571    C        by atmospheric eddy transports.
572    C        J. Climate 7(12), 1870-1882
573    o Will need long integration testing.
574    checkpoint53b_pre
575    
576    o created pkg/mom_common in preparation for variable viscosity code
577    o renamed and split ini_mnc_io
578    o pkg/land: only liquid water can diffuse or run-off
579    o AIM using land: soil wetness availability account for snow
580    o added printMapIncludesZeros to EEPARMS namelist
581    
582    checkpoint53a_post
583    o restore TAMC compatiility for exf_constants.h
584      (no usage of _d in .h)
585    o Additional storing in KPP in an attempt to fix
586      TAMC problem with adjoint of KPP code.
587      (no success...)
588      Identical code works fine with TAF.
589    o Fixed verification/natl_box_adjoint experiments.
590      Initial conditions for Experiments 4 and 5 are now
591      set to "SEAICE_initialHEFF=0", which is more reasonable,
592      and hence gives results closer to that of Experiment 2.
593      Experiment 6 requires "SEAICE_initialHEFF=1" and is working once again.
594    
595    checkpoint53
596    o stable checkpoint with up to date adjoint
597      TLM tests only run for 1st grdchk at the moment. Its unclear if
598      this is a code problem or a problem with resetting code for g_
599      variables that is being autogenerated wrong.
600      Otherwise:
601       o code is testing well on Linux (P4 and AMD), SGI Alitix and Irix.
602       o IBM and HPQ tests could not be run because NCAR (IBM) and NCCS (HPQ)
603         are both still in chaos following their security messes.
604    
605    checkpoint52n_post
606    o various build (mostly genmake2) changes including:
607      - ability to build/run more easily within Cygwin and MacOS X
608      - more complicated Makefiles with a stronger dependence on Gnu Make
609    o pkg/seaice modifications
610      - added SEAICE_deltaTtherm and SEAICE_deltaTdyn in order to reduce cost
611        of pkg/seaice.  For time being, SEAICE_deltaTtherm = deltaTtracer
612        and SEAICE_deltaTdyn must be integer multiple of SEAICE_deltaTtherm
613      - changed proxy for geostrophic velocity and updated
614        verification/lab_sea/results/* accordingly
615      - removed ADI dynamic solver from pkg/seaice
616      - replaced all occurrences of 1/CSTICE or 1/CSUICE
617        they caused a segmentation fault under certain circumstances
618      - removed unused field WATR
619    o bug fix in exf-package:
620      -  u/vwind are defined at tracer/pressure points. The values of u/vwindmask
621         now reflect that.
622    o small fixes in the cal-package:
623      - set reference date to first of January for the 'model'-calendar
624      - failing of cal_CheckDate in cal_FullDate now actually results in a
625        warning
626    o small modifications in exf-package to increase flexibility:
627      - change climtempfreeze from a constant to a runtime parameter that
628        is specified in data.exf_clim (default = -1.9)
629      - introduce runtime parameter windstressmax (default = 2.0) that is
630        specified in data.exf
631    o fix in pkg/thsice:
632      - update Evap (Tsf changes) to be consistent with Latent heat flux
633      - maintain snow precipitation unchanged (even if Tsf >=0).
634     both affect global_ocean.cs32x15/(.thsice) results ; update output.txt.thsice
635    
636    checkpoint52m_post
637    o commented out call to cal_GetMonthsRec in swfrac.F and added some comments;
638      this does not change any of the results, because myIter=0 in swfrac always
639      lead to jwtype=2 (the default anyway), but the call makes little sense
640      as a default
641    o no SHORTWAVE HEATING lost at the bottom; affects lab_sea output; updated
642    o re-work pkg/thsice interface with ocean & bulk_force pkg to allow
643      atmospheric model (AIM) to use it ; add a "slab ocean" component.
644    o pkg/aim_v23 : add interface to call thermodynamic sea-ice (thsice)
645    o update global_ocean.cs32x15/output (.thsice) after changes in pkg/thsice
646    o extending useSingleCpuIO option to work with new exch2 I/O format
647      - old-style, missing-tile I/O is still accessible by defining CPP
648        option MISSING_TILE_IO in pkg/mdsio/MDSIO_OPTIONS.h
649    o monitor package can write to NetCDF files and/or STDOUT
650    o fix re-definition bug in MNC_GRID_INIT_ALL()
651    o doc/api_reference added
652      - a framework for building an API Reference Manual using all
653        the "protex" comments embedded in the code
654      - see pkg/generic_advdiff, pkg/mnc for examples
655    o MNC:
656      - make myThid the last argument for all subroutines
657      - new default for the grid files: NO halos
658      - add useMNC and other run-time flags
659      - separate per-file-group handling of the unlimited dimension
660      - early initialization of MNC w/ flags read at the eeboot stage
661        and lookup tables set in initialize_fixed
662      - write all MNC output to a freshly-created directory
663    o thermodynamics sea-ice (thsice) is now tested with testreport.
664    o test exp. advect_cs: use the standard cs32 grid ; update output.
665    o Added functionality to grdchk:
666      pick global i,j,k position (or nearest wet) where to perform check.
667    
668    checkpoint52l_post
669    o new land formulation:
670      a) use ground enthalpy as prognostic variable to ensure exact
671         energy conservation.
672      b) account for water temperature and for latent heat of freezing
673         in all processes (rain, run-off, ground storage)
674      c) compute surface and ground temperature implicitly.
675    o aim_v23 modification for new land model:
676      a) Treat separately land / ocean / sea-ice surface fluxes
677         to allow implicit computation of land & sea-ice surface temp.
678      b) add snow precipitation.
679    checkpoint52l_pre
680    
681    o fix problem (when 1 tile is empty) in min/max monitor.
682    o fix bug in timeAve (u,vZetatave in common block)
683    o no restoring terms under ice for pkg/seaice
684    o moved SEAICE_DO_DIAGS to S/R do_the_model_io
685    o bug fix for pkg/seaice pickups
686    o change MNC call locations and add coordinate attributes
687    o mask vorticity at the corner ; affects global_ocean.cs output
688    o fix ini_curvilinear_grid ; test CS experiments with mpi
689    o aim.5l_cs: modify sea-ice fraction input file ; update output
690    o exf: Adding simple range check routine for exf fields after read
691    
692    checkpoint52k_post
693    o removed top-level directory "diags"
694      - we migrated these functions to pkg/timeave a long time ago
695    o fixed bug introduced in pkg/monitor
696    o accuracy ctrlprec = 32 insuffient for gradient checks using
697      averaged fields (I/O via cost_averages)
698      -> use ctrl.h in active_file*.F to control I/O precision.
699    o paramter list of CALL ADTHE_MAIN_LOOP
700    o added initial do_the_model_io to ecco
701    o set all ...Freq to zero after cost_final
702      taveFreq, dumpFreq, pChkptFreq
703    o MAIN and ECCO differentiability restored
704    
705    checkpoint52j_post
706    o fix nchklev in tamc.h for global_ocean...
707    o Bringing flow and store directives up-to-date
708    o Really weird thing!
709      call do_the_model_io from the_main_loop with parameters
710      starttime, niter0, instead of mytime, myiter, screws up the adjoint:
711      TAF places a faulty re-call of initialise_varia
712      after cost_final call which should not be there!
713    checkpoint52j_pre
714    
715    o tweaked pkg/mdsio/mdsio_readfield.F to have more useful logic and
716      error messages
717    o fix errorMessageUnit problem. remove old version of mdsio from eesup/src
718    o modified pkg/mdsio/mdsio_readfield.F -- commented out diskspace-eating
719      debug output
720    o modfied mdsio_writefield.F to handle non-global files under exch2
721    o update AIM Cubed-Sphere exp (aim.5l_cs):
722      use land, p*, and Franco.M forcing fields with fractional land-sea mask.
723    o update AIM Equatorial Channel exp:
724      use standard aim_v23 pkg and change the SST fields to be more realistic.
725    o globalFile now works with exch2
726    o merged mom_vi_calc_relvort3.F, mom_vi_del2uv.F and
727      mon_vort3.F to make exch2 friendly with >1 tile per cube face
728    o limit timeave output for hi-res integrations
729    
730    checkpoint52i_post
731    o removed CPP_EEOPTIONS.h from verification/*/code
732    o updated the tiling for global_ocean_pressure
733    checkpoint52i_pre
734    
735    o update aim.5l_LatLon using the standard aim_v23 pkg (instead of
736      the old pkg/aim).
737    
738    checkpoint52h_pre
739    o shortened offending lines in exch2_send_rx?.template
740    o modified W2_EXCH2_TOPOLOGY.h to reflect current topology generator
741    o genmake2: add build & execute test for NetCDF code
742    o modded aim.5l_cs and global_ocean.cs32x15 to use exch2, removed
743      old grid files, moved CPP_EEMACROS.h to pkg/exch2
744    o added useSingleCPUIO flag to monitor.F and mon_out.F
745    o added temporary hack for sea-ice dynamics on the cube:
746      metric terms are ignored and it is assumed that
747      there is no sea-ice on cube faces 1, 2, 4, and 5.
748    
749    checkpoint52f_post
750    o added exch2 cube routines to pkg/exch2
751    o modified hs94.cs-32x32x5 to use it
752    o added new grid files to hs94.cs-32x32x5
753    checkpoint52f_pre
754    
755    checkpoint52e_post
756    o first attempt to solve vertical advection (momemtum & T,S) Implicitly
757      for now, only implemented with Lax-Wendroff (tracer) and VectInv (mom)
758    o step.1 : rewrite (as in MultiDimAdvec) explicit tracer stepping (gad_calc_rhs.F)
759      to work with implicit vertical advection and AB ;
760    o Implic. vert. advect. implemented with flux-limiter, c2, u3 and c4 scheme.
761    o update output (gad_calc_rhs.F changes affect truncation error) of exp:
762     front_relax, global_ocean .cs32x15 & _pressure, internal_wave, lab_sea, natl_box
763    checkpoint52e_pre
764    
765    o fix bug in impldiff when Nr=1 or Nr=2.
766  o change the albedo in pkg/thsice (as GISS model); update output_thsice.txt  o change the albedo in pkg/thsice (as GISS model); update output_thsice.txt
767  o little changes (fftpack, aim/phy_suflux ...) to pass the "checkbounds" test.  o little changes (fftpack, aim/phy_suflux ...) to pass the "checkbounds" test.
768  o fix deltaTfreesurf in update_cg2d.F (left from c45d_post): affects 2  o fix deltaTfreesurf in update_cg2d.F (left from c45d_post): affects 2

Legend:
Removed from v.1.232  
changed lines
  Added in v.1.456

  ViewVC Help
Powered by ViewVC 1.1.22