#!/bin/perl

do "parameters.pl";


# Open the file for reading.
open DATA, "Btoroidal_dp.csv" or die "can't open file $!";
my @data = <DATA>;
close (DATA);

$PlasmaStartRow=`cat PlasmaStartRow`;
$PlasmaEndRow=`cat PlasmaEndRow`;
$PlasmaTimeLength=`cat PlasmaTimeLength`;


$Total=0;
for ($i = $PlasmaStartRow; $i < $PlasmaEndRow; ++$i)  
{
    @row = split(';',$data[$i]);
    $Total+=$row[1]*$DeltaTs;    
}

open RESULT, ">MeanBt";
    printf RESULT "%4.3f", $Total/$PlasmaTimeLength*1000;
close RESULT;


if($Total<0){
    system("echo 1 > ReversedBt;");
}#reverse polarity
else 
{
    system("echo 0 > ReversedBt;");
}


exit;