/[MITgcm]/MITgcm/pkg/mom_fluxform/mom_u_adv_wu.F
ViewVC logotype

Diff of /MITgcm/pkg/mom_fluxform/mom_u_adv_wu.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.4 by jmc, Tue Dec 10 03:02:38 2002 UTC revision 1.5 by jmc, Fri Jan 24 17:15:59 2003 UTC
# Line 47  C  advectiveFluxWU      :: advective flu Line 47  C  advectiveFluxWU      :: advective flu
47  C !LOCAL VARIABLES: ====================================================  C !LOCAL VARIABLES: ====================================================
48  C  i,j                  :: loop indices  C  i,j                  :: loop indices
49        INTEGER I,J        INTEGER I,J
50          _RL tmpFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
51  CEOP  CEOP
52    
53  #ifdef NONLIN_FRSURF  #ifdef NONLIN_FRSURF
# Line 80  C     (linear) Free-surface correction a Line 81  C     (linear) Free-surface correction a
81    
82        DO j=1-Oly,sNy+Oly        DO j=1-Oly,sNy+Oly
83         DO i=1-Olx+1,sNx+Olx         DO i=1-Olx+1,sNx+Olx
84          advectiveFluxWU(i,j) =  c       advectiveFluxWU(i,j) =
85       &     0.5*(  c    &     0.5*(
86    c    &        wFld( i ,j,k,bi,bj)*rA( i ,j,bi,bj)
87    c    &       +wFld(i-1,j,k,bi,bj)*rA(i-1,j,bi,bj)
88    c    &         )*uFld(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)       &        wFld( i ,j,k,bi,bj)*rA( i ,j,bi,bj)
91       &       +wFld(i-1,j,k,bi,bj)*rA(i-1,j,bi,bj)       &       +wFld(i-1,j,k,bi,bj)*rA(i-1,j,bi,bj) )  
92       &         )*uFld(i,j,k,bi,bj)          advectiveFluxWU(i,j) = tmpFld(i,j)*uFld(i,j,k,bi,bj)
93         ENDDO         ENDDO
94        ENDDO        ENDDO
95    
# Line 93  C     (linear) Free-surface correction a Line 98  C     (linear) Free-surface correction a
98  C     Vertical advection - interior ; assume uFld & wFld are masked    C     Vertical advection - interior ; assume uFld & wFld are masked  
99        DO j=1-Oly,sNy+Oly        DO j=1-Oly,sNy+Oly
100         DO i=1-Olx+1,sNx+Olx         DO i=1-Olx+1,sNx+Olx
101          advectiveFluxWU(i,j) =  c       advectiveFluxWU(i,j) =
102       &    0.25*(  c    &    0.25*(
103    c    &        wFld( i ,j,k,bi,bj)*rA( i ,j,bi,bj)
104    c    &       +wFld(i-1,j,k,bi,bj)*rA(i-1,j,bi,bj)
105    c    &         )*( uFld(i,j,k,bi,bj)+uFld(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)       &        wFld( i ,j,k,bi,bj)*rA( i ,j,bi,bj)
108       &       +wFld(i-1,j,k,bi,bj)*rA(i-1,j,bi,bj)       &       +wFld(i-1,j,k,bi,bj)*rA(i-1,j,bi,bj) )
109       &         )*( uFld(i,j,k,bi,bj)+uFld(i,j,k-1,bi,bj) )          advectiveFluxWU(i,j) = tmpFld(i,j)*0.5 _d 0
110         &          *( uFld(i,j,k,bi,bj)+uFld(i,j,k-1,bi,bj) )
111         ENDDO         ENDDO
112        ENDDO        ENDDO
113    
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22