model/src/calc_div_ghat.F

C $Header: /u/gcmpack/models/MITgcmUV/model/src/calc_div_ghat.F,v 1.9 1999/05/18 17:42:01 adcroft Exp $

#include "CPP_OPTIONS.h"

C     /==========================================================\
C     | S/R CALC_DIV_GHAT                                        |
C     | o Form the right hand-side of the surface pressure eqn.  |
C     |==========================================================|
C     \==========================================================/
      SUBROUTINE CALC_DIV_GHAT( 
     I        bi,bj,iMin,iMax,jMin,jMax,
     I        K,
     I        xA,yA,
     I        myThid)

      IMPLICIT NONE

C     == Global variables ==
#include "SIZE.h"
#include "DYNVARS.h"
#include "FFIELDS.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "CG2D.h"
#ifdef ALLOW_NONHYDROSTATIC
#include "CG3D.h"
#endif

C     == Routine arguments ==
C     bi, bj, iMin, iMax, jMin, jMax - Range of points for which calculation
C                                      results will be set.
C     k                              - Index of layer.
C     xA, yA                         - Cell face areas
C     myThid - Instance number for this innvocation of CALC_MOM_RHS
      INTEGER bi,bj,iMin,iMax,jMin,jMax
      INTEGER K
      _RS xA (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
      _RS yA (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
      INTEGER myThid

C     == Local variables ==
      INTEGER i,j
      _RL pf (1-OLx:sNx+OLx,1-OLy:sNy+OLy)

C--   Pressure equation source term
C     Term is the vertical integral of the divergence of the 
C     time tendency terms along with a relaxation term that
C     pulls div(U) + dh/dt back toward zero.

      IF ( k .EQ. Nr ) THEN
C      Initialise source term on first pass
       DO j=1,sNy
        DO i=1,sNx
C Note: The source term containing cg2d_x and cg3d_x (at k=1)
C       has been moved to solve_for_pressure.F for convenience.
         cg2d_b(i,j,bi,bj) =
     &     freeSurfFac*_rA(i,j,bi,bj)*horiVertRatio*( 0.
#ifdef USE_NATURAL_BCS
     &      +EmPmR(I,J,bi,bj)/deltaTMom
#endif
     &     )
        ENDDO
       ENDDO
      ENDIF

      DO j=1,sNy
       DO i=1,sNx+1
        pf(i,j) = xA(i,j)*gUNm1(i,j,k,bi,bj) / deltaTmom
       ENDDO
      ENDDO
      DO j=1,sNy
       DO i=1,sNx
        cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj) +
     &   pf(i+1,j)-pf(i,j)
       ENDDO
      ENDDO

#ifdef ALLOW_NONHYDROSTATIC
      IF (nonHydrostatic) THEN
       DO j=1,sNy
        DO i=1,sNx
         cg3d_b(i,j,k,bi,bj) = pf(i+1,j)-pf(i,j)
        ENDDO
       ENDDO
      ENDIF
#endif

      DO j=1,sNy+1
       DO i=1,sNx
        pf(i,j) = yA(i,j)*gVNm1(i,j,k,bi,bj) / deltatmom
       ENDDO
      ENDDO

      DO j=1,sNy
       DO i=1,sNx
        cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj) +
     &   pf(i,j+1)-pf(i,j)
       ENDDO
      ENDDO

#ifdef ALLOW_NONHYDROSTATIC
      IF (nonHydrostatic) THEN
       DO j=1,sNy
        DO i=1,sNx
         cg3d_b(i,j,k,bi,bj) = cg3d_b(i,j,k,bi,bj) +
     &    pf(i,j+1)-pf(i,j)
        ENDDO
       ENDDO
      ENDIF
#endif

      RETURN
      END
