#!/usr/bin/perl

# This is a perl script to make ps images from snap data 

print("\n");
print (" plottme v. 0.1 ", "\n" );
print(" This is the plotting little program that uses snapplot to make ps images from snapshots ", "\n", "\n");

# Give the input *.snap file
print(" Input file = ");
$input = <STDIN>;
chomp($input);

print (" Starting time: ");
$start = <STDIN>;
chomp($start);

print (" Ending time: ");
$end = <STDIN>;
chomp($end);

print (" Time Step: ");
$tstep = <STDIN>;
chomp($tstep);

print(" Output name: ");
$name = <STDIN>;
chomp($name);

$ii=0;
for($i=$start;$i<=$end;$i+=$tstep)
  {
    $i = sprintf "%.5f", $i;
    if($ii<9)
      {
	$iii = "0000$ii";
      }
    if($ii>=10 && $ii<99)
      {
	$iii = "000$ii";
      }
    if($ii>=100 && $ii<999 )
      {
	$iii = "00$ii";
      }
    if($ii>=1000 && $ii<9999)
      {
	$iii = "0$ii";
      }
    if($ii>=10000)
      {
	$iii = "$ii";
      }
    print(" Printing snapshot for time ", $i, "\n");
    system(`snapplot $input xrange=-75:75 yrange=-80:80 times=$i color=" i<1?0.1:i<1002&i>1000?0.4:i<1000?0.2:i<1500?0.1:0.3" psize="i<1?0.1:i<1002&i>1000?0.1:0.001" fill_circle=f yapp=$i.ps/vcps`);

    #system(`snapplot $input xrange=-3:3 yrange=-3:3 color="i<143&i>141?0.3:i<10334&i>10332?0.4:i<10000?0.1:i>10000?0.2:0.5" psize="i<143&i>141?0.1:i<10334&i>10332?0.1:0.001" fill_circle=f visib="i<143&i>141||i<10334&i>10332" yapp=$i.ps/vcps`);

    #system(`snapplot $input xrange=-4:4 yrange=-4:4 times=$i color="i<143&i>141?0:i<16078&i>16076?0.4:i<10000?0.1:i>=10000?0.2:0.5" psize="i<143&i>141?0.01:i<16079&i>16077?0.1:0.001" fill_circle=f yapp=$i.ps/vcps`);

    print(" Converting into gif format...\n");
    system(" convert -density 196 $i.ps'[0]' -resize 100%  -trim  -bordercolor white -border 3 -size 768x512 $iii.gif ");
    system("rm -f $i.ps & ");
  #  system(" convert -adjoin $i.gif a.gif ");
  #  system("rm -f $i.gif & ");
    $ii += 1;
  }
print(" Merging gifs...");
system(" convert -adjoin *.gif $name.gif ");

system(" mv $name.gif $name.gif.temp &");
print(" Removing temporary files...\n");
system(" rm -f *.gif &");
system(" mv $name.gif.temp $name.gif &");
print(" End \n");