#!/bin/bash

#START=21000
#START=19898
START=8000
END=21097
BASE=/golem/database/operation/shots/;



echo "All comments"
#for i in `seq $START $END`; do  echo -ne $i "\t"; cat /golem/database/operation/shots/$i/basicdiagn/Plasma oddelovac /golem/database/operation/shots/$i/Ub_limit oddelovac /golem/database/operation/shots/$i/Ucd_limit oddelovac /golem/database/operation/shots/$i/pressure_setting_mPa|tr "\n" " "; echo "";done |tee Allshots
#grep "1   800.0  400.0" Allshots |tee SelectedShots
echo "Table generation"
#for i in `awk '{print $1}' SelectedShots`; do cat $BASE/$i/date; echo -ne "\t" $i "\t"; echo `cat  $BASE/$i/basicdiagn/PlasmaTimeLength`*1000|bc|tr "\n" " ";echo -ne "\t"; cat  $BASE/$i/basicdiagn/ElectronTempMax|tr "\n" " ";echo -ne "\t"; echo `cat  $BASE/$i/basicdiagn/StateEqElectronDensity|sed -e 's/[eE]+*/\\*10\\^/'`*0.000000000000000001|bc|tr "\n" " "|sed -e 's/000000000000//';echo -ne "\t"; echo `cat  $BASE/$i/basicdiagn/IplaMean`*0.001|bc|tr "\n" " ";echo -ne "\t";echo ""; done|column -t|tee WholeTablePure

echo "Manualni oprava excesu"
grep -vwFf /GW/DischargeAnalysis/Fails/index.page WholeTablePure|grep -v nan > WholeTable

awk '{print $2 " " $3}' WholeTable > rand_t.dat
gnuplot -e "Titul='Doba vyboje t_{pl} [ms]'" Statistika.gp
convert -rotate 90 statistic.eps DobaVyboje.pdf
convert -rotate 90 statistic.eps DobaVyboje.jpg
sort -n -k3 WholeTable > SortedDobaVyboje

awk '{print $2 " " $4}' WholeTable > rand_t.dat
gnuplot -e "Titul='Elektronova teplota T_e [eV]'" Statistika.gp
convert -rotate 90 statistic.eps ElektronovaTeplota.pdf
convert -rotate 90 statistic.eps ElektronovaTeplota.jpg
sort -n -k4 WholeTable > SortedElTeplota

awk '{print $2 " " $5}' WholeTable > rand_t.dat
gnuplot -e "Titul='Elektronova hustota n_e [x10^{18}/m^3]'" Statistika.gp
convert -rotate 90 statistic.eps ElektronovaHustota.pdf
convert -rotate 90 statistic.eps ElektronovaHustota.jpg
sort -n -k5 WholeTable > SortedElHustota

awk '{print $2 " " $6}' WholeTable > rand_t.dat
gnuplot -e "Titul='Proud Plazmatem I_{pl} [kA]'" Statistika.gp
convert -rotate 90 statistic.eps ProudPlazmatem.pdf
convert -rotate 90 statistic.eps ProudPlazmatem.jpg
sort -n -k6 WholeTable > SortedProudPlazmatem


exit