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

Legend:
Removed from v.1.257  
changed lines
  Added in v.1.427

  ViewVC Help
Powered by ViewVC 1.1.22