anchor.diff:          end if
   
             sampl(ie)    = totals
  -          ename(ie)    = bez(myix)
  -      end do !do j = 1, napx
  -    end if !if(rselect.gt.0 .and. rselect.lt.65) then
  -  end if !if( firstrun ) then
  +      end do
  +    end if
  +  end if
   
   !GRD HERE WE LOOK FOR ADEQUATE DATABASE INFORMATION
     found = .false.
  
positive.code:   call alloc(counteddpop, npart_new, numeffdpop, 0, "counteddpop") !(npart,numeffdpop)
    call alloc(counted2d,   npart_new, numeff, numeffdpop, 0, "counted2d") !(npart,numeff,numeffdpop)
  
    call alloc(ename,    mNameLen, nblz_new, ' ', "ename") !(nblz_new)
    call alloc(nampl,    nblz_new, 0, "nampl") !(nblz_new)
    call alloc(sum_ax,   nblz_new, zero, "sum_ax") !(nblz_new)
    call alloc(sqsum_ax, nblz_new, zero, "sqsum_ax") !(nblz_new)
  
negative.code:     end do
    end if
  
    if(iturn == 1 .and. firstrun) then
      write(outlun,*) ' '
      write(outlun,*)   'Collimator information: '
  