program idealp
c     ideally periodic peptide
      parameter (m0=10)
      character*1 a(m0)
      integer t(m0),b(m0,6),bb(6),nt(m0)
      real x(m0),f(m0),r(m0)
      write(6,*)'(Ala)n'
      write(6,*)'Input n:'
      nt(1)=3
      b(1,1)=2
      b(1,2)=3
      b(1,3)=-5
      nt(2)=1
      b(2,1)=1
      nt(3)=3
      b(3,1)=1
      b(3,2)=4
      b(3,3)=5
      nt(4)=1
      b(4,1)=3
      nt(5)=4
      b(5,1)=3
      b(5,2)=6
      b(5,3)=7
      b(5,4)=11
      nt(6)=1
      b(6,1)=5
      nt(7)=4
      b(7,1)=5
      b(7,2)=8
      b(7,3)=9
      b(7,4)=10
      nt(8)=1
      b(8,1)=7
      nt(9)=1
      b(9,1)=7
      nt(10)=1
      b(10,1)=7
      x( 1)=0.
      x( 2)=0.128
      x( 3)=0.818
      x( 4)=0.708
      x( 5)=1.982
      x( 6)=1.721
      x( 7)=2.416
      x( 8)=3.263
      x( 9)=1.611
      x(10)=2.677
      f( 1)=0.0
      f( 2)=-127.
      f( 3)=63.0
      f( 4)=50.0
      f( 5)=102.
      f( 6)=137.0
      f( 7)=91
      f( 8)=102
      f( 9)=95.
      f(10)=72.
      r( 1)=0.313
      r( 2)=1.168
      r( 3)=1.037
      r( 4)=1.953
      r( 5)=1.287
      r( 6)=1.776
      r( 7)=2.684
      r( 8)=2.996
      r( 9)=3.362
      r(10)=2.956
      a( 1)='C'
      a( 2)='O'
      a( 3)='N'
      a( 4)='H'
      a( 5)='C'
      a( 6)='H'
      a( 7)='C'
      a( 8)='H'
      a( 9)='H'
      a(10)='H'
      t( 1)=98
      t( 2)=99
      t( 3)=96
      t( 4)=4
      t( 5)=97
      t( 6)=100
      t( 7)=101
      t( 8)=102
      t( 9)=102
      t(10)=102
      px=3.14
      m=10
      f0=120.
      read(5,*)n
      open(9,file='i.xyz')
      nat=m*n
      write(9,901)nat
901   format(i6)
      ia=0
      do 1 i=1,n
      do 11 ii=1,m
      ia=ia+1
      xx=x(ii)+px*real(i-1)
      fi=(f(ii)+f0*real(i-1))*4.0*atan(1.)/180.
      yy=r(ii)*sin(fi)
      zz=r(ii)*cos(fi)
      do 2 j=1,nt(ii)
      bb(j)=b(ii,j)+(i-1)*m0
      if(bb(j).lt.0)bb(j)=bb(j)+m*n
2     if(bb(j).gt.nat)bb(j)=1
11    write(9,902)ia,a(ii),xx,yy,zz,t(ii),(bb(j),j=1,nt(ii))
1     continue
902   format(i6,1x,1x,a1,2x,3f12.6,6i6)

      close(9)

      end