anchor.diff:    case(dist_fmtPX)
       yv1(1:napx) = (dist_partCol2(1:napx)/ejfv(1:napx))*c1e3
     case(dist_fmtPXP0)
  -    yv1(1:napx) = (dist_partCol2(1:napx)/e0f)*c1e3
  +    yv1(1:napx) = ((dist_partCol2(1:napx)*e0f)/ejfv(1:napx))*c1e3
     case(dist_fmtPXN)
       doNormal = .true.
     case(dist_fmtPhiX)
  
positive.code:   case(dist_fmtPY)
      yv2(1:napx) = (dist_partCol4(1:napx)/ejfv(1:napx))*c1e3
    case(dist_fmtPYP0)
      yv2(1:napx) = (dist_partCol4(1:napx)/e0f)*c1e3
    case(dist_fmtPYN)
      doNormal = .true.
    case(dist_fmtPhiY)
  
negative.code:   case(dist_fmtNONE)
      sigmv(1:napx) = zero
    case(dist_fmtZETA)
      sigmv(1:napx) = dist_partCol5(1:napx)/rvv(1:napx)
    case(dist_fmtSIGMA)
      sigmv(1:napx) = dist_partCol5(1:napx)
    case(dist_fmtDT)
      sigmv(1:napx) = -beta0*(dist_partCol5(1:napx)*clight)
    case(dist_fmtZN)
      doNormal = .true.
    case(dist_fmtJZ)
  