Read xvg file

read_xvg

import numpy as np
import sys

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

    Returns numpy array of data
    """
    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)