3
# --------------- water ------------------
5
# Create a rhombohedral box of water. (A rectangular box works also.)
7
wat = new SPCE [9].move(3.5526287, 0, 0 )
8
[9].move(1.77631435, 3.3133, 0 )
11
# Optional: Center the water box at the origin. (Not really necessary.)
13
wat[*][*][*].move(-23.9802437, -14.90985, 11.47)
15
# --------------- Note: -----------------
16
# The spacing between water molecules does not matter much as long as it is
17
# reasonable. (I adjusted the spacing try to insure that the waters are spread
18
# uniformly throughout the box. We do not want bubles to form if there are
19
# gaps near the XY periodic boundaries.) We will have to equilibrate it later.