def METHOD_NAME(filename):
    ds = yt.load(filename)
    ad = ds.all_data()
    profile = yt.create_profile(ad, args.binfields, args.fields,
                                n_bins=tuple(args.nbins),
                                weight_field=None)
    print(profile.keys())
    print(profile["enucdot"].shape)
    print(profile["y"].shape)
    for y, e in zip(profile["y"], profile["enucdot"]):
        print('({}, {})'.format(y, e))
    exit()
    return ds.time, emax, elocmax, efwhm