anchor.diff:        yv2(j) = rcyp(j)*c1e3 + torbyp(ie)
         ejv(j) =  rcp(j)*c1e3
   
  -      !  Energy update, as recommended by Frank
  +      ! call part_updatePartEnergy(1,.false.)
         ejfv(j)     = sqrt(ejv(j)**2-nucm(j)**2)
         rvv(j)      = (ejv(j)*e0f)/(e0*ejfv(j))
         dpsv(j)     = (ejfv(j)*(nucm0/nucm(j))-e0f)/e0f
  
positive.code:   use mod_common
    use coll_common
    use mod_settings
    use mod_common_main
    use mod_common_track
    use coll_db
  
negative.code:         sigmv(j)         = zero
          part_abs_pos(j)  = ie
          part_abs_turn(j) = iturn
          numxv(j)         = numx
          nabs_type(j)     = 0
          nhit_stage(j)    = 0
        end if
      end if
    end do
  
  