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

Legend:
Removed from v.1.171  
changed lines
  Added in v.1.369

  ViewVC Help
Powered by ViewVC 1.1.22