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

Legend:
Removed from v.1.163  
changed lines
  Added in v.1.375

  ViewVC Help
Powered by ViewVC 1.1.22