<<Graphics`Polyhedra`

Clear["@"];
(*** ** ** ** ** ********)
statef = "ground.m";
scale = 0.4;
(*** ** ** ** ** ********)
dir = "D:\\sergeev\\bgu\\math\\zilberg";
SetDirectory[dir];
Get[statef];
Print[title];
qpl := Transpose[qatoms + ampl qmodes[[nmode]]];
qmn := Transpose[qatoms - ampl qmodes[[nmode]]];
{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[t_] := Table[
      na = atnumb[[natom]];
      q = qatoms[[natom]] + ampl Sin[t]qmodes[[nmode, 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;
mt = 10;
ampl = 0.8;
rat = 1.6; (* y/x *)
off = 0.3;
size = 200;
SetDirectory["images"];
disp = 0;
Do[
    Print["nmode = ", nmode];
    fnames = " ";
    batchfile = "conv.bat";
    delay = 0.1;
    delay100 = Round[delay 100];
    animfile = 
      "vibr" <> If[nmode < 10, "0", ""] <> ToString[nmode] <> ".gif";
    Do[
      t = 2 Pi nt/mt;
      plt = Show[molec[t],
          PlotRange -> 1.5{{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},
          DisplayFunction -> If[disp === 0, Identity, $DisplayFunction]];
      pictfile = ToString[nt] <> ".tif";
      fnames = fnames <> pictfile <> " ";
      Display[pictfile, plt, "TIFF"];
      , {nt, mt}];
    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];
    , {nmode, nmodes}];