anchor.diff:      call dist_appendFormat(dist_fmtP,        uFac, 6)
     case("DE/E0","DEE0")
       call dist_appendFormat(dist_fmtDEE0,     one,  6)
  -  case("DP/P0","DPP0")
  +  case("DP/P0","DPP0","DELTA")
       call dist_appendFormat(dist_fmtDPP0,     one,  6)
     case("DE/P0","DEP0","PT")
       call dist_appendFormat(dist_fmtPT,       one,  6)
  
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)
  