anchor.diff:    implicit none
   
     real(kind=fPrec), intent(in) :: refEnergy
  -  integer,          intent(in) :: refArray   ! 1 for energy, 2 for momentum
  +  integer,          intent(in) :: refArray
   
     real(kind=fPrec) e0o, e0fo
     integer          j
  
positive.code:       xv(2,j)   = xv(2,j) + clov(2,j)  * real(idz(2),fPrec)
        yv(2,j)   = yv(2,j) + clopv(2,j) * real(idz(2),fPrec)
      end if
      ejfv(j)     = e0f*(one+dpsv(j))
      ejv(j)      = sqrt(ejfv(j)**2 + nucm(j)**2)
      moidpsv(j)  = mtc(j)/(one+dpsv(j))
      omoidpsv(j) = ((one-mtc(j))*oidpsv(j))*c1e3
  
negative.code: 
    ! Modify the Energy Dependent Arrays
    do j=1, napx
      dpsv1(j)    = (dpsv(j)*c1e3)/(one + dpsv(j))
      dpd(j)      = one + dpsv(j)                      ! For thick tracking
      dpsq(j)     = sqrt(dpd(j))                       ! For thick tracking
  