<<Graphics`Polyhedra`
<<Graphics`Animation`

(* Should be executed as Mathematica notebook, pictures should be saved to files *)
Clear["@"];
(*** ** ** ** ** ********)
statef = "ground.m";
scale = 0.4;
(*** ** ** ** ** ********)
dir = "D:\\sergeev\\bgu\\math\\zilberg";
SetDirectory[dir];
Get[statef];
Print[title];
molec = Table[
      na = atnumb[[natom]];
      q = qatoms[[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}];
g = Show[molec,
      Boxed -> False,
      ImageSize -> 400,
      DisplayFunction -> Identity];
SpinShow[g,
    SpinOrigin -> {0, 0, 1.5},
    SpinDistance -> 2,
    SpinTilt -> {-30 Degree, -30 Degree},
    Frames -> 6,
    SpinRange -> {0 Degree, 360 Degree}];

(* Before execution of this part, gif-files should be present! *)
(* gif - animation *)
Clear["@"];
(*** ** ** *** INPUT starts **********)
{nmin, nmax, nstep} = {3, 62, 1};
filename[n_] := "a_gr_" <> ToString[n] <> ".gif";
delay = 0.1;
dir = "D:\\sergeev\\bgu\\math\\zilberg\\html\\images";
(*** ** ** *** INPUT ends **********)
batchfile = "conv.bat";
SetDirectory[dir];
delay100 = Round[delay 100];
fnames = " ";
Do[
    pictfile = filename[n];
    p = ToString[n];
    RenameFile[pictfile, p];
    fnames = fnames <> p <> " ";
    , {n, nmin, nmax, nstep}];
animfile = "animate.gif";
Print[fnames <> "-> " <> animfile];
os = OpenWrite[batchfile, PageWidth -> Infinity, FormatType -> OutputForm];
Write[os, 
    "D:\\programs\\ImageMagick-5.1.1\\convert -adjoin -delay " <> 
      ToString[delay100] <>
			fnames <> animfile];
Close[os];
Run[batchfile];