def get_data_Caffrey1998():
    vd = array([11.0, 2.70, 1.30, 0.14, 0.02, 0.02, 0.01, 0.006, 0.007, 0.009, 0.01])
    dp = array([60.4, 30.0, 21.2, 6.9, 2.4, 1.34, 0.74, 0.40, 0.22, 0.12, 0.07])

    vd_err = array([ [0.04, 0.14, 0.02, 0.01, 0.02, 0.02, 0.01, 0.004, 0.005, 0.005, 0.01] ,
                     [0.04, 0.14, 0.02, 0.01, 0.02, 0.02, 0.01, 0.004, 0.005, 0.005, 0.01] ])
    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, 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_Zufall1998():
    vd = array([0.03, 0.007, 0.01, 0.05, 0.5, 2.5, 6.0])
    dp = array([0.35, 0.80, 1.50, 3.00, 6.00, 15.0, 25.0])

    vd_err = array([ [0.025, 0.0063, 0.0087, 0.041, 0.41, 1.7, 5.2] , [0.13, 0.038, 0.07, 0.15, 1.5, 5.5, 24.0] ])
    dp_err = array([ [0.000, 0.0000, 0.0000, 0.000, 0.00, 0.0, 0.0] , [0.00, 0.000, 0.00, 0.00, 0.0, 0.0, 0.0] ])
 
    return vd, vd_err, dp, dp_err

def get_data_Moller1970():
    vd = array([0.05, 0.045, 0.025, 0.024, 0.018, 0.01, 0.02, 0.01, 0.015, 0.01, 0.016, 0.014, 0.015, 0.03])
    dp = array([0.038, 0.04, 0.052, 0.06, 0.1, 0.14, 0.154, 0.17, 0.3, 0.48, 0.6, 1.4, 1.6, 1.6])

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

def get_data_Zhang2014_15():
    vd = array([2.0, 2.7, 2.2, 2.3, 3.3, 4.2, 5.5, 6.8, 12.8])
    dp = array([1.0, 2.3, 4.0, 7.5, 12.5, 17.5, 22.5, 27.5, 40.0])

    vd_err = array([1.2, 0.99, 0.5, 0.4, 0.8, 0.5, 0.9, 0.8, 0.7])
    dp_err = array([0.5, 0.75, 1.0, 2.5, 2.5, 2.5, 2.5, 2.5, 10.0])

    return vd, vd_err, dp, dp_err

def get_data_Zhang2014_37():
    vd = array([8.0, 5.4, 5.7, 4.4, 6.1, 8.0, 10.1, 11.5, 18.0])
    dp = array([1.0, 2.3, 4.0, 7.5, 12.5, 17.5, 22.5, 27.5, 40.0])

    vd_err = array([7.5, 4.1, 3.3, 0.9, 2.3, 2.3, 2.3, 3.3, 3.2])
    dp_err = array([0.5, 0.75, 1.0, 2.5, 2.5, 2.5, 2.5, 2.5, 10.0])

    return vd, vd_err, dp, dp_err

def get_data_Zhang2014_57():
    vd = array([6.9, 11.5, 10.1, 12.1, 9.8, 12.3, 15.3, 17.7, 24.1])
    dp = array([1.0, 2.3, 4.0, 7.5, 12.5, 17.5, 22.5, 27.5, 40.0])

    vd_err = array([6.5, 10.7, 5.8, 5.0, 3.8, 5.5, 6.9, 6.3, 4.6])
    dp_err = array([0.5, 0.75, 1.0, 2.5, 2.5, 2.5, 2.5, 2.5, 10.0])

    return vd, vd_err, dp, dp_err

def get_data_Sievering1981():
    vd = array([0.4, 0.2, 0.7])
    dp = array([0.15, 0.45, 1.0])

    vd_err = array([ [0.1, 0.1, 0.1], [0.1, 0.1, 0.1] ])
    dp_err = array([ [0.04, 0.15, 0.4], [0.05, 0.15, 0.5] ])