ts=yt.load(["/lunarc/nobackup/users/samvad/FINAL-50-0.5/output/output_00029/info_00029.txt","/lunarc/nobackup/users/samvad/FINAL-50-0.5/output/output_00037/info_00037.txt","/lunarc/nobackup/users/samvad/FINAL-50-0.5/output/output_00042/info_00042.txt"])

for ds in ts:

        ds.define_unit("H", (1.674*10**(-24), "g"))

        ad=ds.all_data()

        dens=ad.cut_region(["obj[('gas','density')].in_units('H/cm**3') > 0"])

      profiles.append(yt.create_profile(ad, ("gas","density"), [("gas","cell_mass")], weight_field=None, fractional=True)) 

My data has highest refinement level of 15 from simulations. However, I'm not sure if this is the same level at which the profiles are being plotted