llcrnrlat=58,
            urcrnrlon=10,
            urcrnrlat=72)

rawdatapath = '../../../DATA/'
figpath = './Figures/'

startYear = 2009
endYear = 2016
numYears = endYear - startYear + 1

xpts_all = []
ypts_all = []
for year in xrange(startYear, endYear + 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)

Bscatdataoutpath = '../../BitBucket/backscatter/Data_output/'

year1 = 2009
year2 = 2015
day1 = 80
day2 = 110
day1str = '%03d' % day1
day2str = '%03d' % day2
lons = load(Bscatdataoutpath + 'lonsA')
lats = load(Bscatdataoutpath + 'latsA')
xpts, ypts = m(lons, lats)
#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/'

startYear=2009
endYear=2016
numYears=endYear-startYear+1

xpts_all=[]
ypts_all=[]
for year in xrange(startYear, endYear+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)


Bscatdataoutpath='../../BitBucket/backscatter/Data_output/'

year1=2009
year2=2015
day1=80
day2=110
day1str='%03d' %day1
day2str='%03d' %day2
lons=load(Bscatdataoutpath+'lonsA')
lats=load(Bscatdataoutpath+'latsA')
            lats.append(lats_t[x])
            lons.append(lons_t[x])

    return lons, lats

#mpl.rc('text', usetex=True)
m=Basemap(projection='stere', lat_0=74, lon_0=-90,llcrnrlon=-110, llcrnrlat=73,urcrnrlon=10, urcrnrlat=72)

rawdatapath='../../../DATA/'
figpath='./Figures/'


xpts_all=[]
ypts_all=[]

lonsT, latsT = ro.calc_icebridge_flights('2015FALL',rawdatapath, 'GR')
xptsT, yptsT=m(lonsT, latsT)
xpts_all.append(xptsT)
ypts_all.append(yptsT)

Bscatdataoutpath='../../BitBucket/backscatter/Data_output/'
year1=2015
year2=2015
day1=280
day2=300
day1str='%03d' %day1
day2str='%03d' %day2
lons=load(Bscatdataoutpath+'lonsA')
lats=load(Bscatdataoutpath+'latsA')
xpts, ypts = m(lons, lats)
image_mean=load(Bscatdataoutpath+'A'+str(year1)+str(year2)+'-'+day1str+'-'+day2str)