Fixed minor (inconsequential to simulation data) bug by setting ener_pe = 0 outside conditional so that the 
cumulative potential energy does not grow due to slow accumulation over many processors. Also added optional diagnostic 
print statements.