#!/bin/bash

N=2400
cat <<eof
<cg_molecule>
  <name>HEX</name>
  <ident>hexane</ident>
  <topology>
    <cg_beads>
eof
for ((i=1;i<=$N;i++)){
  cat <<eof
      <cg_bead>
        <name>SC$i</name>
        <type>SC</type>
        <mapping>UNITY</mapping>
        <beads>1:hex:SC$i</beads>
      </cg_bead>
eof
}
cat <<eof
  </cg_beads>
  <cg_bonded>
    <bond>
      <name>bond</name>
      <beads>
eof
for ((i=1;i<$N;i++)){
  echo "        SC$i SC$(($i+1))"
}
cat <<eof
      </beads>
    </bond>
    <angle>
      <name>angle</name>
      <beads>
eof
for ((i=1;i<$N;i++)){
	echo "        SC$i SC$(($i+1)) SC$(($i+2))"
}
cat <<eof
      </beads>
    </angle>
  </cg_bonded>
  </topology>
  <maps>
    <map>
      <name>UNITY</name>
      <weights>1</weights>
    </map>
  </maps>
</cg_molecule>
eof
