do i=1,nres_molec(molec)
itmp=itmp+1
istype(itmp)=sugarcode(sequence(i,molec)(1:1),i)
- itype(itmp,molec)=rescode(i,sequence(i,molec)(1:2),iscode,molec)
+ sequence(i,molec)=sequence(i,molec)(1:2)
+ itype(itmp,molec)=rescode(i,sequence(i,molec),iscode,molec)
enddo
endif