call etotal(energia(0),fT)
if (constr_homology) energia(0)=energia(0)+
& waga_homology(iset)*energia(22)
+ write (iout,*) "constr_homology",constr_homology," iset",iset,
+ & " waga_homology",waga_homology(iset)
#ifdef DEBUG
write (iout,*) "Conformation",i
call enerprint(energia(0),fT)
c & potE(i,iparm),-entfac(i),rms_nat,iscore
& potE(i,nparmset),-entfac(i),rms_nat,iscore
c endif
- write (iout,'(2i5,3e15.5)') i,me,potE(i,iparm),-entfac(i)
+c write (iout,'(2i5,3e15.5)') i,me,potE(i,iparm),-entfac(i)
#ifndef MPI
if (cxfile) call cxwrite(ixdrf,csingle,potE(i,nparmset),
& -entfac(i),rms_nat,iscore)