& " the value read in: ",energia(0),eini," point",
& iii+1,indstart(me1)+iii," T",
& 1.0d0/(1.987D-3*beta_h(ib,ipar))
+#define DEBUG
#ifdef DEBUG
call enerprint(energia(0),fT)
#endif
+#undef DEBUG
#ifdef DEBUG
write (iout,'(4f10.5,2i5)') 0.0,energia(0),0.0,
& 1.0d0/(beta_h(ib,ipar)*1.987D-3),