def get_linear_fits(method='cifelli'):
    """
    Get coefficients from linear fits to Zdp and Zh
    method = str, Reference for method
        'cifelli': Cifelli et al. (2002) - Amazon
        'cr1995': Carey and Rutledge (1995) - Colorado
        Default: Carey and Rutledge (2000) - Darwin
    """
    if method == 'cifelli':
        fit_a = 0.8178
        fit_b = 12.5088
    elif method == 'cr1995':
        fit_a = 0.9091
        fit_b = 8.5091
    else:
        fit_a = 0.77
        fit_b = 14.0
    return fit_a, fit_b
