def get_data_Allen1991():
    vd = array([0.1])
    dp = array([1.0])

    vd_err = array([ [0.03] , [0.03] ])
    dp_err = array([ [0.9]  , [1.0]  ])
    
    return vd, vd_err, dp, dp_err

def get_data_Vong2004():
    vd = array([0.16, 0.30, 0.44])
    dp = array([0.52, 0.52, 0.52])

    vd_err = array([0.07, 0.07, 0.07])
    dp_err = array([0.0 , 0.0 , 0.0 ])
 
    return vd, vd_err, dp, dp_err

def get_data_Clough1975():
    vd = array([0.024, 0.015])
    dp = array([0.5,   0.5])

    vd_err = array([0.0, 0.0])
    dp_err = array([0.0, 0.0])
    
    return vd, vd_err, dp, dp_err

def get_data_Hummelshoj1994():
    vd = array([0.16, 0.20, 0.28])
    dp = array([0.10, 0.24, 0.60])

    vd_err = array([0.0 , 0.0 , 0.0 ])
    dp_err = array([ [0.034 , 0.09 , 0.2 ], [0.051, 0.16, 0.3] ])
 
    return vd, vd_err, dp, dp_err

def get_data_Gallagher1988():
    vd = array([0.5, 4.0, 2.9, 4.6, 5.8, 5.1, 7.1, 7.9, 8.9, 8.4, 6.6, 6.4, 4.9, 5.1])
    dp = array([5.0, 7.0, 9.0, 11., 13., 15., 17., 19., 21., 23., 25., 27., 29., 31.])

    vd_err = array([0.33 , 0.32 , 0.23, 0.23, 0.12, 0.20, 0.28, 0.32, 0.36, 0.59, 1.1, 2.7, 2.9, 5.0])
    dp_err = array([0.0 , 0.0 , 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0])
 
    return vd, vd_err, dp, dp_err

def get_data_Chamberlain1967():
    vd = array([7.6, 3.9, 1.2, 0.05, 0.08, 4.9, 3.2, 0.8, 0.09, 0.03, 0.03, 3.3, 2.7, 0.16, 0.019, 0.02])
    dp = array([32.0, 19.0, 5.0, 1.0, 0.08, 32.0, 19.0, 5.0, 2.0, 1.0, 0.08, 32.0, 19.0, 5.0, 1.0, 0.08])

    vd_err = zeros(16)
    dp_err = zeros(16)
 
    return vd, vd_err, dp, dp_err

def get_data_Nemitz2002():
    vd = array([0.003, 0.005, 0.01, 0.03, 0.01, 0.01, 0.05, 0.06, 0.004, 0.04, 0.07, 0.24, 0.03, 0.06, 0.12, 0.27])
    dp = array([0.125, 0.155, 0.23, 0.425,0.125,0.155,0.23, 0.425,0.125, 0.155,0.23, 0.425,0.125,0.155,0.23, 0.425])

    vd_err = zeros(16)
    dp_err = array([0.005, 0.005, 0.01, 0.025, 0.005, 0.005, 0.01, 0.025, 0.005, 0.005, 0.01, 0.025, 0.005, 0.005, 0.01, 0.025])
 
    return vd, vd_err, dp, dp_err

def get_data_Wesely1985():
    vd = array([0.22])
    dp = array([1.0])

    vd_err = array([0.06])
    dp_err = array([ [0.9] , [1.0] ])

    return vd, vd_err, dp, dp_err

def get_data_Wesely1977():
    vd = array([0.55])
    dp = array([0.07])

    vd_err = array([0.45])
    dp_err = array([ [0.02], [0.03] ])

    return vd, vd_err, dp, dp_err

def get_data_Sievering1982():
    vd = array([0.38])
    dp = array([0.25])
    
    vd_err = array([0.29])
    dp_err = array([ [0.10], [0.05] ])

    return vd, vd_err, dp, dp_err

def get_data_Neumann1985():
    vd = array([0.05])
    dp = array([0.3])
   
    vd_err = array([0.0])
    dp_err = array([ [0.2], [0.2] ])

    return vd, vd_err, dp, dp_err

def get_data_Sievering1988():
    vd = array([0.37, 1.19])
    dp = array([0.47, 0.47])

    vd_err = array([0.04, 0.18])
    dp_err = array([ [0.0, 0.0], [0.0, 0.0] ])

    return vd, vd_err, dp, dp_err

def get_data_Connan2018():
    vd = array([0.17, 0.073, 0.23, 0.26, 0.059, 0.19, 2.1])
    dp = array([0.24, 0.27, 0.50, 0.60, 1.2, 3., 7.8])

    vd_err = array([0.16, 0.053, 0.13, 0.125, 0.02, 0.105, 0.59])
    dp_err = zeros(7)