def METHOD_NAME(center,profile_times_years,input_prefix):
    def _my_radial_velocity(field,data):
        xv = data["gas","velocity_x"]
        yv = data["gas","velocity_y"]
        zv = data["gas","velocity_z"]
        x_hat = data["enzoe","x"] - center[0]
        y_hat = data["enzoe","y"] - center[1]
        z_hat = data["enzoe","z"] - center[2]
        r = np.sqrt(x_hat*x_hat+y_hat*y_hat+z_hat*z_hat)
        x_hat /= r
        y_hat /= r
        z_hat /= r
        return xv*x_hat + yv*y_hat + zv*z_hat
    def _mass_flux(field,data):
        xv = data["gas","velocity_x"]
        yv = data["gas","velocity_y"]
        zv = data["gas","velocity_z"]
        x_hat = data["enzoe","x"] - center[0]
        y_hat = data["enzoe","y"] - center[1]
        z_hat = data["enzoe","z"] - center[2]
        r = np.sqrt(x_hat*x_hat+y_hat*y_hat+z_hat*z_hat)
        x_hat /= r
        y_hat /= r
        z_hat /= r
        rad_vel = xv*x_hat + yv*y_hat + zv*z_hat
        density = data["gas","density"]
        return density*rad_vel
    ds_pattern = f"{input_prefix}_????/{input_prefix}_????.block_list"
    time_series = yt.load(ds_pattern)
    time_list = [ds.current_time.to("yr") for ds in time_series]
    ind_list = [find_snapshot_index(t,time_list) for t in profile_times_years]
    yt.add_field(("gas","my_radial_velocity"),
             sampling_type = "cell",
             function=_my_radial_velocity,
             units="cm/s",
             take_log=False,
             validators=[ValidateParameter(['center'])])
    yt.add_field(("gas","mass_flux"),
             sampling_type = "cell",
             function=_mass_flux,
             units="g/cm**2/s",
             take_log=False,
             validators=[ValidateParameter(['center'])])
    profiles = [get_profiles(time_series,time_list,ind,center)
                for ind in ind_list if ind >= 0]
    fig,ax = plt.subplots(nrows = 3,sharex = True)
    fig.set_figheight(12)
    fig.set_figwidth(6)
    for profile in profiles:
        ax[0].plot(profile["radius"],profile["density"])
        ax[0].set_ylim(1.0e-19,1.0e-14)
        ax[0].set_ylabel(r"$\rho\,(g\,cm^{-3})$")
        ax[0].set_yscale('log')
        ax[0].axvline(x = 4.0*profile["cell_width"],ls = "--",color = "k",lw = 0.5)
        ax[1].plot(profile["radius"],profile["radial_velocity"])
        ax[1].set_ylim(-5,1)
        ax[1].set_ylabel(r"$v_r \, (km\,s^{-1}$")
        ax[1].axvline(x = 4.0*profile["cell_width"],ls = "--",color = "k",lw = 0.5)
        ax[2].plot(profile["radius"],profile["mass_flux"],
                   label = r"$%1.1f \times 10^3\,\mathrm{yrs}$" %(profile["time"]/1.0e3))
        ax[2].set_yscale('log')
        ax[2].set_ylim(1.0e-6,1.0e-2)
        ax[2].set_ylabel(r"$\dot{M} \, (M_{\odot}/\mathrm{yr})$")
        ax[2].set_xlabel(r"$r\,cm$")
        ax[2].set_xlim(1.0e15,1.0e17)
        ax[2].set_xscale('log')
        ax[2].axvline(x = 4.0*profile["cell_width"],ls = "--",color = "k",lw = 0.5)
        ax[2].axhline(y = 1.5e-4,ls = "--",color = "k")
    ax[2].legend(bbox_to_anchor = (0.0,-0.2),loc = "upper left",ncol = 3)
    fig.savefig("radial_profiles.png")