Read xvg file

read_xvg

import numpy as np
import sys

def read_xvg(fname):
    """Read first two columns of data from xvg file fname

    Returns numpy array of data
    """
    skip = 0
    with open(fname, 'r') as f:
        for i, line in enumerate(f):
            if not line.startswith(('#', '@')):
                skip = i
                break

    return np.genfromtxt(fname, skip_header=skip, usecols=(0, 1))