totfree(i)=energia(0)
c write (iout,'(8f10.5)') ((c(l,k),l=1,3),k=1,nres)
c write (iout,'(8f10.5)') ((c(l,k+nres),l=1,3),k=nnt,nct)
- call enerprint(energia(0),fT)
c call pdbout(totfree(i),16,i)
#ifdef DEBUG
write (iout,*) i," energia",(energia(j),j=0,19)
write (iout,*) "etot", etot
write (iout,*) "ft(6)", ft(6)
+ call enerprint(energia(0),fT)
#endif
do k=1,max_ene
enetb(k,i)=energia(k)