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

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

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

revision 1.3 by adcroft, Wed Sep 26 19:05:21 2001 UTC revision 1.3.6.1 by heimbach, Fri Mar 7 04:46:40 2003 UTC
# Line 46  C  i,j                  :: loop indices Line 46  C  i,j                  :: loop indices
46  CEOP  CEOP
47    
48  C     Zonal flux d/dx V  C     Zonal flux d/dx V
49        DO j=1-Oly,sNy+Oly        DO j=1-Oly+1,sNy+Oly-1
50         DO i=1-Olx+1,sNx+Olx         DO i=1-Olx+1,sNx+Olx
51          fZon(i,j) = drF(k)*hFacZ(i,j)          fZon(i,j) = drF(k)*hFacZ(i,j)
52       &   *_dyU(i,j,bi,bj)       &   *_dyU(i,j,bi,bj)
# Line 60  C     Zonal flux d/dx V Line 60  C     Zonal flux d/dx V
60    
61  C     Meridional flux d/dy V  C     Meridional flux d/dy V
62        DO j=1-Oly,sNy+Oly-1        DO j=1-Oly,sNy+Oly-1
63         DO i=1-Olx,sNx+Olx         DO i=1-Olx+1,sNx+Olx-1
64          fMer(i,j) = drF(k)*_hFacC(i,j,k,bi,bj)          fMer(i,j) = drF(k)*_hFacC(i,j,k,bi,bj)
65       &   *_dxF(i,j,bi,bj)       &   *_dxF(i,j,bi,bj)
66       &   *_recip_dyF(i,j,bi,bj)       &   *_recip_dyF(i,j,bi,bj)
# Line 68  C     Meridional flux d/dy V Line 68  C     Meridional flux d/dy V
68         ENDDO         ENDDO
69        ENDDO        ENDDO
70    
71  C     del^2 U  C     del^2 V
72        DO j=0,sNy+2        DO j=1-Oly+1,sNy+Oly-1
73         DO i=0,sNx+1         DO i=1-Olx+1,sNx+Olx-1
74          del2v(i,j) = recip_drF(k)*_recip_hFacS(i,j,k,bi,bj)          del2v(i,j) = recip_drF(k)*_recip_hFacS(i,j,k,bi,bj)
75       &  *recip_rAs(i,j,bi,bj)       &  *recip_rAs(i,j,bi,bj)
76       &  *( fZon(i+1,j)    - fZon(i, j )       &  *( fZon(i+1,j)    - fZon(i, j )
# Line 81  C     del^2 U Line 81  C     del^2 U
81    
82        IF (no_slip_sides) THEN        IF (no_slip_sides) THEN
83  C-- No-slip BCs impose a drag at walls...  C-- No-slip BCs impose a drag at walls...
84        DO j=0,sNy+2        DO j=1-Oly+1,sNy+Oly-1
85         DO i=0,sNx+1         DO i=1-Olx+1,sNx+Olx-1
86          hFacZClosedW = _hFacS(i,j,k,bi,bj) - hFacZ(i,j)          hFacZClosedW = _hFacS(i,j,k,bi,bj) - hFacZ(i,j)
87          hFacZClosedE = _hFacS(i,j,k,bi,bj) - hFacZ(i+1,j)          hFacZClosedE = _hFacS(i,j,k,bi,bj) - hFacZ(i+1,j)
88          del2v(i,j) = del2v(i,j)          del2v(i,j) = del2v(i,j)
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.3.6.1
  ViewVC Help
Powered by ViewVC 1.1.22