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

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

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

revision 1.6 by baylor, Mon Sep 26 15:27:11 2005 UTC revision 1.7 by jmc, Tue Dec 5 05:30:38 2006 UTC
# Line 12  C !INTERFACE: ========================== Line 12  C !INTERFACE: ==========================
12       I        uFld, del2u,       I        uFld, del2u,
13       O        xViscFluxU,       O        xViscFluxU,
14       I        viscAh_D,viscA4_D,       I        viscAh_D,viscA4_D,
15       I        myThid)       I        myThid )
16    
17  C !DESCRIPTION:  C !DESCRIPTION:
18  C Calculates the area integrated zonal viscous fluxes of U:  C Calculates the area integrated zonal viscous fluxes of U:
19  C \begin{equation*}  C \begin{equation*}
20  C F^x = - \frac{ \Delta y_f \Delta r_f h_c }{\Delta x_f}  C F^x = - \frac{ \Delta y_f \Delta r_f h_c }{\Delta x_f}
21  C  ( A_h \delta_i u - A_4 \delta_i \nabla^2 u )  C  ( A_h \delta_i u - A_4 \delta_i \nabla^2 u )
22  C \end{equation*}  C \end{equation*}
23    
# Line 56  C     - Laplacian  and bi-harmonic terms Line 56  C     - Laplacian  and bi-harmonic terms
56          xViscFluxU(i,j) =          xViscFluxU(i,j) =
57       &    _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)
58       &     *(       &     *(
59       &       -viscAh_D(i,j)*(uFld(i+1,j)-uFld(i,j))       &       -viscAh_D(i,j)*( uFld(i+1,j)-uFld(i,j) )
60       &       *cosFacU(J,bi,bj)       &       *cosFacU(J,bi,bj)
61       &       +viscA4_D(i,j)*(del2u(i+1,j)       -del2u(i,j)       )       &       +viscA4_D(i,j)*(del2u(i+1,j)-del2u(i,j))
62  #ifdef COSINEMETH_III  #ifdef COSINEMETH_III
63       &       *sqcosFacU(J,bi,bj)       &       *sqCosFacU(J,bi,bj)
64  #else  #else
65       &       *cosFacU(J,bi,bj)       &       *cosFacU(J,bi,bj)
66  #endif  #endif
67       &      )*_recip_dxF(i,j,bi,bj)       &      )*_recip_dxF(i,j,bi,bj)
68    c    &       *deepFacC(k)        ! dyF scaling factor
69    c    &       *recip_deepFacC(k)  ! recip_dxF scaling factor
70         ENDDO         ENDDO
71        ENDDO        ENDDO
72    
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22