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

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

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

revision 1.4 by edhill, Thu Oct 9 04:19:20 2003 UTC revision 1.5 by baylor, Fri Sep 16 20:57:09 2005 UTC
# Line 11  C !INTERFACE: ========================== Line 11  C !INTERFACE: ==========================
11       I        bi,bj,k,       I        bi,bj,k,
12       I        vFld, del2v, hFacZ,       I        vFld, del2v, hFacZ,
13       O        xViscFluxV,       O        xViscFluxV,
14         I        viscAh_Z,viscAh_D,viscA4_Z,viscA4_D,
15       I        myThid)       I        myThid)
16    
17  C !DESCRIPTION:  C !DESCRIPTION:
# Line 37  C  myThid               :: thread number Line 38  C  myThid               :: thread number
38        _RL vFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL vFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
39        _RL del2v(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL del2v(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
40        _RS hFacZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RS hFacZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
41          _RL viscAh_Z(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
42          _RL viscAh_D(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
43          _RL viscA4_Z(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
44          _RL viscA4_D(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
45        INTEGER myThid        INTEGER myThid
46    
47  C !OUTPUT PARAMETERS: ==================================================  C !OUTPUT PARAMETERS: ==================================================
# Line 54  C     - Laplacian  and bi-harmonic terms Line 59  C     - Laplacian  and bi-harmonic terms
59          xViscFluxV(i,j) =          xViscFluxV(i,j) =
60       &     _dyU(i,j,bi,bj)*drF(k)*hFacZ(i,j)       &     _dyU(i,j,bi,bj)*drF(k)*hFacZ(i,j)
61       &     *(       &     *(
62       &       -viscAh*(vFld(i,j)-vFld(i-1,j))       &       -viscAh_Z(i,j)*(vFld(i,j)-vFld(i-1,j))
63       &       *cosFacV(J,bi,bj)       &       *cosFacV(J,bi,bj)
64       &       +viscA4*(del2v(i,j)       -del2v(i-1,j)       )       &       +viscA4_Z(i,j)*(del2v(i,j)       -del2v(i-1,j)       )
65  #ifdef COSINEMETH_III  #ifdef COSINEMETH_III
66       &      *sqcosFacV(J,bi,bj)       &      *sqcosFacV(J,bi,bj)
67  #else  #else
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22