looks good now

plot.py

@@ -5,21 +5,31 @@
 mencoder "mf://density0*.png" -mf fps=20 -o density.avi -ovc lavc -lavcopts vcodec=mpeg4
 """
 
+import sys
 from cPickle import load
 
+from numpy import array
 from pylab import pcolor, pcolormesh, show, savefig, clf, colorbar, gca
 
-f = open("../plot217to249")
-#f = open("../plot191to215")
-#f = open("../plot157to189")
-#f = open("../plot1to155")
+C_min = 0.174630617373
+C_max = 2.23176607373
+
+files = [
+        "../plot1to155",
+        "../plot157to189",
+        "../plot191to215",
+        "../plot217to249",
+        ]
+file = files[0]
+f = open(file)
 print "load"
 data = load(f)
-X, Y, C = data[-1]
+X, Y, C = data[150]
 print "pcolor"
 print C.shape
-pcolormesh(X, Y, C)
-#colorbar()
+print "min/max", C.min(), C.max()
+pcolormesh(X, Y, C, vmin=C_min, vmax=C_max)
+colorbar()
 gca().set_aspect("equal")
 print "show"
 show()