pkg/mom_fluxform/mom_v_adv_wv.F

C $Header: /u/gcmpack/MITgcm/pkg/mom_fluxform/mom_v_adv_wv.F,v 1.4 2002/12/10 03:02:38 jmc Exp $
C $Name:  $

#include "CPP_OPTIONS.h"

CBOP
C !ROUTINE: MOM_V_ADV_WV

C !INTERFACE: ==========================================================
      SUBROUTINE MOM_V_ADV_WV(
     I        bi,bj,k,
     I        vFld,wFld,
     O        advectiveFluxWV,
     I        myThid)

C !DESCRIPTION:
C Calculates the vertical advective flux of meridional momentum:
C \begin{equation*}
C F^r = \overline{W}^j \overline{v}^{k}
C \end{equation*}

C !USES: ===============================================================
      IMPLICIT NONE
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#ifdef NONLIN_FRSURF
#include "SURFACE.h"
#endif

C !INPUT PARAMETERS: ===================================================
C  bi,bj                :: tile indices
C  k                    :: vertical level
C  rTrans               :: vertical transport
C  vFld                 :: meridional flow
C  myThid               :: thread number
      INTEGER bi,bj,k
      _RL vFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
      _RL wFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
      INTEGER myThid

C !OUTPUT PARAMETERS: ==================================================
C  advectiveFluxWV      :: advective flux
      _RL advectiveFluxWV(1-OLx:sNx+OLx,1-OLy:sNy+OLy)

C !LOCAL VARIABLES: ====================================================
C  i,j                  :: loop indices
      INTEGER I,J
      _RL tmpFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
CEOP

#ifdef NONLIN_FRSURF
      IF ( k.EQ.Nr+1 .AND.
     &     useRealFreshWaterFlux .AND.
     &     buoyancyRelation .EQ. 'OCEANICP' ) THEN
       DO j=1-Oly+1,sNy+Oly
        DO i=1-Olx,sNx+Olx
         advectiveFluxWV(i,j) = convertEmP2rUnit*
     &      ( PmEpR(i, j ,bi,bj)*rA(i, j ,bi,bj)
     &       +PmEpR(i,j-1,bi,bj)*rA(i,j-1,bi,bj)
     &      )*0.5 _d 0 * vFld(i,j,k-1,bi,bj)
        ENDDO
       ENDDO

      ELSEIF ( k.GT.Nr .OR. (k.EQ.1.AND.rigidLid) ) THEN
#else /* NONLIN_FRSURF */
      IF ( k.GT.Nr .OR.
     &    (k.EQ.1.AND.rigidLid) ) THEN
#endif /* NONLIN_FRSURF */
C     Advective flux = 0  at k=Nr+1 ; = 0 at k=1 if rigid-lid

      DO j=1-Oly,sNy+Oly
       DO i=1-Olx,sNx+Olx
        advectiveFluxWV(i,j) = 0.
       ENDDO
      ENDDO

      ELSEIF (k.EQ.1) THEN
C     (linear) Free-surface correction at k=1

      DO j=1-Oly+1,sNy+Oly
       DO i=1-Olx,sNx+Olx
c       advectiveFluxWV(i,j) =
c    &     0.5*(
c    &        wFld(i, j ,k,bi,bj)*rA(i, j ,bi,bj)
c    &       +wFld(i,j-1,k,bi,bj)*rA(i,j-1,bi,bj)
c    &         )*vFld(i,j,k,bi,bj)
        tmpFld(i,j) = 0.5 _d 0*(
     &        wFld(i, j ,k,bi,bj)*rA(i, j ,bi,bj)
     &       +wFld(i,j-1,k,bi,bj)*rA(i,j-1,bi,bj) )
        advectiveFluxWV(i,j) = tmpFld(i,j)*vFld(i,j,k,bi,bj)
       ENDDO
      ENDDO

      ELSE

C     Vertical advection - interior ; assume vFld & wFld are masked 
      DO j=1-Oly+1,sNy+Oly
       DO i=1-Olx,sNx+Olx
