/[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.128 by heimbach, Thu Feb 6 22:43:09 2003 UTC revision 1.288 by edhill, Thu May 6 02:42:49 2004 UTC
# Line 3  $Header$ Line 3  $Header$
3      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
4      ==============================      ==============================
5    
6    checkpoint52n_post
7    o various build (mostly genmake2) changes including:
8      - ability to build/run more easily within Cygwin and MacOS X
9      - more complicated Makefiles with a stronger dependence on Gnu Make
10    o pkg/seaice modifications
11      - added SEAICE_deltaTtherm and SEAICE_deltaTdyn in order to reduce cost
12        of pkg/seaice.  For time being, SEAICE_deltaTtherm = deltaTtracer
13        and SEAICE_deltaTdyn must be integer multiple of SEAICE_deltaTtherm
14      - changed proxy for geostrophic velocity and updated
15        verification/lab_sea/results/* accordingly
16      - removed ADI dynamic solver from pkg/seaice
17      - replaced all occurrences of 1/CSTICE or 1/CSUICE
18        they caused a segmentation fault under certain circumstances
19      - removed unused field WATR
20    o bug fix in exf-package:
21      -  u/vwind are defined at tracer/pressure points. The values of u/vwindmask
22         now reflect that.
23    o small fixes in the cal-package:
24      - set reference date to first of January for the 'model'-calendar
25      - failing of cal_CheckDate in cal_FullDate now actually results in a
26        warning
27    o small modifications in exf-package to increase flexibility:
28      - change climtempfreeze from a constant to a runtime parameter that
29        is specified in data.exf_clim (default = -1.9)
30      - introduce runtime parameter windstressmax (default = 2.0) that is
31        specified in data.exf
32    o fix in pkg/thsice:
33      - update Evap (Tsf changes) to be consistent with Latent heat flux
34      - maintain snow precipitation unchanged (even if Tsf >=0).
35     both affect global_ocean.cs32x15/(.thsice) results ; update output.txt.thsice
36    
37    checkpoint52m_post
38    o commented out call to cal_GetMonthsRec in swfrac.F and added some comments;
39      this does not change any of the results, because myIter=0 in swfrac always
40      lead to jwtype=2 (the default anyway), but the call makes little sense
41      as a default
42    o no SHORTWAVE HEATING lost at the bottom; affects lab_sea output; updated
43    o re-work pkg/thsice interface with ocean & bulk_force pkg to allow
44      atmospheric model (AIM) to use it ; add a "slab ocean" component.
45    o pkg/aim_v23 : add interface to call thermodynamic sea-ice (thsice)
46    o update global_ocean.cs32x15/output (.thsice) after changes in pkg/thsice
47    o extending useSingleCpuIO option to work with new exch2 I/O format
48      - old-style, missing-tile I/O is still accessible by defining CPP
49        option MISSING_TILE_IO in pkg/mdsio/MDSIO_OPTIONS.h
50    o monitor package can write to NetCDF files and/or STDOUT
51    o fix re-definition bug in MNC_GRID_INIT_ALL()
52    o doc/api_reference added
53      - a framework for building an API Reference Manual using all
54        the "protex" comments embedded in the code
55      - see pkg/generic_advdiff, pkg/mnc for examples
56    o MNC:
57      - make myThid the last argument for all subroutines
58      - new default for the grid files: NO halos
59      - add useMNC and other run-time flags
60      - separate per-file-group handling of the unlimited dimension
61      - early initialization of MNC w/ flags read at the eeboot stage
62        and lookup tables set in initialize_fixed
63      - write all MNC output to a freshly-created directory
64    o thermodynamics sea-ice (thsice) is now tested with testreport.
65    o test exp. advect_cs: use the standard cs32 grid ; update output.
66    o Added functionality to grdchk:
67      pick global i,j,k position (or nearest wet) where to perform check.
68    
69    checkpoint52l_post
70    o new land formulation:
71      a) use ground enthalpy as prognostic variable to ensure exact
72         energy conservation.
73      b) account for water temperature and for latent heat of freezing
74         in all processes (rain, run-off, ground storage)
75      c) compute surface and ground temperature implicitly.
76    o aim_v23 modification for new land model:
77      a) Treat separately land / ocean / sea-ice surface fluxes
78         to allow implicit computation of land & sea-ice surface temp.
79      b) add snow precipitation.
80    checkpoint52l_pre
81    
82    o fix problem (when 1 tile is empty) in min/max monitor.
83    o fix bug in timeAve (u,vZetatave in common block)
84    o no restoring terms under ice for pkg/seaice
85    o moved SEAICE_DO_DIAGS to S/R do_the_model_io
86    o bug fix for pkg/seaice pickups
87    o change MNC call locations and add coordinate attributes
88    o mask vorticity at the corner ; affects global_ocean.cs output
89    o fix ini_curvilinear_grid ; test CS experiments with mpi
90    o aim.5l_cs: modify sea-ice fraction input file ; update output
91    o exf: Adding simple range check routine for exf fields after read
92    
93    checkpoint52k_post
94    o removed top-level directory "diags"
95      - we migrated these functions to pkg/timeave a long time ago
96    o fixed bug introduced in pkg/monitor
97    o accuracy ctrlprec = 32 insuffient for gradient checks using
98      averaged fields (I/O via cost_averages)
99      -> use ctrl.h in active_file*.F to control I/O precision.
100    o paramter list of CALL ADTHE_MAIN_LOOP
101    o added initial do_the_model_io to ecco
102    o set all ...Freq to zero after cost_final
103      taveFreq, dumpFreq, pChkptFreq
104    o MAIN and ECCO differentiability restored
105    
106    checkpoint52j_post
107    o fix nchklev in tamc.h for global_ocean...
108    o Bringing flow and store directives up-to-date
109    o Really weird thing!
110      call do_the_model_io from the_main_loop with parameters
111      starttime, niter0, instead of mytime, myiter, screws up the adjoint:
112      TAF places a faulty re-call of initialise_varia
113      after cost_final call which should not be there!
114    checkpoint52j_pre
115    
116    o tweaked pkg/mdsio/mdsio_readfield.F to have more useful logic and
117      error messages
118    o fix errorMessageUnit problem. remove old version of mdsio from eesup/src
119    o modified pkg/mdsio/mdsio_readfield.F -- commented out diskspace-eating
120      debug output
121    o modfied mdsio_writefield.F to handle non-global files under exch2
122    o update AIM Cubed-Sphere exp (aim.5l_cs):
123      use land, p*, and Franco.M forcing fields with fractional land-sea mask.
124    o update AIM Equatorial Channel exp:
125      use standard aim_v23 pkg and change the SST fields to be more realistic.
126    o globalFile now works with exch2
127    o merged mom_vi_calc_relvort3.F, mom_vi_del2uv.F and
128      mon_vort3.F to make exch2 friendly with >1 tile per cube face
129    o limit timeave output for hi-res integrations
130    
131    checkpoint52i_post
132    o removed CPP_EEOPTIONS.h from verification/*/code
133    o updated the tiling for global_ocean_pressure
134    checkpoint52i_pre
135    
136    o update aim.5l_LatLon using the standard aim_v23 pkg (instead of
137      the old pkg/aim).
138    
139    checkpoint52h_pre
140    o shortened offending lines in exch2_send_rx?.template
141    o modified W2_EXCH2_TOPOLOGY.h to reflect current topology generator
142    o genmake2: add build & execute test for NetCDF code
143    o modded aim.5l_cs and global_ocean.cs32x15 to use exch2, removed
144      old grid files, moved CPP_EEMACROS.h to pkg/exch2
145    o added useSingleCPUIO flag to monitor.F and mon_out.F
146    o added temporary hack for sea-ice dynamics on the cube:
147      metric terms are ignored and it is assumed that
148      there is no sea-ice on cube faces 1, 2, 4, and 5.
149    
150    checkpoint52f_post
151    o added exch2 cube routines to pkg/exch2
152    o modified hs94.cs-32x32x5 to use it
153    o added new grid files to hs94.cs-32x32x5
154    checkpoint52f_pre
155    
156    checkpoint52e_post
157    o first attempt to solve vertical advection (momemtum & T,S) Implicitly
158      for now, only implemented with Lax-Wendroff (tracer) and VectInv (mom)
159    o step.1 : rewrite (as in MultiDimAdvec) explicit tracer stepping (gad_calc_rhs.F)
160      to work with implicit vertical advection and AB ;
161    o Implic. vert. advect. implemented with flux-limiter, c2, u3 and c4 scheme.
162    o update output (gad_calc_rhs.F changes affect truncation error) of exp:
163     front_relax, global_ocean .cs32x15 & _pressure, internal_wave, lab_sea, natl_box
164    checkpoint52e_pre
165    
166    o fix bug in impldiff when Nr=1 or Nr=2.
167    o change the albedo in pkg/thsice (as GISS model); update output_thsice.txt
168    o little changes (fftpack, aim/phy_suflux ...) to pass the "checkbounds" test.
169    o fix deltaTfreesurf in update_cg2d.F (left from c45d_post): affects 2
170      ocean exp. that use NLFS: global_ocean.cs32x15 & global_ocean_pressure
171    o added exf_offset_atemp for Celsius to Kelvin conversion of atemp
172    o eesupp/inc/EEPARAMS.h: MAX_NO_PROCS=1024
173    o pkg/cal/cal_toseconds.F fix for long (>60 years) integrations
174    
175    checkpoint52d_post
176    o re-tiled many expts to allow comparison in multi-processor mode
177      - updated output accordingly
178      - changed hs94.128x64x5 to use intial T from file rather than use the
179        random number generator
180    o add infrastructure in model routines for fizhi and gridalt packages
181    o change global_ocean.cs32x15 because of Pb with GM on CS-grid.
182    o changes to permit hi-res, cubed-sphere, configuration
183      - added useSingleCpuIO capability to mdsio_readfield.F and exf_interp.F
184      - added "#undef ALLOW_USE_MPI" support to eesupp/src/scatter_2d.F
185      - added pkg/exf/exf_set_uv.F for on-the-fly interpolation and rotation
186        of surface winds for the cube
187      (A verification experiment, based on global_ocean.cs32x15, but with
188       pkg/seaice turned on and with on-the-fly interpolation from the NCEP
189       Gaussian grid is described in MITgcm_contrib/high_res_cube/README_ice,
190       complete with example surface forcing files and matlab scripts to look
191       at the output and compare it with that of global_ocean.cs32x15.)
192    o packages pickup: $PKG_write_pickup called from packages_write_pickup ;
193       but call to $PKG_read_pickup remains in $PKG_ini_vars: change CD_code
194       to conform to the common rule.
195    o check-in the first version of the coupler + interfaces for multi-components
196       set-up (Chris's version, adapted to work on cubed-sphere grid); use MPI.
197       - imply modification in eesupp (MPI_COMM_MODEL < MPI_COMM_WORLD).
198       - atmosphere: use pkg aim_v23 ; no sea-ice for now.
199       - ocean and atmos. share the same horizontal grid (no interpolation).
200       - exchanges between component and coupler is done in forward_step
201    o new verification experiment, cpl_aim+ocn:
202       coupled ocean - AIM atmosphere on cs32 grid.
203    
204    checkpoint52d_pre
205    o surface flux time-av diagnostic modified for NLFS & RealFreshWater
206    o add call to AIM pendant S/R in external_forcing & external_fields_load.F
207    
208    checkpoint52c_post
209    o Fixing OPTIONS files for use of exf:
210      - ALLOW_ECCO not needed
211      - ECCO_CPPOPTIONS.h should not be included from
212        EXF_OPTIONS.h/CAL_OPTIONS.h, but from CPP_OPTIONS.h
213        as has been the case in the past
214      - All cal/exf related options are currently set in
215        single header file. Should remain so until we
216        agree on a consistent new policy.
217    o Removed all CPP options for cal package.
218      cal has no more CPP options.
219    o modified STOREs in GAD_ADVECTION
220    o corrected key comp. for passkey
221    
222    checkpoint52b_post
223    o new version: "thsice" of thermodynamic sea-ice pkg:
224     - only 1 call from "forward_step"; conserve energy, fresh-water & salt.
225     - modify fluxes but do not change theta directly;
226     - light cleaning of bulk_force pkg (still needed with thsice).
227     - add forcing fields + data files in global_ocean.cs32x15/inp_thsice
228       + output_thsice.txt to test thsice & bulk_force pkgs.
229    o put back grid & coriolis monitor output.
230    checkpoint52b_pre
231    
232    o call "do_the_model_io" at the beginning (=> dump initial state)
233    o use rhoConstFresh to define convertEmP2rUnit (was done only in ocean-P)
234    o sort out interplay between tamc.h and PTRACERS_OPTIONS.h
235      (need info on NUMBER_OF_PTRACERS)
236    o missing PTRACERS headers in some routines
237    o default tamc.h in pkg/autodiff supposed to crash
238      to ensure proper customization
239    
240    checkpoint52a_post
241    o modifications to make FREEZE (_SURFACE) flux visible to pkg/kpp
242      - moved surfaceTendencyTice from pkg/seaice to main code
243      - subroutine FREEZE_SURFACE now limits only surface temperature
244        and is called in FORWARD_STEP (if allowFreezing=.TRUE.)
245        => new output.txt for global_ocean.cs32x15 & global_with_exf
246      - keep the old version of FREEZE for backward compatibility;
247        global_ocean.90x40x15 left unchanged using useOldFreezing=.TRUE.
248      - corresponding modifs for TAF directives
249    o added surface flux output variables to TIMEAVE_STATVARS
250    o time-averaged output for pkg/ptracers
251    o fix problem with ECCO-related divided adjoint
252      snapshot file was missing some ad-variables,
253      but no warning issued by TAF
254    
255    checkpoint52a_pre
256    o some bug fixes for #undef REAL4_IS_SLOW
257    
258    checkpoint52 = ecco_c52_e35
259    o after testing adjoints of checkpoint51u_post
260    
261    checkpoint51u_post
262    o this is the post-merge tag of the ecco-branch
263      (needs another clean checkout and re-testing)
264    o various natl_box_adjoint configs and results
265      but no input fields
266    o obcs adjoint not yet merged!
267    o adding ECCO global 1x1 and 2x2 degree configs
268    o updated/extended some store directives for seaice
269    o pkg/ecco now containes ecco-specific part of cost function
270    o top level routines the_main_loop, forward_step
271      supersede those in model/src/
272      previous input data.cost now in data.ecco
273      (new namelist ecco_cost_nml)
274    o in ctrl/ cleaned some cross-dependencies and updated CPP options
275    
276    checkpoint51t_post
277    o this is the pre-merge tag of the ecco-branch
278    
279    checkpoint51s_post
280    o cleanup: removed various obsolete #define-s
281    o pkg/grdchk
282      - renamed 'ph-grd 3' to 'grad-res'
283      - merged from ecco-branch
284        (remaining bug fixes for obcs gradient checks)
285      - additional high-precision output for testreport
286        (grep for 'precision_grdchk_result')
287    
288    checkpoint51r_post
289    o cleaned up or removed CPP_OPTIONS.h files
290    o updated results for some of the tests to reflect the new
291      optfiles (-ieee: -O0 instead of -O3)
292    
293    checkpoint51q_post
294    o cd_code:
295      - moved cd_scheme.F -> cd_code_scheme.F
296      - separate read_checkpoint from cd_code_ini_vars.F
297      - separated cd_code part from write_checkpoint
298      - updated AD_SOURCE, generated .flow
299      - added CD_CODE_VARS.h to the_main_loop
300      - Moving cd_code specific store dir.s to pkg/cd_code
301    o ctrl package totally restructured
302      - pack/unpack now optional and decoupled from
303        xx_/adxx_ I/O
304      - ctrl_pack/unpack cleaned
305        (new routines ctrl_init_ctrlvar.F, pkg/ctrl/ctrl_init_wet.F)
306      - confined inclusion of AD_CONFIG.h to where necessary.
307    o genmake2: changed AD dependencies from SRCFILES to AD_FILES
308    
309    checkpoint51p_post
310    o create initial cd_code package
311      - WARNING: the initialization of fields within packages
312        is broken because INI_FIELDS (where READ_CHECKPOINT is
313        called) is called *before* PACKAGES_INIT_VARIABLES.
314        This seems to be backwards since
315        1) its important to initially zero the entire field including
316           the "halo" region
317        2) zeroing needs to be done *before* reading pickup info
318           and the pickup call might or might not happen in any
319           particular run
320      - this needs an overhaul for R2
321    o cleanup verification/*/code/GMREDI_OPTIONS.h
322    o cleanup testreport
323    
324    checkpoint51o_post
325    o bringing AD-related verif's up to date with genmake2
326      - carbon/
327      - global_ocean.90x40x15/
328      - hs.94.1x64x5/
329      - (front_relax/ )
330    o (more) consistent directory structure for AD-related verif.'s
331    o cleaning ALLOW_GRADIENT_CHECK -> ALLOW_GRDCHK
332    o cleaning some ALLOW_TANGENTLINEAR_RUN -> ALLOW_AUTODIFF
333    o bug fix in find_alpha.F for MDJWF:
334      - modif. to alpha = 1/D*( dN/dT - rho*dD/Dt) to account for
335        change rho -> rho-rhoConst
336      - replace call find_rho to find_rhonum
337    o initialisation of rFlx extended to full array (required by TAF)
338      and shifted to thermodynamics
339    o removed PTRACERS.h in ptracers routine
340    o added surfacetendencyPtr to S/R parameter list pracers_forcing
341    o Correcting initialisations for fVerU, fVerV.
342    o helping TAF to recognize dependency on kUp for
343      fVerU, fVerV
344    o Replacing delZ by delR in pkg/ctrl/
345    checkpoint51o_pre
346    
347    o undid all of cp51 pending some cleanups and discussion
348    o added the [#include "AD_CONFIG.h"] statement to all files that need
349       it for adjoint/tl #defines
350     o re-worked the build logic in genmake2 to support AD_CONFIG.h
351     o removed tools/genmake since it no longer works
352    
353    checkpoint51n_post
354    o modifications to make FREEZE flux visible to pkg/kpp
355      - moved surfaceTendencyTice from pkg/seaice to main code
356      - FREEZE & EXTERNAL_FORCING_SURF moved to FORWARD_STEP
357      - subroutine FREEZE now limits only surface temperature
358        (this means new output.txt for global_ocean.90x40x15,
359         global_ocean.cs32x15, and global_with_exf)
360    o added surface flux output variables to TIMEAVE_STATVARS
361    
362    checkpoint51n_pre
363    o added "#undef ALLOW_USE_MPI" support to eesupp/src/gather_2d.F
364    o added [#include "AD_CONFIG.h"] to all files that use the following
365      defines:
366      - ALLOW_ADJOINT_RUN
367      - ALLOW_TANGENTLINEAR_RUN
368      - ALLOW_ECCO_OPTIMIZATION
369    
370    checkpoint51m_post
371    o incorporating changes to get TAMC/TAF working in genmake2
372      - added new TAMC/TAF targets to the Makefile
373      - various cleanups related to package handling
374      - these changes should have no impact on the forward model
375    
376    checkpoint51l_post
377    o changes to pkg/cal and pkg/exf that allow and speed-up long integrations
378     - faster algorithm for cal_addtime.F
379     - higher precision arithmetic, where needed
380     - pkg/exf *startdate variables changed from calendar date to
381       model integration time
382    
383    checkpoint51l_pre
384    o add salt flux to forcing fields (used when salty sea-ice forms or melts)
385    o enable to apply restoring on S & T only in a latitude band.
386    o update output (external_forcing_surf modif. change truncation errors)
387       global_ocean.cs32x15, global_ocean_pressure, lab_sea, natl_box
388    
389    checkpoint51k_post
390    o change default to: useNHMTerms=.FALSE.
391     - new output.txt for atmos. exp (NHMTerms were wrong but now turned off):
392       adjustment.128x64x1, aim.5l_Eq_Ch, aim.5l_LatLon, hs94.128x64x5, hs94.1x64x5
393     - new output.txt for idealized oceanic exp (NHMTerms are now turned off):
394        exp1, ideal_2D_oce
395     - change data file of oceanic exp (leave the results unchanged for now)
396       exp2, global_ocean[90x40x15,pressure,with_exf], lab_sea, natl_box
397    o set gravitySign automatically according to the vertical coordinate used
398      (and change sign: back to "logical" value: +1 when R increases downward)
399    o added an include barrier to the default CPP_OPTIONS.h and
400      removed the #include "PACKAGES_CONFIG.h" from it
401    o fix a bug in mom_fluxform (from chk51j_post)
402    
403    checkpoint51j_post
404    o Added some AD-related initialisations in mom_vecinv/ mom_fluxform/
405    o adjusted some flow directives (dynamics, thermodynamics)
406    o added cubed-sphere case to test cost function
407    o heckpoint_lev?_directives.h now include
408      pkg-specific directives, whenever those are enabled.
409    o Moving pkg-specific flow directives to corresponding pkg
410      Moving pkg-specific differentiation list to corresponding pkg
411      Moving pkg-specific checkpoint-lev storing to corresponding pkg
412      ---> conventions are
413           - pkg_ad.flow
414           - pkg_ad_diff.list
415           - pkg_ad_check_lev?_dir.h
416    
417    checkpoint51i_post
418    o merge of the branch-genmake2
419      - knowingly breaks the older genmake
420      - automatic generation of all ALLOW_${PKG_NAME} defines
421      - new testreport supersedes older testscript
422    o large number of package inclusion and header cleanups
423    checkpoint51i_pre
424    
425    o new output for adjoin of global ocean vector invariant
426    o fix for checkpoint of thermodynamic seaice and bulk_force
427    o Steph gchem various cleanups
428    
429    checkpoint 51h_pre
430    o fixed some single quotes
431    o modify gchem pkg to be more versatile
432    o modifications to clean up dic pkg
433    o add pkg/cfc
434    
435    checkpoint51g_post
436    o provide links for ARPACK and ADM*TLM
437    o start removing verification/vero/
438    o update AD-related stuff
439      (in particular some IF-statements in model/src/ routines)
440    o adding to verification/global_ocean.90x40x15/ adjoint of
441      vector invariant setup
442    
443    checkpoint51f_post
444    o update advect_xz/output.txt (left from chekpoint50f_pre, ini_vel modif)
445    o Mods and bug fixes to pkg/cal, pkg/exf, etc., needed for computation
446      of tracer Green's fucntions for ocean inversion project.
447    
448    checkpoint51f_pre
449    o Added on-the-fly spatial interpolation capability
450        "USE_EXF_INTERPOLATION" to pkg/exf.
451      - This is a temporary Cartesian-grid hack until
452        the super-duper ESMF coupler becomes available.
453      - See verification/global_with_exf/README for usage example.
454      - Removed obsolete EXFwindOnBgrid and SEAICEwindOnCgrid
455        flags and modified pkg/seaice accordingly.
456    o Bug fix to pkg/ptracers, pkg/generic_advdiff/gad_calc_rhs.F,
457        and pkg/kpp/kpp_transport_ptr.F for dealing with tracer
458        non-local transport term.
459    
460    checkpoint51e_post
461    o pkg/mom_vecinv:
462      - changes related to hFac: use S/R from aim.5l_cs/code (coriolis & vort.advec)
463      - add time-ave diagnostic of vorticity advection.
464    o add diagnostic (instantaneous) of Shapiro Filter effect for T,S & UV.
465    o atmospheric geopotential: - include water vapor effect;
466                                - use p^kappa in p* scaling.
467    o land pkg: read initial state from files.
468    o AIM : add stratospheric drag in the upper level
469    
470    checkpoint51d_post
471    o change function ILNBLNK (=> last non-blank char.) to return zero
472      (instead of the string length) when the string. char is empty
473    o fix bug in dxV,dyU definition of cubed-sphere grid (ini_curvilinear_grid)
474    o modif's to restore TAMC compatibility
475      - avoid "_d" constructs in PARAMETER declarations
476      - avoid interrupting lines by comment lines
477    o Use cluster local disks for purely local I/O
478      vs. globally visible disks needed for ctrl stuff
479      (merged from ecco-branch)
480      - new active_file_loc... and mdsio_..._loc routines
481      - affected packages: autodiff, ctrl, grdchk, mdsio
482    
483    checkpoint51c_post
484    o introducing integer flag debugLevel
485      - debLevZero: suppress all output
486      - debLevA: standard
487      - debLevB: equivalent to current debugMode=.TRUE.
488    o introducing pathname variable mdsioLocalDir for mdsio
489      that can be specified at runtime
490      (relevant for cluster I/O to local disk)
491    `
492    >>>checkpoint51b_post
493    adjoint of ptracers
494    o disentangled ALLOW_PTRACERS using new ALLOW_GCHEM
495    o extended state in checkpoint_lev?, etc.
496    o alternative ctrl/cost with tr1 or ptracer
497    o new comlev1_bibj_ptracers
498    o new ptracers.flow
499    o made ptracers_forcing same as external_forcing_s
500    o added verification for adjoint with ptracers
501      - adm. grdchk fail so far
502      - tlm. grdchk OK
503      - checked tapelev/comlev recomputations: OK
504    <<<checkpoint51b_pre
505    
506    checkpoint51a_post
507    o added code for biogeochemistry package
508      - new pkgs: dic, gchem
509      - new subroutine in ptracers
510    o modify code for biogeo packages
511      - in model/src and pkg/ptracers
512    
513    checkpoint51
514    o merged improved packages from ecco-branch (ecco_c51_e34) to MAIN
515      - autodiff, ctrl, exf, grdchk, seaice
516      kept separate packages
517      - cost, ecco, obcs
518    o updating verification/internal_wave, verification/lab_sea
519      (all verif's pass on faulks.lcs.mit.edu RH7.3, g77)
520    
521    checkpoint50i_post
522    o Preparing next differentiable checkpoint and sync
523      of MAIN vs. ecco-branch
524      - 1st step: this one is differentiable
525      - 2nd step: next(! not this) checkpoint will have updated packages
526                  from ecco-branch)
527    
528    checkpoint50h_post
529    o new package "land": provide surface BC for AIM.
530    o pkg/ptracers: apply exchange to all tracers after reading a pickup.
531    o AIM: use true log-P extrapolation for near surface temp (= as in
532        F.Molteni paper) ; update aim.5l_cs output.
533    o New variable in PARM03 for more flexible pickups
534    o somewhat cleaned package initialisation sequence for
535      ctrl/ cost/ ecco/ in packages_... routines
536    
537    checkpoint50g_post
538    o merged with release1_p17 (pkg/seaice and verification/lab_sea)
539      - added SEAICE_MULTILEVEL for 8-category sea-ice thermodynamics
540      - LSR sea-ice dynamic solver moved to SouthWest B-grid location and
541        made the default because of faster convergence than ADI
542    
543    checkpoint50f_post
544    o exf totally restructured
545      - replaced all exf_set_ routines by one generic exf_set_gen.F
546        and made sure it's still differentiable
547      - split exf_getffields.F into new exf_bulkformulae.F
548      - exf_inscal_/exf_outscal_ suffixes consistent with field names
549      - updated adjoint/makefile and relevant verif. accordingly
550    
551    checkpoint50f_pre
552     o AIM: near surface temp evaluated from Pot.Temp: fixed for Part Cell ;
553          + change Evap over land. (update aim.5l_cs output)
554     o added missing .cvsignores and made uniform across expts
555     o changed configuration of tidal_basin_2d expt
556     o added pkg/debug functions for tracing code in fatal crashes
557     o tidied pkg/monitor to use MONITOR_OPTIONS.h
558     o new fn in pkg/monitor to stop model if solution grows too big
559     o removed lots of single quotes from comments
560    
561    checkpoint50e_post
562    o make KPP work with PTRACERS
563     - fix gad_calc_rhs to call new routine kpp_transport_ptr, which is
564       nearly a copy of kpp_transport_s
565     - there is no analogue to SurfaceTendencyS, so I have to use
566       gPtr(of the surface layer) instead
567    o add a new platform SunFire+mpi (SunFire 15000) to genmake
568    checkpoint50e_pre
569    
570    checkpoint50d_post
571    o change kpp output from multiple-record state files to single-record state
572      files analogous to write_state.F
573    o reduce the output frequency of cg3d-related stuff to the monitor frequency,
574      analogous to the cg2d-related output.
575    o fix small problem with in ptracers_write_checkpoint.F: len(suff)=512,
576      so that writing to internal file fn (with length 512) fails.
577    checkpoint50d_pre
578    
579    checkpoint50c_post
580    o fix problem in ini_procs.F when using MPI (was there since checkpoint48f)
581    o add monitor diagnostic of vorticity (mon_vort3.F) and potential energy.
582    o add time-average diagnostic of U*V and vertical diffusive flux for theta.
583    o bug fixed in convective adjustment diagnostic (multi-tile set-up).
584    o 4th Order Adv scheme: return to previous (c40pre7) masking near boundaries.
585    
586    checkpoint50c_pre
587    Merging from release1_p13:
588    o bug fix for pkg/seaice dynamic solver
589    o Added SEAICE_initialHEFF to pkg/seaice
590    
591    checkpoint50b_post
592    o store u*,v* in gU,V instead of in guNm1,gvNm1
593    o in order to put the momForcing out of the Adams-Bashforth:
594      move forcing & CD-scheme calls from mom_fluxform & mom_vecinv
595      to timestep.F
596    o new flag "useCDscheme" (default=F); replace guCD,gvCD by local arrays
597    o re-generate output for highly sensitive test-exp:
598       global_ocean_pressure, ideal_2D_oce, lab_sea, natl_box
599    
600    checkpoint50b_pre
601    o restore the use of momAdvection & useCoriolis flags (mom_vecinv)
602    o corrected CPP option mismatch in exf_Bulk... functions for case
603      #define ALLOW_ATM_WIND && #undef ALLOW_ATM_TEMP
604    
605    checkpoint50a_post
606    o Bug fix for merging between c50 and KPP.
607      ikey was passed from thermodynamics to kpp_calc via
608      common block rather than being recomputed in kpp_calc,
609      in contradiction with new key itdkey.
610      New key ikppkey created, and tamc.h headers updated.
611    
612    checkpoint50
613    - 2nd step of merge with ecco-branch:
614      merge changes in ecco-branch back into main trunk
615    - bug fix in find_alpha
616    - more storing in gmredi_calc_tensor.F, kpp_routines.F
617    - avoid iterative loop in ini_pressure for ALLOW_AUTODIFF_TAMC
618    - merged some sliced AD support routines
619    
620    checkpoint49
621    merged from ecco-branch:
622    o exf:
623      - Enable initialisation of forcing fields to constant
624        (runtime) values.
625      - in exf_getffields.F
626        Reduce i-/j-loop to interior domain, discarding overlaps.
627        That also fixes wrong TAF-key computations for key_1, key_2
628        with bulf formulae.
629      - exf_init.F modify #ifdef for exf_init_evap
630      - exf_getffieldrec.F, ctrl_getrec.F
631        The following INT-usages are not safe:
632          fldsecs  = int(fldsecs/fldperiod)*fldperiod
633          fldcount = int(fldsecs/fldperiod) + 1
634        and were modified.
635    o autodiff: checkpoint_lev...
636      - updated storing for bulk formulae (some arrays were missing)
637      - added state for seaice package
638      - synchronized with main branch
639    o DIVA:
640      Modif.s in
641      - the_main_loop.F
642      - the_model_main.F
643      and added cost_final_restore.F cost_final_store.F
644      o With these changes, it is possible to fully separate
645        the first forward run to compute the cost function (mdthe_main_loop)
646        and adjoint runs over individual DIVA intervals (adthe_main_loop)
647        into separate model executions.
648        This enables better tuning of DIVA to given queue sizes.
649      o Loss of cost function value in consecutive DIVA runs is
650        fixed by storing to 'costfinal' file.
651      o Appropriate call of ctrl_unpack/_pack,grdchk_main at
652        beginning/end of adjoint integration is accounted for.
653      The current changes apply to TAF-1.4.26.
654      TAF has one remaining bug, requiring one manual change
655      (see MITgcm manual).
656    o cal: synchronized between MAIN and ecco-branch
657    
658    checkpoint48i_post
659    o Changes to restore differentiability of code w.r.t. previous tag
660      (mostly adding new routines to make list and replacing
661      pressure by totPhiHyd).
662    o Updated and merged grdchk package
663      - has both ADM and TLM checks
664      - works for single- and multi-proc.
665      - output cleaned
666      - worked successfully for parallel DIVA
667    o genmake:
668      - Adding -platform=SP4
669      - Modif. for DEC Alpha (many don't have makedepend; use mkdep instead).
670    o Modif. for DIVA:
671      ctrl_unpack should not be called in intermediate DIVA interval.
672      (similar thing would apply for active files such as
673      tbar, sbar, psbar, which should not be reinitialised;
674      only relevant for ecco-branch)
675    
676    checkpoint48h_post
677    o improve "ini_p_ground" when using finite difference form for calc_phi_hyd
678    o Added net flux to downward flux conversion to pkg/exf/exf_getffields.F
679    o Added SEAICE_initialHEFF to pkg/seaice
680    o ecco_the_model_main.F: SEAICE_WRITE_PICKUP must precede WRITE_CHECKPOINT
681    
682    checkpoint48g_post
683    --- EOS:
684    o compute locally the pressure for use in EOS: UNESCO, JMD95P, MDJWF
685      store total Potential in totPhyHyd for diagnostics & EOS.
686    o fix restart and overlap Pb when using Z-coord and EOS funct. of P
687    --- Pressure/geopotential gradient:
688    o use Ro_surf & R_low instead of hFac in cal_phi_hyd:
689    o NLFS + p-coord.: eta*Alpha' is not longer lagging 1.time-step behind
690      ( changes results of global_ocean_pressure )
691    o NLFS + z-coord.: add eta*Rho' contribution.
692    o accurate phiHyd when using interface-W at the middle between 2 centers
693     ( changes results of ideal_2D_oce )
694    o includes r* 2nd term in potential gradient.
695    
696    checkpoint48f_post
697    Merging from release1_p12:
698    o Modifications for using pkg/exf with pkg/seaice
699      - improved description of the various forcing configurations
700      - added basic radiation bulk formulae to pkg/exf
701      - units/sign fix for evap computation in exf_getffields.F
702      - updated verification/global_with_exf/results/output.txt
703    o Added pkg/sbo for computing IERS Special Bureau for the Oceans
704      (SBO) core products, including oceanic mass, center-of-mass,
705      angular, and bottom pressure (see pkg/sbo/README.sbo).
706    o Lower bound for viscosity/diffusivity in pkg/kpp/kpp_routines.F
707      to avoid negative values in shallow regions.
708      - updated verification/natl_box/results/output.txt
709      - updated verification/lab_sea/results/output.txt
710    o MPI gather, scatter: eesupp/src/gather_2d.F and scatter_2d.F
711    o Added useSingleCpuIO option (see PARAMS.h).
712    o Updated useSingleCpuIO option in mdsio_writefield.F to
713      work with multi-field files, e.g., for single-file pickup.
714    o pkg/seaice:
715      - bug fix in growth.F: QNET for no shortwave case
716      - added HeffFile for specifying initial sea-ice thickness
717      - changed SEAICE_EXTERNAL_FLUXES wind stress implementation
718    o Added missing /* */ to CPP comments in pkg/seaice, pkg/exf,
719      kpp_transport_t.F, forward_step.F, and the_main_loop.F
720    o pkg/seaice:
721      - adjoint-friendly modifications
722      - added a SEAICE_WRITE_PICKUP at end of the_model_main.F
723    
724    checkpoint48e_post
725    o New S/R for diagnostic of bottom pressure (phi0surf was missing in c48d_post)
726    o dynamics: change definition of computational domain & adapt mom_fluxform
727      accordingly ; when viscA4=0, allows to run the dynamics with Olx=Oly=2.
728    o following exp works with Olx=Oly=2 (and pass through checkbounds option):
729        adjust_nlfs.cs-32x32x1 adjustment.128x64x1 adjustment.cs-32x32x1
730        exp0 exp1 exp2 exp4 exp5 hs94.cs-32x32x5 hs94.128x64x5 hs94.1x64x5
731        global_ocean.90x40x15 global_with_exf global_ocean_pressure
732        inverted_barometer solid-body.cs-32x32x1
733      (high order AdvScheme and/or multiDimAd works with Olx=Oly=2 ? not sure)
734    
735    checkpoint48d_post
736    o in preparation for r*:
737     a) use pre-computed gradient of hydrostatic potential:
738        changes in timestep.F & mom_cdscheme.F affects results of ideal_2D_oce
739     b) move phi0surf from calc_phi_hyd to calc_grad_phi_hyd:
740        => affects results of glob_oce_pressure (different truncation error)  
741      update output.txt of ideal_2D_oce & global_ocean_pressure
742    
743    checkpoint48d_pre
744    o in preparation for r*:
745      new S/R (calc_grad_phi_hyd.F) to compute Hydrostatic potential gradient.
746      pass the 2 comp. of the grad. as arguments to momentum S/R.
747      but for the moment, only used if it does not change the results.
748    o add finite volume form of calc_phi_hyd for buoyancyRelation OCEANIC & OCEANICP
749    
750  checkpoint48c_post  checkpoint48c_post
751  Modifications in treatment of initial flux forcing controls  Modifications in treatment of initial flux forcing controls
752  to avoid extra storing of flux fields.  to avoid extra storing of flux fields.

Legend:
Removed from v.1.128  
changed lines
  Added in v.1.288

  ViewVC Help
Powered by ViewVC 1.1.22