anchor.diff:      ! correct TOT-CSec for energy dependence of qel
       ! TOT CS is here without a Coulomb contribution
       csect(0,ma) = csref(0,ma) + freep(ma) * (pptot - pptref)
  -    bn(ma) = (bnref(ma) * csect(0,ma)) / csref(0,ma)
  +    bn(ma)      = (bnref(ma) * csect(0,ma)) / csref(0,ma)
   
       ! also correct inel-CS
       csect(1,ma) = (csref(1,ma) * csect(0,ma)) / csref(0,ma)
  
positive.code:         linside(j) = .true.
          if(dowrite_impact) then
            if(tiltangle > zero) then
              x_Dump = (x+c_aperture/two+tiltangle*sp)*mirror+c_offset
            else
              x_Dump = (x+c_aperture/two+tiltangle*(sp-c_length))*mirror+c_offset
            end if
            xpDump = (xp+tiltangle)*mirror
            y_Dump = z
            ypDump = zp
            s_Dump = sp+real(j_slices-1,fPrec)*c_length
  
negative.code:         linside(j) = .false.
          if(dowrite_impact) then
            if(tiltangle > zero) then
              x_Dump = (x + c_aperture/two + tiltangle*(s+sp))*mirror + c_offset
            else
              x_Dump = (x + c_aperture/two + tiltangle*(s+sp-c_length))*mirror + c_offset
            end if
            xpDump = (xp+tiltangle)*mirror
            y_Dump = z
  