def METHOD_NAME(shape, calibrated_center, r_list, r_width):
    R = core.radial_grid(calibrated_center, shape)
    II = np.zeros_like(R)
    for r in r_list:
        tmp = 100 * np.exp(-(((R - r) / r_width) ** 2))
        II += tmp
    return II