/[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.75 by jmc, Thu Jun 20 23:52:03 2002 UTC revision 1.184 by edhill, Thu Oct 9 04:19:17 2003 UTC
# Line 2  $Header$ Line 2  $Header$
2    
3      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
4      ==============================      ==============================
5    
6    checkpoint51i_post
7    o merge of the branch-genmake2
8      - knowingly breaks the older genmake
9      - automatic generation of all ALLOW_${PKG_NAME} defines
10      - new testreport supersedes older testscript
11    o large number of package inclusion and header cleanups
12    checkpoint51i_pre
13    
14    o new output for adjoin of global ocean vector invariant
15    o fix for checkpoint of thermodynamic seaice and bulk_force
16    o Steph gchem various cleanups
17    
18    checkpoint 51h_pre
19    o fixed some single quotes
20    o modify gchem pkg to be more versatile
21    o modifications to clean up dic pkg
22    o add pkg/cfc
23    
24    checkpoint51g_post
25    o provide links for ARPACK and ADM*TLM
26    o start removing verification/vero/
27    o update AD-related stuff
28      (in particular some IF-statements in model/src/ routines)
29    o adding to verification/global_ocean.90x40x15/ adjoint of
30      vector invariant setup
31    
32    checkpoint51f_post
33    o update advect_xz/output.txt (left from chekpoint50f_pre, ini_vel modif)
34    o Mods and bug fixes to pkg/cal, pkg/exf, etc., needed for computation
35      of tracer Green's fucntions for ocean inversion project.
36    
37    checkpoint51f_pre
38    o Added on-the-fly spatial interpolation capability
39        "USE_EXF_INTERPOLATION" to pkg/exf.
40      - This is a temporary Cartesian-grid hack until
41        the super-duper ESMF coupler becomes available.
42      - See verification/global_with_exf/README for usage example.
43      - Removed obsolete EXFwindOnBgrid and SEAICEwindOnCgrid
44        flags and modified pkg/seaice accordingly.
45    o Bug fix to pkg/ptracers, pkg/generic_advdiff/gad_calc_rhs.F,
46        and pkg/kpp/kpp_transport_ptr.F for dealing with tracer
47        non-local transport term.
48    
49    checkpoint51e_post
50    o pkg/mom_vecinv:
51      - changes related to hFac: use S/R from aim.5l_cs/code (coriolis & vort.advec)
52      - add time-ave diagnostic of vorticity advection.
53    o add diagnostic (instantaneous) of Shapiro Filter effect for T,S & UV.
54    o atmospheric geopotential: - include water vapor effect;
55                                - use p^kappa in p* scaling.
56    o land pkg: read initial state from files.
57    o AIM : add stratospheric drag in the upper level
58    
59    checkpoint51d_post
60    o change function ILNBLNK (=> last non-blank char.) to return zero
61      (instead of the string length) when the string. char is empty
62    o fix bug in dxV,dyU definition of cubed-sphere grid (ini_curvilinear_grid)
63    o modif's to restore TAMC compatibility
64      - avoid "_d" constructs in PARAMETER declarations
65      - avoid interrupting lines by comment lines
66    o Use cluster local disks for purely local I/O
67      vs. globally visible disks needed for ctrl stuff
68      (merged from ecco-branch)
69      - new active_file_loc... and mdsio_..._loc routines
70      - affected packages: autodiff, ctrl, grdchk, mdsio
71    
72    checkpoint51c_post
73    o introducing integer flag debugLevel
74      - debLevZero: suppress all output
75      - debLevA: standard
76      - debLevB: equivalent to current debugMode=.TRUE.
77    o introducing pathname variable mdsioLocalDir for mdsio
78      that can be specified at runtime
79      (relevant for cluster I/O to local disk)
80    `
81    >>>checkpoint51b_post
82    adjoint of ptracers
83    o disentangled ALLOW_PTRACERS using new ALLOW_GCHEM
84    o extended state in checkpoint_lev?, etc.
85    o alternative ctrl/cost with tr1 or ptracer
86    o new comlev1_bibj_ptracers
87    o new ptracers.flow
88    o made ptracers_forcing same as external_forcing_s
89    o added verification for adjoint with ptracers
90      - adm. grdchk fail so far
91      - tlm. grdchk OK
92      - checked tapelev/comlev recomputations: OK
93    <<<checkpoint51b_pre
94    
95    checkpoint51a_post
96    o added code for biogeochemistry package
97      - new pkgs: dic, gchem
98      - new subroutine in ptracers
99    o modify code for biogeo packages
100      - in model/src and pkg/ptracers
101    
102    checkpoint51
103    o merged improved packages from ecco-branch (ecco_c51_e34) to MAIN
104      - autodiff, ctrl, exf, grdchk, seaice
105      kept separate packages
106      - cost, ecco, obcs
107    o updating verification/internal_wave, verification/lab_sea
108      (all verif's pass on faulks.lcs.mit.edu RH7.3, g77)
109    
110    checkpoint50i_post
111    o Preparing next differentiable checkpoint and sync
112      of MAIN vs. ecco-branch
113      - 1st step: this one is differentiable
114      - 2nd step: next(! not this) checkpoint will have updated packages
115                  from ecco-branch)
116    
117    checkpoint50h_post
118    o new package "land": provide surface BC for AIM.
119    o pkg/ptracers: apply exchange to all tracers after reading a pickup.
120    o AIM: use true log-P extrapolation for near surface temp (= as in
121        F.Molteni paper) ; update aim.5l_cs output.
122    o New variable in PARM03 for more flexible pickups
123    o somewhat cleaned package initialisation sequence for
124      ctrl/ cost/ ecco/ in packages_... routines
125    
126    checkpoint50g_post
127    o merged with release1_p17 (pkg/seaice and verification/lab_sea)
128      - added SEAICE_MULTILEVEL for 8-category sea-ice thermodynamics
129      - LSR sea-ice dynamic solver moved to SouthWest B-grid location and
130        made the default because of faster convergence than ADI
131    
132    checkpoint50f_post
133    o exf totally restructured
134      - replaced all exf_set_ routines by one generic exf_set_gen.F
135        and made sure it's still differentiable
136      - split exf_getffields.F into new exf_bulkformulae.F
137      - exf_inscal_/exf_outscal_ suffixes consistent with field names
138      - updated adjoint/makefile and relevant verif. accordingly
139    
140    checkpoint50f_pre
141     o AIM: near surface temp evaluated from Pot.Temp: fixed for Part Cell ;
142          + change Evap over land. (update aim.5l_cs output)
143     o added missing .cvsignores and made uniform across expts
144     o changed configuration of tidal_basin_2d expt
145     o added pkg/debug functions for tracing code in fatal crashes
146     o tidied pkg/monitor to use MONITOR_OPTIONS.h
147     o new fn in pkg/monitor to stop model if solution grows too big
148     o removed lots of single quotes from comments
149    
150    checkpoint50e_post
151    o make KPP work with PTRACERS
152     - fix gad_calc_rhs to call new routine kpp_transport_ptr, which is
153       nearly a copy of kpp_transport_s
154     - there is no analogue to SurfaceTendencyS, so I have to use
155       gPtr(of the surface layer) instead
156    o add a new platform SunFire+mpi (SunFire 15000) to genmake
157    checkpoint50e_pre
158    
159    checkpoint50d_post
160    o change kpp output from multiple-record state files to single-record state
161      files analogous to write_state.F
162    o reduce the output frequency of cg3d-related stuff to the monitor frequency,
163      analogous to the cg2d-related output.
164    o fix small problem with in ptracers_write_checkpoint.F: len(suff)=512,
165      so that writing to internal file fn (with length 512) fails.
166    checkpoint50d_pre
167    
168    checkpoint50c_post
169    o fix problem in ini_procs.F when using MPI (was there since checkpoint48f)
170    o add monitor diagnostic of vorticity (mon_vort3.F) and potential energy.
171    o add time-average diagnostic of U*V and vertical diffusive flux for theta.
172    o bug fixed in convective adjustment diagnostic (multi-tile set-up).
173    o 4th Order Adv scheme: return to previous (c40pre7) masking near boundaries.
174    
175    checkpoint50c_pre
176    Merging from release1_p13:
177    o bug fix for pkg/seaice dynamic solver
178    o Added SEAICE_initialHEFF to pkg/seaice
179    
180    checkpoint50b_post
181    o store u*,v* in gU,V instead of in guNm1,gvNm1
182    o in order to put the momForcing out of the Adams-Bashforth:
183      move forcing & CD-scheme calls from mom_fluxform & mom_vecinv
184      to timestep.F
185    o new flag "useCDscheme" (default=F); replace guCD,gvCD by local arrays
186    o re-generate output for highly sensitive test-exp:
187       global_ocean_pressure, ideal_2D_oce, lab_sea, natl_box
188    
189    checkpoint50b_pre
190    o restore the use of momAdvection & useCoriolis flags (mom_vecinv)
191    o corrected CPP option mismatch in exf_Bulk... functions for case
192      #define ALLOW_ATM_WIND && #undef ALLOW_ATM_TEMP
193    
194    checkpoint50a_post
195    o Bug fix for merging between c50 and KPP.
196      ikey was passed from thermodynamics to kpp_calc via
197      common block rather than being recomputed in kpp_calc,
198      in contradiction with new key itdkey.
199      New key ikppkey created, and tamc.h headers updated.
200    
201    checkpoint50
202    - 2nd step of merge with ecco-branch:
203      merge changes in ecco-branch back into main trunk
204    - bug fix in find_alpha
205    - more storing in gmredi_calc_tensor.F, kpp_routines.F
206    - avoid iterative loop in ini_pressure for ALLOW_AUTODIFF_TAMC
207    - merged some sliced AD support routines
208    
209    checkpoint49
210    merged from ecco-branch:
211    o exf:
212      - Enable initialisation of forcing fields to constant
213        (runtime) values.
214      - in exf_getffields.F
215        Reduce i-/j-loop to interior domain, discarding overlaps.
216        That also fixes wrong TAF-key computations for key_1, key_2
217        with bulf formulae.
218      - exf_init.F modify #ifdef for exf_init_evap
219      - exf_getffieldrec.F, ctrl_getrec.F
220        The following INT-usages are not safe:
221          fldsecs  = int(fldsecs/fldperiod)*fldperiod
222          fldcount = int(fldsecs/fldperiod) + 1
223        and were modified.
224    o autodiff: checkpoint_lev...
225      - updated storing for bulk formulae (some arrays were missing)
226      - added state for seaice package
227      - synchronized with main branch
228    o DIVA:
229      Modif.s in
230      - the_main_loop.F
231      - the_model_main.F
232      and added cost_final_restore.F cost_final_store.F
233      o With these changes, it is possible to fully separate
234        the first forward run to compute the cost function (mdthe_main_loop)
235        and adjoint runs over individual DIVA intervals (adthe_main_loop)
236        into separate model executions.
237        This enables better tuning of DIVA to given queue sizes.
238      o Loss of cost function value in consecutive DIVA runs is
239        fixed by storing to 'costfinal' file.
240      o Appropriate call of ctrl_unpack/_pack,grdchk_main at
241        beginning/end of adjoint integration is accounted for.
242      The current changes apply to TAF-1.4.26.
243      TAF has one remaining bug, requiring one manual change
244      (see MITgcm manual).
245    o cal: synchronized between MAIN and ecco-branch
246    
247    checkpoint48i_post
248    o Changes to restore differentiability of code w.r.t. previous tag
249      (mostly adding new routines to make list and replacing
250      pressure by totPhiHyd).
251    o Updated and merged grdchk package
252      - has both ADM and TLM checks
253      - works for single- and multi-proc.
254      - output cleaned
255      - worked successfully for parallel DIVA
256    o genmake:
257      - Adding -platform=SP4
258      - Modif. for DEC Alpha (many don't have makedepend; use mkdep instead).
259    o Modif. for DIVA:
260      ctrl_unpack should not be called in intermediate DIVA interval.
261      (similar thing would apply for active files such as
262      tbar, sbar, psbar, which should not be reinitialised;
263      only relevant for ecco-branch)
264    
265    checkpoint48h_post
266    o improve "ini_p_ground" when using finite difference form for calc_phi_hyd
267    o Added net flux to downward flux conversion to pkg/exf/exf_getffields.F
268    o Added SEAICE_initialHEFF to pkg/seaice
269    o ecco_the_model_main.F: SEAICE_WRITE_PICKUP must precede WRITE_CHECKPOINT
270    
271    checkpoint48g_post
272    --- EOS:
273    o compute locally the pressure for use in EOS: UNESCO, JMD95P, MDJWF
274      store total Potential in totPhyHyd for diagnostics & EOS.
275    o fix restart and overlap Pb when using Z-coord and EOS funct. of P
276    --- Pressure/geopotential gradient:
277    o use Ro_surf & R_low instead of hFac in cal_phi_hyd:
278    o NLFS + p-coord.: eta*Alpha' is not longer lagging 1.time-step behind
279      ( changes results of global_ocean_pressure )
280    o NLFS + z-coord.: add eta*Rho' contribution.
281    o accurate phiHyd when using interface-W at the middle between 2 centers
282     ( changes results of ideal_2D_oce )
283    o includes r* 2nd term in potential gradient.
284    
285    checkpoint48f_post
286    Merging from release1_p12:
287    o Modifications for using pkg/exf with pkg/seaice
288      - improved description of the various forcing configurations
289      - added basic radiation bulk formulae to pkg/exf
290      - units/sign fix for evap computation in exf_getffields.F
291      - updated verification/global_with_exf/results/output.txt
292    o Added pkg/sbo for computing IERS Special Bureau for the Oceans
293      (SBO) core products, including oceanic mass, center-of-mass,
294      angular, and bottom pressure (see pkg/sbo/README.sbo).
295    o Lower bound for viscosity/diffusivity in pkg/kpp/kpp_routines.F
296      to avoid negative values in shallow regions.
297      - updated verification/natl_box/results/output.txt
298      - updated verification/lab_sea/results/output.txt
299    o MPI gather, scatter: eesupp/src/gather_2d.F and scatter_2d.F
300    o Added useSingleCpuIO option (see PARAMS.h).
301    o Updated useSingleCpuIO option in mdsio_writefield.F to
302      work with multi-field files, e.g., for single-file pickup.
303    o pkg/seaice:
304      - bug fix in growth.F: QNET for no shortwave case
305      - added HeffFile for specifying initial sea-ice thickness
306      - changed SEAICE_EXTERNAL_FLUXES wind stress implementation
307    o Added missing /* */ to CPP comments in pkg/seaice, pkg/exf,
308      kpp_transport_t.F, forward_step.F, and the_main_loop.F
309    o pkg/seaice:
310      - adjoint-friendly modifications
311      - added a SEAICE_WRITE_PICKUP at end of the_model_main.F
312    
313    checkpoint48e_post
314    o New S/R for diagnostic of bottom pressure (phi0surf was missing in c48d_post)
315    o dynamics: change definition of computational domain & adapt mom_fluxform
316      accordingly ; when viscA4=0, allows to run the dynamics with Olx=Oly=2.
317    o following exp works with Olx=Oly=2 (and pass through checkbounds option):
318        adjust_nlfs.cs-32x32x1 adjustment.128x64x1 adjustment.cs-32x32x1
319        exp0 exp1 exp2 exp4 exp5 hs94.cs-32x32x5 hs94.128x64x5 hs94.1x64x5
320        global_ocean.90x40x15 global_with_exf global_ocean_pressure
321        inverted_barometer solid-body.cs-32x32x1
322      (high order AdvScheme and/or multiDimAd works with Olx=Oly=2 ? not sure)
323    
324    checkpoint48d_post
325    o in preparation for r*:
326     a) use pre-computed gradient of hydrostatic potential:
327        changes in timestep.F & mom_cdscheme.F affects results of ideal_2D_oce
328     b) move phi0surf from calc_phi_hyd to calc_grad_phi_hyd:
329        => affects results of glob_oce_pressure (different truncation error)  
330      update output.txt of ideal_2D_oce & global_ocean_pressure
331    
332    checkpoint48d_pre
333    o in preparation for r*:
334      new S/R (calc_grad_phi_hyd.F) to compute Hydrostatic potential gradient.
335      pass the 2 comp. of the grad. as arguments to momentum S/R.
336      but for the moment, only used if it does not change the results.
337    o add finite volume form of calc_phi_hyd for buoyancyRelation OCEANIC & OCEANICP
338    
339    checkpoint48c_post
340    Modifications in treatment of initial flux forcing controls
341    to avoid extra storing of flux fields.
342    
343    checkpoint48c_pre
344    o Added tidal_basin_2d experiment (currently 1D!)
345    o can use r* Coordinate with OBC (but Orlanski Rad. BC not yet implemented)
346    
347    checkpoint48b_post
348    o r* coordinate added in #ifdef NONLIN_FRSURF block.
349     (modification to pressure gradient not yet implemented)
350    o mom_fluxform: new S/R mom_calc_rtrans.F computes vertical transports
351      (for U & V) that are used for vertical advection of momentum.
352    
353    checkpoint48a_post
354    o In preparation for r*, use temp. array to compute vertical transport
355       in flux-form momentum pkg. => change results of few test-exp.
356      Update output of: exp4, exp5, front_relax, global_ocean_pressure,
357                        lab_sea, natl_box & plume_on_slope
358    
359    checkpoint48
360    o fix small bug that was causing lot of warning (TIMER_STOP 'UPDATE_CG2D')
361    
362    checkpoint47j_post
363    o add new exp. (ideal_2D_oce) to test "exotic" parameters and options.
364    o moved adexch from forward_step to addummy_in_stepping
365    o added exch in ctrl_map_...
366    o aveFreq.=0. at end of cost_final to avoid overwrite in
367      reverse checkpointing loop
368    o Yet more changes in gmredi:
369      * adgmredi_calc_tensor
370        avoiding all recomputation of gmredi_slope_limit
371      * adgmredi_x/y/rtransport
372        added flag for excessive storing to avoid recomp. of
373        u/v/rtans, dTdx/y/z
374        -> this is not really necessary and very memory-consuming
375      * adgmredi_slope_psi:
376        consistency with gmredi_slope_limit in treatment of GM_slopeSqCutoff
377      * gmredi_slope_limit
378        re-activated full calculation of taperfct for case 'ac02'
379    o updating verification/ adjoints of carbon, global_ocean.90x40x15
380    
381    checkpoint47i_post
382    o GM_EXCLUDE_() replace GM_TAPER_(): by default, all the code is compiled.
383    o add GM_Small_Number & GM_slopeSqCutoff as run-time parameters (data.gmredi)
384      default values = 1.D-12 & 1.D+48 (=> give same results as checkpoint47f_post)
385    o cg2dTargetResWunit: define tolerance using deltaTmom (and not deltaTfreesurf)
386    o gmredi: fix few bugs.
387    
388    checkpoint47h_post (duplicate tag for PH's c48 tag, in order to move c48)
389    This will be a good checkpoint to plug into the ecco-branch.
390    o few modif.'s to get latest version adjointed
391      (mainly kick out code in ini_linear_phisurf)
392    o modif's to run adjoint with exactConserv
393    o case GM_BOLUS_ADVEC should be cleaned
394      S/R gmredi_slope_psi should be cleaned
395    o verification/carbon now has exactConserv=.TRUE.
396    
397    checkpoint47g_post
398    o time-average AIM output: becomes independent from statvars time-average
399      (can be used with #undef ALLOW_TIMEAVE) ; add new variables (e.g., RH)
400    o After TAF bug fix:
401      Tested adjoint of gmredi package.
402      Tests were done with checkpoint47 rather than with latest
403      checkpoint47f_post since the latter is broken for the adjoint
404      verification/carbon/ contains both ADM and TLM gradient checks.
405    Here they are, before they get changed/lost/stolen.
406    Mostly modif.'s to fix numerical sensitivities.
407    Gradient checks OK for
408    - GM_taper_scheme:
409      * clipping
410      * ac02
411      * linear
412      * glw91
413      * dm95
414      * ldd97
415    - GMREDI_OPTIONS:
416      * GM_VISBECK_VARIABLE_K
417      * GM_NON_UNITY_DIAGONAL
418      * GM_EXTRA_DIAGONAL
419      * GM_BOLUS_ADVEC
420    in conjunction with data.gmredi parameters to be checked in
421    in a few minutes under verification/carbon/code/
422    
423    checkpoint47f_post
424    Merging from release1_p10:
425    o modifications for using pkg/exf with pkg/seaice
426      - pkg/seaice CPP options SEAICE_EXTERNAL_FORCING
427        and SEAICE_EXTERNAL_FLUXES
428      - pkg/exf CPP options EXF_READ_EVAP and
429        EXF_NO_BULK_COMPUTATIONS
430      - usage examples are Experiments 8 and 9 in
431        verification/lab_sea/README
432      - verification/lab_sea default experiment now uses
433        pkg/gmredi, pkg/kpp, pkg/seaice, and pkg/exf
434    
435    checkpoint47e_post
436    o created a tag prior to creating a branch
437    
438    checkpoint47d_post
439    o fixed bug in the oceanic pressure coordinates code: vertical viscosity
440      at the bottom boundary had an erroneous half slip boundary condition
441    o added bottom drag and no slip boundary condition capabilities to
442      oceanic pressure coordinates code
443    o added new verification experiment global_ocean_pressure, this experiment is
444      also described in the tutorial section of the manual
445    checkpoint47d_pre
446    
447    checkpoint47c_post
448    o Added a verification experiment for a global cubed-sphere ocean run
449      -  verification/global_ocean.cubed32x32x30
450      -  uses pkg/therm_seaice and pkg/bulk_forcing
451    
452    o allow to run AIM physics with SPEEDY input files (from Franco Molteni)
453    o allow a more accurate definition of Ro_Surf (selectFindRoSurf=1)
454      when using P-coordinate; only implemented for atmospheric config.
455    o OCEANICP & realFreshWater: include P-E direct effect on wVel ;
456      NOTES: requires option NONLIN_FRSURF to be "#define".
457    o update advect_xz/results/output.txt (left from checkpoint44g_post)
458    
459    checkpoint47b_post
460    Merging from release1_p9:
461    o pkg/seaice
462      - removed GOTO's and added taf directives
463      - double precision constants to reduce the g77 (Linux)
464        to F77 (SGI) differences reported in release1_p8
465    o tools/genmake
466      - added SGI options
467    o verification/testscript
468      - updated to that of checkpoint47a_post
469    o verification/global_ocean.90x40x15/input/eedata
470      - modified for SGI f77 compatibility
471    o verification/lab_sea
472      - added description of sea-ice model
473      - added missing matlab routines
474      - added test of thermodynamics parallelization
475    
476    checkpoint47a_post
477    o new pkg: aim_v23
478     = F. Molteni atmos.physics (SPEEDY, ver23) adapted to MITgcm.
479       for now, keep the same surface forcing as before.
480     - Part-Cell implemented into AIM; check that Heat & Water are conserved
481     - aim.5l_cs: use new aim pkg (run 10yr & get better results than before)
482     - aim.5l_LatLon & aim.5l_Equatorial_Channel are still using old aim pkg
483    o add diagnostic of surface correction term in monitor output
484    
485    o bulk_force and therm_seaice
486    - Two packages:  bulk_force (Bulk forcing)
487      and therm_seaice (thermodynamic_seaice) - adopted from LANL CICE.v2.0.2
488    - Earlier integration from Stephaine Dutkiewicz
489      and Patrick Heimbach.
490    - Two ifdef statements for compile time,
491       ALLOW_THERM_SEAICE and ALLOW_BULK_FORCE
492    
493    checkpoint47
494    differentiable version of checkpoint46n_post
495    o bug fix in quasihydrostaticterms
496    o * "clean" adjoint code (in terms of extensive recomputations)
497        can now be obtained for all GMREDI options (i.e. for
498        - GM_VISBECK_VARIABLE_K
499        - GM_NON_UNITY_DIAGONAL
500        - GM_EXTRA_DIAGONAL
501        - GM_BOLUS_ADVEC )
502      * However, wrong gradient check problem remains unsolved.
503      * New CPP options have been introduced for different
504        tapering schemes
505    o external_fields_load now part of differentiation list
506    o pressure needs multiple storing;
507      would be nice to have store_pressure at beginning or
508      end of forward_step, e.g. by having phiHyd global (5-dim.)
509      (NB: pressure is needed for certain cases in find_rho,
510      which is also invoked through convective_adjustment).
511    o recomputations in find_rho for cases
512     'JMD95'/'UNESCO' or 'MDJWF' are OK.
513    o #define ATMOSPHERIC_LOADING should be differentiable
514    o ini_forcing shifted to begining of initialise_varia
515    o Incorporating QNVS line search routines into MITgcm
516      (this is separate code, and therefore not under pkg)
517      lsopt/, optim/
518    o Updated verification/carbon/
519      - #define GM_VISBECK_VARIABLE_K
520      - GM_taper_scheme    = 'ldd97'    
521    
522    checkpoint46n_post
523    Merging from release1_p8:
524    o verification/natl_box:
525      updating new external_fields_load routine
526    o New package: pkg/seaice
527      Sea ice model by D. Menemenlis (JPL) and Jinlun Zhang (Seattle).
528      The sea-ice code is based on Hibler (1979-1980).
529      Two sea-ice dynamic solvers, ADI and LSR, are included.
530      In addition to computing prognostic sea-ice variables and diagnosing
531      the forcing/external data fields that drive the ocean model,
532      SEAICE_MODEL also sets theta to the freezing point under sea-ice.
533      The implied surface heat flux is then stored in variable
534      surfaceTendencyTice, which is needed by KPP package (kpp_calc.F and
535      kpp_transport_t.F) to diagnose surface buoyancy fluxes and for the
536      non-local transport term.  Because this call precedes model
537      thermodynamics, temperature under sea-ice may not be "exactly" at
538      the freezing point by the time theta is dumped or time-averaged.
539      N.B.: fairly large differences in accuracy occur across
540            different platforms/compilers; comparison between
541            g77 (Linux) and F77 (SGI Origin 2000/MIPSpro)
542            generated output gives:
543       T           S           U           V
544      C D M    c        m  s        m  s        m  s        m  s
545      n p a R  g  m  m  e  .  m  m  e  .  m  m  e  .  m  m  e  .
546      f n k u  2  i  a  a  d  i  a  a  d  i  a  a  d  i  a  a  d
547      g d e n  d  n  x  n  .  n  x  n  .  n  x  n  .  n  x  n  .
548    
549      Y Y Y Y  5  5  7  7  7  8 10  9  6  6  6  6  7  5  7  5  7 FAIL  lab_sea
550    o GAD:
551      - generated new common blocks to account for call of
552        same gad routines with differing traceridentities
553        (needed to modify tracerIdentity indices in GAD.h)
554      - generated separate common blocks for case useCubedSphereExchange
555        (Department of Futurology)
556      - parameter lists to gmredi_?transport: added tracerIdentity
557      - added new key indices to tamc.h
558    o external_fields_load:
559      - added this routine to TAF list
560      - needed to make some common blocks global and additional storing
561        along the same lines as exf package (checkpoint_lev?_directives.h)
562    o exf:
563      updated external forcing package
564      - bug fixes carried over from ecco-branch
565        (missing OBCS_OPTIONS.h in two routines)
566      - enable easy to use "no forcing".
567      - added exf I/O for atmospheric loading
568      - added exf I/O for runoff data
569      - transfered scaling between exf <-> MITgcm to exf namelist
570    o Adding new verification experiment global_with_exf.
571      This verification experiment is almost identical to the
572      global_ocean.90x40x15 experiment, except that it uses
573      the exf (and cal) package instead of the
574      external_fields_load routine (cf README file).
575    
576    checkpoint46m_post
577    o modified testscript to return -- when there are no lines of "%MON" output
578    o added new flag "quasiHydrostatic" for QH mode
579      and new routine quasihydrostaticterms()
580    o added 2*Omega*Cos(phi)*W to u equations (non-hydrostatic Coriolis term)
581      New routine: mom_u_coriolis_nh()
582    o fixed sign of non-hydrostatic metric terms in u and v equations.
583      Depend on correct sign of gravitySign.
584    o added work.pc* to .cvsignore in aim directories using ifc
585    o introduced flag for controllin non-hydrostatic metric terms
586    o removed gravitySign from PARM04
587    o switched to g77/ifc from pgf77 in aim experiments
588    
589    checkpoint46l_post
590    o made convective adjustment work with pressure coordinates:
591     - changed the direction of k-loop in convective_adjustment.F for the
592       case of pressure coordinates (OCEANICP,ATMOSPHERIC buoyancyRelation)
593     - adjusted the reference pressure k-index in convective_adjustment.F
594     - adjusted the convection condition in convect.F (in analogy to
595       calc_ivdc.F)
596     - convective_adjustment no longer computes anything on the halos
597     - removed the warnings about negative salinity from find_rho.F and
598       find_alpha.F; instead the new routine look_for_neg_salinity, called
599       at the beginning of find_rho, find_alpha, and find_beta, does a
600       check of the entire slice, if CPP-option
601       CHECK_SALINITY_FOR_NEGATIVE_VALUES is defined
602    checkpoint46l_pre
603    
604    checkpoint46k_post
605    o fixed the verification/global_ocean.90x40x15 experiment:
606     - new bathymetry (the world according to A., JMC, and M.)
607     - new initial fields and forcing fields (*.bin files)
608     - new POLY3.COEFFS (for the next release one should switch to a full
609       equation of state: JMD95P or MDJWF)
610     - fixed several errors and redundancies in the data file
611     - experiment uses looped cells
612     - added matlab directory with diagnostic scripts for plotting of output
613    
614    o S/R aim_initialise.F replace S/R aim_init from file aim_do_inphys.F:
615     - read AIM physics parameters from a file (data.aimphys)
616     - set defaults values = F.Molteni paper (Clim.Dyn., 2002)
617    
618    checkpoint46j_post
619    o split calc_exact_eta in 2 S/R : integr_continuity & update_etaH
620    o move wVel computation at the end of the time step, in S/R integr_continuity
621    o create specific S/R to exchange T,S before DYNAMICS (for stagger time step)
622    o update timeave pkg for wVel diagnostic.
623    checkpoint46j_pre
624    
625    checkpoint46i_post
626    o Clean up AIM package (and keep the results unchanged):
627      include CPP_OPTION and use IMPLICT NONE in all routines ;
628      declare all the variables _RL ; use _d 0 for all numerical constants.
629      use ifdef ALLOW_AIM everywhere. And now AIM can be used with g77 !
630    
631    checkpoint46h_post
632    o cleaned up the use of rhoNil and rhoConst.
633      - rhoNil should only appear in the LINEAR equation of state, everywhere
634        else rhoNil is replaced by rhoConst, e.g. find_rho computes rho-rhoConst
635        and the dynamical equations are all divided by rhoConst
636    o introduced new parameter rhoConstFresh, a reference density of fresh
637      water, to remove the fresh water flux's dependence on rhoNil. The default
638      value is 999.8 kg/m^3
639    o cleanup up external_forcing.F and external_forcing_surf.F
640      - can now be used by both OCEANIC and OCEANICP
641    checkpoint46h_pre
642    
643     Added code to convert surface volume flux (fresh water) into
644     a mass flux when using P coordinates in the ocean (OCEANICP).
645     Note: It assumes you have set rho0=rhoConst=density of fresh water.
646    
647    checkpoint46g_post
648    o Include a new diagnostic variable phiHydLow for the ocean model
649      - in z-coordinates, it is the bottom pressure anomaly
650      - in p-coordinates, it is the sea surface elevation
651      - in both cases, these variable have global drift, reflecting the mass
652        drift in z-coordinates and the volume drift in p-coordinates
653      - included time averaging for phiHydLow, be aware of the drift!
654    o depth-dependent computation of Bo_surf for pressure coordinates
655      in the ocean (buoyancyRelation='OCEANICP')
656      - requires a new routine (FIND_RHO_SCALAR) to compute density with only
657        Theta, Salinity, and Pressure in the parameter list. This routine is
658        presently contained in find_rho.F. This routine does not give the
659        correct density for 'POLY3', which would be a z-dependent reference
660        density.
661    o cleaned up find_rho
662      - removed obsolete 'eqn' from the parameter list.
663    o added two new verification experiments: gop and goz
664      (4x4 degree global ocean, 15 layers in pressure and height coordinates)
665    checkpoint46g_pre
666    
667    checkpoint46f_post
668    o Enable tangent linear (forward mode) gradient checks:
669      - extended active file handling to g_... files
670      - added TANGENT_SIMULATION to theSimulationMode
671      - extended grdchk package accordingly
672    o added tangent linear setup in verification/carbon/code_ftl
673    o added adjoint setup in verification/front_relax/code_ad
674    
675    checkpoint46e_post
676    o Added new equation of state -> MDJWF
677      - EOS of McDougall et al., 2002, JAOT, submitted
678      - caveat: the equation of state is only valid for a smaller (more
679        realistic?) range of values than JMD95P/Z and UNESCO
680      - added masks to the calculation of pressure in store_pressure
681      - added more check values for density in check_eos (ini_eos.F), some of
682        the old check values are out of the range of the MDJWF-eos, so don't
683        expect perfect matches for those
684    checkpoint46e_pre
685    
686    checkpoint46d_post
687    o fixed store_pressure to work with both buoyancy relation = 'OCEANIC' and
688      'OCEANICP', also initialised field pressure correctly in ini_eos in the
689      case of pressure coordinates. eosType='JMD95Z' in
690      combination with buoyancyRelation='OCEANICP' now causes an error.
691    o Changed p = pressure(i,j,k,bi,bj)  to  p = pressure(i,j,kRef,bi,bj)
692      in find_alpha/beta.
693    checkpoint46d_pre
694    
695    o Changed p = pressure(i,j,k,bi,bj)  to  p = pressure(i,j,kRef,bi,bj)
696      so that JMD95Z and JMD95P give approptiate static stability in find_rho.
697    
698    checkpoint46c_post
699    o Changes necessary for ocean in p-coordinates
700     -  Added new buoyancy relation = 'OCEANICP'
701     -  Added new parameters = gravitySign (this used to be contained inside
702        the factor dRdZ which I added when we first switched to R coordinates).
703     X GM/Redi is not compatible (yet)
704     X bottom drag and no-slip need to be debugged.
705    checkpoint46c_pre
706    
707    checkpoint46b_post
708    o Added new equation of state -> JMD95Z and JMD95P
709      - EOS of Jackett and McDougall, 1995, JAOT
710      - moved all EOS parameters into EOS.h
711      - new routines ini_eos.F, store_pressure.F
712    o Added UNESCO EOS, but not recommended because it requires
713      in-situ temperature (see JMD95)
714    o Modified formatting for knudsen2.f in utils/knudsen2 and added
715      unesco.f to be used with POLY3
716    checkpoint46b_pre
717    
718    checkpoint46a_post
719    o Added (atmospheric) pressure loading
720      - new field in FFIELDS.h, etc...
721      - new cpp flag ATMOSPHERIC_LOADING
722    o Changed hFacC to _hFacC in calc_phi_hyd.F
723    o Added SHORTWAVE_HEATING to some files for consistency
724    checkpoint46a_pre
725    
726    checkpoint46
727    further merges from release1_p5
728    o new ctrl package
729      - adopted from ECCO environment to enable optimization
730      - added Eliassen Palm fluxes to controls
731    o added Eliassen Palm flux controls to gradient check package
732    o cost package
733      - Compatible with new ctrl package
734      - added Eliassen Palm cost hooks
735      - modif's of existing cost functions
736    o Adjoint-related bug fixes in kpp:
737      - kpp_calc: sore of kpphbl avoids recomputation/call to S/R kppmix
738      - kpp_routines: store of Rib avoids partial recomputation bug of TAF.
739    o autodiff package
740      - flow directives:
741        + new for S/R do_the_model_io
742        + enhanced for S/R checkpoint
743        + new gmredi.flow
744      - Introduced CPP option ALLOW_AUTODIFF_MONITOR to
745        disable adcommon.h/g_common.h by default.
746      - fixed problem for adjoint of global_max
747    o modified verification/carbon/
748    o New setup for adjoint of global_ocean
749    o added hooks for Eliassen Palm fluxes (dynamics.F, ini_fields.F)
750    o cleaned TAF keys iikey, idkey (dynamics.F, thermodynamics.F)
751    o incorporated hooks for sponge layer code (external_forcing.F)
752    
753    checkpoint45d_post
754    o import 2 fixes from release1_p5 (CG2D.h & packages_readparms.F)
755    o use recip_dx*,recip_dy* instead of /dx*,/dy* in orlanski_E,W,N,S
756      (affects plume_on_slope "cg2d" output)
757    o add new parameter deltaTfreesurf for free-surface time stepping
758    o exchange gT,gS when using staggerTimeStep & GM & Oly<4
759  o NONLIN_FRSURF: rescale gNm1 to get a better conservation with A-B  o NONLIN_FRSURF: rescale gNm1 to get a better conservation with A-B
760  o change phiMin in exp1 (to agree with documentation)  o change phiMin in exp1 (to agree with documentation)
761  o add new flags:  o add new flags:

Legend:
Removed from v.1.75  
changed lines
  Added in v.1.184

  ViewVC Help
Powered by ViewVC 1.1.22