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

Annotation of /MITgcm_contrib/natl_12/code/mom_vi_u_xviscflux.F

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


Revision 1.1 - (hide annotations) (download)
Tue Aug 5 21:22:44 2003 UTC (21 years, 11 months ago) by cnh
Branch: MAIN
CVS Tags: HEAD
Adding set of files for 1/12 Atlantic configuration

1 cnh 1.1 C $Header: /u/u0/gcmpack/MITgcm/pkg/mom_vecinv/mom_vi_u_xviscflux.F,v 1.2 2001/05/29 14:01:39 adcroft Exp $
2     C $Name: $
3    
4     #include "CPP_OPTIONS.h"
5    
6     SUBROUTINE MOM_VI_U_XVISCFLUX(
7     I bi,bj,k,
8     I hDiv,del2u,
9     O xViscFluxU,
10     I myThid)
11     IMPLICIT NONE
12     C
13     C Calculate viscous flux in X direction for U eqn
14     C - is not vector invariant ...
15     C
16    
17     C == Global variables ==
18     #include "SIZE.h"
19     #include "EEPARAMS.h"
20     #include "PARAMS.h"
21     #include "GRID.h"
22    
23     C == Routine arguments ==
24     INTEGER bi,bj,k
25     _RL hDiv(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
26     _RL del2u(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
27     _RL xViscFluxU(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
28     INTEGER myThid
29    
30     C == Local variables ==
31     INTEGER I,J
32    
33     C - Laplacian and bi-harmonic terms
34     DO j=1-Oly,sNy+Oly-1
35     DO i=1-Olx,sNx+Olx-1
36     xViscFluxU(i,j) =
37     & _dyF(i,j,bi,bj)*drF(k)*_hFacC(i,j,k,bi,bj)
38     & *(-viscAh*hDiv(i,j)*cosFacU(J,bi,bj)
39     #ifdef COSINEMETH_III
40     & +viscA4*(del2u(i+1,j)-del2u(i,j))*sqcosFacU(J,bi,bj)
41     #else
42     & +viscA4*(del2u(i+1,j)-del2u(i,j))*cosFacU(J,bi,bj)
43     #endif
44     & *recip_dxf(i,j,bi,bj) )
45     ENDDO
46     ENDDO
47    
48     RETURN
49     END

  ViewVC Help
Powered by ViewVC 1.1.22