virial function gives T,W,Pi Need Debuggit add getinfo.F90 main.F90 virial.F90 makefile