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

Legend:
Removed from v.1.107  
changed lines
  Added in v.1.246

  ViewVC Help
Powered by ViewVC 1.1.22