integer ilen
external ilen
integer iperm,tperm
- integer i,j,ii,k,l,itrial,itmp,i1,i2,it1,it2,nres_temp
+ integer i,j,ii,k,l,itrial,itmp,i1,i2,it1,it2,nres_temp,itemp
double precision sumv
C
C Read PDB structure if applicable
335 continue
unres_pdb=.false.
nres_temp=nres
-c call readpdb
- call readpdb_template(nmodel_start+1)
+ call readpdb
close(ipdbin)
if (nres.ge.nres_temp) then
nmodel_start=nmodel_start+1
enddo
enddo
else
- if (me.eq.king .or. .not. out1file)
+c itemp=nres
+c nres=nres_temp
+c call gen_rand_conf(itemp,*115)
+c nmodel_start=nmodel_start+1
+c do i=1,2*nres
+c do j=1,3
+c chomo(j,i,nmodel_start)=c(j,i)
+c enddo
+c enddo
+c goto 116
+ 115 if (me.eq.king .or. .not. out1file)
& write (iout,'(a,2i5,1x,a)')
& "Different number of residues",nres_temp,nres,
& " model skipped."
endif
+ 116 continue
nres=nres_temp
enddo
332 continue
endif
sigma_odl(k,ii)=1.0d0/(sigma_odl(k,ii)*sigma_odl(k,ii))
else
- ii=ii+1
- l_homo(k,ii)=.false.
+ct ii=ii+1
+ct l_homo(k,ii)=.false.
endif
enddo
enddo
if(me.eq.king .or. .not. out1file .and. fg_rank.eq.0) then
write (iout,*) "Distance restraints from templates"
do ii=1,lim_odl
- write(iout,'(3i5,100(2f8.2,1x,l1,4x))')
+ write(iout,'(i8,2i5,100(f8.2,e12.4,1x,l1,4x))')
& ii,ires_homo(ii),jres_homo(ii),
& (odl(ki,ii),1.0d0/dsqrt(sigma_odl(ki,ii)),l_homo(ki,ii),
& ki=1,constr_homology)