[MITgcm] Log of /MITgcm/model/src/ini

o compute locally the pressure for use in EOS : UNESCO, JMD95P or MDJWF
o store total Potential in totPhyHyd for diagnostic (DYNVARS.h)
o fix restart and overlap Pb when using Z-coord and EOS funct. of P

Revision 1.9 - (view) (download) (annotate) - [select for diffs]
Mon Feb 3 21:58:09 2003 UTC (21 years, 3 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint48d_post, checkpoint48d_pre, checkpoint48e_post, checkpoint48f_post
Changes since 1.8: +23 -5 lines
Diff to previous 1.8

Wrong Restart when using Z-coord with EOS function of P. => add a warning

Revision 1.7 - (view) (download) (annotate) - [select for diffs]
Tue Oct 29 20:16:28 2002 UTC (21 years, 6 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint46l_post, checkpoint46l_pre
Changes since 1.6: +78 -57 lines
Diff to previous 1.6

o fixed a bug in ini_eos: changed e to _d for all MDJWF-coefficients

Revision 1.6 - (view) (download) (annotate) - [select for diffs]
Wed Sep 25 19:36:50 2002 UTC (21 years, 7 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint46h_post, checkpoint46i_post, checkpoint46j_post, checkpoint46j_pre, checkpoint46k_post
Changes since 1.5: +4 -4 lines
Diff to previous 1.5

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

Revision 1.5 - (view) (download) (annotate) - [select for diffs]
Wed Sep 18 16:38:02 2002 UTC (21 years, 8 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint46g_post, checkpoint46h_pre
Changes since 1.4: +15 -7 lines
Diff to previous 1.4

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)

Revision 1.4 - (view) (download) (annotate) - [select for diffs]
Thu Sep 5 20:49:33 2002 UTC (21 years, 8 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint46e_post, checkpoint46f_post, checkpoint46g_pre
Changes since 1.3: +85 -26 lines
Diff to previous 1.3

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

Revision 1.3 - (view) (download) (annotate) - [select for diffs]
Tue Aug 20 14:45:18 2002 UTC (21 years, 9 months ago) by mlosch
Branch: MAIN
Changes since 1.2: +2 -2 lines
Diff to previous 1.2

o bug-fix in ini_eos.F: check_eos still called find_rhop0 with the wrong
  number of arguments

Revision 1.2 - (view) (download) (annotate) - [select for diffs]
Mon Aug 19 14:21:30 2002 UTC (21 years, 9 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint46d_post, checkpoint46e_pre
Changes since 1.1: +47 -12 lines
Diff to previous 1.1

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.

Revision 1.1 - (view) (download) (annotate) - [select for diffs]
Wed Aug 7 16:55:52 2002 UTC (21 years, 9 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint46b_post, checkpoint46c_post, checkpoint46c_pre, checkpoint46d_pre

o Added new equation of state -> JMD95Z and JMD95P
  - EOS of Jackett and McDougall, 1995, JPO
  - 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
  in-situ temperature (see JMD95)
o Modified formatting for knudsen2.f in utils/knudsen2 and added
  unesco.f to be used with POLY3

This form allows you to request diffs between any two revisions of this file. For each of the two "sides" of the diff, select a symbolic revision name using the selection box, or choose 'Use Text Field' and enter a numeric revision.

Links to HEAD:	(view) (download) (annotate)
Links to checkpoint55a_post:	(view) (download) (annotate)
Sticky Tag:

	ViewVC Help
Powered by ViewVC 1.1.22

Log of /MITgcm/model/src/ini_eos.F