/[MITgcm]/MITgcm/pkg/mom_common/mom_hdissip.F
ViewVC logotype

Diff of /MITgcm/pkg/mom_common/mom_hdissip.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 Feb 11 20:22:51 2014 UTC
# Line 4  C $Name$ Line 4  C $Name$
4  #include "MOM_COMMON_OPTIONS.h"  #include "MOM_COMMON_OPTIONS.h"
5    
6        SUBROUTINE MOM_HDISSIP(        SUBROUTINE MOM_HDISSIP(
7       I        bi,bj,k,       I        bi, bj, k,
8       I        hDiv,vort3,tension,strain,KE,       I        tension, strain, hFacZ,
9       I        hFacZ,       I        viscAh_s, viscAh_t, viscA4_s, viscA4_t,
10       I        viscAh_s,viscAh_t,viscA4_s,viscA4_t,       I        harmonic, biharmonic, useVariableViscosity,
11       I        harmonic,biharmonic,useVariableViscosity,       O        uDissip, vDissip,
12       O        uDissip,vDissip,       I        myThid )
13       I        myThid)  
14        IMPLICIT NONE        IMPLICIT NONE
15  C  
16  C     Calculate horizontal dissipation terms in terms of tension and strain  C     Calculate horizontal dissipation terms in terms of tension and strain
 C  
17  C       Du = d/dx At Tension + d/dy As Strain  C       Du = d/dx At Tension + d/dy As Strain
18  C       Dv = d/dx As Strain  - d/dy At Tension  C       Dv = d/dx As Strain  - d/dy At Tension
19    
20  C     == Global variables ==  C     == Global variables ==
21  #include "SIZE.h"  #include "SIZE.h"
22  #include "EEPARAMS.h"  #include "EEPARAMS.h"
 #include "GRID.h"  
23  #include "PARAMS.h"  #include "PARAMS.h"
24    #include "GRID.h"
25    
26  C     == Routine arguments ==  C     == Routine arguments ==
27        INTEGER bi,bj,k        INTEGER bi,bj,k
       _RL hDiv(1-OLx:sNx+OLx,1-OLy:sNy+OLy)  
       _RL vort3(1-OLx:sNx+OLx,1-OLy:sNy+OLy)  
