print(date)
dat = ECMWF('FULL-EA',date,exp='VOZ')
dat._get_var('T')
dat._get_var('VO')
dat._get_var('O3')
dat._get_var('U')
dat._get_var('V')
dat._mkp()
dat._mkz()
dat._mkthet()
dat._mkpv()
if date >= datetime(2017,10,9,0):
datr = dat.shift2west(-180)
dats[i] = datr.interpolPT(pts,varList=varList,lonRange=(-140,-20),latRange=(30,60))
elif date >= datetime(2017,9,22,0):
dats[i] = dat.interpolPT(pts,varList=varList,lonRange=(140,260),latRange=(35,65))
elif date >= datetime(2017,9,10,0):
dats[i] = dat.interpolPT(pts,varList=varList,lonRange=(60,180),latRange=(35,65))
else:
datr = dat.shift2west(-180)
dats[i] = datr.interpolPT(pts,varList=varList,lonRange=(-20,100),latRange=(35,65))
dat.close()
date += timedelta(hours=6)
i += 1
#%%
if predates:
outfile = 'ERA5-extract-pre.pkl'
else:
outfile = 'ERA5-extract-B1-PT5.pkl'
with gzip.open(outfile,'wb') as f:
while date < day2:
print(date)
dat = ECMWF('FULL-EA', date, exp='VOZ')
dat._get_var('T')
dat._get_var('VO')
dat._get_var('O3')
dat._get_var('U')
dat._get_var('V')
dat._mkp()
dat._mkz()
dat._mkthet()
dat._mkpv()
if date >= datetime(2017, 10, 12, 0):
dats[i] = dat.interpolPT(pts,
varList=varList,
latRange=(5, 35),
lonRange=(100, 220))
elif date >= datetime(2017, 10, 3, 0):
datr = dat.shift2west(-180)
dats[i] = datr.interpolPT(pts,
varList=varList,
latRange=(20, 50),
lonRange=(-180, -60))
elif date >= datetime(2017, 9, 23, 0):
datr = dat.shift2west(-180)
dats[i] = datr.interpolPT(pts,
varList=varList,
latRange=(20, 50),
lonRange=(-120, 0))
elif date >= datetime(2017, 9, 11, 0):
datr = dat.shift2west(-180)
while date < day2:
print(date)
dat = ECMWF('FULL-EA', date, exp='VOZ')
dat._get_var('T')
dat._get_var('VO')
dat._get_var('O3')
dat._get_var('U')
dat._get_var('V')
dat._mkp()
dat._mkz()
dat._mkthet()
dat._mkpv()
if date >= datetime(2017, 10, 9, 0):
dats[i] = dat.interpolPT(pts,
varList=varList,
lonRange=(30, 150),
latRange=(40, 70))
elif date >= datetime(2017, 10, 1, 0):
datr = dat.shift2west(-180)
dats[i] = datr.interpolPT(pts,
varList=varList,
lonRange=(-50, 70),
latRange=(40, 70))
elif date >= datetime(2017, 9, 26, 0):
datr = dat.shift2west(-180)
dats[i] = datr.interpolPT(pts,
varList=varList,
lonRange=(-130, -10),
latRange=(40, 70))
elif date >= datetime(2017, 9, 20, 0):
dats[i] = dat.interpolPT(pts,
