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	C $Header: /u/gcmpack/models/MITgcmUV/model/src/calc_phi_hyd.F,v 1.3 1998/09/05 17:52:13 cnh Exp $
2
3	#include "CPP_EEOPTIONS.h"
4
5	SUBROUTINE CALC_PHI_HYD( bi, bj, iMin, iMax, jMin, jMax, K,
6	I buoyKM1, buoyKP1, phiHyd, myThid)
7	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	_RL phiHyd(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
24	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	phiHyd(i,j,K)=phiHyd(i,j,Km1)-rhoConst*halfLayer
42	& 0.5 _d 0( drF(Km1)+drF(K) )*recip_HoriVertRatio
43	& 0.5 _d 0( buoyKM1(i,j)+buoyKP1(i,j) )
44	C & *rkFac
45	ENDDO
46	ENDDO
47
48	! ------------------------------------------------------------------------------
49	return
50	end
51	! ==============================================================================