#!/bin/bash
grep ucvol= $1 | awk '{print $5}' > tmp1
grep -a ETOT $1 | awk '{if ($2==1 && FNR!=1){print prevL;} prevL=$1 "  " $2 "  " $3;}END{print prevL;}' | awk '{print $3}' > tmp2
paste tmp1 tmp2 > EvsUcvol

E0=`cat EvsUcvol | awk '{if(FNR == 0){min=$2}if($2<min){min=$2}}END{print min}'`
V0=`cat EvsUcvol | awk '{if(FNR == 0){min=$2;minV=$1}if($2<min){min=$2;minV=$1}}END{print minV}'`


cat > EvsUcvol.plt << EOF
f(x)=E0+9.*V0*B0/16.*(((V0/x)**(2./3.)-1)**3.*B0p+((V0/x)**(2./3.)-1)**2.*(6.-4.*(V0/x)**(2./3.)))
B0p=3.
B0=0.0004
V0=$V0
E0=$E0
fit f(x) "EvsUcvol" via E0,V0,B0,B0p
plot "EvsUcvol" w lp, f(x)
print B0*29421.0119374
EOF

gnuplot -persist EvsUcvol.plt

rm -f tmp1
rm -f tmp2