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

Diff of /MITgcm_contrib/natl_12/code/mom_v_yviscflux.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 = (_dyF(i,j,bi,bj)**4)/(5000.**4)          gridScalingA4 = (_dyF(i,j,bi,bj)**4)/(5000.**4)
56            gridScalingAh = (_dyF(i,j,bi,bj)**2)/(5000.**2)
57          yViscFluxV(i,j) =          yViscFluxV(i,j) =
58       & _dxF(i,j,bi,bj)*drF(k)*_hFacC(i,j,k,bi,bj)       & _dxF(i,j,bi,bj)*drF(k)*_hFacC(i,j,k,bi,bj)
59       &     *(       &     *(
60       &       -viscAh*(vFld(i,j+1)-vFld(i,j))       &       -viscAh*(vFld(i,j+1)-vFld(i,j))*gridScalingAh
61  #ifdef ISOTROPIC_COS_SCALING  #ifdef ISOTROPIC_COS_SCALING
62       &       *cosFacU(J,bi,bj)       &       *cosFacU(J,bi,bj)
63  #endif  #endif
64       &       +viscA4*(del2v(i,j+1)       -del2v(i,j)       )       &       +viscA4*(del2v(i,j+1)       -del2v(i,j)       )
65         &              *gridScalingA4
66  #ifdef ISOTROPIC_COS_SCALING  #ifdef ISOTROPIC_COS_SCALING
67  #ifdef COSINEMETH_III  #ifdef COSINEMETH_III
68       &       *sqCosFacU(J,bi,bj)       &       *sqCosFacU(J,bi,bj)
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22