/[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.3 by adcroft, Wed Sep 26 19:05:21 2001 UTC revision 1.7 by jmc, Tue Dec 5 05:30:38 2006 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
4  #include "CPP_OPTIONS.h"  #include "MOM_FLUXFORM_OPTIONS.h"
5    
6  CBOP  CBOP
7  C !ROUTINE: MOM_U_XVISCFLUX  C !ROUTINE: MOM_U_XVISCFLUX
# Line 11  C !INTERFACE: ========================== Line 11  C !INTERFACE: ==========================
11       I        bi,bj,k,       I        bi,bj,k,
12       I        uFld, del2u,       I        uFld, del2u,
13       O        xViscFluxU,       O        xViscFluxU,
14       I        myThid)       I        viscAh_D,viscA4_D,
15         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 36  C  myThid               :: thread number Line 37  C  myThid               :: thread number
37        INTEGER bi,bj,k        INTEGER bi,bj,k
38        _RL uFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL uFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
39        _RL del2u(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL del2u(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
40          _RL viscAh_D(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
41          _RL viscA4_D(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
42        INTEGER myThid        INTEGER myThid
43    
44  C !OUTPUT PARAMETERS: ==================================================  C !OUTPUT PARAMETERS: ==================================================
# Line 53  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*(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*(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.3  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22