<<Graphics`Polyhedra`

(* Compare ground vs. excited state *)
 Clear["@"];
scale = 0.4;
dir = "D:\\sergeev\\bgu\\math\\zilberg";
SetDirectory[dir];
<< ground.m;
qat[1] = qatoms;
<< excit.m;
qat[2] = qatoms;
qpl := Transpose[qat[1]];
qmn := Transpose[qat[2]];
{xmin, ymin, zmin} := {Min[qpl[[1]], qmn[[1]]], Min[qpl[[2]], qmn[[2]]], 
      Min[qpl[[3]], qmn[[3]]]};
{xmax, ymax, zmax} := {Max[qpl[[1]], qmn[[1]]], Max[qpl[[2]], qmn[[2]]], 
      Max[qpl[[3]], qmn[[3]]]};
molec[qa_] := Table[
      na = atnumb[[natom]];
      q = qa[[natom]];
      If[na === 1,
        Polyhedron[Dodecahedron, q, 0.6 scale],
        If[na === 6, Polyhedron[Icosahedron, q, scale],
          Truncate[Polyhedron[Dodecahedron, q, 1.0 scale], 0.4]
          ]]
      , {natom, matoms}];
rxy = 2.7;
phi = 2 Pi 13/22;
ampl = 0.8;
rat = 1.6; (* y/x *)
off = 0.2;
size = 400;
SetDirectory["images"];
disp = 0;
fnames = " ";
batchfile = "conv.bat";
delay = 0.3;
delay100 = Round[delay 100];
animfile = "compare.gif";
Do[
    plt = Show[molec[qat[nstate]],
        PlotRange -> 1.2{{xmin, xmax}, {ymin, ymax}, {zmin, zmax}},
        ViewPoint -> {rxy Cos[phi], rxy Sin[phi], 2.},
        Boxed -> False,
        PlotRegion -> {{-off, 1 + off}, {-off, 1 + off}},
        ImageSize -> size{1/rat, 1},
        Background -> 
          If[nstate == 1, RGBColor[1, 1, 1], RGBColor[1, 0.95, 0.95]],
        DisplayFunction -> If[disp === 0, Identity, $DisplayFunction]];
    pictfile = ToString[nstate] <> ".tif";
    fnames = fnames <> pictfile <> " ";
    Display[pictfile, plt, "TIFF"];
    , {nstate, 2}];
os = OpenWrite[batchfile, PageWidth -> Infinity, FormatType -> OutputForm];
Write[os, 
    "D:\\programs\\ImageMagick-5.1.1\\convert -adjoin -loop 4096 -delay " <> 
      ToString[delay100] <>
			fnames <> animfile];
Close[os];
Run[batchfile];