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

Diff of /MITgcm_contrib/natl_12/code/mom_u_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 46  C !LOCAL VARIABLES: ==================== Line 46  C !LOCAL VARIABLES: ====================
46  C  i,j                  :: loop indices  C  i,j                  :: loop indices
47        INTEGER I,J        INTEGER I,J
48        _RS gridScalingA4        _RS gridScalingA4
49          _RS gridScalingAh
50  CEOP  CEOP
51    
52  C     - Laplacian  and bi-harmonic terms  C     - Laplacian  and bi-harmonic terms
53        DO j=1-Oly,sNy+Oly-1        DO j=1-Oly,sNy+Oly-1
54         DO i=1-Olx,sNx+Olx-1         DO i=1-Olx,sNx+Olx-1
55          gridScalingA4 = (_dxF(i,j,bi,bj)**4)/(5000.**4)          gridScalingA4 = (_dxF(i,j,bi,bj)**4)/(5000.**4)
56            gridScalingAh = (_dxF(i,j,bi,bj)**2)/(5000.**2)
57          xViscFluxU(i,j) =          xViscFluxU(i,j) =
58       &    _dyF(i,j,bi,bj)*drF(k)*_hFacC(i,j,k,bi,bj)       &    _dyF(i,j,bi,bj)*drF(k)*_hFacC(i,j,k,bi,bj)
59       &     *(       &     *(
60       &       -viscAh*(uFld(i+1,j)-uFld(i,j))       &       -viscAh*(uFld(i+1,j)-uFld(i,j))*gridScalingAh
61       &       *cosFacU(J,bi,bj)       &       *cosFacU(J,bi,bj)
62       &       +viscA4*(del2u(i+1,j)       -del2u(i,j)       )       &       +viscA4*(del2u(i+1,j)       -del2u(i,j)       )
63       &              *gridScalingA4       &              *gridScalingA4
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22