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

Legend:
Removed from v.1.157  
changed lines
  Added in v.1.356

  ViewVC Help
Powered by ViewVC 1.1.22