for i in be vr nh an sd
do
	grompp -f ${i}.mdp -c et.gro -p et.top -o et_${i}.tpr
	mdrun -deffnm et_${i} -v -nt 1
	echo 2 | trjconv -f et_${i}.trr -s et_${i}.tpr -o et_${i}.pdb
	echo 1 11 | g_energy -f et_${i}.edr -o et_${i}_en.xvg
	echo -e 'set datafile commentschars "#@&"'"\n"'set term png'"\n"'set output "'${i}'.png"'"\n"'plot "./et_'${i}'_en.xvg" using 1:2,  "./et_'${i}'_en.xvg" using 1:3' |gnuplot
	g_bond -f et_${i}.trr -s et_${i}.tpr -o bond_${i}.xvg -n b.ndx
	echo -e 'set datafile commentschars "#@&"'"\n"'set term png'"\n"'set output "'${i}'_bond.png"'"\n"'plot "./bond_'${i}'.xvg" with boxes' |gnuplot
done
