add min and max data limits

Author: zingale

analysis/plotcompact.py

@@ -7,7 +7,8 @@
 
 
 # plot an output file using the solver's dovis script
-def makeplot(plotfile, variable, outfile):
+
+def makeplot(plotfile, variable, outfile, vmin=None, vmax=None):
 
     sim = io.read(plotfile)
     myd = sim.cc_data
@@ -28,11 +29,18 @@ def makeplot(plotfile, variable, outfile):
     else:
         v = myd.get_var(variable)
 
+    if vmin is None:
+        vmin = v.v().min()
+
+    if vmax is None:
+        vmax = v.v().max()
+
     plt.figure(num=1, figsize=(6.0, 6.0), dpi=100, facecolor='w')
 
     plt.imshow(np.transpose(v[myg.ilo:myg.ihi+1, myg.jlo:myg.jhi+1]),
                interpolation="nearest", origin="lower",
-               extent=[myg.xmin, myg.xmax, myg.ymin, myg.ymax])
+               extent=[myg.xmin, myg.xmax, myg.ymin, myg.ymax],
+               vmin=vmin, vmax=vmax)
 
     plt.axis("off")
     plt.subplots_adjust(bottom=0.0, top=1.0, left=0.0, right=1.0)
@@ -45,6 +53,10 @@ def get_args():
 
     parser.add_argument("-o", type=str, default="plot.png",
                         metavar="plot.png", help="output file name")
+    parser.add_argument("-m", type=float, default=None,
+                        metavar="data-min", help="minimum data value to show")
+    parser.add_argument("-M", type=float, default=None,
+                        metavar="data-max", help="maximum data value to show")
     parser.add_argument("plotfile", type=str, nargs=1,
                         help="the plotfile you wish to plot")
     parser.add_argument("variable", type=str, nargs=1,
@@ -59,4 +71,4 @@ def get_args():
 
     args = get_args()
 
-    makeplot(args.plotfile[0], args.variable[0], args.o)
+    makeplot(args.plotfile[0], args.variable[0], args.o, vmin=args.m, vmax=args.M)