#!/usr/bin/perl

$pi=3.1415926535;

$n=180;

for ($i=0; $i<$n; $i++) {
  $az = 70*sin(2.*$pi*$i/$n);
  print "az = $az\n";

  open(IN,">sh.in");
  print IN <<EOF
50:10:38	/* zeměpisná šířka (phi) stanoviště [+DD:MM:SS] */
15:50:21	/* zeměpisná délka (lambda) stanoviště [+DD:MM:SS, kladně na V] */
$az 0		/* azimut a výška normály ke stěně [deg] */
-150 -250	/* levý dolní roh číselníku x, y [cm] */
150 10		/* pravý horní roh číselníku x, y [cm] */
75		/* vzdálenost nodu od paty ukazatele (délka ukazatele) [cm]; pro gnómon vzdálenost od roviny číselníku */
1:0:0		/* interval hodinových rysek dt [HH:MM:SS] */
1:0:0		/* interval hodin od východu/západu slunce a temporálních dvz [HH:MM:SS] */
0:0:10		/* interval datových křivek dd [HH:MM:SS] */
10		/* interval azimutálních úseček daz [deg] */
10		/* interval výškových křivek dht [deg] */
2		/* krok pro výpočet výškových křivek dah [deg] */
7		/* počet datových křivek (0, 3 nebo 7) */
0		/* polos nebo gnómon? [0|1] */
0		/* zobrazit i neskutečné čáry pod horizontem? [0|1] */
1		/* zobrazit hodinové rysky? [0|1] */
0		/* pravý sluneční čas nebo korigovaný o délku lambda na SEČ? [0|1] */
0		/* "zimní" nebo letní čas? [0|1] */
1		/* prodloužení rysek: žádné, k patě polosu, k okraji číselníku, obojí? [0|1|2|3] */
0		/* datové čáry pro celočíselné délky trvání dne? [0|1] */
0		/* počítat rysky od východu slunce? [0|1] */
0		/* počítat rysky od západu slunce předchozího dne? [0|1] */
0		/* rysky temporálních hodin? [0|1] */
1		/* počítat analemu, pro 12 h nebo všechny hodinové rysky? [0|1|2] */
0		/* část analemy: celá, letní nebo zimní? [0|1|2] */
0		/* azimutální úsečky? [0|1] */
0		/* výškové křivky? [0|1] */

EOF
;
  close(IN);

  open(IN,">napis.plt");
  printf(IN "set label \"a = %5.1f°\" at graph 0.85, graph 0.9 font \"Luxi, 32\" front", $az);
  close(IN);

  open(IN,">view.plt");
  printf(IN "set view 75, 360-20, 1.25, 0.75");
  close(IN);

  `./sh.sh`;
  $tmp = sprintf("%04d", $i);
  `mv sh.png png/sh_azimut_j_3d$tmp.png`;
}

exit;

