###########################################################################
### Read in model data
polarave = []
high = []
for i in range(len(runnames)):
polaraveq = CAP.PolarCap(runnamesdata[i], variable, level, latpolar,
period)
highq = SH.SHI(runnamesdata[i], period)
polarave.append(polaraveq)
high.append(highq)
###########################################################################
### Read in reanalysis data
years = np.arange(1979, 2019 + 1, 1)
late, lone, leve, thicke = REAN.readOBS('ERA5', variable, level, period)
late, lone, leve, slpe = REAN.readOBS('ERA5', 'SLP', level, period)
she = REAN.calcOBS_SHI(slpe, late, lone)
vare = REAN.calcOBS_PolarCap(thicke, late, lone, latpolar)
latrr, lonrr, levr, thickr = REAN.readOBS('NCEP1', variable, level, period)
latr, lonr, levr, slpr = REAN.readOBS('NCEP1', 'SLP', level, period)
shr = REAN.calcOBS_SHI(slpr, latr, lonr)
varr = REAN.calcOBS_PolarCap(thickr, latrr, lonrr, latpolar)
###############################################################################
###############################################################################
###############################################################################
### Calculate ensemble means
meanPOL = np.empty((len(polarave)))
meanSHI = np.empty((len(high)))
for i in range(len(runnames)):
meanPOL[i] = np.nanmean(polarave[i])
high.append(highq)
###########################################################################
### Read in model data for E3SM
polarave_E3SM = []
high_E3SM = []
for i in range(len(runnames_E3SM)):
polaraveq_E3SM = CAP.PolarCap(runnamesdata_E3SM[i],variable,level,latpolar,period)
highq_E3SM = SH.SHI(runnamesdata_E3SM[i],period)
polarave_E3SM.append(polaraveq_E3SM)
high_E3SM.append(highq_E3SM)
###########################################################################
### Read in reanalysis data
years = np.arange(1979,2019+1,1)
late,lone,leve,thicke = REAN.readOBS('ERA5',variable,level,period)
late,lone,leve,slpe = REAN.readOBS('ERA5','SLP',level,period)
she = REAN.calcOBS_SHI(slpe,late,lone)
vare= REAN.calcOBS_PolarCap(thicke,late,lone,latpolar)
latrr,lonrr,levr,thickr = REAN.readOBS('NCEP1',variable,level,period)
latr,lonr,levr,slpr = REAN.readOBS('NCEP1','SLP',level,period)
shr = REAN.calcOBS_SHI(slpr,latr,lonr)
varr= REAN.calcOBS_PolarCap(thickr,latrr,lonrr,latpolar)
###########################################################################
### Read in AMIP data
yearsAA = np.arange(1979,2016+1,1)
lataa,lonaa,timeaa,levaa,t_AA = AMAA.readAMIPAA(variable,'AMIP-AA',level,False,True,period)
lataa,lonaa,timeaa,levaa,slp_AA = AMAA.readAMIPAA('SLP','AMIP-AA',level,False,True,period)
t_AAe = np.nanmean(t_AA,axis=0)
slp_AAe = np.nanmean(slp_AA,axis=0)
sh_AA = REAN.calcOBS_SHI(slp_AAe,lataa,lonaa)