Geomagnetic Indices Plotting

In [1]:

    

%pylab inline
import datetime as dt
import numpy as np
import matplotlib.pyplot as mp

from davitpy import pydarn
from davitpy import gme


    

    




Populating the interactive namespace from numpy and matplotlib

In [2]:

    

#Let's set a start and end date...
sTime = dt.datetime(2012,11,19)
eTime = dt.datetime(2012,11,20)


    

In [3]:

    

#Now we read in AE Data...
ae=gme.ind.readAe(sTime,eTime)


    

    




returning a list with 25 records of ae data

In [4]:

    

#You can plot AE in one easy command...
gme.plotting.plotGME(ae)


    

In [5]:

    

#Or make things a little more fancy...
plotsTime = dt.datetime(2010,11,19,8,0)
ploteTime = dt.datetime(2010,11,19,12,0)

myFig = mp.figure(figsize=(9,10))
ax=myFig.add_subplot(311)
gme.plotting.plotGME(ae)
ax=myFig.add_subplot(312)
gme.plotting.plotGME(ae,parameter=('al','au'),NoLegend=True)
ax=myFig.add_subplot(313)
gme.plotting.plotGME(ae,parameter=('al','au','ao'))


    

In [6]:

    

#Now we do the SYM and ASYM indices...
symasy=gme.ind.symasy.readSymAsy(sTime=sTime, eTime=eTime, symh=None, symd=None, asyh=None, asyd=None)


    

    




returning a list with 1441 records of sym/asy data

In [7]:

    

myFig = mp.figure(figsize=(9,10))
ax=myFig.add_subplot(211)
gme.plotting.plotGME(symasy,parameter=['symh','symd'])
ax=myFig.add_subplot(212)
gme.plotting.plotGME(symasy,parameter=['asyh','asyd'])


    

In [8]:

    

#And finally DST...
dst=gme.ind.readDst(sTime=sTime, eTime=eTime)


    

    




returning a list with 25 records of dst data

In [9]:

    

myFig = mp.figure(figsize=(9,5))
gme.plotting.plotGME(dst)


    

In [ ]: