#!/usr/bin/env python

#
# $Id: pts.py,v 1.2 2009/04/28 12:04:43 dintrans Exp $
# Plot the data/time_series.dat file: bmx, bmy, bmz and brms
#

import numpy as N
import pylab as P
import pencil as pc

par2=pc.read_param(quiet=True, param2=True)

ts=pc.read_ts(plot_data=False)
P.subplot(211)
P.plot(ts.t, ts.uxbm, label='uxbm')
P.xlabel('time')

P.subplot(212)
P.plot(ts.t, ts.brms, label='brms')
P.xlabel('time')

P.show()