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

Legend:
Removed from v.1.104  
changed lines
  Added in v.1.251

  ViewVC Help
Powered by ViewVC 1.1.22