168 |
DO i=1-Olx,sNx+Olx |
DO i=1-Olx,sNx+Olx |
169 |
C First compute the fraction of open water for the w-control volume |
C First compute the fraction of open water for the w-control volume |
170 |
C at the southern face |
C at the southern face |
171 |
hFacCtmp=max(hFacC(I,J,K-1,bi,bj)-Half,0. _d 0) |
hFacCtmp=max( _hFacC(I,J,K-1,bi,bj)-Half,0. _d 0 ) |
172 |
& + min(hFacC(I,J,K ,bi,bj),Half) |
& + min( _hFacC(I,J,K ,bi,bj) ,Half ) |
173 |
IF (hFacCtmp .GT. 0.) THEN |
IF (hFacCtmp .GT. 0.) THEN |
174 |
recip_hFacCtmp = 1./hFacCtmp |
recip_hFacCtmp = 1./hFacCtmp |
175 |
ELSE |
ELSE |
202 |
DO I=iMin,iMax |
DO I=iMin,iMax |
203 |
C First compute the fraction of open water for the w-control volume |
C First compute the fraction of open water for the w-control volume |
204 |
C at the southern face |
C at the southern face |
205 |
hFacStmp=max(hFacS(I,J,K-1,bi,bj)-Half,0. _d 0) |
hFacStmp=max(_hFacS(I,J,K-1,bi,bj)-Half,0. _d 0) |
206 |
& + min(hFacS(I,J,K ,bi,bj),Half) |
& + min(_hFacS(I,J,K ,bi,bj),Half) |
207 |
tmp_VbarZ=Half*( |
tmp_VbarZ=Half*( |
208 |
& _hFacS(I,J,K-1,bi,bj)*vVel( I ,J,K-1,bi,bj) |
& _hFacS(I,J,K-1,bi,bj)*vVel( I ,J,K-1,bi,bj) |
209 |
& +_hFacS(I,J, K ,bi,bj)*vVel( I ,J, K ,bi,bj)) |
& +_hFacS(I,J, K ,bi,bj)*vVel( I ,J, K ,bi,bj)) |
234 |
DO I=iMin,iMax+1 |
DO I=iMin,iMax+1 |
235 |
C First compute the fraction of open water for the w-control volume |
C First compute the fraction of open water for the w-control volume |
236 |
C at the western face |
C at the western face |
237 |
hFacWtmp=max(hFacW(I,J,K-1,bi,bj)-Half,0. _d 0) |
hFacWtmp=max(_hFacW(I,J,K-1,bi,bj)-Half,0. _d 0) |
238 |
& + min(hFacW(I,J,K ,bi,bj),Half) |
& + min(_hFacW(I,J,K ,bi,bj),Half) |
239 |
tmp_UbarZ=Half*( |
tmp_UbarZ=Half*( |
240 |
& _hFacW(I,J,K-1,bi,bj)*uVel( I ,J,K-1,bi,bj) |
& _hFacW(I,J,K-1,bi,bj)*uVel( I ,J,K-1,bi,bj) |
241 |
& +_hFacW(I,J, K ,bi,bj)*uVel( I ,J, K ,bi,bj)) |
& +_hFacW(I,J, K ,bi,bj)*uVel( I ,J, K ,bi,bj)) |