& wtor_d,wsccor,wbond
#endif
call etotal(energia(0),fT)
+#define DEBUG
#ifdef DEBUG
write (iout,*) "Conformation",i
call enerprint(energia(0),fT)
c write (iout,*) "ftors",ftors
c call intout
#endif
+#undef DEBUG
if (energia(0).ge.1.0d20) then
write (iout,*) "NaNs detected in some of the energy",
& " components for conformation",ii+1
do k=1,21
enetb(k,iii+1,iparm)=energia(k)
enddo
-#define DEBUG
#ifdef DEBUG
write (iout,'(2i5,f10.1,3e15.5)') i,iii,
& 1.0d0/(beta_h(ib,ipar)*1.987D-3),energia(0),eini,efree
call pdbout(ii+1,beta_h(ib,ipar),efree,energia(0),0.0d0,rmsdev)
close(ipdb)
#endif
-#undef DEBUG
endif
enddo ! iparm