/[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.178 by dimitri, Wed Sep 24 06:50:26 2003 UTC revision 1.370 by edhill, Mon Oct 18 16:04:19 2004 UTC
# Line 3  $Header$ Line 3  $Header$
3      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
4      ==============================      ==============================
5    
6  untagged  o remove all the *_ioinc flags and replace with the single global
7      outputTypesInclusive flag
8    o allow to set a vertical profile of vertical diffusivity for T & S
9    
10    checkpoint55g_post
11    o Renaming carbon AD experiments to comply with testreport conventions
12    o Updating adjoint of Held-Suarez hs94.1x64x5
13    o Fixing directives for non-ECCO forward_step, after ECCO directives
14      have been fixed to avoid recomp. of do_oceanic_phys
15      Works well for global_ocean.90x40x15,
16      but recomp. issue remains for carbon
17    o Updating adjoint output for natl_box_adjoint after
18      pkg/ecco has been modified (see taga c55d_post)
19    
20    checkpoint55f_post
21    o rename to timeave_init_fixed
22    
23    checkpoint55e_post
24    o some delicate re-shuffle of store directives to avoid one
25      extra call of do_oceanic_physics
26      affects ecco_the_main_loop.F, do_oceanic_physics.F
27      checkpoint_lev1_directives.h
28    o NB: this may break global_ocean adjoint temporarily,
29      but it is clear how to fix it. Will do later, need this now.
30    
31    checkpoint55d_post
32    o ECCO specific cost function terms (up-to-date with 1x1 runs)
33    o ecco_cost_weights is modified to 1x1 runs
34    o modifs to allow observations to be read in as
35      single file or yearly files
36    o enable to read exf forcing fields as either
37      single file or yearly files (flag useExfYearlyFields)
38    checkpoint55d_pre
39    
40    o switch between 2-level and 3-level checkpointing
41    o add ability of MNC to write local and "partial" (eg. 2D slices where
42      the full 3D field is never actually stored) arrays to NetCDF files
43      with the correct (that is, the complete multi-dimensional) set of
44      array indicies
45      - used in mom_vecinv() to write the diagFreq output
46      - tested (demonstrated) in verification/aim.5l_cs
47    o further mnc-ification
48      - add verification/lab_sea to the list with mnc setups
49      - fix undefined mnc v-names in cd_code
50      - add kpp to mnc output (both snapshot and timeave)
51    o non-hydrostatic code
52      - added horizontal viscosity parameter viscAhW for w, that replaces viscAh
53        in calc_gw.F
54      - defaults to viscAh
55      - this is useful, when I want to use viscAh=0 (e.g., only biharmonic
56        friction or Leith) in non-hydrostatic mode, as long as there is
57        no biharmonic friction and/or Leigth or equivalent for vertical velocity.
58    o pkg/cost
59      extend simple cost functions to treat time mean quantities
60      (T,U,V) via cost_accumulate_state from cost_tile
61    o multidimAdvection on CS-grid: implement the full 3 pass ; overlap needs to
62        be 2 times larger than the minimum size the advection scheme requires.
63        now get a perfect conservation of salt in advect_cs test-exp (updated)
64    o GGL90 (Gaspar's TKE mixed layer)
65      - add horizontal diffusion of TKE, requires exchanges in
66        do_fields_blocking_exchanges, horizontal diffusivity is zero by
67        default. In OPA there is no horizontal diffusion of TKE but the
68        mixing coefficients are computed from a horizontal average of TKE of
69        6 points or so. I think that diffusion has a little more physical
70        justification.
71      - clean up ggl90_calc in the hope of reducing memory usage (this hope
72        was in vain)
73      - mask tke-variable in ggl90_init
74    o fix mask Pb in multidimAdvection on CS-grid (still a small Pb with ipass)
75    o fix biharmonic diffusion on CS-grid.
76    o more obcs stuff
77      - added obcs_external_fields_load as an alternative to the EXF-pkg option
78        for reading boundary values. This routine is called from
79        obcs_prescribe_read
80      - updated obcs_check
81      - changed sponge layer relaxation parameter units for days to seconds
82      - turned off sponge layer by default
83      - changed argument list of OBCS_CALC to myTime+deltaTclock
84      - fixed about 35 bugs in obcs_external_fields_load
85    
86    checkpoint55c_post
87    o obcs stuff
88      - new runtime flags useOBCSprescribe, useOBCSbalance
89      - re-arranged obcs_calc
90      - moved exf_getobcs to obcs_prescibe_read
91        obcs_prescibe_read is called after ORLANSKI stuff
92        and independent of it (executable can contain all).
93      - updated obcs_check
94      - added exchanges for U,V,gU,gV for useOBCSprescribe=.TRUE.
95        These may be necessary only if tangential velocities are
96        prescribed at OB, not sure, needs testing.
97    o make the MNC_CW_ADD_VATTR_* subroutines all take scalar arguments
98      which simplifies the calls and passes the ifc check-bounds test
99    
100    checkpoint55b_post
101    o add a 2D scalar, Z location, exchange_2 routine.
102    o EXCH_UV for CS-grid: add one u,v in the corner-halo region, allows
103      to compute vorticity on a wider stencil [e.g., vort3(0,1)&(1,0)]
104      affects advect_cs results (for the wrong reason): update output.txt
105    o use a more standard discretization for biharmonic viscosity ;
106     (original version still available with #define MOM_VI_ORIGINAL_VISCA4 )
107      affects results of dome set-up: update output.txt
108    o fix gad_advection.F (CubedSphere part), back to version 1.23
109    o fix bug in mom_vi_del2uv.F (CubedSphere part).
110    
111    checkpoint55a_post
112    o merged code to
113      * prescribe/read time-dependent open boundaries
114        (works in conjunction with exf, cal)
115      * sponge layer code for open boundaries
116      * each boundary N/S/E/W now has its own CPP option
117        (healthy for the adjoint)
118    
119    checkpoint55
120    o up-to-date adjoint of previous breaks
121      - natl_box_adjoint: problem for exp 12, 13: cost = 0. needs fix
122      - global_ocean_ebm: adjoint still needs fix (TLM ok)
123      - global_ocean.90x40x15: extended GMREDI_OPTIONS.h
124      - carbon: * remove tr1-related adjoint verification
125                * add ptracer as control variable
126    o remove all tr1-related code (ALLOW_PASSIVE_TRACER)
127      (adjoint stuff still has some tr1 'names', but all use ptracer arrays)
128    
129    checkpoint54f_post
130    o in verification/carbin/ remove experiment related to
131      ALLOW_PASSIVE_TRACER (array tr1 etc)
132      code_ad/ input_ad/ results_ad
133    
134    o add two new packages
135      - ggl90, TKE-model of Gaspar et al. (1990), Blanke+Delecluse (1993)
136        + subject to optimization (memory)
137      - opps, OPPS (ocean penetrative plume scheme) of Paluszkiewicz+Romea (1997)
138        + subject to reordering (many GOTOs in the original code)
139        + so far opps is called for every (i,j) point, and then operates in
140          the vertical; one may want to reorder the scheme in order operate in
141          horizontal slabs at the cost of (a lot of) memory, at least ten new 3D
142          fields)
143    o overhaul of IO so that we now have flags for MDSIO and/or MNC
144      - all verification tests compile and run with linux_ia32_g77
145      - defaults are compatible with current input files--nothing
146        should change if you were not previously using MNC
147      - MNC output has been added in numerous places (eg. timeave)
148        but there are still a few writes not yet do-able with MNC
149        (this is in progress)
150      - flags now allow for either/or/both use of MDSIO and MNC and
151        documentation will soon follow
152      - numerous small formatting cleanups for ProTeX
153    o add version info (via BUILD_INFO.h) to MNC and stdout
154    
155    checkpoint54e_post
156    o add two new packages
157      - pp81 (Packanowski and Philander, 1981), Richardson number and
158        stratification dependent mixing
159      - my82 (Mellor and Yamada, 1982) level 2 turbulence closure scheme
160      - these scheme operate in a similar manner as KPP (without ghat, of
161        course), however there is only one diffusivity that is used for both
162        temperature and salinity (and passive tracers)ç
163    o add offline package
164    o add "verification/dic_example" which contains a dic-gchem-ptracer setup
165      provided by Stephanie Dutkiewicz
166    o genmake2: add ability to search for C compilers and cleanup test output
167      - fixes recent SX-6 problems and should help with others (eg. HP-UX)
168    
169    checkpoint54d_post
170    o data.diagnostics: if levels() is unspecified then ALL levels are diagnosed
171    o bug fix in pkg/cal/cal_addtime.F for TheCalendar='model' day count is wrong,
172      leading to additional shift by 5 days each year of integration
173       (i.e. leading to considerable offset over the years)
174    o update global_ocean.cs32x15/thsice: use z* + RealFW + staggerTimeStep
175    o start with no AB and no initial filtering (staggerTimeStep):
176      change the results of test-exp: hs94.128x64x5, hs94.1x64x5 & hs94.cs-32x32x5
177    o move ocean specific S/R calls from forward_step.F to do_oceanic_phys.F
178      (external loading & cpl import/export stay in forward_step.F)
179    o update coupling interface to work with new staggerTimeStep
180    o change MDSIO ptracers_pickup.
181    o fix a bug in Flux deriv. vs Ts (pkg/aim_v23/phy_suflux_sice.F);
182    o genmake2 fixes so netcdf now works on halem ("osf1_alpha_*")
183    o bugs fixed in kpp_calc.F (ustarY) and external_forcing.F (Qsw & Partial-cell)
184    
185    checkpoint54c_post
186    o replace surfaceTendency U,V,S,T,Tice,Ptr by surfaceForcing U,V,S,T,Tice,Ptr
187      this affects (truncation error) output.txt of test-exp:
188       global_ocean_ebm,_pressure,.cs32x15(&.thsice), lab_sea & natl_box
189    o pkg/kpp: replace pMask by maskC to reduce memory size (one 3.D array)
190    o change Qnet to always be the net heat flux, (+upward).
191      affects (truncation error) output of natl_box exp.: updated
192    o define the max number of passive tracers in PTRACERS_SIZE.h
193    o start to test mnc routinely in exp0
194    o add sea-ice loading (sIceLoad) as a new forcing field.
195    
196    checkpoint54b_post
197    o add testreport "-j JOBS" option for parallel builds (tested and pretty
198       impressive on Altix and Opteron systems)
199    o diagnostics can now write NetCDF files through the mnc package
200    o staggerTimeStep & NLFS: update pkg aim_v23 & thsice; update aim.5l_cs output
201    o staggerTimeStep & multi-Dim-Advec: Exch(U,V) added
202    o staggerTimeStep & NLFS: rescale T & S surface tendencies.
203      affects (truncation error) output of lab_sea & natl_box exp.: updated
204    o change pickup _cd & _ptracers (old pickup restart: usePickupBeforeC54=.T.)
205    
206    checkpoint54a_post
207    o re-write staggerTimeStep: step fwd U,V and then T,S (reverse the order).
208    o change pickup file (to restart from old pickup, set usePickupBeforeC54=.T.)
209    o get a perfect restart when using RealFreshWaterFlux + NLFS (not finished)
210    o update output of all exp. that use staggerTimeStep (except hs94.cs & 2D)
211      + internal_wave (shift in writing cg2d_ini_res)
212    checkpoint54a_pre
213    
214    o require forcing_In_AB=.FALSE. and remove surfaceTendencyTice from pkg/seaice
215      - modified verification/lab_sea and natl_box_adjoint accordingly
216    
217    checkpoint54
218    checkpoint53g_post
219    o prepare splitting of thermodynamics: store convect.counter in common block
220    o store d.etaH/dt (instead of Div.hV) in common ; affects (truncation error)
221      output of: global_ocean.cs32x15, global_ocean_pressure & internal_wave
222    o cnh's modifs to gad_advection.F for cube-sphere multi-dim advection
223    o new parameter "nShapS": allow to treat differently S from T (nShapT)
224    o separate Vert.Advec.Scheme from horizontal Advec.Scheme (T & S)
225    o restore adjointability
226    
227    checkpoint53f_post
228    o modifications that affect AIM experiment results (update output.txt):
229      - no advection of S at top interface Nr (pkg/generic_advdiff)
230      - air surf. temp: use lower value between T0 & T1 (phy_suflux_prep.F)
231    o AIM pkg:
232      - include stability function into surf.Flux derivative relative to Tsurf
233      - calculate clear-sky radiation & surface temp. change (for diagnostics)
234    
235    checkpoint53e_post
236    o merged cylindrical coord configuration
237    o added rotating_tank experiment
238    checkpoint53e_pre
239    
240    checkpoint53d_post
241    o faulks upgraded: re-generate output.txt that have changed:
242       global_ocean.cs32x15(.thsice), global_ocean_ebm, lab_sea, natl_box
243    o Added Bryan and Lewis, 1979, vertical diffusivity profile
244    o pkg/land: fix problem in land_impl_temp.F when Ts arround freezing.
245    o added Sadourny discretization of Coriolis
246      - moved some PARAMETERS from mom_*_coriolis.F to PARAMS.h
247      - re-enabled use of omega3 (useAbsVorticity)
248    o added Leith (Phys. Fluids, 1968, 10, 1409-1416) variable viscosity for
249      vector-invariant equations
250      - note: does *NOT* yet use same variable visc. for no-slip BC's.
251    o recoded the lateral boundary conditions for non-hydrostatic w, this time
252      hopefully without new bugs
253      - still to do: biharmonic friction for non-hydrostatic w
254    o bug fix for mom_*_sidedrag: rA -> rAs,rAw
255    o regenerated output for global_ocean_ebm (using -ieee)
256    o fix mom_common/mom_*_sidedrag.F to recover numeric output prior to creation
257      of mom_common/
258      - simple re-arrangement of instructions within computing loop
259    checkpoint53d_pre
260    
261    checkpoint53c_post
262    o new version of coupler for AIM+thSIce and Ocean models
263     - new fields (SSS,h_MixLayer ...) are exchanged through the coupler
264     - land: change Evaporation of snow
265     - aim with thSIce: fix bug in snow-precip.
266    o changes to calc_gw.F
267     - bug fix for w^2 term near the bottom boundary
268     - (hopefully) improve the lateral slip boundary condtions for use with
269        partial/looped cells
270      - because the bug fix changes two verification experiments anyway (exp5
271        and plume on slope), change the lateral boundary condition from half slip
272        to the value of no_slip_sides
273    o ebm package wasnt quite up to the parallel task. Fixed that.
274    o Adjoint of EBM package + verification
275      * TLM OK
276      * ADM not sufficient accuracy yet (needs further testing)
277    
278    checkpoint53b_post
279    o Commiting new energy balance model to repository
280      * package is pkg/ebm
281      * verif. is verification/global_ocean_ebm
282      * references are in ebm_driver.F
283    C      * X. Wang, P. Stone and J. Marotzke, 1999:
284    C        Global thermohaline circulation. Part I:
285    C        Sensitivity to atmospheric moisture transport.
286    C        J. Climate 12(1), 71-82
287    C      * X. Wang, P. Stone and J. Marotzke, 1999:
288    C        Global thermohaline circulation. Part II:
289    C        Sensitivity with interactive transport.
290    C        J. Climate 12(1), 83-91
291    C      * M. Nakamura, P. Stone and J. Marotzke, 1994:
292    C        Destabilization of the thermohaline circulation
293    C        by atmospheric eddy transports.
294    C        J. Climate 7(12), 1870-1882
295    o Will need long integration testing.
296    checkpoint53b_pre
297    
298    o created pkg/mom_common in preparation for variable viscosity code
299    o renamed and split ini_mnc_io
300    o pkg/land: only liquid water can diffuse or run-off
301    o AIM using land: soil wetness availability account for snow
302    o added printMapIncludesZeros to EEPARMS namelist
303    
304    checkpoint53a_post
305    o restore TAMC compatiility for exf_constants.h
306      (no usage of _d in .h)
307    o Additional storing in KPP in an attempt to fix
308      TAMC problem with adjoint of KPP code.
309      (no success...)
310      Identical code works fine with TAF.
311    o Fixed verification/natl_box_adjoint experiments.
312      Initial conditions for Experiments 4 and 5 are now
313      set to "SEAICE_initialHEFF=0", which is more reasonable,
314      and hence gives results closer to that of Experiment 2.
315      Experiment 6 requires "SEAICE_initialHEFF=1" and is working once again.
316    
317    checkpoint53
318    o stable checkpoint with up to date adjoint
319      TLM tests only run for 1st grdchk at the moment. Its unclear if
320      this is a code problem or a problem with resetting code for g_
321      variables that is being autogenerated wrong.
322      Otherwise:
323       o code is testing well on Linux (P4 and AMD), SGI Alitix and Irix.
324       o IBM and HPQ tests could not be run because NCAR (IBM) and NCCS (HPQ)
325         are both still in chaos following their security messes.
326    
327    checkpoint52n_post
328    o various build (mostly genmake2) changes including:
329      - ability to build/run more easily within Cygwin and MacOS X
330      - more complicated Makefiles with a stronger dependence on Gnu Make
331    o pkg/seaice modifications
332      - added SEAICE_deltaTtherm and SEAICE_deltaTdyn in order to reduce cost
333        of pkg/seaice.  For time being, SEAICE_deltaTtherm = deltaTtracer
334        and SEAICE_deltaTdyn must be integer multiple of SEAICE_deltaTtherm
335      - changed proxy for geostrophic velocity and updated
336        verification/lab_sea/results/* accordingly
337      - removed ADI dynamic solver from pkg/seaice
338      - replaced all occurrences of 1/CSTICE or 1/CSUICE
339        they caused a segmentation fault under certain circumstances
340      - removed unused field WATR
341    o bug fix in exf-package:
342      -  u/vwind are defined at tracer/pressure points. The values of u/vwindmask
343         now reflect that.
344    o small fixes in the cal-package:
345      - set reference date to first of January for the 'model'-calendar
346      - failing of cal_CheckDate in cal_FullDate now actually results in a
347        warning
348    o small modifications in exf-package to increase flexibility:
349      - change climtempfreeze from a constant to a runtime parameter that
350        is specified in data.exf_clim (default = -1.9)
351      - introduce runtime parameter windstressmax (default = 2.0) that is
352        specified in data.exf
353    o fix in pkg/thsice:
354      - update Evap (Tsf changes) to be consistent with Latent heat flux
355      - maintain snow precipitation unchanged (even if Tsf >=0).
356     both affect global_ocean.cs32x15/(.thsice) results ; update output.txt.thsice
357    
358    checkpoint52m_post
359    o commented out call to cal_GetMonthsRec in swfrac.F and added some comments;
360      this does not change any of the results, because myIter=0 in swfrac always
361      lead to jwtype=2 (the default anyway), but the call makes little sense
362      as a default
363    o no SHORTWAVE HEATING lost at the bottom; affects lab_sea output; updated
364    o re-work pkg/thsice interface with ocean & bulk_force pkg to allow
365      atmospheric model (AIM) to use it ; add a "slab ocean" component.
366    o pkg/aim_v23 : add interface to call thermodynamic sea-ice (thsice)
367    o update global_ocean.cs32x15/output (.thsice) after changes in pkg/thsice
368    o extending useSingleCpuIO option to work with new exch2 I/O format
369      - old-style, missing-tile I/O is still accessible by defining CPP
370        option MISSING_TILE_IO in pkg/mdsio/MDSIO_OPTIONS.h
371    o monitor package can write to NetCDF files and/or STDOUT
372    o fix re-definition bug in MNC_GRID_INIT_ALL()
373    o doc/api_reference added
374      - a framework for building an API Reference Manual using all
375        the "protex" comments embedded in the code
376      - see pkg/generic_advdiff, pkg/mnc for examples
377    o MNC:
378      - make myThid the last argument for all subroutines
379      - new default for the grid files: NO halos
380      - add useMNC and other run-time flags
381      - separate per-file-group handling of the unlimited dimension
382      - early initialization of MNC w/ flags read at the eeboot stage
383        and lookup tables set in initialize_fixed
384      - write all MNC output to a freshly-created directory
385    o thermodynamics sea-ice (thsice) is now tested with testreport.
386    o test exp. advect_cs: use the standard cs32 grid ; update output.
387    o Added functionality to grdchk:
388      pick global i,j,k position (or nearest wet) where to perform check.
389    
390    checkpoint52l_post
391    o new land formulation:
392      a) use ground enthalpy as prognostic variable to ensure exact
393         energy conservation.
394      b) account for water temperature and for latent heat of freezing
395         in all processes (rain, run-off, ground storage)
396      c) compute surface and ground temperature implicitly.
397    o aim_v23 modification for new land model:
398      a) Treat separately land / ocean / sea-ice surface fluxes
399         to allow implicit computation of land & sea-ice surface temp.
400      b) add snow precipitation.
401    checkpoint52l_pre
402    
403    o fix problem (when 1 tile is empty) in min/max monitor.
404    o fix bug in timeAve (u,vZetatave in common block)
405    o no restoring terms under ice for pkg/seaice
406    o moved SEAICE_DO_DIAGS to S/R do_the_model_io
407    o bug fix for pkg/seaice pickups
408    o change MNC call locations and add coordinate attributes
409    o mask vorticity at the corner ; affects global_ocean.cs output
410    o fix ini_curvilinear_grid ; test CS experiments with mpi
411    o aim.5l_cs: modify sea-ice fraction input file ; update output
412    o exf: Adding simple range check routine for exf fields after read
413    
414    checkpoint52k_post
415    o removed top-level directory "diags"
416      - we migrated these functions to pkg/timeave a long time ago
417    o fixed bug introduced in pkg/monitor
418    o accuracy ctrlprec = 32 insuffient for gradient checks using
419      averaged fields (I/O via cost_averages)
420      -> use ctrl.h in active_file*.F to control I/O precision.
421    o paramter list of CALL ADTHE_MAIN_LOOP
422    o added initial do_the_model_io to ecco
423    o set all ...Freq to zero after cost_final
424      taveFreq, dumpFreq, pChkptFreq
425    o MAIN and ECCO differentiability restored
426    
427    checkpoint52j_post
428    o fix nchklev in tamc.h for global_ocean...
429    o Bringing flow and store directives up-to-date
430    o Really weird thing!
431      call do_the_model_io from the_main_loop with parameters
432      starttime, niter0, instead of mytime, myiter, screws up the adjoint:
433      TAF places a faulty re-call of initialise_varia
434      after cost_final call which should not be there!
435    checkpoint52j_pre
436    
437    o tweaked pkg/mdsio/mdsio_readfield.F to have more useful logic and
438      error messages
439    o fix errorMessageUnit problem. remove old version of mdsio from eesup/src
440    o modified pkg/mdsio/mdsio_readfield.F -- commented out diskspace-eating
441      debug output
442    o modfied mdsio_writefield.F to handle non-global files under exch2
443    o update AIM Cubed-Sphere exp (aim.5l_cs):
444      use land, p*, and Franco.M forcing fields with fractional land-sea mask.
445    o update AIM Equatorial Channel exp:
446      use standard aim_v23 pkg and change the SST fields to be more realistic.
447    o globalFile now works with exch2
448    o merged mom_vi_calc_relvort3.F, mom_vi_del2uv.F and
449      mon_vort3.F to make exch2 friendly with >1 tile per cube face
450    o limit timeave output for hi-res integrations
451    
452    checkpoint52i_post
453    o removed CPP_EEOPTIONS.h from verification/*/code
454    o updated the tiling for global_ocean_pressure
455    checkpoint52i_pre
456    
457    o update aim.5l_LatLon using the standard aim_v23 pkg (instead of
458      the old pkg/aim).
459    
460    checkpoint52h_pre
461    o shortened offending lines in exch2_send_rx?.template
462    o modified W2_EXCH2_TOPOLOGY.h to reflect current topology generator
463    o genmake2: add build & execute test for NetCDF code
464    o modded aim.5l_cs and global_ocean.cs32x15 to use exch2, removed
465      old grid files, moved CPP_EEMACROS.h to pkg/exch2
466    o added useSingleCPUIO flag to monitor.F and mon_out.F
467    o added temporary hack for sea-ice dynamics on the cube:
468      metric terms are ignored and it is assumed that
469      there is no sea-ice on cube faces 1, 2, 4, and 5.
470    
471    checkpoint52f_post
472    o added exch2 cube routines to pkg/exch2
473    o modified hs94.cs-32x32x5 to use it
474    o added new grid files to hs94.cs-32x32x5
475    checkpoint52f_pre
476    
477    checkpoint52e_post
478    o first attempt to solve vertical advection (momemtum & T,S) Implicitly
479      for now, only implemented with Lax-Wendroff (tracer) and VectInv (mom)
480    o step.1 : rewrite (as in MultiDimAdvec) explicit tracer stepping (gad_calc_rhs.F)
481      to work with implicit vertical advection and AB ;
482    o Implic. vert. advect. implemented with flux-limiter, c2, u3 and c4 scheme.
483    o update output (gad_calc_rhs.F changes affect truncation error) of exp:
484     front_relax, global_ocean .cs32x15 & _pressure, internal_wave, lab_sea, natl_box
485    checkpoint52e_pre
486    
487    o fix bug in impldiff when Nr=1 or Nr=2.
488    o change the albedo in pkg/thsice (as GISS model); update output_thsice.txt
489    o little changes (fftpack, aim/phy_suflux ...) to pass the "checkbounds" test.
490    o fix deltaTfreesurf in update_cg2d.F (left from c45d_post): affects 2
491      ocean exp. that use NLFS: global_ocean.cs32x15 & global_ocean_pressure
492    o added exf_offset_atemp for Celsius to Kelvin conversion of atemp
493    o eesupp/inc/EEPARAMS.h: MAX_NO_PROCS=1024
494    o pkg/cal/cal_toseconds.F fix for long (>60 years) integrations
495    
496    checkpoint52d_post
497    o re-tiled many expts to allow comparison in multi-processor mode
498      - updated output accordingly
499      - changed hs94.128x64x5 to use intial T from file rather than use the
500        random number generator
501    o add infrastructure in model routines for fizhi and gridalt packages
502    o change global_ocean.cs32x15 because of Pb with GM on CS-grid.
503    o changes to permit hi-res, cubed-sphere, configuration
504      - added useSingleCpuIO capability to mdsio_readfield.F and exf_interp.F
505      - added "#undef ALLOW_USE_MPI" support to eesupp/src/scatter_2d.F
506      - added pkg/exf/exf_set_uv.F for on-the-fly interpolation and rotation
507        of surface winds for the cube
508      (A verification experiment, based on global_ocean.cs32x15, but with
509       pkg/seaice turned on and with on-the-fly interpolation from the NCEP
510       Gaussian grid is described in MITgcm_contrib/high_res_cube/README_ice,
511       complete with example surface forcing files and matlab scripts to look
512       at the output and compare it with that of global_ocean.cs32x15.)
513    o packages pickup: $PKG_write_pickup called from packages_write_pickup ;
514       but call to $PKG_read_pickup remains in $PKG_ini_vars: change CD_code
515       to conform to the common rule.
516    o check-in the first version of the coupler + interfaces for multi-components
517       set-up (Chris's version, adapted to work on cubed-sphere grid); use MPI.
518       - imply modification in eesupp (MPI_COMM_MODEL < MPI_COMM_WORLD).
519       - atmosphere: use pkg aim_v23 ; no sea-ice for now.
520       - ocean and atmos. share the same horizontal grid (no interpolation).
521       - exchanges between component and coupler is done in forward_step
522    o new verification experiment, cpl_aim+ocn:
523       coupled ocean - AIM atmosphere on cs32 grid.
524    
525    checkpoint52d_pre
526    o surface flux time-av diagnostic modified for NLFS & RealFreshWater
527    o add call to AIM pendant S/R in external_forcing & external_fields_load.F
528    
529    checkpoint52c_post
530    o Fixing OPTIONS files for use of exf:
531      - ALLOW_ECCO not needed
532      - ECCO_CPPOPTIONS.h should not be included from
533        EXF_OPTIONS.h/CAL_OPTIONS.h, but from CPP_OPTIONS.h
534        as has been the case in the past
535      - All cal/exf related options are currently set in
536        single header file. Should remain so until we
537        agree on a consistent new policy.
538    o Removed all CPP options for cal package.
539      cal has no more CPP options.
540    o modified STOREs in GAD_ADVECTION
541    o corrected key comp. for passkey
542    
543    checkpoint52b_post
544    o new version: "thsice" of thermodynamic sea-ice pkg:
545     - only 1 call from "forward_step"; conserve energy, fresh-water & salt.
546     - modify fluxes but do not change theta directly;
547     - light cleaning of bulk_force pkg (still needed with thsice).
548     - add forcing fields + data files in global_ocean.cs32x15/inp_thsice
549       + output_thsice.txt to test thsice & bulk_force pkgs.
550    o put back grid & coriolis monitor output.
551    checkpoint52b_pre
552    
553    o call "do_the_model_io" at the beginning (=> dump initial state)
554    o use rhoConstFresh to define convertEmP2rUnit (was done only in ocean-P)
555    o sort out interplay between tamc.h and PTRACERS_OPTIONS.h
556      (need info on NUMBER_OF_PTRACERS)
557    o missing PTRACERS headers in some routines
558    o default tamc.h in pkg/autodiff supposed to crash
559      to ensure proper customization
560    
561    checkpoint52a_post
562    o modifications to make FREEZE (_SURFACE) flux visible to pkg/kpp
563      - moved surfaceTendencyTice from pkg/seaice to main code
564      - subroutine FREEZE_SURFACE now limits only surface temperature
565        and is called in FORWARD_STEP (if allowFreezing=.TRUE.)
566        => new output.txt for global_ocean.cs32x15 & global_with_exf
567      - keep the old version of FREEZE for backward compatibility;
568        global_ocean.90x40x15 left unchanged using useOldFreezing=.TRUE.
569      - corresponding modifs for TAF directives
570    o added surface flux output variables to TIMEAVE_STATVARS
571    o time-averaged output for pkg/ptracers
572    o fix problem with ECCO-related divided adjoint
573      snapshot file was missing some ad-variables,
574      but no warning issued by TAF
575    
576    checkpoint52a_pre
577    o some bug fixes for #undef REAL4_IS_SLOW
578    
579    checkpoint52 = ecco_c52_e35
580    o after testing adjoints of checkpoint51u_post
581    
582    checkpoint51u_post
583    o this is the post-merge tag of the ecco-branch
584      (needs another clean checkout and re-testing)
585    o various natl_box_adjoint configs and results
586      but no input fields
587    o obcs adjoint not yet merged!
588    o adding ECCO global 1x1 and 2x2 degree configs
589    o updated/extended some store directives for seaice
590    o pkg/ecco now containes ecco-specific part of cost function
591    o top level routines the_main_loop, forward_step
592      supersede those in model/src/
593      previous input data.cost now in data.ecco
594      (new namelist ecco_cost_nml)
595    o in ctrl/ cleaned some cross-dependencies and updated CPP options
596    
597    checkpoint51t_post
598    o this is the pre-merge tag of the ecco-branch
599    
600    checkpoint51s_post
601    o cleanup: removed various obsolete #define-s
602    o pkg/grdchk
603      - renamed 'ph-grd 3' to 'grad-res'
604      - merged from ecco-branch
605        (remaining bug fixes for obcs gradient checks)
606      - additional high-precision output for testreport
607        (grep for 'precision_grdchk_result')
608    
609    checkpoint51r_post
610    o cleaned up or removed CPP_OPTIONS.h files
611    o updated results for some of the tests to reflect the new
612      optfiles (-ieee: -O0 instead of -O3)
613    
614    checkpoint51q_post
615    o cd_code:
616      - moved cd_scheme.F -> cd_code_scheme.F
617      - separate read_checkpoint from cd_code_ini_vars.F
618      - separated cd_code part from write_checkpoint
619      - updated AD_SOURCE, generated .flow
620      - added CD_CODE_VARS.h to the_main_loop
621      - Moving cd_code specific store dir.s to pkg/cd_code
622    o ctrl package totally restructured
623      - pack/unpack now optional and decoupled from
624        xx_/adxx_ I/O
625      - ctrl_pack/unpack cleaned
626        (new routines ctrl_init_ctrlvar.F, pkg/ctrl/ctrl_init_wet.F)
627      - confined inclusion of AD_CONFIG.h to where necessary.
628    o genmake2: changed AD dependencies from SRCFILES to AD_FILES
629    
630    checkpoint51p_post
631    o create initial cd_code package
632      - WARNING: the initialization of fields within packages
633        is broken because INI_FIELDS (where READ_CHECKPOINT is
634        called) is called *before* PACKAGES_INIT_VARIABLES.
635        This seems to be backwards since
636        1) its important to initially zero the entire field including
637           the "halo" region
638        2) zeroing needs to be done *before* reading pickup info
639           and the pickup call might or might not happen in any
640           particular run
641      - this needs an overhaul for R2
642    o cleanup verification/*/code/GMREDI_OPTIONS.h
643    o cleanup testreport
644    
645    checkpoint51o_post
646    o bringing AD-related verif's up to date with genmake2
647      - carbon/
648      - global_ocean.90x40x15/
649      - hs.94.1x64x5/
650      - (front_relax/ )
651    o (more) consistent directory structure for AD-related verif.'s
652    o cleaning ALLOW_GRADIENT_CHECK -> ALLOW_GRDCHK
653    o cleaning some ALLOW_TANGENTLINEAR_RUN -> ALLOW_AUTODIFF
654    o bug fix in find_alpha.F for MDJWF:
655      - modif. to alpha = 1/D*( dN/dT - rho*dD/Dt) to account for
656        change rho -> rho-rhoConst
657      - replace call find_rho to find_rhonum
658    o initialisation of rFlx extended to full array (required by TAF)
659      and shifted to thermodynamics
660    o removed PTRACERS.h in ptracers routine
661    o added surfacetendencyPtr to S/R parameter list pracers_forcing
662    o Correcting initialisations for fVerU, fVerV.
663    o helping TAF to recognize dependency on kUp for
664      fVerU, fVerV
665    o Replacing delZ by delR in pkg/ctrl/
666    checkpoint51o_pre
667    
668    o undid all of cp51 pending some cleanups and discussion
669    o added the [#include "AD_CONFIG.h"] statement to all files that need
670       it for adjoint/tl #defines
671     o re-worked the build logic in genmake2 to support AD_CONFIG.h
672     o removed tools/genmake since it no longer works
673    
674    checkpoint51n_post
675    o modifications to make FREEZE flux visible to pkg/kpp
676      - moved surfaceTendencyTice from pkg/seaice to main code
677      - FREEZE & EXTERNAL_FORCING_SURF moved to FORWARD_STEP
678      - subroutine FREEZE now limits only surface temperature
679        (this means new output.txt for global_ocean.90x40x15,
680         global_ocean.cs32x15, and global_with_exf)
681    o added surface flux output variables to TIMEAVE_STATVARS
682    
683    checkpoint51n_pre
684    o added "#undef ALLOW_USE_MPI" support to eesupp/src/gather_2d.F
685    o added [#include "AD_CONFIG.h"] to all files that use the following
686      defines:
687      - ALLOW_ADJOINT_RUN
688      - ALLOW_TANGENTLINEAR_RUN
689      - ALLOW_ECCO_OPTIMIZATION
690    
691    checkpoint51m_post
692    o incorporating changes to get TAMC/TAF working in genmake2
693      - added new TAMC/TAF targets to the Makefile
694      - various cleanups related to package handling
695      - these changes should have no impact on the forward model
696    
697    checkpoint51l_post
698    o changes to pkg/cal and pkg/exf that allow and speed-up long integrations
699     - faster algorithm for cal_addtime.F
700     - higher precision arithmetic, where needed
701     - pkg/exf *startdate variables changed from calendar date to
702       model integration time
703    
704    checkpoint51l_pre
705    o add salt flux to forcing fields (used when salty sea-ice forms or melts)
706    o enable to apply restoring on S & T only in a latitude band.
707    o update output (external_forcing_surf modif. change truncation errors)
708       global_ocean.cs32x15, global_ocean_pressure, lab_sea, natl_box
709    
710    checkpoint51k_post
711    o change default to: useNHMTerms=.FALSE.
712     - new output.txt for atmos. exp (NHMTerms were wrong but now turned off):
713       adjustment.128x64x1, aim.5l_Eq_Ch, aim.5l_LatLon, hs94.128x64x5, hs94.1x64x5
714     - new output.txt for idealized oceanic exp (NHMTerms are now turned off):
715        exp1, ideal_2D_oce
716     - change data file of oceanic exp (leave the results unchanged for now)
717       exp2, global_ocean[90x40x15,pressure,with_exf], lab_sea, natl_box
718    o set gravitySign automatically according to the vertical coordinate used
719      (and change sign: back to "logical" value: +1 when R increases downward)
720    o added an include barrier to the default CPP_OPTIONS.h and
721      removed the #include "PACKAGES_CONFIG.h" from it
722    o fix a bug in mom_fluxform (from chk51j_post)
723    
724    checkpoint51j_post
725    o Added some AD-related initialisations in mom_vecinv/ mom_fluxform/
726    o adjusted some flow directives (dynamics, thermodynamics)
727    o added cubed-sphere case to test cost function
728    o heckpoint_lev?_directives.h now include
729      pkg-specific directives, whenever those are enabled.
730    o Moving pkg-specific flow directives to corresponding pkg
731      Moving pkg-specific differentiation list to corresponding pkg
732      Moving pkg-specific checkpoint-lev storing to corresponding pkg
733      ---> conventions are
734           - pkg_ad.flow
735           - pkg_ad_diff.list
736           - pkg_ad_check_lev?_dir.h
737    
738    checkpoint51i_post
739    o merge of the branch-genmake2
740      - knowingly breaks the older genmake
741      - automatic generation of all ALLOW_${PKG_NAME} defines
742      - new testreport supersedes older testscript
743    o large number of package inclusion and header cleanups
744    checkpoint51i_pre
745    
746    o new output for adjoin of global ocean vector invariant
747    o fix for checkpoint of thermodynamic seaice and bulk_force
748    o Steph gchem various cleanups
749    
750    checkpoint 51h_pre
751    o fixed some single quotes
752    o modify gchem pkg to be more versatile
753    o modifications to clean up dic pkg
754    o add pkg/cfc
755    
756    checkpoint51g_post
757    o provide links for ARPACK and ADM*TLM
758    o start removing verification/vero/
759    o update AD-related stuff
760      (in particular some IF-statements in model/src/ routines)
761    o adding to verification/global_ocean.90x40x15/ adjoint of
762      vector invariant setup
763    
764    checkpoint51f_post
765    o update advect_xz/output.txt (left from chekpoint50f_pre, ini_vel modif)
766  o Mods and bug fixes to pkg/cal, pkg/exf, etc., needed for computation  o Mods and bug fixes to pkg/cal, pkg/exf, etc., needed for computation
767    of tracer Green's fucntions for ocean inversion project.    of tracer Green's fucntions for ocean inversion project.
768    

Legend:
Removed from v.1.178  
changed lines
  Added in v.1.370

  ViewVC Help
Powered by ViewVC 1.1.22