c       advectiveFluxWV(i,j) =
c    &    0.25*(
c    &        wFld(i, j ,k,bi,bj)*rA(i, j ,bi,bj)
c    &       +wFld(i,j-1,k,bi,bj)*rA(i,j-1,bi,bj)
c    &         )*( vFld(i,j,k,bi,bj)+vFld(i,j,k-1,bi,bj) )
        tmpFld(i,j) = 0.5 _d 0*(
     &        wFld(i, j ,k,bi,bj)*rA(i, j ,bi,bj)
     &       +wFld(i,j-1,k,bi,bj)*rA(i,j-1,bi,bj) )
        advectiveFluxWV(i,j) = tmpFld(i,j)*0.5 _d 0
     &          *( vFld(i,j,k,bi,bj)+vFld(i,j,k-1,bi,bj) )
       ENDDO
      ENDDO

      IF (.NOT. rigidLid) THEN
C     (linear) Free-surface correction at k>1
        DO j=1-Oly+1,sNy+Oly
         DO i=1-Olx,sNx+Olx
          advectiveFluxWV(i,j) = advectiveFluxWV(i,j)
     &     +0.25*(
     &          wFld(i, j ,k,bi,bj)*rA(i, j ,bi,bj)*
     &          (maskC(i,j,k,bi,bj) - maskC(i,j,k-1,bi,bj))
     &         +wFld(i,j-1,k,bi,bj)*rA(i,j-1,bi,bj)*
     &          (maskC(i,j-1,k,bi,bj)-maskC(i,j-1,k-1,bi,bj))
     &           )*vFld(i,j,k,bi,bj)
         ENDDO
        ENDDO
C- endif NOT rigidLid
      ENDIF

      ENDIF

      RETURN
      END
