202 |
mask13c(i,j,bi,bj) = 0. _d 0 |
mask13c(i,j,bi,bj) = 0. _d 0 |
203 |
endif |
endif |
204 |
|
|
|
www(i,j,bi,bj) = cosphi(i,j,bi,bj) |
|
205 |
enddo |
enddo |
206 |
enddo |
enddo |
207 |
|
|
299 |
do j = jmin,jmax-1,2 |
do j = jmin,jmax-1,2 |
300 |
do i = i6min,imax-5, 6 |
do i = i6min,imax-5, 6 |
301 |
fctile_drift = fctile_drift |
fctile_drift = fctile_drift |
302 |
& + (wud(i,j,bi,bj)*www(i,j,bi,bj)*mask6c(i,j,bi,bj)* |
& + (wud(i,j,bi,bj)*cosphi(i,j,bi,bj)* |
303 |
|
& mask6c(i,j,bi,bj)* |
304 |
& (u6bar(i,j,bi,bj) - uddat(i,j,bi,bj))* |
& (u6bar(i,j,bi,bj) - uddat(i,j,bi,bj))* |
305 |
& (u6bar(i,j,bi,bj) - uddat(i,j,bi,bj)) ) |
& (u6bar(i,j,bi,bj) - uddat(i,j,bi,bj)) ) |
306 |
& + (wvd(i,j,bi,bj)*www(i,j,bi,bj)*mask6c(i,j,bi,bj)* |
& + (wvd(i,j,bi,bj)*cosphi(i,j,bi,bj)* |
307 |
|
& mask6c(i,j,bi,bj)* |
308 |
& (v6bar(i,j,bi,bj) - vddat(i,j,bi,bj))* |
& (v6bar(i,j,bi,bj) - vddat(i,j,bi,bj))* |
309 |
& (v6bar(i,j,bi,bj) - vddat(i,j,bi,bj)) ) |
& (v6bar(i,j,bi,bj) - vddat(i,j,bi,bj)) ) |
310 |
|
if ( cosphi(i,j,bi,bj)*mask6c(i,j,bi,bj) .ne. 0. ) then |
311 |
|
if ( wud(i,j,bi,bj) .ne. 0. ) |
312 |
|
& num_drift(bi,bj) = num_drift(bi,bj) + 1. _d 0 |
313 |
|
if ( wvd(i,j,bi,bj) .ne. 0. ) |
314 |
|
& num_drift(bi,bj) = num_drift(bi,bj) + 1. _d 0 |
315 |
|
endif |
316 |
enddo |
enddo |
317 |
enddo |
enddo |
318 |
|
|