pkg/mom_fluxform/mom_u_adv_wu.F

C $Header: /u/gcmpack/MITgcm/pkg/mom_fluxform/mom_u_adv_wu.F,v 1.3 2001/09/26 19:05:21 adcroft Exp $
C $Name:  $

#include "CPP_OPTIONS.h"

CBOP
C !ROUTINE: MOM_U_ADV_WU

C !INTERFACE: ==========================================================
      SUBROUTINE MOM_U_ADV_WU(
     I        bi,bj,k,
     I        uFld,wFld,
     O        advectiveFluxWU,
     I        myThid)

C !DESCRIPTION:
C Calculates the vertical advective flux of zonal momentum:
C \begin{equation*}
C F^r = \overline{W}^i \overline{u}^{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  uFld                 :: zonal flow
C  myThid               :: thread number
      INTEGER bi,bj,k
      _RL uFld(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  advectiveFluxWU      :: advective flux
      _RL advectiveFluxWU(1-OLx:sNx+OLx,1-OLy:sNy+OLy)

C !LOCAL VARIABLES: ====================================================
C  i,j                  :: loop indices
      INTEGER I,J
CEOP

#ifdef NONLIN_FRSURF
      IF ( k.EQ.Nr+1 .AND.
     &     useRealFreshWaterFlux .AND.
     &     buoyancyRelation .EQ. 'OCEANICP' ) THEN 
       DO j=1-Oly,sNy+Oly
        DO i=1-Olx+1,sNx+Olx
         advectiveFluxWU(i,j) = convertEmP2rUnit*
     &      ( PmEpR( i ,j,bi,bj)*rA( i ,j,bi,bj)
     &       +PmEpR(i-1,j,bi,bj)*rA(i-1,j,bi,bj)
     &      )*0.5 _d 0 * uFld(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
        advectiveFluxWU(i,j) = 0.
       ENDDO
      ENDDO

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

      DO j=1-Oly,sNy+Oly
       DO i=1-Olx+1,sNx+Olx
        advectiveFluxWU(i,j) =
     &     0.5*(
     &        wFld( i ,j,k,bi,bj)*rA( i ,j,bi,bj)
     &       +wFld(i-1,j,k,bi,bj)*rA(i-1,j,bi,bj)
     &         )*uFld(i,j,k,bi,bj)
       ENDDO
      ENDDO

      ELSE

C     Vertical advection - interior ; assume uFld & wFld are masked  
      DO j=1-Oly,sNy+Oly
       DO i=1-Olx+1,sNx+Olx
        advectiveFluxWU(i,j) =
     &    0.25*(
     &        wFld( i ,j,k,bi,bj)*rA( i ,j,bi,bj)
     &       +wFld(i-1,j,k,bi,bj)*rA(i-1,j,bi,bj)
     &         )*( uFld(i,j,k,bi,bj)+uFld(i,j,k-1,bi,bj) )
       ENDDO
      ENDDO

      IF (.NOT. rigidLid) THEN
C     (linear) Free-surface correction at k>1
        DO j=1-Oly,sNy+Oly
         DO i=1-Olx+1,sNx+Olx
          advectiveFluxWU(i,j) = advectiveFluxWU(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-1,j,k,bi,bj)*rA(i-1,j,bi,bj)*
     &          (maskC(i-1,j,k,bi,bj)-maskC(i-1,j,k-1,bi,bj))
     &           )*uFld(i,j,k,bi,bj)
         ENDDO
        ENDDO
C- endif NOT rigidLid
      ENDIF

      ENDIF

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