anchor.diff:  !  logical firstrun
   
     integer, save :: nsurvive, nsurvive_end, num_selhit, n_impact
  -
  -  integer, save :: napx00
   
     real(kind=fPrec), allocatable, save :: db_tilt(:,:) !(max_ncoll,2)
   
  
positive.code:   write(lout,"(a,e15.8)") 'COLL> Info: SIGSECUT3           = ', sigsecut3
    write(lout,"(a)")
  
    napx00 = napx
  
    write(lout,"(a,i0)")    'COLL> Info: NAPX                = ', napx
    write(lout,"(a,e15.8)") 'COLL> Info: Sigma_x0            = ', sqrt(mybetax*myemitx0_dist)
    write(lout,"(a,e15.8)") 'COLL> Info: Sigma_y0            = ', sqrt(mybetay*myemity0_dist)
  
negative.code:   firstcoll = .true.
  
  !GRD HERE WE NEED TO INITIALIZE SOME COLLIMATION PARAMETERS
  
    do j = 1, napx
      part_hit_pos(j)   = 0
  