model/src/calc_phi_hyd.F

C $Header: /u/gcmpack/models/MITgcmUV/model/src/calc_phi_hyd.F,v 1.3 1998/09/05 17:52:13 cnh Exp $

#include "CPP_EEOPTIONS.h"

      SUBROUTINE CALC_PHI_HYD( bi, bj, iMin, iMax, jMin, jMax, K,
     I    buoyKM1, buoyKP1, phiHyd, myThid)
C     /==========================================================\
C     | SUBROUTINE CALC_PHI_HYD                                  |
C     | o Integrate the hydrostatic relation to find phiHyd.     |
C     |                                                          |
C     \==========================================================/
      IMPLICIT NONE
C     == Global variables ==
#include "SIZE.h"
#include "DYNVARS.h"
#include "GRID.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
C     == Routine arguments ==
      INTEGER bi,bj,iMin,iMax,jMin,jMax,K
      _RL buoyKM1(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
      _RL buoyKP1(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
      _RL phiHyd(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
      integer myThid
C     == Local variables ==
      INTEGER i,j,Km1
      _RL halfLayer

      if (K.eq.1) then
       Km1=1
       halfLayer=0.5 _d 0
      else
       Km1=K-1
       halfLayer=1.0 _d 0
      endif

C--   Contribution to phiHyd(:,:,K) from buoy(:,:,K-1) + buoy(:,:,K)
C     (This is now the actual hydrostatic pressure|height at the T/S points)
      DO j=jMin,jMax
       DO i=iMin,iMax
        phiHyd(i,j,K)=phiHyd(i,j,Km1)-rhoConst*halfLayer
     &  *0.5 _d 0*( drF(Km1)+drF(K) )*recip_HoriVertRatio
     &  *0.5 _d 0*( buoyKM1(i,j)+buoyKP1(i,j) )
C    &  *rkFac
       ENDDO
      ENDDO

! ------------------------------------------------------------------------------
      return
      end
! ==============================================================================
1	cnh	1.4	C $Header: /u/gcmpack/models/MITgcmUV/model/src/calc_phi_hyd.F,v 1.3 1998/09/05 17:52:13 cnh Exp $
2	cnh	1.1
3			#include "CPP_EEOPTIONS.h"
4
5			SUBROUTINE CALC_PHI_HYD( bi, bj, iMin, iMax, jMin, jMax, K,
6	cnh	1.2	I buoyKM1, buoyKP1, phiHyd, myThid)
7	cnh	1.1	C /==========================================================\
8			C \| SUBROUTINE CALC_PHI_HYD \|
9			C \| o Integrate the hydrostatic relation to find phiHyd. \|
10			C \| \|
11			C \==========================================================/
12			IMPLICIT NONE
13			C == Global variables ==
14			#include "SIZE.h"
15			#include "DYNVARS.h"
16			#include "GRID.h"
17			#include "EEPARAMS.h"
18			#include "PARAMS.h"
19			C == Routine arguments ==
20			INTEGER bi,bj,iMin,iMax,jMin,jMax,K
21			_RL buoyKM1(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
22			_RL buoyKP1(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
23	cnh	1.2	_RL phiHyd(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
24	cnh	1.1	integer myThid
25			C == Local variables ==
26			INTEGER i,j,Km1
27			_RL halfLayer
28
29			if (K.eq.1) then
30			Km1=1
31			halfLayer=0.5 _d 0
32			else
33			Km1=K-1
34			halfLayer=1.0 _d 0
35			endif
36
37			C-- Contribution to phiHyd(:,:,K) from buoy(:,:,K-1) + buoy(:,:,K)
38			C (This is now the actual hydrostatic pressure\|height at the T/S points)
39			DO j=jMin,jMax
40			DO i=iMin,iMax
41	cnh	1.3	phiHyd(i,j,K)=phiHyd(i,j,Km1)-rhoConst*halfLayer
42	cnh	1.4	& 0.5 _d 0( drF(Km1)+drF(K) )*recip_HoriVertRatio
43	cnh	1.1	& 0.5 _d 0( buoyKM1(i,j)+buoyKP1(i,j) )
44	cnh	1.4	C & *rkFac
45	cnh	1.1	ENDDO
46			ENDDO
47
48			! ------------------------------------------------------------------------------
49			return
50			end
51			! ==============================================================================