SUBROUTINE energies
!
!  Compute different energies
!
  USE basic
  IMPLICIT NONE
  REAL(kind=db) :: add
  INTEGER :: i, J1, J2, J3, J4
!
  ekin=0._db
  epot=0._db
  etot=0._db
!  msim/2*(VR(i)*VR(i)+VZ(i)*VZ(i)+VTHET(i)*VTHET(i))*vnorm*vnorm!
  DO i=1,nplasma
     ekin=ekin+msim*vlight*vlight*(1/sqrt(1-VR(i)*VR(i)-VZ(i)*VZ(i)-VTHET(i)*VTHET(i))-1)
  END DO

  DO i=1,nplasma
     J1=(rindex(i)-1)*(nz+1)+zindex(i)
     J2=J1+1
     J3=rindex(i)*(nz+1)+zindex(i)
     J4=J3+1
!
! Interpolation for electric field
     add = (1-WZ(i))*(1-WR(i))*pot(J1)+WZ(i)*(1-WR(i))*pot(J2)+(1-WZ(i))*WR(i)*pot(J3)+ &
          &  WZ(i)*WR(i)*pot(J4)
!       
!add=potinn+(potout-potinn)/(log(rout/rinn))*log(R(i)/rinn*rnorm)

     epot=epot+qsim*add*phinorm
  END DO
!DO i=1,(nrank(1)-1)*nrank(2)
!epot=epot+sol(i)*rhs(i)
!END DO
  etot=epot+ekin
!  
END SUBROUTINE energies