/[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.8 by adcroft, Wed Jul 1 22:49:29 1998 UTC revision 1.111 by mlosch, Mon Dec 16 21:51:50 2002 UTC
# Line 3  $Header$ Line 3  $Header$
3      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
4      ==============================      ==============================
5    
6    checkpoint47d_post
7    o fixed bug in the oceanic pressure coordinates code: vertical viscosity
8      at the bottom boundary had an erroneous half slip boundary condition
9    o added bottom drag and no slip boundary condition capabilities to
10      oceanic pressure coordinates code
11    o added new verification experiment global_ocean_pressure, this experiment is
12      also described in the tutorial section of the manual
13    checkpoint47d_pre
14    
15    checkpoint47c_post
16    o Added a verification experiment for a global cubed-sphere ocean run
17      -  verification/global_ocean.cubed32x32x30
18      -  uses pkg/therm_seaice and pkg/bulk_forcing
19    
20    
21    o allow to run AIM physics with SPEEDY input files (from Franco Molteni)
22    o allow a more accurate definition of Ro_Surf (selectFindRoSurf=1)
23      when using P-coordinate; only implemented for atmospheric config.
24    o OCEANICP & realFreshWater: include P-E direct effect on wVel ;
25      NOTES: requires option NONLIN_FRSURF to be "#define".
26    o update advect_xz/results/output.txt (left from checkpoint44g_post)
27    
28    checkpoint47b_post
29    Merging from release1_p9:
30    o pkg/seaice
31      - removed GOTO's and added taf directives
32      - double precision constants to reduce the g77 (Linux)
33        to F77 (SGI) differences reported in release1_p8
34    o tools/genmake
35      - added SGI options
36    o verification/testscript
37      - updated to that of checkpoint47a_post
38    o verification/global_ocean.90x40x15/input/eedata
39      - modified for SGI f77 compatibility
40    o verification/lab_sea
41      - added description of sea-ice model
42      - added missing matlab routines
43      - added test of thermodynamics parallelization
44    
45    checkpoint47a_post
46    o new pkg: aim_v23
47     = F. Molteni atmos.physics (SPEEDY, ver23) adapted to MITgcm.
48       for now, keep the same surface forcing as before.
49     - Part-Cell implemented into AIM; check that Heat & Water are conserved
50     - aim.5l_cs: use new aim pkg (run 10yr & get better results than before)
51     - aim.5l_LatLon & aim.5l_Equatorial_Channel are still using old aim pkg
52    o add diagnostic of surface correction term in monitor output
53    
54    o bulk_force and therm_seaice
55    - Two packages:  bulk_force (Bulk forcing)
56      and therm_seaice (thermodynamic_seaice) - adopted from LANL CICE.v2.0.2
57    - Earlier integration from Stephaine Dutkiewicz
58      and Patrick Heimbach.
59    - Two ifdef statements for compile time,
60       ALLOW_THERM_SEAICE and ALLOW_BULK_FORCE
61    
62    checkpoint47
63    differentiable version of checkpoint46n_post
64    o bug fix in quasihydrostaticterms
65    o * "clean" adjoint code (in terms of extensive recomputations)
66        can now be obtained for all GMREDI options (i.e. for
67        - GM_VISBECK_VARIABLE_K
68        - GM_NON_UNITY_DIAGONAL
69        - GM_EXTRA_DIAGONAL
70        - GM_BOLUS_ADVEC )
71      * However, wrong gradient check problem remains unsolved.
72      * New CPP options have been introduced for different
73        tapering schemes
74    o external_fields_load now part of differentiation list
75    o pressure needs multiple storing;
76      would be nice to have store_pressure at beginning or
77      end of forward_step, e.g. by having phiHyd global (5-dim.)
78      (NB: pressure is needed for certain cases in find_rho,
79      which is also invoked through convective_adjustment).
80    o recomputations in find_rho for cases
81     'JMD95'/'UNESCO' or 'MDJWF' are OK.
82    o #define ATMOSPHERIC_LOADING should be differentiable
83    o ini_forcing shifted to begining of initialise_varia
84    o Incorporating QNVS line search routines into MITgcm
85      (this is separate code, and therefore not under pkg)
86      lsopt/, optim/
87    o Updated verification/carbon/
88      - #define GM_VISBECK_VARIABLE_K
89      - GM_taper_scheme    = 'ldd97'    
90    
91    checkpoint46n_post
92    Merging from release1_p8:
93    o verification/natl_box:
94      updating new external_fields_load routine
95    o New package: pkg/seaice
96      Sea ice model by D. Menemenlis (JPL) and Jinlun Zhang (Seattle).
97      The sea-ice code is based on Hibler (1979-1980).
98      Two sea-ice dynamic solvers, ADI and LSR, are included.
99      In addition to computing prognostic sea-ice variables and diagnosing
100      the forcing/external data fields that drive the ocean model,
101      SEAICE_MODEL also sets theta to the freezing point under sea-ice.
102      The implied surface heat flux is then stored in variable
103      surfaceTendencyTice, which is needed by KPP package (kpp_calc.F and
104      kpp_transport_t.F) to diagnose surface buoyancy fluxes and for the
105      non-local transport term.  Because this call precedes model
106      thermodynamics, temperature under sea-ice may not be "exactly" at
107      the freezing point by the time theta is dumped or time-averaged.
108      N.B.: fairly large differences in accuracy occur across
109            different platforms/compilers; comparison between
110            g77 (Linux) and F77 (SGI Origin 2000/MIPSpro)
111            generated output gives:
112       T           S           U           V
113      C D M    c        m  s        m  s        m  s        m  s
114      n p a R  g  m  m  e  .  m  m  e  .  m  m  e  .  m  m  e  .
115      f n k u  2  i  a  a  d  i  a  a  d  i  a  a  d  i  a  a  d
116      g d e n  d  n  x  n  .  n  x  n  .  n  x  n  .  n  x  n  .
117    
118      Y Y Y Y  5  5  7  7  7  8 10  9  6  6  6  6  7  5  7  5  7 FAIL  lab_sea
119    o GAD:
120      - generated new common blocks to account for call of
121        same gad routines with differing traceridentities
122        (needed to modify tracerIdentity indices in GAD.h)
123      - generated separate common blocks for case useCubedSphereExchange
124        (Department of Futurology)
125      - parameter lists to gmredi_?transport: added tracerIdentity
126      - added new key indices to tamc.h
127    o external_fields_load:
128      - added this routine to TAF list
129      - needed to make some common blocks global and additional storing
130        along the same lines as exf package (checkpoint_lev?_directives.h)
131    o exf:
132      updated external forcing package
133      - bug fixes carried over from ecco-branch
134        (missing OBCS_OPTIONS.h in two routines)
135      - enable easy to use "no forcing".
136      - added exf I/O for atmospheric loading
137      - added exf I/O for runoff data
138      - transfered scaling between exf <-> MITgcm to exf namelist
139    o Adding new verification experiment global_with_exf.
140      This verification experiment is almost identical to the
141      global_ocean.90x40x15 experiment, except that it uses
142      the exf (and cal) package instead of the
143      external_fields_load routine (cf README file).
144    
145    checkpoint46m_post
146    o modified testscript to return -- when there are no lines of "%MON" output
147    o added new flag "quasiHydrostatic" for QH mode
148      and new routine quasihydrostaticterms()
149    o added 2*Omega*Cos(phi)*W to u equations (non-hydrostatic Coriolis term)
150      New routine: mom_u_coriolis_nh()
151    o fixed sign of non-hydrostatic metric terms in u and v equations.
152      Depend on correct sign of gravitySign.
153    o added work.pc* to .cvsignore in aim directories using ifc
154    o introduced flag for controllin non-hydrostatic metric terms
155    o removed gravitySign from PARM04
156    o switched to g77/ifc from pgf77 in aim experiments
157    
158    checkpoint46l_post
159    o made convective adjustment work with pressure coordinates:
160     - changed the direction of k-loop in convective_adjustment.F for the
161       case of pressure coordinates (OCEANICP,ATMOSPHERIC buoyancyRelation)
162     - adjusted the reference pressure k-index in convective_adjustment.F
163     - adjusted the convection condition in convect.F (in analogy to
164       calc_ivdc.F)
165     - convective_adjustment no longer computes anything on the halos
166     - removed the warnings about negative salinity from find_rho.F and
167       find_alpha.F; instead the new routine look_for_neg_salinity, called
168       at the beginning of find_rho, find_alpha, and find_beta, does a
169       check of the entire slice, if CPP-option
170       CHECK_SALINITY_FOR_NEGATIVE_VALUES is defined
171    checkpoint46l_pre
172    
173    checkpoint46k_post
174    o fixed the verification/global_ocean.90x40x15 experiment:
175     - new bathymetry (the world according to A., JMC, and M.)
176     - new initial fields and forcing fields (*.bin files)
177     - new POLY3.COEFFS (for the next release one should switch to a full
178       equation of state: JMD95P or MDJWF)
179     - fixed several errors and redundancies in the data file
180     - experiment uses looped cells
181     - added matlab directory with diagnostic scripts for plotting of output
182    
183    o S/R aim_initialise.F replace S/R aim_init from file aim_do_inphys.F:
184     - read AIM physics parameters from a file (data.aimphys)
185     - set defaults values = F.Molteni paper (Clim.Dyn., 2002)
186    
187    checkpoint46j_post
188    o split calc_exact_eta in 2 S/R : integr_continuity & update_etaH
189    o move wVel computation at the end of the time step, in S/R integr_continuity
190    o create specific S/R to exchange T,S before DYNAMICS (for stagger time step)
191    o update timeave pkg for wVel diagnostic.
192    checkpoint46j_pre
193    
194    checkpoint46i_post
195    o Clean up AIM package (and keep the results unchanged):
196      include CPP_OPTION and use IMPLICT NONE in all routines ;
197      declare all the variables _RL ; use _d 0 for all numerical constants.
198      use ifdef ALLOW_AIM everywhere. And now AIM can be used with g77 !
199    
200    checkpoint46h_post
201    o cleaned up the use of rhoNil and rhoConst.
202      - rhoNil should only appear in the LINEAR equation of state, everywhere
203        else rhoNil is replaced by rhoConst, e.g. find_rho computes rho-rhoConst
204        and the dynamical equations are all divided by rhoConst
205    o introduced new parameter rhoConstFresh, a reference density of fresh
206      water, to remove the fresh water flux's dependence on rhoNil. The default
207      value is 999.8 kg/m^3
208    o cleanup up external_forcing.F and external_forcing_surf.F
209      - can now be used by both OCEANIC and OCEANICP
210    checkpoint46h_pre
211    
212     Added code to convert surface volume flux (fresh water) into
213     a mass flux when using P coordinates in the ocean (OCEANICP).
214     Note: It assumes you have set rho0=rhoConst=density of fresh water.
215    
216    checkpoint46g_post
217    o Include a new diagnostic variable phiHydLow for the ocean model
218      - in z-coordinates, it is the bottom pressure anomaly
219      - in p-coordinates, it is the sea surface elevation
220      - in both cases, these variable have global drift, reflecting the mass
221        drift in z-coordinates and the volume drift in p-coordinates
222      - included time averaging for phiHydLow, be aware of the drift!
223    o depth-dependent computation of Bo_surf for pressure coordinates
224      in the ocean (buoyancyRelation='OCEANICP')
225      - requires a new routine (FIND_RHO_SCALAR) to compute density with only
226        Theta, Salinity, and Pressure in the parameter list. This routine is
227        presently contained in find_rho.F. This routine does not give the
228        correct density for 'POLY3', which would be a z-dependent reference
229        density.
230    o cleaned up find_rho
231      - removed obsolete 'eqn' from the parameter list.
232    o added two new verification experiments: gop and goz
233      (4x4 degree global ocean, 15 layers in pressure and height coordinates)
234    checkpoint46g_pre
235    
236    checkpoint46f_post
237    o Enable tangent linear (forward mode) gradient checks:
238      - extended active file handling to g_... files
239      - added TANGENT_SIMULATION to theSimulationMode
240      - extended grdchk package accordingly
241    o added tangent linear setup in verification/carbon/code_ftl
242    o added adjoint setup in verification/front_relax/code_ad
243    
244    checkpoint46e_post
245    o Added new equation of state -> MDJWF
246      - EOS of McDougall et al., 2002, JAOT, submitted
247      - caveat: the equation of state is only valid for a smaller (more
248        realistic?) range of values than JMD95P/Z and UNESCO
249      - added masks to the calculation of pressure in store_pressure
250      - added more check values for density in check_eos (ini_eos.F), some of
251        the old check values are out of the range of the MDJWF-eos, so don't
252        expect perfect matches for those
253    checkpoint46e_pre
254    
255    checkpoint46d_post
256    o fixed store_pressure to work with both buoyancy relation = 'OCEANIC' and
257      'OCEANICP', also initialised field pressure correctly in ini_eos in the
258      case of pressure coordinates. eosType='JMD95Z' in
259      combination with buoyancyRelation='OCEANICP' now causes an error.
260    o Changed p = pressure(i,j,k,bi,bj)  to  p = pressure(i,j,kRef,bi,bj)
261      in find_alpha/beta.
262    checkpoint46d_pre
263    
264    o Changed p = pressure(i,j,k,bi,bj)  to  p = pressure(i,j,kRef,bi,bj)
265      so that JMD95Z and JMD95P give approptiate static stability in find_rho.
266    
267    checkpoint46c_post
268    o Changes necessary for ocean in p-coordinates
269     -  Added new buoyancy relation = 'OCEANICP'
270     -  Added new parameters = gravitySign (this used to be contained inside
271        the factor dRdZ which I added when we first switched to R coordinates).
272     X GM/Redi is not compatible (yet)
273     X bottom drag and no-slip need to be debugged.
274    checkpoint46c_pre
275    
276    checkpoint46b_post
277    o Added new equation of state -> JMD95Z and JMD95P
278      - EOS of Jackett and McDougall, 1995, JAOT
279      - moved all EOS parameters into EOS.h
280      - new routines ini_eos.F, store_pressure.F
281    o Added UNESCO EOS, but not recommended because it requires
282      in-situ temperature (see JMD95)
283    o Modified formatting for knudsen2.f in utils/knudsen2 and added
284      unesco.f to be used with POLY3
285    checkpoint46b_pre
286    
287    checkpoint46a_post
288    o Added (atmospheric) pressure loading
289      - new field in FFIELDS.h, etc...
290      - new cpp flag ATMOSPHERIC_LOADING
291    o Changed hFacC to _hFacC in calc_phi_hyd.F
292    o Added SHORTWAVE_HEATING to some files for consistency
293    checkpoint46a_pre
294    
295    checkpoint46
296    further merges from release1_p5
297    o new ctrl package
298      - adopted from ECCO environment to enable optimization
299      - added Eliassen Palm fluxes to controls
300    o added Eliassen Palm flux controls to gradient check package
301    o cost package
302      - Compatible with new ctrl package
303      - added Eliassen Palm cost hooks
304      - modif's of existing cost functions
305    o Adjoint-related bug fixes in kpp:
306      - kpp_calc: sore of kpphbl avoids recomputation/call to S/R kppmix
307      - kpp_routines: store of Rib avoids partial recomputation bug of TAF.
308    o autodiff package
309      - flow directives:
310        + new for S/R do_the_model_io
311        + enhanced for S/R checkpoint
312        + new gmredi.flow
313      - Introduced CPP option ALLOW_AUTODIFF_MONITOR to
314        disable adcommon.h/g_common.h by default.
315      - fixed problem for adjoint of global_max
316    o modified verification/carbon/
317    o New setup for adjoint of global_ocean
318    o added hooks for Eliassen Palm fluxes (dynamics.F, ini_fields.F)
319    o cleaned TAF keys iikey, idkey (dynamics.F, thermodynamics.F)
320    o incorporated hooks for sponge layer code (external_forcing.F)
321    
322    checkpoint45d_post
323    o import 2 fixes from release1_p5 (CG2D.h & packages_readparms.F)
324    o use recip_dx*,recip_dy* instead of /dx*,/dy* in orlanski_E,W,N,S
325      (affects plume_on_slope "cg2d" output)
326    o add new parameter deltaTfreesurf for free-surface time stepping
327    o exchange gT,gS when using staggerTimeStep & GM & Oly<4
328    o NONLIN_FRSURF: rescale gNm1 to get a better conservation with A-B
329    o change phiMin in exp1 (to agree with documentation)
330    o add new flags:
331     - T,S forcing inside/outside Adams-Bashforth
332     - turn off temp,salt advection and/or forcing
333     - for each tracer, define internal flag for multiDimAdvection & A-B
334    o monitor: change definition of KE to check conservation of Energy
335    
336    checkpoint45c_post
337    o Removed f90 comments ("!")
338    o Modified initialisations to enable exf package for MITgcm
339      (without invoking ECCO_PACKAGE).
340    o modifications to gradient check package (Martin Losch)
341      - enable centered differences
342      - modified format of standard output
343    
344    checkpoint45b_post
345    o cleaned exchange calls in cg2d
346    o Included CPP option ALLOW_AUTODIFF_MONITOR in forward_step
347    o included CPP option SINGLE_LAYER_MODE in thermodynamics
348      to configure barotropic setup (Martin Losch)
349    o moved some initialisations within bi,bj loop in dynamics
350      (TAF flow dependency)
351    o in initialise_varia, moved call packages_init_variables
352      before convective_adjustment_init
353      (corresponding adjustment of TAF tape initialisation).
354    o in calc_diffusivity/calc_viscosity extended domain to full overlap.
355    
356    checkpoint45a_post
357    o compute Rho over the whole domain (extend iMin,jMin indices)
358    o specific exchange for X-slice domain (case Ny=1)
359    o argument futureIter of S/R obcs_calc missing in exp4 & internal_wave /code
360    o solve_for_pressure : set cg2d_x to zero where OBCS are applied
361    
362    checkpoint44h_post = checkpoint45
363    o Fix to restore differentiability:
364      place at which packages_init_variables is called in
365      initialise_varia is crucial for reverse control flow.
366    o modifications in GMRedi to enable stable adjoint
367      (corresponding modif's in thermodynamics, dynamics)
368    o added missing hook for sponge layer code in external_forcing
369    o modified test cost functions
370    o storing in gad_advection (plus removal of write(0,*)!
371    o missing headers TR1.h added in convective_adjustment routines
372    o errorMessageUnit set non-equal zero in eeboot
373      (conflict on some platforms, e.g. IBM SP3)
374    o modified carbon verif.
375    
376    checkpoint44h_pre
377    o vertical grid option: allow to put the Interface at the middle between
378      2 cell-centers ; replace delZ in pkg/kpp by drF.
379    o GM Advective form: Tracers are advected using the residual transport (=
380      Euler+GM-bolus); set param GM_AdvSeparate=T to return to previous form.
381    
382    checkpoint44g_post
383    o fix surface correction term in multiDimAdvection (affects plume_on_slope
384      results) ; ifdef missing in thermodynamics.F
385    
386    checkpoint44f_pre,checkpoint44f_post
387    o added PTRACERS package
388      This allows an arbitrary number of passive tracers to be integrated
389      forward simultaneously with the dynamicaly model.
390      + Implemented so far:
391        - basic forward algorithm (time-stepping, advection, diffusion, convection)
392        - I/O and checkpointing
393        - GM/Redi  *but*  using the GM/Redi coefficient of Salt
394      + Not implemented so far:
395        - KPP
396        - OBCS
397      + No specific example supplied (yet) but global_ocean.90x40x15 has the
398        necessary data.ptracer file. Simply use -enable=ptracers and uncomment
399        line in data.pkg.
400      + This package is disabled by default since it increases storage.
401    
402    
403    o testscript: added -cleanup option
404    o Shap_filt: a) new shap_filt S/R to use no-slip BC with S2 filter
405                 b) enable to filter 2D fields.
406    
407    checkpoint44e_pre,post
408    o re-wrote convective adjustment to use pre-calculated weights
409      Affects: model/src/convective_adjustment.F, convective_adjustment_ini.F
410      Adds:    model/src/convective_weights.F, convectively_mixtracer.F
411      Original calls to convect() still in place.
412    
413    chkpt44d_pre,post
414    o added missing EXCLUDE_MONITOR flags
415    o changed "e" to "_d" in gmredi_slope_limit, gmredi_slope_psi
416      (incompatible typ in MIN/MAX expressions caused problems
417      on IBM SP3)
418    o in genmake added variable MAKEDEPEND
419      plus resetting for case SunOS
420    o added timer_stats.c routine for IBM SP3
421    o removed variables in dynamics
422    !!! I forgot to create tag chkpt44d_pre !!!
423    !!! But chkpt44d_pre = chkpt44c_post    !!!
424    
425    o real fresh water flux implemented with non-linear free-surface.
426    
427    chkpt44c_pre,post
428    o few fix (mask in shap_s2, EmPmR in external_field_load,
429      USE_NATURAL_BCS in solve_for_P);
430    o add arguments myIter & myTime to S/R obcs_calc & solve_for_P
431    
432    checkpoint44b_pre/_post
433    o merge of relevant stuff from the ecco-branch:
434      - genmake: removed $S64 overwrite for case SunOS
435      - pkg/exf: update and corrections for field swapping and obcs
436      - pkg/ecco: parameter lists for the_model_main, the_main_loop
437                  harmonized between ECCO and MITgcm
438      - pkg/autodiff: added flow directives for obcs, mdsio_gl_slice
439                      updated checkpointing_lev... lists for obcs
440      - model/src: minor changes in forward_step, plot_field
441                   added directive for divided adjoint in the_main_loop
442      - pkg/mdsio: added mdsio_gl_slice
443    
444    o check parameters & config (chkpt44a_pre,post)
445    o OBC and NonLin_FrSurf.
446    
447    checkpoint44
448    o fix bug in mom_vi_del2uv
449    o select when filters are applied ; add options to zonal_filter (data.zonfilt)
450    o gmredi: fix Pb in the adiabatic form ; add options (.e.g. Bolus advection)
451    o update AIM experiments (NCEP input files)
452    o improve and extend diagnostics (Monitor, TimeAve with NonLin-FrSurf)
453    o added some stuff for AD
454    
455    checkpoint43
456    o added two new advection tests
457    o added global 4x4 degree ocean test
458    o added test of H&S94 with free surface on cube
459    o some final mods for gradient check, adjoint,...
460    
461    checkpoint42
462      Fixed AD-related problems:
463    o Store directives up-to-date with re-arranged Adams-Bashforth
464      (mainly thermodynamics.F)
465    o New store directives for multi-dim. advection schemes
466      * new CPP flag ALLOW_MULTI_DIM_ADVECTION
467      * new common block and key passkey
468      (mainly gad_advection.F)
469    o Modified store directives for split of dynamics/thermodynamics
470      for the case ALLOW_KPP
471    o Cleaned argument list for timestep_tracer.F
472    
473    checkpoint41
474    o re-formatted/added comments in prologues of key subroutines
475      to work with protex
476    
477    checkpoint40
478    o moved calc_gs, calc_gt and calc_tr1 back to model/src
479    o added FLT (floats) package from Arne along with demo. experiment
480    o re-arranged Adams-Bashforth scheme for tracers to allow easier
481      implementation of forward-in-time forcing terms
482    o more mods for non-linear free-surface including new variable
483    o modified multi-dim method to work on cube
484    o bug fix in KPP (from DM)
485    
486    checkpoint40pre9
487    o added DST advection schemes
488    o fix MPI bug on SGI's
489    o modifications for non-linear free-surface
490    o fixed conditionals for fixed flow/passive tracer only mode
491    o added tension/strain form of stress tensor
492    o changes to commenting of arguments for TAF/TAMC
493    o added/updated many matlab scripts
494    o added multi-dimensional advection method for tracers
495      
496    checkpoint40pre8
497    o new package rw/
498      Added method for dumping intermediate local arrays:
499      mdsio_writetile - same as mdsio_writefield except works from inside
500                        bi,bj loop
501      mdsio_writelocal - same as mdsio_writetile except works for local
502                         arrays
503      write_local_r? - higher-level wrapper for mdsio_writelocal
504      Controlled by diagFreq. Defaults to zero (ie. no dumps)
505      Example given at end of mom_vecinv.F that dumps some local arrays.
506    o replaced calc_mom_rhs.F by mom_fluxform.F/mom_vecinc.F to
507      control flux-form versus vector invariant code.
508      Switch is runtime flag vectorInvariantMomentum (defaults to false)
509      or genmake -DISABLE ...
510    o TAMC-related update to cope with mom_fluxform.F/mom_vecinc.F
511      (additional storing is now needed for case vecinv;
512      currently nothing implemented to switch this storing off
513      in case fluxform is used since no CPP flag available).
514    
515    checkpoint40pre7
516    o bug fix in dynamics.F
517      calc_viscosity called before k=1,Nr loop.
518    o added some #ifdef's ALLOW_PASSIVE_TRACER
519    o cleaned initialisations in ctrl_map_ini
520      related to TAF
521    
522    checkpoint40pre6
523    o split dynamics into thermodynamics & dynamics
524    o corresponding split of calc_diffusivity into
525      calc_diffusivity, calc_viscosity
526      (kpp_calc_diff, kpp_calc_visc).
527    o Added diffkr, kapgm to control variables
528    o bug fix in gmredi_slope_limit for ldd97
529    
530    checkpoint40pre5
531    o proper initialization of Adams-Bashforth time stepping
532      + uses forward step for very first step of model (iter=0)
533    o re-generated *all* output since almost all output is modified by this
534      change (sorry but we've been putting this one off for about 2 years)
535      + natl_box, aim.5l_Equatorial_Channel and aim.5l_zon-ave were
536        *not* regenerated since they are already out of date and I don't
537        want to create the impression that they are suddenly working.
538    o removed KAP stuff for DEC's
539    o modified/added many .cvsignore files in */code and */input
540    o added new expt "solid-body.cs-32x32x1"
541      + this uses an alternative form of ini_curvilinear_grid.F which
542        is not necessarily any better than the current one...
543    
544    checkpoint40pre4
545    o added re-initialisations of utrans,vtrans,rtrans in
546      calc_common_factors which are crucial for TAF.
547    o extended iMin,jMin range for calc_common_factors,
548      calc_diffusivity
549    o added option ALLOW_PASSIVE_TRACER
550    o a few minor updates and fixes
551    
552    checkpoint40pre3
553    o created calc_gtr1.F for gad package
554    o split the_main_loop and re-introduced forward_step
555    o some recomputations of nIter0 for TAMC
556    
557    checkpoint40pre2
558    o Updated to latest branch of AJA, JMC
559    o Added gradient check package
560    o Added routines for passive tracer handling
561    o Added carbon sequestration sensitivity experiment
562    
563    checkpoint39
564     o Modifications and fixes to support TAMC
565       (mainly missing or wrong directives)
566     o Corrected common blocks when fooling TAMC to produce adjoint I/O
567       in addummy_in_stepping
568     o Updated verification/vero/ (bug fixes in code/, correct pickups in input/)
569     o Enabled "sliced" (xz-,yz-) exchanges, needed for obcs.
570     o Included following packages:
571       cal: calendar
572       exf: external forcing
573       ecco: essential routines to control use of cal, exf
574             (NB: this is not the full ECCO environment)
575     o disabled some packages by default in genmake
576    
577    checkpoint38
578     o see notes_c37_adj.txt
579     o corrected wrong layer index in external_forcing_t
580     o replaced the_main_loop/forward_step by the_main_loop
581     o removed forceoutput from write_state
582       (replaced criteria for initial/final state)
583     o added runtime flag for pre-C35 pickups: usePickupBeforeC35
584    
585    checkpoint37
586     o mdsio.F routine split and packaged in pkg/mdsio
587       (mdsio.F and mdsio_gl.F still exist in eesupp/src but are not
588        used unless genmake -disable=mdsio)
589     o allow a more precise surface-geopotential gradient for the atmosphere
590       taking into account the change in space of Alpha (=1/rho) at the surface.
591     o variable "etaN" (replace cg2d_x) is the surface R-anomaly [r unit]
592       (ocean: sea surface height [m] / atmos: surface pressure anomaly [Pa])
593     o Change units: Phi_Hydrostatic (=phiHyd), Phi_Surface(=cg2d_x) and
594       Phi_NonHydrostatic (=cg3d_x), all have now the dimension of a potential
595       [m2/s2] (pressure/rhoConst in the ocean ; Geo-Potential in the atmosphere).
596     o Time-Average diagnostics defined as a pseudo package "timeave".
597     o pickup files : old pickup (before c35) did not work with new code (c35
598       and after) unless activating 2 commented lines in initialise-varia.
599     o Non-Hydrostatic pickup file Pb fixed.
600     o bug with Rigid-lid version (since c35) fixed.  
601    
602    checkpoint36
603     o implement Crank-Nickelson barotropic time-stepping
604     o recover Time-Average diagnostics
605     o fix Pb in checkpoint (PR#19)
606     o AIM become a standard package ; clean unused CPP_kees
607    
608    checkpoint35
609     o subtantial rewrite of dynamics.f
610     o allows staggered time-stepping
611     o packaged obcs, aim (Atmospheric Physics), zonal_filt, shap_filt
612     o added bottom drag (linear + quadratic)
613    
614    checkpoint34
615     o Modified exchange routines for correct reverse simulation.
616     o Added verification output to comply with bug fix in impldiff.F
617       as of c33 (exp2, natl_box).
618     o Corrected store directives in dynamics.F
619     o Corrected array boundaries in impldiff.F
620     o Corrected array dimensions in gmredi.
621       Added temporary scalars to avoid storage by TAMC (gmredi, kpp).
622     o Added routine GLOBAL_SUM_INT to global_sum.F
623     o Added mdsio_gl.F for ECCO purposes.
624    
625    checkpoint33
626     o Fixed bug in lower BC in impldiff.F
627     o Fixed confusion about units of forcing arrays in FFIELDS.h
628       namely Fu,Fv,Qnet,Qsw,EmPmR:
629       - Removed verification/natl_box/code/external_fields_scale.F
630         (did not differ from that in model/src)
631       - Changed units of fu,fv,Qnet,Qsw,EmPmR back to proper units
632         (see FFIELDS.h for description)
633       - Scale fu,fv,Qnet,Qsw,EmPmR when used in external_forcing_surf.F,
634         kpp_calc.F and kpp_transport_t.F
635       - Removed model/src/external_fields_scale.F and calls to it
636       - verification/natl_box uses flux data with "atmospheric" sign so
637         a special version of external_fields_load.F is used to
638         change the data as it's read in. This way, the arrays
639         have the right units and signs at all times tha a user could
640         possibly use them.
641     o Corrected genmake to delete a temporary file during configuration
642    
643    checkpoint32
644     o Scaling of forcing fields moved from external_forcing.F to
645       external_fields_scale.F, called directly after loading fields.
646     o Surface relaxation terms added to surface tendency fields in
647       extra subroutine.
648     o impldiff.F changed to avoid storing by TAMC of huge intermediate
649       fields.
650     o Changes in GMRedi to comply with packaging structure.
651     o Modifications and bug fixes in KPP package.
652     o Verification experiment for KPP added (verification/nat_box/)
653    
654    checkpoint31
655     o Packages interface brought in line with defined standard
656       -> new routines packages_...F
657       -> changes in initialise_fixed.F, initialise_varia.F
658       -> enables initialisation interface for ECCO package
659     o Rescaling of external forcing fields immediately after read
660       (removed from external_forcing.F)
661       -> new: surface_fields_scale.F
662       -> affects ini_foorcing.F, external_fields_load.F
663     o Computation of surface tendencies before tendency updates.
664       -> new: external_forcing_surf.F
665     o Change of variable names useGMRedi, useKPP, useECCO
666     o Further changes in dynamics.F, impldiff.F for TAMC
667     o Tested for experiments 0, 2, 4.
668       - not yet tested for GMRedi.
669       - not yet tested for KPP.
670    
671    checkpoint30
672     o Updates for OBCs
673     o New experiment for test OBCs (internal_wave)
674     o New output.txt in verification/exp4/results/
675     o Updates for KPP: calculation of viscosity corrected, new variables
676       kpp_dumpFreq and kpp_taveFreq
677     o CPP'd calls to system routines for TARGET_CRAY_VECTOR (P.H.)
678       and entry for T90 in genmake
679    
680    checkpoint29
681     o packaged KPP and GM/Redi
682    
683    checkpoint28
684     o Corrected initialization of hFacMin for difference vertical coordinates
685     o Modified calculation of partial cell thickness to more robust algorithm
686     o config_summary.F: replaced write_1D_... by write_0D_...
687       read_write.F: added subroutines write_0D_...
688     o routines calc_common_factors, calc_gs, calc_gt, calc_phi_hyd, find_rho:
689       included some initialisations required for TAMC.
690     o routine calc_isoslopes: modified sqrt computation to make it
691       differentiable for zero argument.
692     o routines dynamics, impldiff:
693       included store directives and tape key computations required for TAMC.
694     o ALL TAMC related changes are between ifdef ALLOW_AUTODIFF_TAMC
695    
696    checkpoint27
697     o fixed calc_isoslopes() to so that cancellation of terms in tensor
698       formulation works properly with variable resolution.
699     o restructured the calling tree between the_model_main() and dynamics()
700     o split initialise() into initialise_fixed() and initialise_varia()
701     o introduced initialization file pSurfInitFile for free surface (cg2d_x)
702       in new routine ini_psurf()
703    
704    checkpoint26
705     o moved some #ifdef in optional routines to encompass everything
706       except the argument declarations. This was using up some memory.
707     o fixed macros like _ADT to work with both versions of Linux cpp
708     o cleaned up some unused and uninitialized variables
709       (helps when debugging with strict compile options)
710     o split up CPP_EEOPTIONS.h into options and macros
711       (macros are now in eesupp/inc/CPP_EEMACROS.h)
712     o patch for 2 processor JAM mode
713    
714    checkpoint25
715     o updates for OBCs and NH
716     o dramatically reduced memory usage
717     o added "JAM" routines for use on Artic network
718     o parameterization of convection by implicit vertical diffusion
719     o added a diagnostic of convective index
720     o moved call to calc_divg_hat() from dynamics() to solve_for_pressure()
721       (this is was partly for OBCs and NH but also in anticipation of
722        an explicit free-surface option)
723     o initial state input files for U and V
724    
725    checkpoint24
726     o Removed some spurious "REAL"s
727     o Updated KPP to latest version from JPL
728     o Added a COS(latitude) in horizontal viscosity
729    
730    checkpoint23
731     o Added del^4 diffusion for salt and temperature
732     o Multiple minor fixes (implicit none's, arguments to global_sum, etc.)
733    
734    checkpoint22
735     o Made many small changes: implicit none, format(i5), etc...
736     o Introduced KPP mixing scheme (compiles and runs but *not* tested)
737    
738    checkpoint21
739     o Introduced new I/O package that supports direct-access records (mdsio.F)
740     o Split-up mid-level I/O routines
741     o Updated input data for experiments 0,1,2,4 and 5
742     o Finally removed the genmake.* scripts (everyone seems happy with genmake)
743    
744    checkpoint20
745     o Tidied up open-boundaries some more
746     o Added non-hydrostatic code
747     o Fixed some input data
748     o Modified MPI start-up to work in "coupled" context.
749    
750    checkpoint19
751     o Tidied up directory structure for verification experiments.
752     o Added new code for open boundaries.
753     o Corrected advective terms near boundaries.
754     o Added CPP flag for "old UV" geometry and corrected spherical geometry.
755     o Implemented free-slip/no-slip BCs and del^4 dissipation.
756     o Split default parameters out of ini_parms.F into set_defaults.F.
757    
758    checkpoint18
759     o Further changes for optional code generation with
760       CPP - ugh!
761    
762    checkpoint17
763     o Beginning to incorporating modularity
764       for supporting atmos. and ocean configs.
765     o Small changes for TAMC complinace
766     o New ultra simple test case. Barotropic box
767       configuration "exp0".
768    
769    checkpoint16
770     o Minor changes to exchange routines to support
771       accumulation in reverse mode. Adjoint of an
772       assignment is an addition.
773     o Changes to support compilation by TAMC (LR(1)
774       syntax) and by g77 under Linux ( 72 character
775       limit ).
776    
777    checkpoint15
778     Rewrote exchange routines to allow
779     o General tile <-> tile connectivity.
780     o TAMC reverse mode flag ( is this really needed? )
781     o Variable width overlap and exchange region widths.
782     o Hooks for sharing data via DMA, Arctic, Memory Channel,
783       shmput/shmget, VIA, SYSV shmem and every other cool
784       communication method.
785    
786    checkpoint14
787     Consistent isomorphism chages made. These allow dynamical core
788     to support both pressure and meters as a vertical coordinate.
789    
790    checkpoint13
791     Release which corrects global_max problem and a couple
792     of KGM time-averaging diagnostic errors. Also contains
793     a double-gyre experiment setup as well as the baseline
794     4 degree global template. The double gyre template is
795     organised to allow it to be applied automatically to
796     the baseline configuration.
797    
798    branch-point-rdot
799     A few tidy-ups have been made bt the real point of this
800     check-point (excuse the pun) is to allow a branch to be
801     made for the purposes of adding the rDot alterations.
802     I don't want to commit the rDot to the main-trunk until
803     we've actually tested it...
804    
805    checkpoint12
806     Tidy up and added generic genmake script (formerly
807     configure script).
808     Released to testers.
809    
810    checkpoint11
811     Re-arranged initialise() and introduced ini_vertical_grid()
812     and ini_masks_etc().
813    
814  checkpoint10  checkpoint10
815   Introduced the configure script.   Introduced the configure script.
816   Separated out the diags package.   Separated out the diags package.

Legend:
Removed from v.1.8  
changed lines
  Added in v.1.111

  ViewVC Help
Powered by ViewVC 1.1.22