model/src/quasihydrostaticterms.F

C $Header: /u/gcmpack/MITgcm/model/src/quasihydrostaticterms.F,v 1.3 2005/04/06 18:20:22 jmc Exp $
C $Name:  $

#include "CPP_OPTIONS.h"

      SUBROUTINE QUASIHYDROSTATICTERMS(
     I      bi,bj,k,
     O      effectiveBuoy,
     I      myThid )

C     /==========================================================\
C     | o SUBROUTINE BUOYANCY                                    |
C     |   Calculates buoyancy for an XY slice of a tile.         |
C     |==========================================================|
C     |                                                          |
C     | k - is the density level                                 |
C     |                                                          |
C     \==========================================================/
      IMPLICIT NONE

C--   == Global data ==
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "DYNVARS.h"

C--   == Routine arguments ==
      INTEGER bi,bj,k
      _RL effectiveBuoy(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
      INTEGER myThid

C--   == Local variables ==
      INTEGER i,j
      _RL scalingFactor

      scalingFactor=rhoConst/gravity
      scalingFactor=1./gravity

      DO j=1-Oly,sNy+Oly
       DO i=1-Olx,sNx+Olx-1
        effectiveBuoy(i,j)=effectiveBuoy(i,j)
     &    -scalingFactor*
     &     fCoriCos(i,j,bi,bj)*
     &     0.5*( uVel(i,j,k,bi,bj)+uVel(i+1,j,k,bi,bj) )
       ENDDO
      ENDDO

      IF (useNHMTerms) THEN
       DO j=1-Oly,sNy+Oly-1
        DO i=1-Olx,sNx+Olx-1
         effectiveBuoy(i,j)=effectiveBuoy(i,j)
     &    -scalingFactor*
     &     (
     &       0.5*( uVel(i,j,k,bi,bj)**2+uVel(i+1,j,k,bi,bj)**2 )
     &      +0.5*( vVel(i,j,k,bi,bj)**2+vVel(i,j+1,k,bi,bj)**2 )
     &     )*recip_RSphere
        ENDDO
       ENDDO
      ENDIF

      RETURN
      END
1	C $Header: /u/gcmpack/MITgcm/model/src/quasihydrostaticterms.F,v 1.3 2005/04/06 18:20:22 jmc Exp $
2	C $Name: $
3
4	#include "CPP_OPTIONS.h"
5
6	SUBROUTINE QUASIHYDROSTATICTERMS(
7	I bi,bj,k,
8	O effectiveBuoy,
9	I myThid )
10
11	C /==========================================================\
12	C \| o SUBROUTINE BUOYANCY \|
13	C \| Calculates buoyancy for an XY slice of a tile. \|
14	C \|==========================================================\|
15	C \| \|
16	C \| k - is the density level \|
17	C \| \|
18	C \==========================================================/
19	IMPLICIT NONE
20
21	C-- == Global data ==
22	#include "SIZE.h"
23	#include "EEPARAMS.h"
24	#include "PARAMS.h"
25	#include "GRID.h"
26	#include "DYNVARS.h"
27
28	C-- == Routine arguments ==
29	INTEGER bi,bj,k
30	_RL effectiveBuoy(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
31	INTEGER myThid
32
33	C-- == Local variables ==
34	INTEGER i,j
35	_RL scalingFactor
36
37	scalingFactor=rhoConst/gravity
38	scalingFactor=1./gravity
39
40	DO j=1-Oly,sNy+Oly
41	DO i=1-Olx,sNx+Olx-1
42	effectiveBuoy(i,j)=effectiveBuoy(i,j)
43	& -scalingFactor*
44	& fCoriCos(i,j,bi,bj)*
45	& 0.5*( uVel(i,j,k,bi,bj)+uVel(i+1,j,k,bi,bj) )
46	ENDDO
47	ENDDO
48
49	IF (useNHMTerms) THEN
50	DO j=1-Oly,sNy+Oly-1
51	DO i=1-Olx,sNx+Olx-1
52	effectiveBuoy(i,j)=effectiveBuoy(i,j)
53	& -scalingFactor*
54	& (
55	& 0.5( uVel(i,j,k,bi,bj)2+uVel(i+1,j,k,bi,bj)*2 )
56	& +0.5( vVel(i,j,k,bi,bj)2+vVel(i,j+1,k,bi,bj)*2 )
57	& )*recip_RSphere
58	ENDDO
59	ENDDO
60	ENDIF
61
62	RETURN
63	END