--- MITgcm/pkg/mom_fluxform/mom_u_adv_wu.F 2002/12/10 03:02:38 1.4 +++ MITgcm/pkg/mom_fluxform/mom_u_adv_wu.F 2003/01/24 17:15:59 1.5 @@ -1,4 +1,4 @@ -C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/mom_fluxform/mom_u_adv_wu.F,v 1.4 2002/12/10 03:02:38 jmc Exp $ +C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/mom_fluxform/mom_u_adv_wu.F,v 1.5 2003/01/24 17:15:59 jmc Exp $ C $Name: $ #include "CPP_OPTIONS.h" @@ -47,6 +47,7 @@ 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 @@ -80,11 +81,15 @@ DO j=1-Oly,sNy+Oly DO i=1-Olx+1,sNx+Olx - advectiveFluxWU(i,j) = - & 0.5*( +c advectiveFluxWU(i,j) = +c & 0.5*( +c & wFld( i ,j,k,bi,bj)*rA( i ,j,bi,bj) +c & +wFld(i-1,j,k,bi,bj)*rA(i-1,j,bi,bj) +c & )*uFld(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-1,j,k,bi,bj)*rA(i-1,j,bi,bj) - & )*uFld(i,j,k,bi,bj) + & +wFld(i-1,j,k,bi,bj)*rA(i-1,j,bi,bj) ) + advectiveFluxWU(i,j) = tmpFld(i,j)*uFld(i,j,k,bi,bj) ENDDO ENDDO @@ -93,11 +98,16 @@ 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*( +c advectiveFluxWU(i,j) = +c & 0.25*( +c & wFld( i ,j,k,bi,bj)*rA( i ,j,bi,bj) +c & +wFld(i-1,j,k,bi,bj)*rA(i-1,j,bi,bj) +c & )*( uFld(i,j,k,bi,bj)+uFld(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-1,j,k,bi,bj)*rA(i-1,j,bi,bj) - & )*( uFld(i,j,k,bi,bj)+uFld(i,j,k-1,bi,bj) ) + & +wFld(i-1,j,k,bi,bj)*rA(i-1,j,bi,bj) ) + advectiveFluxWU(i,j) = tmpFld(i,j)*0.5 _d 0 + & *( uFld(i,j,k,bi,bj)+uFld(i,j,k-1,bi,bj) ) ENDDO ENDDO