def load_rv_data(filename, indep, dep, indweight=None, dir='./'):

    """
    load dictionary with rv data.
    """

    if '/' in filename:
        path, filename = os.path.split(filename)
    else:
        path = dir

    load_file = os.path.join(path, filename)
    rvdata = np.loadtxt(load_file)

    d ={}
    d['phoebe_rv_time'] = rvdata[:,0]
    d['phoebe_rv_vel'] = rvdata[:,1]
    ncol = len(rvdata[0])

    if indweight=="Standard deviation":

        if ncol >= 3:
            d['phoebe_rv_sigmarv'] = rvdata[:,2]
        else:
            logger.warning('A sigma column is mentioned in the .phoebe file but is not present in the rv data file')
    elif indweight =="Standard weight":
                if ncol >= 3:
                    sigma = np.sqrt(1/rvdata[:,2])
                    d['phoebe_rv_sigmarv'] = sigma
                    logger.warning('Standard weight has been converted to Standard deviation.')
    else:
        logger.warning('Phoebe 2 currently only supports standard deviaton')

    return d
def load_rv_data ( filename , indep , dep , indweight = None , dir = './' ) : if '/' in filename : path , filename = os . path . split ( filename ) else : path = dir load_file = os . path . join ( path , filename ) rvdata = np . loadtxt ( load_file ) d = { } d [ 'phoebe_rv_time' ] = rvdata [ : , 0 ] d [ 'phoebe_rv_vel' ] = rvdata [ : , 1 ] ncol = len ( rvdata [ 0 ] ) if indweight == "Standard deviation" : if ncol >= 3 : d [ 'phoebe_rv_sigmarv' ] = rvdata [ : , 2 ] else : logger . warning ( 'A sigma column is mentioned in the .phoebe file but is not present in the rv data file' ) elif indweight == "Standard weight" : if ncol >= 3 : sigma = np . sqrt ( 1 / rvdata [ : , 2 ] ) d [ 'phoebe_rv_sigmarv' ] = sigma logger . warning ( 'Standard weight has been converted to Standard deviation.' ) else : logger . warning ( 'Phoebe 2 currently only supports standard deviaton' ) return d
load dictionary with rv data.