# Quick script to generate an ice Ih crystal

newModel("Ice Ih");
cell(4.5019, 7.7978, 7.3280, 90, 90, 90);
newAtomFrac(O, 0, 0.6648, 0.0631);
newAtomFrac(O, 0.5, 0.8255, -0.0631);
newAtomFrac(H, 0, 0.6636, 0.1963);
newAtomFrac(H, 0, 0.5363, 0.0183);
newAtomFrac(H, 0.6766, -0.2252, -0.0183);
spacegroup("Cmc21");
pack();
replicate(0, 0, 0, 4, 4, 4);
rebond();
