#!/bin/bash

#python np.py $@ >data
#gnuplot Ef.graph && mv graph.png Ef.png && gnuplot cond_part.graph && mv graph.png cond_part.png
#gwenview Ef.png >/dev/null 2>&1
#gwenview cond_part.png >/dev/null 2>&1
#rm Ef.png
#rm cond_part.png

rm -rf images/
mkdir images

A="donor acceptor"
VALUES="1e16 5e16 1e17 5e17 1e18"
for y in $A ; do
    for x in $VALUES ; do
        python np.py $y $x >data
    gnuplot <<EOF
#!/usr/bin/gnuplot -persist
set title "Ef(T) $y $x"
set xlabel "1/kT, kT in eV"
set ylabel "Ef(T)"
set grid
set terminal png
set output "images/$y-$x-Ef.png"
set term png size 960, 640
plot "data" u 1:3 w l
EOF
        gnuplot <<EOF
#!/usr/bin/gnuplot -persist
set title "conductive_particles(T) $y $x"
set xlabel "1/kT, kT in eV"
set ylabel "log(conductive_particles)"
set grid
set terminal png
set output "images/$y-$x-cond_part.png"
set term png size 960, 640
plot "data" u 1:2 w l
EOF
    done
done
rm data