/[MITgcm]/MITgcm/pkg/ecco/cost_drifter.F
ViewVC logotype

Diff of /MITgcm/pkg/ecco/cost_drifter.F

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

revision 1.2 by heimbach, Mon Oct 11 16:38:53 2004 UTC revision 1.3 by heimbach, Mon Mar 28 23:49:49 2005 UTC
# Line 202  c                 set mask13c=0 in areas Line 202  c                 set mask13c=0 in areas
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    
# Line 300  c             drifters. Line 299  c             drifters.
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    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22