Hello, I am making profile plots as a function of height from the galactic disk with several variables for a FLASH simulation of a galaxy. For the low values of the height, the profile plot gives me nonsense for several variables. After loading in the data and setting some units, these are the few lines I am using: dd = ds.all_data() d = yt.ProfilePlot(dd, "z", ["magnetic_field_magnitude"], weight_field='ones') d.set_log("magnetic_field_magnitude", False) I then change some labels and save the plot, but these two lines are what's most important. For the first ~1 kpc on the height scale, the plot just has a "zig-zag" pattern for the y-axis, magnetic field magnitude from this code, running from the lower to upper bounds, which is complete nonsense for my data. It then follows a nice curve for the remainder of the plot. Does anyone have any suggestions on how to fix the "zig-zag" at the lower height values? Thanks, Jordan Gardipee