from glob import glob
rcParams['axes.labelsize'] =10
rcParams['xtick.labelsize']=10
rcParams['ytick.labelsize']=10
rcParams['legend.fontsize']=10
rcParams['font.size']=10
rc('font',**{'family':'sans-serif','sans-serif':['Arial']})
#mpl.rc('text', usetex=True)
m=Basemap(projection='stere', lat_0=74, lon_0=-90,llcrnrlon=-150, llcrnrlat=58,urcrnrlon=10, urcrnrlat=72)
rawdatapath='../../../DATA/'
figpath='./Figures/'
xptsW, yptsW, xvel, yvel, wind_speed = ro.get_era_winds(m, rawdatapath, 2009, 2014, 0, 3, 75)
xpts_all=[]
ypts_all=[]
for year in xrange(2009, 2014+1, 1):
print year
lonsT, latsT = ro.calc_icebridge_flights(year,rawdatapath, 'GR')
xptsT, yptsT=m(lonsT, latsT)
xpts_all.append(xptsT)
ypts_all.append(yptsT)
#xpts_all, ypts_all = ro.calc_icebridge_flights_years(m,2009,2014, 'GR')
#xpts_type, ypts_type, ice_type_mean, ice_type, latsT = get_mean_icetype(m, 1)
