/[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.95 by adcroft, Tue Nov 5 19:14:26 2002 UTC revision 1.151 by jmc, Thu May 1 22:43:25 2003 UTC
# Line 3  $Header$ Line 3  $Header$
3      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
4      ==============================      ==============================
5    
6    o add time-average diagnostic of U*V and vertical diffusive flux for theta.
7    o bug fixed in convective adjustment diagnostic (multi-tile set-up).
8    o 4th Order Adv scheme: return to previous (c40pre7) masking near boundaries.
9    
10    checkpoint50c_pre
11    Merging from release1_p13:
12    o bug fix for pkg/seaice dynamic solver
13    o Added SEAICE_initialHEFF to pkg/seaice
14    
15    checkpoint50b_post
16    o store u*,v* in gU,V instead of in guNm1,gvNm1
17    o in order to put the momForcing out of the Adams-Bashforth:
18      move forcing & CD-scheme calls from mom_fluxform & mom_vecinv
19      to timestep.F
20    o new flag "useCDscheme" (default=F); replace guCD,gvCD by local arrays
21    o re-generate output for highly sensitive test-exp:
22       global_ocean_pressure, ideal_2D_oce, lab_sea, natl_box
23    
24    checkpoint50b_pre
25    o restore the use of momAdvection & useCoriolis flags (mom_vecinv)
26    o corrected CPP option mismatch in exf_Bulk... functions for case
27      #define ALLOW_ATM_WIND && #undef ALLOW_ATM_TEMP
28    
29    checkpoint50a_post
30    o Bug fix for merging between c50 and KPP.
31      ikey was passed from thermodynamics to kpp_calc via
32      common block rather than being recomputed in kpp_calc,
33      in contradiction with new key itdkey.
34      New key ikppkey created, and tamc.h headers updated.
35    
36    checkpoint50
37    - 2nd step of merge with ecco-branch:
38      merge changes in ecco-branch back into main trunk
39    - bug fix in find_alpha
40    - more storing in gmredi_calc_tensor.F, kpp_routines.F
41    - avoid iterative loop in ini_pressure for ALLOW_AUTODIFF_TAMC
42    - merged some sliced AD support routines
43    
44    checkpoint49
45    merged from ecco-branch:
46    o exf:
47      - Enable initialisation of forcing fields to constant
48        (runtime) values.
49      - in exf_getffields.F
50        Reduce i-/j-loop to interior domain, discarding overlaps.
51        That also fixes wrong TAF-key computations for key_1, key_2
52        with bulf formulae.
53      - exf_init.F modify #ifdef for exf_init_evap
54      - exf_getffieldrec.F, ctrl_getrec.F
55        The following INT-usages are not safe:
56          fldsecs  = int(fldsecs/fldperiod)*fldperiod
57          fldcount = int(fldsecs/fldperiod) + 1
58        and were modified.
59    o autodiff: checkpoint_lev...
60      - updated storing for bulk formulae (some arrays were missing)
61      - added state for seaice package
62      - synchronized with main branch
63    o DIVA:
64      Modif.s in
65      - the_main_loop.F
66      - the_model_main.F
67      and added cost_final_restore.F cost_final_store.F
68      o With these changes, it is possible to fully separate
69        the first forward run to compute the cost function (mdthe_main_loop)
70        and adjoint runs over individual DIVA intervals (adthe_main_loop)
71        into separate model executions.
72        This enables better tuning of DIVA to given queue sizes.
73      o Loss of cost function value in consecutive DIVA runs is
74        fixed by storing to 'costfinal' file.
75      o Appropriate call of ctrl_unpack/_pack,grdchk_main at
76        beginning/end of adjoint integration is accounted for.
77      The current changes apply to TAF-1.4.26.
78      TAF has one remaining bug, requiring one manual change
79      (see MITgcm manual).
80    o cal: synchronized between MAIN and ecco-branch
81    
82    checkpoint48i_post
83    o Changes to restore differentiability of code w.r.t. previous tag
84      (mostly adding new routines to make list and replacing
85      pressure by totPhiHyd).
86    o Updated and merged grdchk package
87      - has both ADM and TLM checks
88      - works for single- and multi-proc.
89      - output cleaned
90      - worked successfully for parallel DIVA
91    o genmake:
92      - Adding -platform=SP4
93      - Modif. for DEC Alpha (many don't have makedepend; use mkdep instead).
94    o Modif. for DIVA:
95      ctrl_unpack should not be called in intermediate DIVA interval.
96      (similar thing would apply for active files such as
97      tbar, sbar, psbar, which should not be reinitialised;
98      only relevant for ecco-branch)
99    
100    checkpoint48h_post
101    o improve "ini_p_ground" when using finite difference form for calc_phi_hyd
102    o Added net flux to downward flux conversion to pkg/exf/exf_getffields.F
103    o Added SEAICE_initialHEFF to pkg/seaice
104    o ecco_the_model_main.F: SEAICE_WRITE_PICKUP must precede WRITE_CHECKPOINT
105    
106    checkpoint48g_post
107    --- EOS:
108    o compute locally the pressure for use in EOS: UNESCO, JMD95P, MDJWF
109      store total Potential in totPhyHyd for diagnostics & EOS.
110    o fix restart and overlap Pb when using Z-coord and EOS funct. of P
111    --- Pressure/geopotential gradient:
112    o use Ro_surf & R_low instead of hFac in cal_phi_hyd:
113    o NLFS + p-coord.: eta*Alpha' is not longer lagging 1.time-step behind
114      ( changes results of global_ocean_pressure )
115    o NLFS + z-coord.: add eta*Rho' contribution.
116    o accurate phiHyd when using interface-W at the middle between 2 centers
117     ( changes results of ideal_2D_oce )
118    o includes r* 2nd term in potential gradient.
119    
120    checkpoint48f_post
121    Merging from release1_p12:
122    o Modifications for using pkg/exf with pkg/seaice
123      - improved description of the various forcing configurations
124      - added basic radiation bulk formulae to pkg/exf
125      - units/sign fix for evap computation in exf_getffields.F
126      - updated verification/global_with_exf/results/output.txt
127    o Added pkg/sbo for computing IERS Special Bureau for the Oceans
128      (SBO) core products, including oceanic mass, center-of-mass,
129      angular, and bottom pressure (see pkg/sbo/README.sbo).
130    o Lower bound for viscosity/diffusivity in pkg/kpp/kpp_routines.F
131      to avoid negative values in shallow regions.
132      - updated verification/natl_box/results/output.txt
133      - updated verification/lab_sea/results/output.txt
134    o MPI gather, scatter: eesupp/src/gather_2d.F and scatter_2d.F
135    o Added useSingleCpuIO option (see PARAMS.h).
136    o Updated useSingleCpuIO option in mdsio_writefield.F to
137      work with multi-field files, e.g., for single-file pickup.
138    o pkg/seaice:
139      - bug fix in growth.F: QNET for no shortwave case
140      - added HeffFile for specifying initial sea-ice thickness
141      - changed SEAICE_EXTERNAL_FLUXES wind stress implementation
142    o Added missing /* */ to CPP comments in pkg/seaice, pkg/exf,
143      kpp_transport_t.F, forward_step.F, and the_main_loop.F
144    o pkg/seaice:
145      - adjoint-friendly modifications
146      - added a SEAICE_WRITE_PICKUP at end of the_model_main.F
147    
148    checkpoint48e_post
149    o New S/R for diagnostic of bottom pressure (phi0surf was missing in c48d_post)
150    o dynamics: change definition of computational domain & adapt mom_fluxform
151      accordingly ; when viscA4=0, allows to run the dynamics with Olx=Oly=2.
152    o following exp works with Olx=Oly=2 (and pass through checkbounds option):
153        adjust_nlfs.cs-32x32x1 adjustment.128x64x1 adjustment.cs-32x32x1
154        exp0 exp1 exp2 exp4 exp5 hs94.cs-32x32x5 hs94.128x64x5 hs94.1x64x5
155        global_ocean.90x40x15 global_with_exf global_ocean_pressure
156        inverted_barometer solid-body.cs-32x32x1
157      (high order AdvScheme and/or multiDimAd works with Olx=Oly=2 ? not sure)
158    
159    checkpoint48d_post
160    o in preparation for r*:
161     a) use pre-computed gradient of hydrostatic potential:
162        changes in timestep.F & mom_cdscheme.F affects results of ideal_2D_oce
163     b) move phi0surf from calc_phi_hyd to calc_grad_phi_hyd:
164        => affects results of glob_oce_pressure (different truncation error)  
165      update output.txt of ideal_2D_oce & global_ocean_pressure
166    
167    checkpoint48d_pre
168    o in preparation for r*:
169      new S/R (calc_grad_phi_hyd.F) to compute Hydrostatic potential gradient.
170      pass the 2 comp. of the grad. as arguments to momentum S/R.
171      but for the moment, only used if it does not change the results.
172    o add finite volume form of calc_phi_hyd for buoyancyRelation OCEANIC & OCEANICP
173    
174    checkpoint48c_post
175    Modifications in treatment of initial flux forcing controls
176    to avoid extra storing of flux fields.
177    
178    checkpoint48c_pre
179    o Added tidal_basin_2d experiment (currently 1D!)
180    o can use r* Coordinate with OBC (but Orlanski Rad. BC not yet implemented)
181    
182    checkpoint48b_post
183    o r* coordinate added in #ifdef NONLIN_FRSURF block.
184     (modification to pressure gradient not yet implemented)
185    o mom_fluxform: new S/R mom_calc_rtrans.F computes vertical transports
186      (for U & V) that are used for vertical advection of momentum.
187    
188    checkpoint48a_post
189    o In preparation for r*, use temp. array to compute vertical transport
190       in flux-form momentum pkg. => change results of few test-exp.
191      Update output of: exp4, exp5, front_relax, global_ocean_pressure,
192                        lab_sea, natl_box & plume_on_slope
193    
194    checkpoint48
195    o fix small bug that was causing lot of warning (TIMER_STOP 'UPDATE_CG2D')
196    
197    checkpoint47j_post
198    o add new exp. (ideal_2D_oce) to test "exotic" parameters and options.
199    o moved adexch from forward_step to addummy_in_stepping
200    o added exch in ctrl_map_...
201    o aveFreq.=0. at end of cost_final to avoid overwrite in
202      reverse checkpointing loop
203    o Yet more changes in gmredi:
204      * adgmredi_calc_tensor
205        avoiding all recomputation of gmredi_slope_limit
206      * adgmredi_x/y/rtransport
207        added flag for excessive storing to avoid recomp. of
208        u/v/rtans, dTdx/y/z
209        -> this is not really necessary and very memory-consuming
210      * adgmredi_slope_psi:
211        consistency with gmredi_slope_limit in treatment of GM_slopeSqCutoff
212      * gmredi_slope_limit
213        re-activated full calculation of taperfct for case 'ac02'
214    o updating verification/ adjoints of carbon, global_ocean.90x40x15
215    
216    checkpoint47i_post
217    o GM_EXCLUDE_() replace GM_TAPER_(): by default, all the code is compiled.
218    o add GM_Small_Number & GM_slopeSqCutoff as run-time parameters (data.gmredi)
219      default values = 1.D-12 & 1.D+48 (=> give same results as checkpoint47f_post)
220    o cg2dTargetResWunit: define tolerance using deltaTmom (and not deltaTfreesurf)
221    o gmredi: fix few bugs.
222    
223    checkpoint47h_post (duplicate tag for PH's c48 tag, in order to move c48)
224    This will be a good checkpoint to plug into the ecco-branch.
225    o few modif.'s to get latest version adjointed
226      (mainly kick out code in ini_linear_phisurf)
227    o modif's to run adjoint with exactConserv
228    o case GM_BOLUS_ADVEC should be cleaned
229      S/R gmredi_slope_psi should be cleaned
230    o verification/carbon now has exactConserv=.TRUE.
231    
232    checkpoint47g_post
233    o time-average AIM output: becomes independent from statvars time-average
234      (can be used with #undef ALLOW_TIMEAVE) ; add new variables (e.g., RH)
235    o After TAF bug fix:
236      Tested adjoint of gmredi package.
237      Tests were done with checkpoint47 rather than with latest
238      checkpoint47f_post since the latter is broken for the adjoint
239      verification/carbon/ contains both ADM and TLM gradient checks.
240    Here they are, before they get changed/lost/stolen.
241    Mostly modif.'s to fix numerical sensitivities.
242    Gradient checks OK for
243    - GM_taper_scheme:
244      * clipping
245      * ac02
246      * linear
247      * glw91
248      * dm95
249      * ldd97
250    - GMREDI_OPTIONS:
251      * GM_VISBECK_VARIABLE_K
252      * GM_NON_UNITY_DIAGONAL
253      * GM_EXTRA_DIAGONAL
254      * GM_BOLUS_ADVEC
255    in conjunction with data.gmredi parameters to be checked in
256    in a few minutes under verification/carbon/code/
257    
258    checkpoint47f_post
259    Merging from release1_p10:
260    o modifications for using pkg/exf with pkg/seaice
261      - pkg/seaice CPP options SEAICE_EXTERNAL_FORCING
262        and SEAICE_EXTERNAL_FLUXES
263      - pkg/exf CPP options EXF_READ_EVAP and
264        EXF_NO_BULK_COMPUTATIONS
265      - usage examples are Experiments 8 and 9 in
266        verification/lab_sea/README
267      - verification/lab_sea default experiment now uses
268        pkg/gmredi, pkg/kpp, pkg/seaice, and pkg/exf
269    
270    checkpoint47e_post
271    o created a tag prior to creating a branch
272    
273    checkpoint47d_post
274    o fixed bug in the oceanic pressure coordinates code: vertical viscosity
275      at the bottom boundary had an erroneous half slip boundary condition
276    o added bottom drag and no slip boundary condition capabilities to
277      oceanic pressure coordinates code
278    o added new verification experiment global_ocean_pressure, this experiment is
279      also described in the tutorial section of the manual
280    checkpoint47d_pre
281    
282    checkpoint47c_post
283    o Added a verification experiment for a global cubed-sphere ocean run
284      -  verification/global_ocean.cubed32x32x30
285      -  uses pkg/therm_seaice and pkg/bulk_forcing
286    
287    o allow to run AIM physics with SPEEDY input files (from Franco Molteni)
288    o allow a more accurate definition of Ro_Surf (selectFindRoSurf=1)
289      when using P-coordinate; only implemented for atmospheric config.
290    o OCEANICP & realFreshWater: include P-E direct effect on wVel ;
291      NOTES: requires option NONLIN_FRSURF to be "#define".
292    o update advect_xz/results/output.txt (left from checkpoint44g_post)
293    
294    checkpoint47b_post
295    Merging from release1_p9:
296    o pkg/seaice
297      - removed GOTO's and added taf directives
298      - double precision constants to reduce the g77 (Linux)
299        to F77 (SGI) differences reported in release1_p8
300    o tools/genmake
301      - added SGI options
302    o verification/testscript
303      - updated to that of checkpoint47a_post
304    o verification/global_ocean.90x40x15/input/eedata
305      - modified for SGI f77 compatibility
306    o verification/lab_sea
307      - added description of sea-ice model
308      - added missing matlab routines
309      - added test of thermodynamics parallelization
310    
311    checkpoint47a_post
312    o new pkg: aim_v23
313     = F. Molteni atmos.physics (SPEEDY, ver23) adapted to MITgcm.
314       for now, keep the same surface forcing as before.
315     - Part-Cell implemented into AIM; check that Heat & Water are conserved
316     - aim.5l_cs: use new aim pkg (run 10yr & get better results than before)
317     - aim.5l_LatLon & aim.5l_Equatorial_Channel are still using old aim pkg
318    o add diagnostic of surface correction term in monitor output
319    
320    o bulk_force and therm_seaice
321    - Two packages:  bulk_force (Bulk forcing)
322      and therm_seaice (thermodynamic_seaice) - adopted from LANL CICE.v2.0.2
323    - Earlier integration from Stephaine Dutkiewicz
324      and Patrick Heimbach.
325    - Two ifdef statements for compile time,
326       ALLOW_THERM_SEAICE and ALLOW_BULK_FORCE
327    
328    checkpoint47
329    differentiable version of checkpoint46n_post
330    o bug fix in quasihydrostaticterms
331    o * "clean" adjoint code (in terms of extensive recomputations)
332        can now be obtained for all GMREDI options (i.e. for
333        - GM_VISBECK_VARIABLE_K
334        - GM_NON_UNITY_DIAGONAL
335        - GM_EXTRA_DIAGONAL
336        - GM_BOLUS_ADVEC )
337      * However, wrong gradient check problem remains unsolved.
338      * New CPP options have been introduced for different
339        tapering schemes
340    o external_fields_load now part of differentiation list
341    o pressure needs multiple storing;
342      would be nice to have store_pressure at beginning or
343      end of forward_step, e.g. by having phiHyd global (5-dim.)
344      (NB: pressure is needed for certain cases in find_rho,
345      which is also invoked through convective_adjustment).
346    o recomputations in find_rho for cases
347     'JMD95'/'UNESCO' or 'MDJWF' are OK.
348    o #define ATMOSPHERIC_LOADING should be differentiable
349    o ini_forcing shifted to begining of initialise_varia
350    o Incorporating QNVS line search routines into MITgcm
351      (this is separate code, and therefore not under pkg)
352      lsopt/, optim/
353    o Updated verification/carbon/
354      - #define GM_VISBECK_VARIABLE_K
355      - GM_taper_scheme    = 'ldd97'    
356    
357    checkpoint46n_post
358    Merging from release1_p8:
359    o verification/natl_box:
360      updating new external_fields_load routine
361    o New package: pkg/seaice
362      Sea ice model by D. Menemenlis (JPL) and Jinlun Zhang (Seattle).
363      The sea-ice code is based on Hibler (1979-1980).
364      Two sea-ice dynamic solvers, ADI and LSR, are included.
365      In addition to computing prognostic sea-ice variables and diagnosing
366      the forcing/external data fields that drive the ocean model,
367      SEAICE_MODEL also sets theta to the freezing point under sea-ice.
368      The implied surface heat flux is then stored in variable
369      surfaceTendencyTice, which is needed by KPP package (kpp_calc.F and
370      kpp_transport_t.F) to diagnose surface buoyancy fluxes and for the
371      non-local transport term.  Because this call precedes model
372      thermodynamics, temperature under sea-ice may not be "exactly" at
373      the freezing point by the time theta is dumped or time-averaged.
374      N.B.: fairly large differences in accuracy occur across
375            different platforms/compilers; comparison between
376            g77 (Linux) and F77 (SGI Origin 2000/MIPSpro)
377            generated output gives:
378       T           S           U           V
379      C D M    c        m  s        m  s        m  s        m  s
380      n p a R  g  m  m  e  .  m  m  e  .  m  m  e  .  m  m  e  .
381      f n k u  2  i  a  a  d  i  a  a  d  i  a  a  d  i  a  a  d
382      g d e n  d  n  x  n  .  n  x  n  .  n  x  n  .  n  x  n  .
383    
384      Y Y Y Y  5  5  7  7  7  8 10  9  6  6  6  6  7  5  7  5  7 FAIL  lab_sea
385    o GAD:
386      - generated new common blocks to account for call of
387        same gad routines with differing traceridentities
388        (needed to modify tracerIdentity indices in GAD.h)
389      - generated separate common blocks for case useCubedSphereExchange
390        (Department of Futurology)
391      - parameter lists to gmredi_?transport: added tracerIdentity
392      - added new key indices to tamc.h
393    o external_fields_load:
394      - added this routine to TAF list
395      - needed to make some common blocks global and additional storing
396        along the same lines as exf package (checkpoint_lev?_directives.h)
397    o exf:
398      updated external forcing package
399      - bug fixes carried over from ecco-branch
400        (missing OBCS_OPTIONS.h in two routines)
401      - enable easy to use "no forcing".
402      - added exf I/O for atmospheric loading
403      - added exf I/O for runoff data
404      - transfered scaling between exf <-> MITgcm to exf namelist
405    o Adding new verification experiment global_with_exf.
406      This verification experiment is almost identical to the
407      global_ocean.90x40x15 experiment, except that it uses
408      the exf (and cal) package instead of the
409      external_fields_load routine (cf README file).
410    
411    checkpoint46m_post
412    o modified testscript to return -- when there are no lines of "%MON" output
413    o added new flag "quasiHydrostatic" for QH mode
414      and new routine quasihydrostaticterms()
415    o added 2*Omega*Cos(phi)*W to u equations (non-hydrostatic Coriolis term)
416      New routine: mom_u_coriolis_nh()
417    o fixed sign of non-hydrostatic metric terms in u and v equations.
418      Depend on correct sign of gravitySign.
419    o added work.pc* to .cvsignore in aim directories using ifc
420    o introduced flag for controllin non-hydrostatic metric terms
421  o removed gravitySign from PARM04  o removed gravitySign from PARM04
422  o switched to g77/ifc from pgf77 in aim experiments  o switched to g77/ifc from pgf77 in aim experiments
423    

Legend:
Removed from v.1.95  
changed lines
  Added in v.1.151

  ViewVC Help
Powered by ViewVC 1.1.22