energies.f90
6bc6183feba0geom_testpardiso

energies.f90
View Options

			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))vnormvnorm!
			DO i=1,nplasma
			ekin=ekin+msimvlightvlight(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)/rinnrnorm)

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