5
# Make a chain of monomers
6
monomers = new Monomer [12].rot(180, 1,0,0).move(2.0, 0, 0)
9
# Now, link the monomers together this way:
11
$bond:bb1 @bond:Monomer/bb $atom:monomers[0]/CA $atom:monomers[1]/CA
12
$bond:bb2 @bond:Monomer/bb $atom:monomers[1]/CA $atom:monomers[2]/CA
13
$bond:bb3 @bond:Monomer/bb $atom:monomers[2]/CA $atom:monomers[3]/CA
14
$bond:bb4 @bond:Monomer/bb $atom:monomers[3]/CA $atom:monomers[4]/CA
15
$bond:bb5 @bond:Monomer/bb $atom:monomers[4]/CA $atom:monomers[5]/CA
16
$bond:bb6 @bond:Monomer/bb $atom:monomers[5]/CA $atom:monomers[6]/CA
17
$bond:bb7 @bond:Monomer/bb $atom:monomers[6]/CA $atom:monomers[7]/CA
18
$bond:bb8 @bond:Monomer/bb $atom:monomers[7]/CA $atom:monomers[8]/CA
19
$bond:bb9 @bond:Monomer/bb $atom:monomers[8]/CA $atom:monomers[9]/CA
20
$bond:bb10 @bond:Monomer/bb $atom:monomers[9]/CA $atom:monomers[10]/CA
21
$bond:bb11 @bond:Monomer/bb $atom:monomers[10]/CA $atom:monomers[11]/CA
24
create_var { $mol } # Create a molecule ID number for this polymer
26
# This causes monomer[0], monomer[1], ... to share the same molecule counter
27
# because in the 2bead.lt file, the "..." in "$mol:..." preferentially looks
28
# for a counter of that type in a parent molecule or earlier ancestor.
34
# Angle, dihedral and improper interactions will be generated
35
# automatically according to the instructions in "monomer.lt"