1	C $Header: /u/gcmpack/MITgcm/pkg/mom_fluxform/mom_v_adv_wv.F,v 1.4 2002/12/10 03:02:38 jmc Exp $
2	C $Name: $
3
4	#include "CPP_OPTIONS.h"
5
6	CBOP
7	C !ROUTINE: MOM_V_ADV_WV
8
9	C !INTERFACE: ==========================================================
10	SUBROUTINE MOM_V_ADV_WV(
11	I bi,bj,k,
12	I vFld,wFld,
13	O advectiveFluxWV,
14	I myThid)
15
16	C !DESCRIPTION:
17	C Calculates the vertical advective flux of meridional momentum:
18	C \begin{equation*}
19	C F^r = \overline{W}^j \overline{v}^{k}
20	C \end{equation*}
21
22	C !USES: ===============================================================
23	IMPLICIT NONE
24	#include "SIZE.h"
25	#include "EEPARAMS.h"
26	#include "PARAMS.h"
27	#include "GRID.h"
28	#ifdef NONLIN_FRSURF
29	#include "SURFACE.h"
30	#endif
31
32	C !INPUT PARAMETERS: ===================================================
33	C bi,bj :: tile indices
34	C k :: vertical level
35	C rTrans :: vertical transport
36	C vFld :: meridional flow
37	C myThid :: thread number
38	INTEGER bi,bj,k
39	_RL vFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
40	_RL wFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
41	INTEGER myThid
42
43	C !OUTPUT PARAMETERS: ==================================================
44	C advectiveFluxWV :: advective flux
45	_RL advectiveFluxWV(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
46
47	C !LOCAL VARIABLES: ====================================================
48	C i,j :: loop indices
49	INTEGER I,J
50	_RL tmpFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
51	CEOP
52
53	#ifdef NONLIN_FRSURF
54	IF ( k.EQ.Nr+1 .AND.
55	& useRealFreshWaterFlux .AND.
56	& buoyancyRelation .EQ. 'OCEANICP' ) THEN
57	DO j=1-Oly+1,sNy+Oly
58	DO i=1-Olx,sNx+Olx
59	advectiveFluxWV(i,j) = convertEmP2rUnit*
60	& ( PmEpR(i, j ,bi,bj)*rA(i, j ,bi,bj)
61	& +PmEpR(i,j-1,bi,bj)*rA(i,j-1,bi,bj)
62	& )0.5 _d 0 vFld(i,j,k-1,bi,bj)
63	ENDDO
64	ENDDO
65
66	ELSEIF ( k.GT.Nr .OR. (k.EQ.1.AND.rigidLid) ) THEN
67	#else /* NONLIN_FRSURF */
68	IF ( k.GT.Nr .OR.
69	& (k.EQ.1.AND.rigidLid) ) THEN
70	#endif /* NONLIN_FRSURF */
71	C Advective flux = 0 at k=Nr+1 ; = 0 at k=1 if rigid-lid
72
73	DO j=1-Oly,sNy+Oly
74	DO i=1-Olx,sNx+Olx
75	advectiveFluxWV(i,j) = 0.
76	ENDDO
77	ENDDO
78
79	ELSEIF (k.EQ.1) THEN
80	C (linear) Free-surface correction at k=1
81
82	DO j=1-Oly+1,sNy+Oly
83	DO i=1-Olx,sNx+Olx
84	c advectiveFluxWV(i,j) =
85	c & 0.5*(
86	c & wFld(i, j ,k,bi,bj)*rA(i, j ,bi,bj)
87	c & +wFld(i,j-1,k,bi,bj)*rA(i,j-1,bi,bj)
88	c & )*vFld(i,j,k,bi,bj)
89	tmpFld(i,j) = 0.5 _d 0*(
90	& wFld(i, j ,k,bi,bj)*rA(i, j ,bi,bj)
91	& +wFld(i,j-1,k,bi,bj)*rA(i,j-1,bi,bj) )
92	advectiveFluxWV(i,j) = tmpFld(i,j)*vFld(i,j,k,bi,bj)
93	ENDDO
94	ENDDO
95
96	ELSE
97
98	C Vertical advection - interior ; assume vFld & wFld are masked
99	DO j=1-Oly+1,sNy+Oly
100	DO i=1-Olx,sNx+Olx
101	c advectiveFluxWV(i,j) =
102	c & 0.25*(
103	c & wFld(i, j ,k,bi,bj)*rA(i, j ,bi,bj)
104	c & +wFld(i,j-1,k,bi,bj)*rA(i,j-1,bi,bj)
105	c & )*( vFld(i,j,k,bi,bj)+vFld(i,j,k-1,bi,bj) )
106	tmpFld(i,j) = 0.5 _d 0*(
107	& wFld(i, j ,k,bi,bj)*rA(i, j ,bi,bj)
108	& +wFld(i,j-1,k,bi,bj)*rA(i,j-1,bi,bj) )
109	advectiveFluxWV(i,j) = tmpFld(i,j)*0.5 _d 0
110	& *( vFld(i,j,k,bi,bj)+vFld(i,j,k-1,bi,bj) )
111	ENDDO
112	ENDDO
113
114	IF (.NOT. rigidLid) THEN
115	C (linear) Free-surface correction at k>1
116	DO j=1-Oly+1,sNy+Oly
117	DO i=1-Olx,sNx+Olx
118	advectiveFluxWV(i,j) = advectiveFluxWV(i,j)
119	& +0.25*(
120	& wFld(i, j ,k,bi,bj)rA(i, j ,bi,bj)
121	& (maskC(i,j,k,bi,bj) - maskC(i,j,k-1,bi,bj))
122	& +wFld(i,j-1,k,bi,bj)rA(i,j-1,bi,bj)
123	& (maskC(i,j-1,k,bi,bj)-maskC(i,j-1,k-1,bi,bj))
124	& )*vFld(i,j,k,bi,bj)
125	ENDDO
126	ENDDO
127	C- endif NOT rigidLid
128	ENDIF
129
130	ENDIF
131
132	RETURN
133	END