28        _RL tension(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL tension(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
29        _RL strain(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL strain (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
30        _RL KE(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RS hFacZ  (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
       _RS hFacZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy)  
       _RL uDissip(1-OLx:sNx+OLx,1-OLy:sNy+OLy)  
       _RL vDissip(1-OLx:sNx+OLx,1-OLy:sNy+OLy)  
31        _RL viscAh_s(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL viscAh_s(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
32        _RL viscAh_t(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL viscAh_t(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
33        _RL viscA4_s(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL viscA4_s(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
34        _RL viscA4_t(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL viscA4_t(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
35        LOGICAL harmonic, biharmonic, useVariableViscosity        LOGICAL harmonic, biharmonic, useVariableViscosity
36          _RL uDissip(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
37          _RL vDissip(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
38        INTEGER myThid        INTEGER myThid
39    
40  C     == Local variables ==  C     == Local variables ==
41        INTEGER I,J        INTEGER i, j
42    
43  C     - Laplacian and bi-harmonic terms  C     - Laplacian and bi-harmonic terms
44        IF (harmonic) THEN        IF (harmonic) THEN
45         DO j=2-Oly,sNy+Oly-1         DO j=2-OLy,sNy+OLy-1
46          DO i=2-Olx,sNx+Olx-1          DO i=2-OLx,sNx+OLx-1
47    
48           uDissip(i,j) =           uDissip(i,j) =
49       &    recip_dyg(i,j,bi,bj)*recip_dyg(i,j,bi,bj)       &    recip_dyG(i,j,bi,bj)*recip_dyG(i,j,bi,bj)
50       &    *recip_dxc(i,j,bi,bj)       &    *recip_dxC(i,j,bi,bj)
51       &    *(       &    *(
52       &       dyf( i , j ,bi,bj)*dyf( i , j ,bi,bj)       &       dyF( i , j ,bi,bj)*dyF( i , j ,bi,bj)
53       &         *viscAh_t( i , j )*tension( i , j )       &         *viscAh_t( i , j )*tension( i , j )
54       &      -dyf(i-1, j ,bi,bj)*dyf(i-1, j ,bi,bj)       &      -dyF(i-1, j ,bi,bj)*dyF(i-1, j ,bi,bj)
55       &         *viscAh_t(i-1, j )*tension(i-1, j )       &         *viscAh_t(i-1, j )*tension(i-1, j )
56       &     )       &     )
57       &    +recip_dxc(i,j,bi,bj)*recip_dxc(i,j,bi,bj)       &    +recip_dxC(i,j,bi,bj)*recip_dxC(i,j,bi,bj)
58       &    *recip_dyg(i,j,bi,bj)       &    *recip_dyG(i,j,bi,bj)
59       &    *(       &    *(
60       &       dxv( i ,j+1,bi,bj)*dxv( i ,j+1,bi,bj)       &       dxV( i ,j+1,bi,bj)*dxV( i ,j+1,bi,bj)
61       &         *viscAh_s(i,j+1)*strain( i ,j+1)       &         *viscAh_s(i,j+1)*strain( i ,j+1)
62       &      -dxv( i , j ,bi,bj)*dxv( i , j ,bi,bj)       &      -dxV( i , j ,bi,bj)*dxV( i , j ,bi,bj)
63       &         *viscAh_s(i, j )*strain( i , j )       &         *viscAh_s(i, j )*strain( i , j )
64       &     )       &     )
65    
66           vDissip(i,j) =           vDissip(i,j) =
67       &    recip_dyc(i,j,bi,bj)*recip_dyc(i,j,bi,bj)       &    recip_dyC(i,j,bi,bj)*recip_dyC(i,j,bi,bj)
68       &    *recip_dxg(i,j,bi,bj)       &    *recip_dxG(i,j,bi,bj)
69       &    *(       &    *(
70       &       dyu(i+1, j ,bi,bj)*dyu(i+1, j ,bi,bj)       &       dyU(i+1, j ,bi,bj)*dyU(i+1, j ,bi,bj)
71       &         *viscAh_s(i+1,j)*strain(i+1,j)       &         *viscAh_s(i+1,j)*strain(i+1,j)
72       &      -dyu( i , j ,bi,bj)*dyu( i , j ,bi,bj)       &      -dyU( i , j ,bi,bj)*dyU( i , j ,bi,bj)
73       &         *viscAh_s( i ,j)*strain( i ,j)       &         *viscAh_s( i ,j)*strain( i ,j)
74       &     )       &     )
75       &    -recip_dxg(i,j,bi,bj)*recip_dxg(i,j,bi,bj)       &    -recip_dxG(i,j,bi,bj)*recip_dxG(i,j,bi,bj)
76       &    *recip_dyc(i,j,bi,bj)       &    *recip_dyC(i,j,bi,bj)
77       &    *(       &    *(
78       &       dxf( i , j ,bi,bj)*dxf( i , j ,bi,bj)       &       dxF( i , j ,bi,bj)*dxF( i , j ,bi,bj)
79       &         *viscAh_t(i, j )*tension(i, j )       &         *viscAh_t(i, j )*tension(i, j )
80       &      -dxf( i ,j-1,bi,bj)*dxf( i ,j-1,bi,bj)       &      -dxF( i ,j-1,bi,bj)*dxF( i ,j-1,bi,bj)
81       &         *viscAh_t(i,j-1)*tension(i,j-1)       &         *viscAh_t(i,j-1)*tension(i,j-1)
82       &     )       &     )
83    
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22