pkg/mom_common/mom_quasihydrostatic.F

C $Header: /u/gcmpack/MITgcm/pkg/mom_common/mom_quasihydrostatic.F,v 1.2 2006/07/13 03:02:48 jmc Exp $
C $Name:  $

#include "MOM_COMMON_OPTIONS.h"

CBOP
C !ROUTINE: MOM_QUASIHYDROSTATIC

C !INTERFACE: ==========================================================
      SUBROUTINE MOM_QUASIHYDROSTATIC(
     I                bi,bj,k,
     I                uFld, vFld,
     O                effectiveBuoy,
     I                myThid )

C !DESCRIPTION:
C     *==========================================================*
C     | o SUBROUTINE MOM_QUASIHYDROSTATIC
C     |   Add Quasi-Hydrostatic Terms to buoyancy
C     *==========================================================*

C !USES: ===============================================================
      IMPLICIT NONE

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

C--   == Routine arguments ==
C !INPUT VARIABLES: ====================================================
C  bi,bj         :: tile indices
C  k             :: vertical level
C  uFld          :: zonal flow
C  vFld          :: meridional flow
C  myThid        :: my Thread Id number
      INTEGER bi,bj,k
      _RL uFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
      _RL vFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
      INTEGER myThid

C !OUTPUT/MODIFIED VARIABLES: ==========================================
C  effectiveBuoy :: Density anomaly
      _RL effectiveBuoy(1-Olx:sNx+Olx,1-Oly:sNy+Oly)

C !LOCAL VARIABLES: ====================================================
C  i,j           :: loop indices
C  scalingFactor :: scaling factor (from acceleration to density)
      INTEGER i,j
      _RL scalingFactor
CEOP

      scalingFactor=rhoConst*gravitySign*recip_gravity
C-    to reproduce old (wrong) results:
c     scalingFactor=gravitySign*recip_gravity

      IF ( use3dCoriolis ) THEN
       DO j=1-Oly,sNy+Oly-1
        DO i=1-Olx,sNx+Olx-1
         effectiveBuoy(i,j)=effectiveBuoy(i,j)
     &    +scalingFactor*
     &     fCoriCos(i,j,bi,bj)*
     &       ( angleCosC(i,j,bi,bj)*0.5 _d 0 *
     &                (uFld(i,j,k,bi,bj)+uFld(i+1,j,k,bi,bj))
     &        -angleSinC(i,j,bi,bj)*0.5 _d 0 *
     &                (vFld(i,j,k,bi,bj)+vFld(i,j+1,k,bi,bj))
     &       )
        ENDDO
       ENDDO
      ENDIF

      IF ( useNHMTerms ) THEN
       DO j=1-Oly,sNy+Oly-1
        DO i=1-Olx,sNx+Olx-1
         effectiveBuoy(i,j)=effectiveBuoy(i,j)
     &    +scalingFactor*
     &     (   (uFld( i ,j,k,bi,bj)*uFld( i ,j,k,bi,bj)
     &         +uFld(i+1,j,k,bi,bj)*uFld(i+1,j,k,bi,bj))
     &       + (vFld(i, j ,k,bi,bj)*vFld(i, j ,k,bi,bj)
     &         +vFld(i,j+1,k,bi,bj)*vFld(i,j+1,k,bi,bj))
     &     )* 0.5 _d 0 * recip_rSphere*recip_deepFacC(k)
        ENDDO
       ENDDO
      ENDIF

      RETURN
      END
1	C $Header: /u/gcmpack/MITgcm/pkg/mom_common/mom_quasihydrostatic.F,v 1.2 2006/07/13 03:02:48 jmc Exp $
2	C $Name: $
3
4	#include "MOM_COMMON_OPTIONS.h"
5
6	CBOP
7	C !ROUTINE: MOM_QUASIHYDROSTATIC
8
9	C !INTERFACE: ==========================================================
10	SUBROUTINE MOM_QUASIHYDROSTATIC(
11	I bi,bj,k,
12	I uFld, vFld,
13	O effectiveBuoy,
14	I myThid )
15
16	C !DESCRIPTION:
17	C ==========================================================
18	C \| o SUBROUTINE MOM_QUASIHYDROSTATIC
19	C \| Add Quasi-Hydrostatic Terms to buoyancy
20	C ==========================================================
21
22	C !USES: ===============================================================
23	IMPLICIT NONE
24
25	C-- == Global data ==
26	#include "SIZE.h"
27	#include "EEPARAMS.h"
28	#include "PARAMS.h"
29	#include "GRID.h"
30
31	C-- == Routine arguments ==
32	C !INPUT VARIABLES: ====================================================
33	C bi,bj :: tile indices
34	C k :: vertical level
35	C uFld :: zonal flow
36	C vFld :: meridional flow
37	C myThid :: my Thread Id number
38	INTEGER bi,bj,k
39	_RL uFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
40	_RL vFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
41	INTEGER myThid
42
43	C !OUTPUT/MODIFIED VARIABLES: ==========================================
44	C effectiveBuoy :: Density anomaly
45	_RL effectiveBuoy(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
46
47	C !LOCAL VARIABLES: ====================================================
48	C i,j :: loop indices
49	C scalingFactor :: scaling factor (from acceleration to density)
50	INTEGER i,j
51	_RL scalingFactor
52	CEOP
53
54	scalingFactor=rhoConstgravitySignrecip_gravity
55	C- to reproduce old (wrong) results:
56	c scalingFactor=gravitySign*recip_gravity
57
58	IF ( use3dCoriolis ) THEN
59	DO j=1-Oly,sNy+Oly-1
60	DO i=1-Olx,sNx+Olx-1
61	effectiveBuoy(i,j)=effectiveBuoy(i,j)
62	& +scalingFactor*
63	& fCoriCos(i,j,bi,bj)*
64	& ( angleCosC(i,j,bi,bj)0.5 _d 0
65	& (uFld(i,j,k,bi,bj)+uFld(i+1,j,k,bi,bj))
66	& -angleSinC(i,j,bi,bj)0.5 _d 0
67	& (vFld(i,j,k,bi,bj)+vFld(i,j+1,k,bi,bj))
68	& )
69	ENDDO
70	ENDDO
71	ENDIF
72
73	IF ( useNHMTerms ) THEN
74	DO j=1-Oly,sNy+Oly-1
75	DO i=1-Olx,sNx+Olx-1
76	effectiveBuoy(i,j)=effectiveBuoy(i,j)
77	& +scalingFactor*
78	& ( (uFld( i ,j,k,bi,bj)*uFld( i ,j,k,bi,bj)
79	& +uFld(i+1,j,k,bi,bj)*uFld(i+1,j,k,bi,bj))
80	& + (vFld(i, j ,k,bi,bj)*vFld(i, j ,k,bi,bj)
81	& +vFld(i,j+1,k,bi,bj)*vFld(i,j+1,k,bi,bj))
82	& )* 0.5 _d 0 * recip_rSphere*recip_deepFacC(k)
83	ENDDO
84	ENDDO
85	ENDIF
86
87	RETURN
88	END