X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Fwham%2Fsrc%2Fmolread_zs.F;h=9c999ccad6f02dac468b8de514be01143a9d7777;hb=7e6e06667b2bab90b607bd8adecfbba5da26b0c9;hp=a98f47c91526792238456cddd5d03578c186b227;hpb=bee3b66ba871f85ef86d5328a081dcf2e974f6f5;p=unres.git diff --git a/source/wham/src/molread_zs.F b/source/wham/src/molread_zs.F index a98f47c..9c999cc 100644 --- a/source/wham/src/molread_zs.F +++ b/source/wham/src/molread_zs.F @@ -57,6 +57,19 @@ C Convert sequence to numeric code do i=1,nres itype(i)=rescode(i,sequence(i),iscode) enddo + if (itype(2).eq.10) then + write (iout,*) + & "Glycine is the first full residue, initial dummy deleted" + do i=1,nres + itype(i)=itype(i+1) + enddo + nres=nres-1 + endif + if (itype(nres).eq.10) then + write (iout,*) + & "Glycine is the last full residue, terminal dummy deleted" + nres=nres-1 + endif write (iout,*) "Numeric code:" write (iout,'(20i4)') (itype(i),i=1,nres) do i=1,nres-1