123 |
vShearTerm(I,J)=-0.5* |
vShearTerm(I,J)=-0.5* |
124 |
& ( (wBarYp*vZp+wBarYm*vZm) |
& ( (wBarYp*vZp+wBarYm*vZm) |
125 |
& +(ABS(wBarYp)*vZp-ABS(wBarYm)*vZm) |
& +(ABS(wBarYp)*vZp-ABS(wBarYm)*vZm) |
126 |
& )*recip_hFacS(i,j,k,bi,bj) |
& )*_recip_hFacS(i,j,k,bi,bj) |
127 |
& *recip_drF(K) |
& * recip_drF(K) |
128 |
c3 & *recip_rAs(i,j,bi,bj) |
c3 & * recip_rAs(i,j,bi,bj) |
129 |
ELSE |
ELSE |
130 |
vShearTerm(I,J)=-0.5*(wBarYp*vZp+wBarYm*vZm) |
vShearTerm(I,J)=-0.5*(wBarYp*vZp+wBarYm*vZm) |
131 |
& *recip_hFacS(i,j,k,bi,bj) |
& *_recip_hFacS(i,j,k,bi,bj) |
132 |
& *recip_drF(K) |
& * recip_drF(K) |
133 |
c3 & *recip_rAs(i,j,bi,bj) |
c3 & * recip_rAs(i,j,bi,bj) |
134 |
ENDIF |
ENDIF |
135 |
ENDDO |
ENDDO |
136 |
ENDDO |
ENDDO |