anchor.diff:                  *(((crabfreq*two)*pi)/clight)*c1m6)*sin_mb(kcrab))*e0f
           end do
           call part_updatePartEnergy(1,.true.)
  -        if(ithick == 1) call envarsv
           goto 640
         case (60) ! JBG RF CC Multipoles
           xory=1
  
positive.code:                 ((three*crkve**2)*cikve)))*(((crabfreq*two)*pi)/clight))*c1m6)*(sin_mb(kcrab)*e0f)
          end do
          call part_updatePartEnergy(1,.true.)
          if(ithick == 1) call envarsv
          goto 640
        case (61) ! JBG RF CC Multipoles
          xory=1
  
negative.code:           ejv(j)=ejv(j) - ((((crabamp4)*((crkve**3*cikve)-(cikve**3*crkve)))*(((crabfreq*two)*pi)/clight))*c1m9)*(sin_mb(kcrab)*e0f)
          end do
          call part_updatePartEnergy(1,.true.)
          goto 640
        case (63) ! Elens
          do j=1,napx
  