/[MITgcm]/MITgcm_contrib/natl_12/code/mom_v_xviscflux.F
ViewVC logotype

Diff of /MITgcm_contrib/natl_12/code/mom_v_xviscflux.F

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

revision 1.1 by cnh, Tue Aug 5 21:22:44 2003 UTC revision 1.2 by cnh, Thu Aug 7 13:03:31 2003 UTC
# Line 47  C !LOCAL VARIABLES: ==================== Line 47  C !LOCAL VARIABLES: ====================
47  C  i,j                  :: loop indices  C  i,j                  :: loop indices
48        INTEGER I,J        INTEGER I,J
49        _RS     gridScalingA4        _RS     gridScalingA4
50          _RS     gridScalingAh
51  CEOP  CEOP
52    
53  C     - Laplacian  and bi-harmonic terms  C     - Laplacian  and bi-harmonic terms
54        DO j=1-Oly,sNy+Oly        DO j=1-Oly,sNy+Oly
55         DO i=1-Olx+1,sNx+Olx         DO i=1-Olx+1,sNx+Olx
56          gridScalingA4 = (_dxV(i,j,bi,bj)**4)/(5000.**4)          gridScalingA4 = (_dxV(i,j,bi,bj)**4)/(5000.**4)
57            gridScalingAh = (_dxV(i,j,bi,bj)**2)/(5000.**2)
58          xViscFluxV(i,j) =          xViscFluxV(i,j) =
59       &     _dyU(i,j,bi,bj)*drF(k)*hFacZ(i,j)       &     _dyU(i,j,bi,bj)*drF(k)*hFacZ(i,j)
60       &     *(       &     *(
61       &       -viscAh*(vFld(i,j)-vFld(i-1,j))       &       -viscAh*(vFld(i,j)-vFld(i-1,j))*gridScalingAh
62       &       *cosFacV(J,bi,bj)       &       *cosFacV(J,bi,bj)
63       &       +viscA4*(del2v(i,j)       -del2v(i-1,j)       )       &       +viscA4*(del2v(i,j)       -del2v(i-1,j)       )
64       &              *gridScalingA4       &              *gridScalingA4
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22