unres_pdb=.false.
if (read2sigma) then
call readpdb_template(k)
+ close(ipdbin)
else
call readpdb
+ close(ipdbin)
endif
c call readpdb
write (iout,*) "read_constr_homology: after reading pdb file"
call flush(iout)
#endif
-
c
c Distance restraints
c
do irec=nnt,nct ! loop for reading res sim
if (read2sigma) then
read (ientin,*,end=1401) i_tmp,rescore2_tmp,rescore_tmp,
- & idomain_tmp,
& rescore3_tmp,idomain_tmp
i_tmp=i_tmp+nnt-1
idomain(k,i_tmp)=idomain_tmp
rescore(k,i_tmp)=rescore_tmp
rescore2(k,i_tmp)=rescore2_tmp
rescore3(k,i_tmp)=rescore3_tmp
- write(iout,'(a7,i5,2f10.5,i5)') "rescore",
+ write(iout,'(a7,i5,3f10.5,i5)') "rescore",
& i_tmp,rescore2_tmp,rescore_tmp,
& rescore3_tmp,idomain_tmp
else