- write (iunit,10) iatom,restyp(iti,1),chainid(ichain),&
+ if (molnum(i).eq.1) then
+
+ write (iunit,10) iatom,restyp(iti,molnum(i)),chainid(ichain),&
+ ires,(c(j,i),j=1,3),vtot(i)
+ elseif(molnum(i).eq.2) then
+ if (istype(i).eq.0) istype(i)=1
+ write (iunit,40) iatom,sugartyp(istype(i)),restyp(iti,2), &
+ chainid(ichain),ires,(c(j,i),j=1,3),vtot(i)
+ else
+ write (iunit,60) iatom,restyp(iti,molnum(i)),chainid(ichain),&