49 |
DO K=1,Nr |
DO K=1,Nr |
50 |
DO J=1,sNy |
DO J=1,sNy |
51 |
DO I=1,sNx |
DO I=1,sNx |
52 |
volTile=volTile+rA(i,j,bi,bj)*drF(k)*hFacC(i,j,k,bi,bj) |
volTile=volTile+rA(i,j,bi,bj)*drF(k)*_hFacC(i,j,k,bi,bj) |
53 |
|
|
54 |
C- Vector Invariant form (like in pkg/mom_vecinv/mom_vi_calc_ke.F) |
C- Vector Invariant form (like in pkg/mom_vecinv/mom_vi_calc_ke.F) |
55 |
c tmpVal=0.25*( uVel( I , J ,K,bi,bj)*uVel( I , J ,K,bi,bj) |
c tmpVal=0.25*( uVel( I , J ,K,bi,bj)*uVel( I , J ,K,bi,bj) |
64 |
C with no assumption on how grid spacing & area are defined. |
C with no assumption on how grid spacing & area are defined. |
65 |
tmpVal=0.25*( |
tmpVal=0.25*( |
66 |
& uVel( i ,j,k,bi,bj)*uVel( i ,j,k,bi,bj) |
& uVel( i ,j,k,bi,bj)*uVel( i ,j,k,bi,bj) |
67 |
& *dyG( i ,j,bi,bj)*dxC( i ,j,bi,bj)*hFacW( i ,j,k,bi,bj) |
& *dyG( i ,j,bi,bj)*dxC( i ,j,bi,bj)*_hFacW( i ,j,k,bi,bj) |
68 |
& +uVel(i+1,j,k,bi,bj)*uVel(i+1,j,k,bi,bj) |
& +uVel(i+1,j,k,bi,bj)*uVel(i+1,j,k,bi,bj) |
69 |
& *dyG(i+1,j,bi,bj)*dxC(i+1,j,bi,bj)*hFacW(i+1,j,k,bi,bj) |
& *dyG(i+1,j,bi,bj)*dxC(i+1,j,bi,bj)*_hFacW(i+1,j,k,bi,bj) |
70 |
& +vVel(i, j ,k,bi,bj)*vVel(i, j ,k,bi,bj) |
& +vVel(i, j ,k,bi,bj)*vVel(i, j ,k,bi,bj) |
71 |
& *dxG(i, j ,bi,bj)*dyC(i, j ,bi,bj)*hFacS(i, j ,k,bi,bj) |
& *dxG(i, j ,bi,bj)*dyC(i, j ,bi,bj)*_hFacS(i, j ,k,bi,bj) |
72 |
& +vVel(i,j+1,k,bi,bj)*vVel(i,j+1,k,bi,bj) |
& +vVel(i,j+1,k,bi,bj)*vVel(i,j+1,k,bi,bj) |
73 |
& *dxG(i,j+1,bi,bj)*dyC(i,j+1,bi,bj)*hFacS(i,j+1,k,bi,bj) |
& *dxG(i,j+1,bi,bj)*dyC(i,j+1,bi,bj)*_hFacS(i,j+1,k,bi,bj) |
74 |
& ) |
& ) |
75 |
volMeanTile= volMeanTile + tmpVal*drF(k) |
volMeanTile= volMeanTile + tmpVal*drF(k) |
76 |
tmpVal= tmpVal*recip_hFacC(i,j,k,bi,bj)*recip_rA(i,j,bi,bj) |
tmpVal= tmpVal*_recip_hFacC(i,j,k,bi,bj)*recip_rA(i,j,bi,bj) |
77 |
|
|
78 |
theMax=max(theMax,tmpVal) |
theMax=max(theMax,tmpVal) |
79 |
IF (tmpVal.NE.0.) THEN |
IF (tmpVal.NE.0.) THEN |