anchor.diff:                  *(((crabfreq*two)*pi)/clight))*c1m9)*(sin_mb(kcrab)*e0f)
           end do
           call part_updatePartEnergy(1,.true.)
  -        if(ithick == 1) call envarsv
           goto 640
         case (62) ! JBG RF CC Multipoles
           xory=1
  
positive.code:           ejv(j)=ejv(j) - ((((crabamp2)*(cikve*crkve))*(((crabfreq*two)*pi)/clight))*c1m3)*(sin_mb(kcrab)*e0f)
          end do
          call part_updatePartEnergy(1,.true.)
          if(ithick == 1) call envarsv
          goto 640
        case (59) ! JBG RF CC Multipoles
          xory=1
  
negative.code:           ejv(j)=ejv(j) - ((((half*(crabamp2))*(crkve**2-cikve**2))*(((crabfreq*two)*pi)/clight))*c1m3)*(sin_mb(kcrab)*e0f)
          end do
          call part_updatePartEnergy(1,.true.)
          goto 640
        case (58) ! JBG RF CC Multipoles
          xory=1
  