- call contact(.true.,ncont_ref,icont_ref,co)
+ call contact(.false.,ncont_ref,icont_ref,co)
! do k=1,nres
! write(icsa_native_int,200) (crefjlee(i,k),i=1,3)
! do k=1,nres
! write(icsa_native_int,200) (crefjlee(i,k),i=1,3)
- call contact(.true.,ncont_ref,icont_ref,co)
+ call contact(.false.,ncont_ref,icont_ref,co)
endif
! write (iout,*) "constr_dist",constr_dist,nstart_sup,nsup
! call flush(iout)
endif
! write (iout,*) "constr_dist",constr_dist,nstart_sup,nsup
! call flush(iout)
if (pdbref) then
if(me.eq.king.or..not.out1file) &
write (2,*) 'Shifting contacts:',nstart_seq,nstart_sup
if (pdbref) then
if(me.eq.king.or..not.out1file) &
write (2,*) 'Shifting contacts:',nstart_seq,nstart_sup
restyp(itype(icont_ref(2,i),1),1),' ',icont_ref(2,i)
enddo
endif
restyp(itype(icont_ref(2,i),1),1),' ',icont_ref(2,i)
enddo
endif
if (constr_homology.gt.0) then
! write (iout,*) "Calling read_constr_homology"
! call flush(iout)
if (constr_homology.gt.0) then
! write (iout,*) "Calling read_constr_homology"
! call flush(iout)
if (itype(chain_border(1,i)+k,mnum).ne.&
itype(chain_border(1,j)+k,mnum)) exit
enddo
if (itype(chain_border(1,i)+k,mnum).ne.&
itype(chain_border(1,j)+k,mnum)) exit
enddo
write(iout,'(20i4)') (tabpermchain(j,i),j=1,nchain)
enddo
return
write(iout,'(20i4)') (tabpermchain(j,i),j=1,nchain)
enddo
return
!c---------------------------------------------------------------------
integer function tperm(i,iperm,tabpermchain)
! implicit none
!c---------------------------------------------------------------------
integer function tperm(i,iperm,tabpermchain)
! implicit none