New routin added for energy conservation --Cettg ligne, et les suivantes ci-dessous, seront ignorées--
A phylmd/ener_conserv.F90