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,