diff --git a/util/econst.plt b/util/econst.plt
index bd34395..3462b8f 100755
--- a/util/econst.plt
+++ b/util/econst.plt
@@ -7,6 +7,8 @@ fi
 gnuplot -persist <<EOF
 fit $range a*x+b "<grep -h econst $*" u 0:2 via a,b
 fit $range c "<grep -h econst $*" u 0:2 via c
+set xlabel "step"
+set ylabel "energy (Ha)"
 p $range "<grep -h econst $*" u 2 w l, a*x+b, c
 print "Econst_avg=",c,"   dE/dt=",a
 EOF
diff --git a/util/econst_cmp.plt b/util/econst_cmp.plt
index 79894a2..2146682 100755
--- a/util/econst_cmp.plt
+++ b/util/econst_cmp.plt
@@ -5,5 +5,7 @@ then
   shift 2
 fi
 gnuplot -persist <<EOF
+set xlabel "step"
+set ylabel "energy (Ha)"
 plot $range "<grep econst $1" u 2 w l, "<grep econst $2" u 2 w l
 EOF
diff --git a/util/econste.plt b/util/econste.plt
index 5498394..2750621 100755
--- a/util/econste.plt
+++ b/util/econste.plt
@@ -5,5 +5,7 @@ then
   shift 2
 fi
 gnuplot -persist <<EOF
+set xlabel "step"
+set ylabel "energy (Ha)"
 p $range "<grep -h econst $*" u 2 w l, "<grep -h '<etotal>' $*" u 2 w l
 EOF
diff --git a/util/econste_cmp.plt b/util/econste_cmp.plt
index 922cebf..340fa08 100755
--- a/util/econste_cmp.plt
+++ b/util/econste_cmp.plt
@@ -7,6 +7,8 @@ then
   shift 2
 fi
 gnuplot -persist <<EOF
+set xlabel "step"
+set ylabel "energy (Ha)"
 p $range "<grep -h econst $1" u 2 w l, "<grep -h '<etotal>' $1" u 2 w l, \
   "<grep -h econst $2" u 2 w l, "<grep -h '<etotal>' $2" u 2 w l
 
diff --git a/util/econste_cp.plt b/util/econste_cp.plt
index 4bf04f2..b3a8b40 100755
--- a/util/econste_cp.plt
+++ b/util/econste_cp.plt
@@ -5,5 +5,7 @@ then
   shift 2
 fi
 gnuplot -persist <<EOF
+set xlabel "step"
+set ylabel "energy (Ha)"
 p $range  "<grep -h econst $*" u 2 w l, "<grep -h '<etotal>' $*" u 2 w l, "<grep -h '<ekin_ec>' $*" u 2 w l
 EOF
diff --git a/util/etotal.plt b/util/etotal.plt
index c78bb8c..053f4bf 100755
--- a/util/etotal.plt
+++ b/util/etotal.plt
@@ -5,5 +5,7 @@ then
   shift 2
 fi
 gnuplot -persist <<EOF
+set xlabel "step"
+set ylabel "energy (Ha)"
 plot $range "<grep -h '<etotal>' $*" u 2 w l
 EOF
diff --git a/util/etotal_cmp.plt b/util/etotal_cmp.plt
index f489b0c..4d349db 100755
--- a/util/etotal_cmp.plt
+++ b/util/etotal_cmp.plt
@@ -5,5 +5,7 @@ then
   shift 2
 fi
 gnuplot -persist <<EOF
+set xlabel "step"
+set ylabel "energy (Ha)"
 plot $range "<grep '<etotal>' $1" u 2 w l, "<grep '<etotal>' $2" u 2 w l
 EOF
diff --git a/util/etotal_int.plt b/util/etotal_int.plt
index b8a409c..ef25715 100755
--- a/util/etotal_int.plt
+++ b/util/etotal_int.plt
@@ -6,5 +6,7 @@ then
 fi
 gnuplot -persist <<EOF
 set grid
+set xlabel "step"
+set ylabel "energy (Ha)"
 plot $range "<grep -h etotal_int $*" u 2 w l
 EOF
diff --git a/util/etotal_int_cmp.plt b/util/etotal_int_cmp.plt
index 8869c52..3886316 100755
--- a/util/etotal_int_cmp.plt
+++ b/util/etotal_int_cmp.plt
@@ -5,5 +5,7 @@ then
   shift 2
 fi
 gnuplot -persist <<EOF
+set xlabel "step"
+set ylabel "energy (Ha)"
 plot $range "<grep etotal $1" u 2 w l, "<grep etotal $2" u 2 w l
 EOF