do i=1,nres
itype(i)=rescode(i,sequence(i),iscode)
enddo
- if (itype(2).eq.10) then
+ if (itype(2).eq.10.and.itype(1).eq.ntyp1) then
write (iout,*)
& "Glycine is the first full residue, initial dummy deleted"
do i=1,nres
enddo
nres=nres-1
endif
- if (itype(nres).eq.10) then
+ if (itype(nres-1).eq.10.and.itype(nres).eq.ntyp1) then
write (iout,*)
& "Glycine is the last full residue, terminal dummy deleted"
nres=nres-1