write(iout,*) 'NNT=',NNT,' NCT=',NCT
c Read distance restraints
if (constr_dist.gt.0) then
+ if (refstr) call read_ref_structure(*11)
call read_dist_constr
call hpb_partition
endif
endif
write (iout,'(a)')
return
+ 11 stop "Error reading reference structure"
end
c-----------------------------------------------------------------------------
logical function seq_comp(itypea,itypeb,length)