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

Legend:
Removed from v.1.223  
changed lines
  Added in v.1.428

  ViewVC Help
Powered by ViewVC 1.1.22