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

Diff of /MITgcm/pkg/mom_common/mom_calc_ke.F

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

revision 1.2 by jmc, Wed Aug 24 15:33:59 2005 UTC revision 1.3 by jmc, Fri Sep 23 19:16:25 2005 UTC
# Line 53  C                 1  /  ___I    ___J   \ Line 53  C                 1  /  ___I    ___J   \
53  C           KE = --- |   U    +  V     |  C           KE = --- |   U    +  V     |
54  C                 2  \                 /  C                 2  \                 /
55  C  C
56        IF (KEscheme.EQ.1) THEN        IF (KEscheme.EQ.-1) THEN
57         DO J=1-OLy,sNy+OLy-1         DO J=1-OLy,sNy+OLy-1
58          DO I=1-OLx,sNx+OLx-1          DO I=1-OLx,sNx+OLx-1
59           KE(i,j) = 0.125*(           KE(i,j) = 0.125*(
# Line 62  C Line 62  C
62         ENDDO         ENDDO
63        ENDDO        ENDDO
64    
65        ELSEIF (KEscheme.EQ.2) THEN        ELSEIF (KEscheme.EQ.0) THEN
66  C       This defn of KE should be used for the vector invariant equations.  C       This defn of KE should be used for the vector invariant equations.
67  C                       _____I  _____J  C                       _____I  _____J
68  C                 1  /     2       2   \  C                 1  /     2       2   \
# Line 72  C Line 72  C
72         DO J=1-OLy,sNy+OLy-1         DO J=1-OLy,sNy+OLy-1
73          DO I=1-OLx,sNx+OLx-1          DO I=1-OLx,sNx+OLx-1
74           KE(i,j) = 0.25*(           KE(i,j) = 0.25*(
75       &             uFld( i , j )*uFld( i , j )       &           ( uFld( i , j )*uFld( i , j )
76       &            +uFld(i+1, j )*uFld(i+1, j )       &            +uFld(i+1, j )*uFld(i+1, j ) )
77       &            +vFld( i , j )*vFld( i , j )       &         + ( vFld( i , j )*vFld( i , j )
78       &            +vFld( i ,j+1)*vFld( i ,j+1) )       &            +vFld( i ,j+1)*vFld( i ,j+1) )
79         &                  )
80          ENDDO          ENDDO
81         ENDDO         ENDDO
82    
83        ELSEIF (KEscheme.EQ.3) THEN        ELSEIF (KEscheme.EQ.1) THEN
84  C       As above but including the lopping factors and should be used  C       As above but including the area
85           DO J=1-OLy,sNy+OLy-1
86            DO I=1-OLx,sNx+OLx-1
87             KE(i,j) = 0.25*(
88         &           ( uFld(i, j )*uFld(i, j )*rAw(i ,j, bi,bj)
89         &            +uFld(i+1,j)*uFld(i+1,j)*rAw(i+1,j,bi,bj) )
90         &         + ( vFld(i, j )*vFld(i, j )*rAs(i ,j, bi,bj)
91         &            +vFld(i,j+1)*vFld(i,j+1)*rAs(i,j+1,bi,bj) )
92         &                  )*recip_rA(i,j,bi,bj)
93            ENDDO
94           ENDDO
95    
96          ELSEIF (KEscheme.EQ.2) THEN
97    C       As KEscheme=0 but including the lopping factors and should be used
98  C       for the conservative form of the momentum equations.  C       for the conservative form of the momentum equations.
99         DO J=1-OLy,sNy+OLy-1         DO J=1-OLy,sNy+OLy-1
100          DO I=1-OLx,sNx+OLx-1          DO I=1-OLx,sNx+OLx-1
101           KE(i,j) = 0.25*(           KE(i,j) = 0.25*(
102       &             uFld( i , j )*uFld( i , j )*hFacW(i,j,k,bi,bj)       &           ( uFld( i , j )*uFld( i , j )*hFacW(i,j,k,bi,bj)
103       &            +uFld(i+1, j )*uFld(i+1, j )*hFacW(i+1,j,k,bi,bj)       &            +uFld(i+1, j )*uFld(i+1, j )*hFacW(i+1,j,k,bi,bj) )
104       &            +vFld( i , j )*vFld( i , j )*hFacS(i,j,k,bi,bj)       &         + ( vFld( i , j )*vFld( i , j )*hFacS(i,j,k,bi,bj)
105       &            +vFld( i ,j+1)*vFld( i ,j+1)*hFacS(i,j+1,k,bi,bj)       &            +vFld( i ,j+1)*vFld( i ,j+1)*hFacS(i,j+1,k,bi,bj) )
106       &            )*recip_hFacC(i,j,k,bi,bj)       &            )*recip_hFacC(i,j,k,bi,bj)
107          ENDDO          ENDDO
108         ENDDO         ENDDO
109    
110        ELSEIF (KEscheme.EQ.4) THEN        ELSEIF (KEscheme.EQ.3) THEN
111  C       As above but including the area and should be used  C       As above but including the area
 C       for the conservative form of the momentum equations.  
112         DO J=1-OLy,sNy+OLy-1         DO J=1-OLy,sNy+OLy-1
113          DO I=1-OLx,sNx+OLx-1          DO I=1-OLx,sNx+OLx-1
114           KE(i,j) = 0.25*(           KE(i,j) = 0.25*(
115         &           (
116       &    uFld(i, j )*uFld(i, j )*hFacW(i ,j, k,bi,bj)*rAw(i ,j, bi,bj)       &    uFld(i, j )*uFld(i, j )*hFacW(i ,j, k,bi,bj)*rAw(i ,j, bi,bj)
117       &   +uFld(i+1,j)*uFld(i+1,j)*hFacW(i+1,j,k,bi,bj)*rAw(i+1,j,bi,bj)       &   +uFld(i+1,j)*uFld(i+1,j)*hFacW(i+1,j,k,bi,bj)*rAw(i+1,j,bi,bj)
118       &   +vFld(i, j )*vFld(i, j )*hFacS(i, j, k,bi,bj)*rAs(i ,j, bi,bj)       &           )  
119         &         + (
120         &    vFld(i, j )*vFld(i, j )*hFacS(i, j, k,bi,bj)*rAs(i ,j, bi,bj)
121       &   +vFld(i,j+1)*vFld(i,j+1)*hFacS(i,j+1,k,bi,bj)*rAs(i,j+1,bi,bj)       &   +vFld(i,j+1)*vFld(i,j+1)*hFacS(i,j+1,k,bi,bj)*rAs(i,j+1,bi,bj)
122       &                  )*recip_hFacC(i,j,k,bi,bj)*recip_rA(i,j,bi,bj)       &           )      )*recip_hFacC(i,j,k,bi,bj)*recip_rA(i,j,bi,bj)
123          ENDDO          ENDDO
124         ENDDO         ENDDO
125    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22