startYear = 2000
endYear = 2016
numYears = endYear - startYear + 1
startMon = 0
endMon = 11
numMons = endMon - startMon + 1
alg = 1 #0 is NT, 1 is BT
areaF = reshape(
fromfile(file=open(dataPath + '/OTHER/psn25area_v3.dat', 'rb'),
dtype='<i4') / 1000., [448, 304]) / 1e6
lats, lons = cF.get_psnlatslons(dataPath)
xpts, ypts = m(lons, lats)
#iceExtMon=np.zeros((numYears, numMons))
#iceAreaMon=np.zeros((numYears, numMons))
for month in xrange(startMon, endMon + 1):
mstr = '%02d' % (month + 1)
iceExtMon = []
iceAreaMon = []
for year in xrange(startYear, endYear + 1):
dateStr = str(year) + mstr
pmask = cF.get_pmask(year, month)
#pmask=89.2
rcParams['axes.linewidth'] = .5
rcParams['lines.linewidth'] = .5
rcParams['patch.linewidth'] = .5
rcParams['ytick.labelsize'] = 8
rcParams['legend.fontsize'] = 8
rcParams['font.size'] = 8
rc('font', **{'family': 'sans-serif', 'sans-serif': ['Arial']})
dataPath = '../../../../DATA'
figpath = '../Figures/'
dataOutPath = '../DataOutput/Intense/'
dataOutPathConc = '../DataOutput/Conc/'
m = Basemap(projection='npstere', boundinglat=53, lon_0=0, resolution='l')
#m = Basemap(projection='npstere',boundinglat=30.52,lon_0=0, resolution='l' )
latsI, lonsI = cF.get_psnlatslons(dataPath)
xptsI, yptsI = m(lonsI, latsI)
alg = 0
year = 2016
Intmonth = []
for month in xrange(0, 12):
mstr = '%02d' % (month + 1)
dateStr = str(year) + mstr
Intmonth.append(
load(dataOutPath + 'iceIntense' + dateStr + 'Alg' + str(alg)))
Concmonth = []
for month in xrange(0, 12):