Hi, I've encountered a problem when I tried to use yt to visualize simple 1D detonation simulation data generated by Castro. The figure I expected to get is a simple step function just as attachment "1d_detonation_expected.png". However, if I don't set a customized axis limit by hand, it will produce something weird at the origin, which can be seen in the attachment "1d_detonation_weird.png". Following is my simple script to make plots. I wonder where I did wrong thus make this happened. Thank you for answering.

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)


Best, 
Mu-Hung