Modifying all the routines to call polarisation energy. Cannot be compiled at the moment.