def initial_guess(initial_amplitude,initial_translational_offset,initial_stddev,initial_vertical_offset):
    """
    Create array with amplitude, phase and initial offset to be used in the curve fit
    
    Args:
        initial_amplitude (array): guess for the initial values of the amplitudes of the gaussians
        initial_translational_offset (array): guess for the initial values of the translational offsets of the gaussians
        initial_stddev (array): guess for the initial values of the standard deviations of the gaussians 
        initial_vertical_offset (float): guess for the initial values of the vertical offset 
    
    Returns:
        p0 (array): lists the initial_amplitude, initial_translational_offset, initial_stddev, initial_vertical_offset in the correct format for the curve fit.
    """
    #p0=[]
    #for a,mu,stddev in zip(initial_amplitude,initial_translational_offset,initial_stddev):
    #    p0.append([a,mu,stddev])
    #p0.append(initial_vertical_offset)
    
    p0 = [i for i in initial_amplitude]\
         + [i for i in initial_translational_offset]\
         + [i for i in initial_stddev]\
         + [initial_vertical_offset]
    
    return p0