+#ifdef DEBUG
+ write (iout,*) "itime",itime," temperature",t_bath,
+ & " potential energy",potE,potEcomp(0)
+ call enerprint(potEcomp)
+#endif
+ if(itime.eq.0) then
+ write(format,'(a1,a4,a1,a4,a10)') "(",format1,",",format2,
+ & ",100a12)"
+ write (istat,format) "#","",
+ & (ename(print_order(i)),i=1,nprint_ene)
+ endif