With the hydrostatic assumption
The first term on the right hand side in (17) and (18) is the external pressure-gradient due to surface slopes, and the second the internal pressure-gradient due to the density gradient. The internal pressure-gradient will only be established by gradients of mean potential temperature and mean salinity . Sediment concentration is assumed to be horizontally homogeneous.
In this subroutine, first, the horizontal buoyancy gradients, and , are calculated from the prescribed gradients of salinity, and , and temperature, and , according to the finite-difference expression
The horizontal salinity and temperature gradients have to supplied by the user, either as constant values or as profiles given in a file (see obs.nml).
use meanflow, only: T,S use meanflow, only: gravity,rho_0,h use observations, only: dsdx,dsdy,dtdx,dtdy use observations, only: idpdx,idpdy,int_press_method use eqstate, only: eqstate1 IMPLICIT NONEINPUT PARAMETERS:
number of vertical layers integer, intent(in) :: nlevREVISION HISTORY:
Original author(s): Hans Burchard & Karsten Bolding $Log: intpressure.F90,v $ Revision 1.8 2007-01-06 11:49:15 kbk namelist file extension changed .inp --> .nml Revision 1.7 2005/08/11 12:32:50 lars corrected error in Latex referencing Revision 1.6 2005/06/27 13:44:07 kbk modified + removed traling blanks Revision 1.5 2004/08/18 11:43:51 lars updated documentation Revision 1.4 2003/03/28 09:20:35 kbk added new copyright to files Revision 1.3 2003/03/28 08:56:56 kbk removed tabs Revision 1.2 2003/03/10 08:50:06 gotm Improved documentation and cleaned up code Revision 18.104.22.168 2001/02/12 15:55:57 gotm initial import into CVS
Karsten Bolding 2012-01-24