#!/usr/bin/env python
# coding: utf-8

# Geomagnetic Indices Plotting
# ==
# ***

# You can easily load data and make stackplots of Kyoto Geomagnetic indices.

# In[1]:


get_ipython().run_line_magic('pylab', 'inline')
import datetime as dt
import numpy as np
import matplotlib.pyplot as mp

from davitpy import pydarn
from davitpy import gme


# In[2]:


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


# ## AE/AL/AU/AO

# In[3]:


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


# 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'))


# ## SYMH/SYMD/ASYH/ASYD

# 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)


# 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'])


# ## DST

# In[8]:


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


# In[9]:


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


# In[ ]: