diff --git a/plot_scripts/testPlot.R b/plot_scripts/testPlot.R
index c78d5e958..244e2a1bb 100755
--- a/plot_scripts/testPlot.R
+++ b/plot_scripts/testPlot.R
@@ -12,7 +12,17 @@ a = read.table(args[1], header=T);
 pkgLimit = sort(unique(a$pkg.clamp.0), decreasing = T)[1]
 
 a$avgNodePower = rowSums(a[, c('avg.pkg.0','avg.pkg.1','avg.dram.0','avg.dram.1')])
+a$coresPerNode = a$cores / a$nodes
 
 sel = which(a$pkg.clamp.0 == pkgLimit)
-
-plot(a$cores[sel], a$avgNodePower[sel])
+coresPerNode = unique(a[sel, 'coresPerNode'])
+sel1 = intersect(sel, which(a$coresPerNode == coresPerNode[1]))
+
+totalEnergy = a$avgNodePower[sel]*a$nodes[sel]*a$avg.time[sel]
+plot(a$cores[sel], totalEnergy, 
+		    main='total energy vs nodes & cores',
+#		    sub=paste(coresPerNode[1], 'cores per node'),
+		    xlab='cores',
+		    ylab='log total energy (joules)',
+		    log='y',
+		    ylim=c(1,10*max(totalEnergy)))