--- MITgcm/doc/tag-index	2002/08/15 17:25:26	1.81
+++ MITgcm/doc/tag-index	2002/10/07 16:35:31	1.89
@@ -1,8 +1,88 @@
-$Header: /home/ubuntu/mnt/e9_copy/MITgcm/doc/tag-index,v 1.81 2002/08/15 17:25:26 adcroft Exp $
+$Header: /home/ubuntu/mnt/e9_copy/MITgcm/doc/tag-index,v 1.89 2002/10/07 16:35:31 jmc Exp $
 
     Notes on tags used in MITgcmUV
     ==============================
 
+checkpoint46j_post
+o split calc_exact_eta in 2 S/R : integr_continuity & update_etaH
+o move wVel computation at the end of the time step, in S/R integr_continuity
+o create specific S/R to exchange T,S before DYNAMICS (for stagger time step)
+o update timeave pkg for wVel diagnostic.
+checkpoint46j_pre
+
+checkpoint46i_post
+o Clean up AIM package (and keep the results unchanged):
+  include CPP_OPTION and use IMPLICT NONE in all routines ; 
+  declare all the variables _RL ; use _d 0 for all numerical constants.
+  use ifdef ALLOW_AIM everywhere. And now AIM can be used with g77 !
+
+checkpoint46h_post
+o cleaned up the use of rhoNil and rhoConst. 
+  - rhoNil should only appear in the LINEAR equation of state, everywhere 
+    else rhoNil is replaced by rhoConst, e.g. find_rho computes rho-rhoConst
+    and the dynamical equations are all divided by rhoConst
+o introduced new parameter rhoConstFresh, a reference density of fresh 
+  water, to remove the fresh water flux's dependence on rhoNil. The default
+  value is 999.8 kg/m^3
+o cleanup up external_forcing.F and external_forcing_surf.F
+  - can now be used by both OCEANIC and OCEANICP
+checkpoint46h_pre
+
+ Added code to convert surface volume flux (fresh water) into
+ a mass flux when using P coordinates in the ocean (OCEANICP).
+ Note: It assumes you have set rho0=rhoConst=density of fresh water.
+
+checkpoint46g_post
+o Include a new diagnostic variable phiHydLow for the ocean model
+  - in z-coordinates, it is the bottom pressure anomaly
+  - in p-coordinates, it is the sea surface elevation
+  - in both cases, these variable have global drift, reflecting the mass
+    drift in z-coordinates and the volume drift in p-coordinates
+  - included time averaging for phiHydLow, be aware of the drift!
+o depth-dependent computation of Bo_surf for pressure coordinates 
+  in the ocean (buoyancyRelation='OCEANICP')
+  - requires a new routine (FIND_RHO_SCALAR) to compute density with only 
+    Theta, Salinity, and Pressure in the parameter list. This routine is 
+    presently contained in find_rho.F. This routine does not give the
+    correct density for 'POLY3', which would be a z-dependent reference
+    density.
+o cleaned up find_rho
+  - removed obsolete 'eqn' from the parameter list.
+o added two new verification experiments: gop and goz
+  (4x4 degree global ocean, 15 layers in pressure and height coordinates)
+checkpoint46g_pre
+
+checkpoint46f_post
+o Enable tangent linear (forward mode) gradient checks:
+  - extended active file handling to g_... files
+  - added TANGENT_SIMULATION to theSimulationMode
+  - extended grdchk package accordingly
+o added tangent linear setup in verification/carbon/code_ftl
+o added adjoint setup in verification/front_relax/code_ad
+
+checkpoint46e_post
+o Added new equation of state -> MDJWF
+  - EOS of McDougall et al., 2002, JAOT, submitted
+  - caveat: the equation of state is only valid for a smaller (more
+    realistic?) range of values than JMD95P/Z and UNESCO
+  - added masks to the calculation of pressure in store_pressure
+  - added more check values for density in check_eos (ini_eos.F), some of 
+    the old check values are out of the range of the MDJWF-eos, so don't
+    expect perfect matches for those 
+checkpoint46e_pre
+
+checkpoint46d_post
+o fixed store_pressure to work with both buoyancy relation = 'OCEANIC' and
+  'OCEANICP', also initialised field pressure correctly in ini_eos in the 
+  case of pressure coordinates. eosType='JMD95Z' in 
+  combination with buoyancyRelation='OCEANICP' now causes an error.
+o Changed p = pressure(i,j,k,bi,bj)  to  p = pressure(i,j,kRef,bi,bj)
+  in find_alpha/beta.
+checkpoint46d_pre
+
+o Changed p = pressure(i,j,k,bi,bj)  to  p = pressure(i,j,kRef,bi,bj)
+  so that JMD95Z and JMD95P give approptiate static stability in find_rho. 
+
 checkpoint46c_post
 o Changes necessary for ocean in p-coordinates
  -  Added new buoyancy relation = 'OCEANICP'
@@ -14,7 +94,7 @@
 
 checkpoint46b_post
 o Added new equation of state -> JMD95Z and JMD95P
-  - EOS of Jackett and McDougall, 1995, JPO
+  - EOS of Jackett and McDougall, 1995, JAOT
   - moved all EOS parameters into EOS.h
   - new routines ini_eos.F, store_pressure.F
 o Added UNESCO EOS, but not recommended because it requires