def bound_maker_erfrecterf(amplitude_bounds,translational_offset_bounds,stddev_bounds,vertical_offset_bounds,number_erfs):
    """
    Create tuple with lower and upper bounds to be used in the curve fit
    
    Args:
        amplitude_bounds (tuple): bounds on the amplitudes of the gaussians
        translational_offset_bounds (tuple): bounds on the translational offsets of the gaussians
        stddev_bounds (tuple): bounds on the standard deviations of the gaussians
        vertical_offset_bounds (tuple): bounds on the vertical offset of the gaussians
        number_erfs (int): the number of erf-rect-erf features in the fit
    
    Returns:
        bounds (tuple): lists the bounds on the parameters used in the erf-rect-erf fits
    """
    lower = [amplitude_bounds[0]]*number_erfs + [translational_offset_bounds[0]]*number_erfs + [stddev_bounds[0]]*number_erfs + [translational_offset_bounds[0]]*number_erfs + [stddev_bounds[0]]*number_erfs + [vertical_offset_bounds[0]]
    upper = [amplitude_bounds[1]]*number_erfs + [translational_offset_bounds[1]]*number_erfs + [stddev_bounds[1]]*number_erfs + [translational_offset_bounds[1]]*number_erfs + [stddev_bounds[1]]*number_erfs + [vertical_offset_bounds[1]]       
    bounds = (lower, upper)
    return bounds