1	adcroft	1.10	C $Header: /u/gcmpack/models/MITgcmUV/model/src/calc_div_ghat.F,v 1.9 1999/05/18 17:42:01 adcroft Exp $
2	cnh	1.1
3	cnh	1.6	#include "CPP_OPTIONS.h"
4	cnh	1.1
5			C /==========================================================\
6			C \| S/R CALC_DIV_GHAT \|
7			C \| o Form the right hand-side of the surface pressure eqn. \|
8	adcroft	1.9	C \|==========================================================\|
9	cnh	1.1	C \==========================================================/
10			SUBROUTINE CALC_DIV_GHAT(
11			I bi,bj,iMin,iMax,jMin,jMax,
12			I K,
13			I xA,yA,
14			I myThid)
15
16			IMPLICIT NONE
17
18			C == Global variables ==
19			#include "SIZE.h"
20			#include "DYNVARS.h"
21			#include "FFIELDS.h"
22			#include "EEPARAMS.h"
23			#include "PARAMS.h"
24			#include "GRID.h"
25			#include "CG2D.h"
26	adcroft	1.8	#ifdef ALLOW_NONHYDROSTATIC
27			#include "CG3D.h"
28			#endif
29	cnh	1.1
30			C == Routine arguments ==
31			C bi, bj, iMin, iMax, jMin, jMax - Range of points for which calculation
32			C results will be set.
33	adcroft	1.9	C k - Index of layer.
34			C xA, yA - Cell face areas
35	cnh	1.1	C myThid - Instance number for this innvocation of CALC_MOM_RHS
36			INTEGER bi,bj,iMin,iMax,jMin,jMax
37			INTEGER K
38			_RS xA (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
39			_RS yA (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
40			INTEGER myThid
41
42			C == Local variables ==
43			INTEGER i,j
44			_RL pf (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
45
46			C-- Pressure equation source term
47			C Term is the vertical integral of the divergence of the
48			C time tendency terms along with a relaxation term that
49			C pulls div(U) + dh/dt back toward zero.
50
51	cnh	1.2	IF ( k .EQ. Nr ) THEN
52	cnh	1.1	C Initialise source term on first pass
53	adcroft	1.10	DO j=1,sNy
54			DO i=1,sNx
55			C Note: The source term containing cg2d_x and cg3d_x (at k=1)
56			C has been moved to solve_for_pressure.F for convenience.
57	cnh	1.1	cg2d_b(i,j,bi,bj) =
58	adcroft	1.10	& freeSurfFac_rA(i,j,bi,bj)horiVertRatio*( 0.
59	adcroft	1.7	#ifdef USE_NATURAL_BCS
60			& +EmPmR(I,J,bi,bj)/deltaTMom
61			#endif
62			& )
63	cnh	1.1	ENDDO
64			ENDDO
65			ENDIF
66
67	adcroft	1.10	DO j=1,sNy
68			DO i=1,sNx+1
69	cnh	1.1	pf(i,j) = xA(i,j)*gUNm1(i,j,k,bi,bj) / deltaTmom
70			ENDDO
71			ENDDO
72			DO j=1,sNy
73			DO i=1,sNx
74	adcroft	1.8	cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj) +
75	cnh	1.1	& pf(i+1,j)-pf(i,j)
76			ENDDO
77			ENDDO
78
79	adcroft	1.8	#ifdef ALLOW_NONHYDROSTATIC
80			IF (nonHydrostatic) THEN
81			DO j=1,sNy
82			DO i=1,sNx
83			cg3d_b(i,j,k,bi,bj) = pf(i+1,j)-pf(i,j)
84			ENDDO
85			ENDDO
86			ENDIF
87			#endif
88
89	adcroft	1.10	DO j=1,sNy+1
90			DO i=1,sNx
91	cnh	1.1	pf(i,j) = yA(i,j)*gVNm1(i,j,k,bi,bj) / deltatmom
92			ENDDO
93			ENDDO
94
95			DO j=1,sNy
96			DO i=1,sNx
97	adcroft	1.8	cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj) +
98	cnh	1.1	& pf(i,j+1)-pf(i,j)
99			ENDDO
100			ENDDO
101	cnh	1.4
102	adcroft	1.8	#ifdef ALLOW_NONHYDROSTATIC
103			IF (nonHydrostatic) THEN
104			DO j=1,sNy
105			DO i=1,sNx
106			cg3d_b(i,j,k,bi,bj) = cg3d_b(i,j,k,bi,bj) +
107			& pf(i,j+1)-pf(i,j)
108			ENDDO
109			ENDDO
110			ENDIF
111			#endif
112	cnh	1.1
113			RETURN
114			END