This routine computes energy, forces and the virial. Derivatives by James Kermode <jrk33@cam.ac.uk>.
e = System total energy, local_e = energy of each atom, vector dimensioned as at%N.
e = System total energy, local_e = energy of each atom, vector dimensioned as at%N.
Forces, dimensioned as f(3,at%N)
Virial