In [1]:

    

import pandas as pd
import pmagpy.ipmag as ipmag
import numpy as np
from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
plt.rcParams['image.cmap'] = 'viridis'
%matplotlib inline


    

In [2]:

    

GPMDB = pd.read_excel('gpmdb_beta_20150811.xlsx')


    

In [3]:

    

GPMDB_Australia = GPMDB.loc[GPMDB['CONTINENT']=='Australia']
GPMDB_Australia_50Ma = GPMDB_Australia.ix[(GPMDB_Australia['HIGHAGE']<65) &#Only choose Cenozoic poles
                                           (GPMDB_Australia['LOWAGE']>1) & # Don't want archaeomagnetic stuff
                                           (GPMDB_Australia['YEAR']>1984) &# Choose more recent results (relatively)
                                           (GPMDB_Australia['N'] >= 20) & # Choose results with greater N
                                           (GPMDB_Australia['AUTHORS'] != 'Smith,M.L., Pillans,B.J., McQueen,K.G.') &#This has no age control
                                           (GPMDB_Australia['AUTHORS'] != 'Hill, M.J., Gratton, M.N., Shaw, J.')#Something is fishy about this study...
                                           ]
GPMDB_Australia_50Ma.to_csv('GPMDB_Australia_50Ma.csv')
GPMDB_Australia_50Ma[['AUTHORS', 'YEAR','PLAT','PLONG','HIGHAGE','LOWAGE', 'N']]


    

    
Out[3]:






  

    

      

      
AUTHORS
      
YEAR
      
PLAT
      
PLONG
      
HIGHAGE
      
LOWAGE
      
N
    
  
  

    

      
139
      
Idnurm,M.
      
1985
      
68.4
      
298.7
      
28
      
23
      
34.0
    
    

      
140
      
Idnurm,M.
      
1985
      
77.2
      
303.5
      
14
      
10
      
48.0
    
    

      
141
      
Idnurm,M.
      
1985
      
83.2
      
283.6
      
4
      
2
      
58.0
    
    

      
142
      
Idnurm,M.
      
1985
      
61.7
      
298.4
      
62
      
58
      
20.0
    
    

      
7103
      
Idnurm,M.
      
1994
      
65.5
      
292.5
      
39
      
34
      
66.0
    
  

In [4]:

    

m = Basemap(projection='ortho',lat_0=-45,lon_0=260-180.,resolution='c',area_thresh=50000)
plt.figure(figsize=(6, 6))
m.drawcoastlines(linewidth=0.25)
m.fillcontinents(color='bisque',lake_color='white',zorder=1)
m.drawmapboundary(fill_color='white')
m.drawmeridians(np.arange(0,360,30))
m.drawparallels(np.arange(-90,90,30))

for n in GPMDB_Australia_50Ma.index:
    pole = ipmag.plot_pole_colorbar(m,GPMDB_Australia_50Ma['PLONG'][n]-180.,
                             -GPMDB_Australia_50Ma['PLAT'][n],
                             GPMDB_Australia_50Ma['DP'][n],
                             GPMDB_Australia_50Ma['DP'][n]
                             -GPMDB_Australia_50Ma['HIGHAGE'][n],
                             -65,0,markersize=80, color="k", alpha=1)
m.colorbar(pole,location='bottom',pad="5%",label='Age of magnetization (Ma)')
plt.show()


    

In [ ]: