& dlog(dexp(quotl)+dexp(-quotl))
enddo
c write (iout,*) 1.0d0/(beta_h(ib)*1.987D-3),ft
- call flush(iout)
+c#ifdef AIX
+c call flush_(iout)
+c#else
+c call flush(iout)
+c#endif
endif
#ifdef MPI
ii=i
#endif
c write (iout,*) "i",i," ii",ii
- call flush(iout)
+c call flush(iout)
if (ib.eq.1) then
do j=1,nres
do k=1,3
call etotal(energia(0),fT)
totfree(i)=energia(0)
#ifdef DEBUG
- write (iout,*) i," energia",(energia(j),j=0,20)
+ write (iout,*) "Conformation",i
call enerprint(energia(0),ft)
write (iout,'(a,3(3x,a5),5x,3(3x,a5))')
& "Residue","X(CA)","Y(CA)","Z(CA)","X(SC)","Y(SC)","Z(SC)"
& (c(j,ires+nres),j=1,3)
enddo
- call intout
- call flush(iout)
+c call intout
+c call flush(iout)
#endif
do k=1,max_ene
enetb(k,i)=energia(k)