num_kms = 10
fadd = ''
ftype = '1km_xyres2m_' + str(thresh) + 'cm' + fadd
figpath = '../Figures/Paper/'
datapath2D = '../Data_output/' + ftype + '/2D/ATMO/'
rawdatapath = '../../../../DATA/'
BSCATdatapath = '../../backscatter/Data_output/'
startYear = 2009
endYear = 2015
numYears = endYear - startYear + 1
ascatOIB = ma.masked_all((numYears, 1530, 1530))
for year in xrange(startYear, endYear + 1):
xptsB, yptsB, ascatOIB[year - startYear] = ro.get_ASCATOIB(
BSCATdatapath, mplot, year)
xpts = []
ypts = []
Cds = []
for year in xrange(startYear, endYear + 1):
xptsT = load(datapath2D + 'mean_xptsB' + str(year) + str(num_kms) +
'km_2D.txt')
yptsT = load(datapath2D + 'mean_yptsB' + str(year) + str(num_kms) +
'km_2D.txt')
CdT = load(datapath2D + 'mean_Cd' + str(year) + str(num_kms) + 'km_2D.txt')
xpts.append(xptsT)
ypts.append(yptsT)
Cds.append(CdT)
aspect = mplot.ymax / mplot.xmax
