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

Legend:
Removed from v.1.123  
changed lines
  Added in v.1.283

  ViewVC Help
Powered by ViewVC 1.1.22