From d8cd5af60b627992bf201bd16854ad24397cf8e7 Mon Sep 17 00:00:00 2001 From: mer0m Date: Fri, 20 Jan 2017 10:51:13 +0100 Subject: [PATCH] - --- allanplot.py | 21 +++++++++++---------- 1 file changed, 11 insertions(+), 10 deletions(-) diff --git a/allanplot.py b/allanplot.py index ad7e4b1..33bfc15 100644 --- a/allanplot.py +++ b/allanplot.py @@ -24,14 +24,14 @@ def parse(): dest='filename', default=FILENAME, help='File(s) to import (default '+FILENAME+')') - + parser.add_argument('-c', nargs = '+', action='store', dest='columns', default=COLUMNS, help='Columns to import (default '+str(COLUMNS)+')') - + parser.add_argument('-r', action='store', dest='rate', @@ -61,30 +61,31 @@ def main(): filename = filename.split() for name in filename: list_files.extend(sorted(glob.glob(name))) - + data = [] - + for f in list_files: with open(f, 'r') as dest_f: data_iter = csv.reader(dest_f, delimiter = '\t', quotechar = '"') temp_data = [value for value in data_iter] data.extend(temp_data) - + data = numpy.asarray(data, dtype = float) - + del(temp_data, list_files, value, f) - - g = Gnuplot.Gnuplot(persist = 1) + + g = Gnuplot.Gnuplot(persist = 1) g('set logscale x') g('set logscale y') g('set grid') + g('set format y "%0.1e"') g.xlabel('Tau (s)') g.ylabel('Adev') - + for i in columns: (tau2, ad, ade, adn) = allantools.adev(data[:,i], rate=rate, data_type="freq", taus='decade') g.replot(Gnuplot.Data(tau2, ad, ade, with_='yerrorbars', title='%s : %s'%(name, str(i)))) - + except Exception as ex: print 'Oups '+str(ex) print 'Program ending\n' -- 2.16.4