#!/usr/bin/perl

$p = $ARGV[3];

open(F1, "$ARGV[0]");

while(<F1>){
    split;
    next if ($_[$p] eq "");
    push(@lag, $_[$p]);
    shift(@lag) if ($#lag > 2);
    if($#lag==2){
	$proj = &proj(@lag);
	$r = sqrt($lag[1]**2 + $lag[0]**2 +$lag[2]**2);
	$theta = atan2($_[1], $_[2]);
	$theta *= 12;
	print  "$proj $r $theta";
	print " red\n";
    }
}

open(F2, "$ARGV[1]");

while(<F2>){
    split;
    next if ($_[$p] eq "");
    push(@lag, $_[$p]);
    shift(@lag) if ($#lag > 2);
    if($#lag==2){
	$proj = &proj(@lag);
	$r = sqrt($lag[1]**2 + $lag[0]**2 +$lag[2]**2);
	$theta = atan2($_[1], $_[2]);
	$theta *= 12;
	print  "$proj $r $theta";
	print " blue\n";
    }
}

open(F3, "$ARGV[2]");

while(<F3>){
    split;
    next if ($_[$p] eq "");
    push(@lag, $_[$p]);
    shift(@lag) if ($#lag > 2);
    if($#lag==2){
	$proj = &proj(@lag);
	$r = sqrt($lag[1]**2 + $lag[0]**2 +$lag[2]**2);
	$theta = atan2($_[1], $_[2]);
	$theta *= 12;
	print  "$proj $r $theta";
	print " black\n";
    }
}


sub proj {
    ($x, $y, $z) = @_;

    $mag = sqrt(1/3);
    $proj = 0;
    $proj = ($x+$y+$z)*$mag*$mag;
    $dmag = (($x-$proj)*($x-$proj))+(($y-$proj)*($y-$proj))+
	(($z-$proj)*($z-$proj));
    return sqrt($dmag);
}