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

Legend:
Removed from v.1.116  
changed lines
  Added in v.1.274

  ViewVC Help
Powered by ViewVC 1.1.22