/[PAMELA software]/DarthVader/TrackerLevel2/src/F77/reductionflight.f
ViewVC logotype

Diff of /DarthVader/TrackerLevel2/src/F77/reductionflight.f

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

revision 1.20 by pam-fi, Thu May 24 16:45:48 2007 UTC revision 1.21 by pam-fi, Tue Aug 7 13:56:29 2007 UTC
# Line 454  c     add strips exceeding inclusion cut Line 454  c     add strips exceeding inclusion cut
454  c------------------------------------------------------------------------  c------------------------------------------------------------------------
455                    if( (rmax-lmax+1).ge.nclstrp )goto 210   !exits inclusion loop                    if( (rmax-lmax+1).ge.nclstrp )goto 210   !exits inclusion loop
456    
457                    if(rstop.eq.0) then !if right cluster morder has not been reached                    if(rstop.eq.0) then !if right cluster border has not been reached
458                       if(value(ir).gt.clinclcut(ir)) then                       if(value(ir).gt.clinclcut(ir)) then
459                          rmax=ir !include a strip on the right                          rmax=ir !include a strip on the right
460                       else                       else
# Line 464  c--------------------------------------- Line 464  c---------------------------------------
464    
465                    if( (rmax-lmax+1).ge.nclstrp )goto 210   !exits inclusion loop                    if( (rmax-lmax+1).ge.nclstrp )goto 210   !exits inclusion loop
466    
467                    if(lstop.eq.0) then !if left cluster morder has not been reached                    if(lstop.eq.0) then !if left cluster border has not been reached
468                       if(value(il).gt.clinclcut(il)) then                       if(value(il).gt.clinclcut(il)) then
469                          lmax=il !include a strip on the left                          lmax=il !include a strip on the left
470                       else                       else
# Line 472  c--------------------------------------- Line 472  c---------------------------------------
472                       endif                       endif
473                    endif                    endif
474    
475    c                  if( (rmax-lmax+1).ge.nclstrp )goto 210   !exits inclusion loop
476    
477                 enddo            !ends strip inclusion loop                 enddo            !ends strip inclusion loop
478                 goto 211                 goto 211
479   210           continue         !jumps here if more than nclstrp have been included   210           continue         !jumps here if more than nclstrp have been included
# Line 492  c--------------------------------------- Line 494  c---------------------------------------
494                    rmax = rmax+1                    rmax = rmax+1
495                    if( (rmax-lmax+1).gt.nclstrp )lmax=lmax+1                    if( (rmax-lmax+1).gt.nclstrp )lmax=lmax+1
496                 endif                 endif
497    c-------------------------------------------------------------------------------
498    c     adjust the cluster in order to have at least ANOTHER strip around the seed
499    c-------------------------------------------------------------------------------
500    c$$$               if(iseed-1.eq.lmax.and.lmax.ne.first)then
501    c$$$                  lmax = lmax-1
502    c$$$                  if( (rmax-lmax+1).gt.nclstrp )rmax=rmax-1
503    c$$$               endif
504    c$$$               if(iseed+1.eq.rmax.and.rmax.ne.last )then
505    c$$$                  rmax = rmax+1
506    c$$$                  if( (rmax-lmax+1).gt.nclstrp )lmax=lmax+1
507    c$$$               endif
508    c---------------------------------------------------
509    c     now we have 5 stored-strips around the maximum
510    c---------------------------------------------------
511    
512  c------------------------------------------------------------------------  c------------------------------------------------------------------------
513  c     adjust the cluster in order to store a minimum number of strips  c     adjust the cluster in order to store a minimum number of strips
Legend:
Removed from v.1.20  
changed lines
  Added in v.1.21
  ViewVC Help
Powered by ViewVC 1.1.23