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

Legend:
Removed from v.1.1  
changed lines
  Added in v.1.114

  ViewVC Help
Powered by ViewVC 1.1.22