import pylab, scipy

array = raw_input('Select the array you wish to plot: ')

dataarray = scipy.load('%(array)s'%vars())

pylab.hist(dataarray, bins = 40)
pylab.show()

vMin = input('Select the lower limit: ')
vMax = input('Select the upper limit: ')

pylab.title('%(array)s'%vars())
pylab.pcolormesh(dataarray, vmin = vMin, vmax = vMax)
pylab.colorbar()
pylab.show()