call etotal(energia(0),fT)
totfree(i)=energia(0)
#ifdef DEBUG
- write (iout,*) i," energia",(energia(j),j=0,19)
+ write (iout,*) i," energia",(energia(j),j=0,21)
+ call enerprint(energia(0),ft)
+ call flush(iout)
#endif
do k=1,max_ene
enetb(k,i)=energia(k)