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)