X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Fwham%2Fsrc-HCD%2Freadpdb.F;fp=source%2Fwham%2Fsrc-HCD%2Freadpdb.F;h=55dcc04de22bb65a7372fb8bad28a9aa37742be9;hb=225ad1b2239a64aee5d46570c477e74c19cf6ef7;hp=83c63abcede45c2f7544f0ee1ddf1dcda2712735;hpb=29dca221ada246dbf1a26c9e0e4fe12ee1e6b978;p=unres.git diff --git a/source/wham/src-HCD/readpdb.F b/source/wham/src-HCD/readpdb.F index 83c63ab..55dcc04 100644 --- a/source/wham/src-HCD/readpdb.F +++ b/source/wham/src-HCD/readpdb.F @@ -27,10 +27,9 @@ C geometry. ibeg=1 ishift1=0 sccalc=.false. - ires=0 do read (ipdbin,'(a80)',end=10) card -c write (iout,'(a)') card +! write (iout,'(a)') card if (card(:5).eq.'HELIX') then nhfrag=nhfrag+1 lsecondary=.true. @@ -59,7 +58,7 @@ c write (iout,'(a)') card iterter(ires_old)=1 ishift1=ishift1+1 ibeg=2 - write (iout,*) "Chain ended",ires,ishift,ires_old +! write (iout,*) "Chain ended",ires,ishift,ires_old if (unres_pdb) then do j=1,3 dc(j,ires)=sccor(j,iii) @@ -120,7 +119,7 @@ c write (iout,*) "BEG ires",ires ! Start a new chain ishift=-ires_old+ires-1 !!!!! ishift1=ishift1-1 !!!!! - write (iout,*) "New chain started",ires,ishift,ishift1,"!" +! write (iout,*) "New chain started",ires,ishift,ishift1,"!" ires=ires-ishift+ishift1 ires_old=ires ibeg=0 @@ -326,7 +325,7 @@ c write (iout,*) i,i-1,(c(j,i),j=1,3),(c(j,i-1),j=1,3),dist(i,i-1) if (itype(i-1).ne.ntyp1 .and. itype(i).ne.ntyp1 .and. & (dist(i,i-1).lt.1.0D0 .or. dist(i,i-1).gt.6.0D0)) then write (iout,'(a,i4)') 'Bad Cartesians for residue',i - stop +c stop endif vbld(i)=dist(i-1,i) vbld_inv(i)=1.0d0/vbld(i)