%matplotlib inline import yt ds = yt.load("A2052_merged_0.3-2_match-core_tmap_bgecorr.fits", slave_files=["A2052_core_tmap_b1_m2000_.fits"], override_fields=[("flux","counts/s/pixel"),("projected_temperature","keV")]) ds.field_list slc = yt.SlicePlot(ds, "z", ["flux","projected_temperature"], origin="native") slc.set_log("flux",True) slc.set_width(250.) slc.show() dd = ds.all_data() dd["flux"].units dd["flux"].units.latex_representation()