import sys
import yt
import yt.units as u
def doit(plotfile):
ds = yt.load(plotfile)
ds.periodicity = (False, False, False)
ad = ds.all_data()
plot = yt.ProfilePlot(ad, "x", "temperature")
plot.set_unit("temperature", "K")
plot.set_unit("x", "m")
plot.set_log('x', False)
plot.set_xlim(1e-20, 450.0)
pid = plotfile.split("plt")[1]
plot.save("1d_detonation_{}.png".format(pid))
if __name__ == "__main__":
# Choose a field
plotfile = ""
try:
plotfile = sys.argv[1]
except:
sys.exit("ERROR: no plotfile specified")
for plt in sys.argv[1:]:
plotfile = plt
doit(plotfile)