anchor.diff:    write(outlun,*) 'parameters. Distribution is flat in the band.'
     write(outlun,*) 'X and Y are fully uncorrelated.'
     write(outlun,*)
  -  write(outlun,*) 'INFO>  Number of particles   = ', mynp
  +  write(outlun,*) 'INFO>  Number of particles   = ', napx
     write(outlun,*) 'INFO>  Av number of x sigmas = ', mynex
     write(outlun,*) 'INFO>  +- spread in x sigmas = ', mdex
     write(outlun,*) 'INFO>  Av number of y sigmas = ', myney
  
positive.code:   iturn = 1
    ie    = 1
    n_tot_absorbed = 0
  
    if(int(mynp/napx00) .eq. 0) then
      write (lout,"(a)")    "COLL> ERROR Setting up collimation tracking: Number of samples is zero!"
      write (lout,"(a)")    "COLL> ERROR Did you forget the COLL block in fort.3?"
      write (lout,"(a,l1)") "COLL> * Value of do_coll = ", do_coll
      write (lout,"(a,i0)") "COLL> * Value of mynp    = ", mynp
      write (lout,"(a,i0)") "COLL> * Value of napx00  = ", napx00
      call prror(-1)
    end if
  
  
    call funit_requestUnit('CollPositions.dat', CollPositions_unit)
    open(unit=CollPositions_unit, file='CollPositions.dat')
  
negative.code:   write(outlun,*) 'INFO>  Alpha y  = ', myalphay
    write(outlun,*)
  
    do while (j.lt.napx)
      j = j + 1
      myemitx = myemitx0*(mynex + ((two*real(rndm4()-half,fPrec))*mdex) )**2
      xsigmax = sqrt(mybetax*myemitx)
  