ERA5数据meteoinfo一页多图

lixiliang2008

f = addfile("F:/1221-22.nc")
gh= f['z'][2,"500","30":"50","90":"130"]/100
udata = f['u'][2,"500","30":"50","90":"130"]
vdata = f['v'][2,"500","30":"50","90":"130"]
rmdata=f['r'][2,"500","30":"50","90":"130"]
tmdata=f['t'][2,"500","30":"50","90":"130"]-273.15
lon = f['longitude']["90":"130"]
lat = f['latitude']["30":"50"]
udata=udata[::9,::9]
vdata=vdata[::9,::9]
speed = sqrt(udata*udata+vdata*vdata)
lon, lat = meshgrid(lon, lat)

figure(figsize=[900, 520], newfig=False)
axesm(position=[0.04,0.52,0.42,0.44],tickfontsize=12,tickbold=True)

levgh=arange(510,588,6)
layergh = contourm(gh,levgh,size=2,color='black')#绘制高度场等值线
clabel(layergh,drawshadow=True,fontsize=12,bold=True)#标注值
#风场
layerw = barbsm(udata,vdata,size=14,color="black")
#相对湿度阴影
levsrh = arange(60,100,10)
colsrh = [(255,255,255),(0,229,238),(0,238,0),(0,205,0),(0,160,0),(0,130,0)]
layerrh = contourfm(rmdata,levsrh,colors=colsrh)
#colorbar(layerrh,fontsize=12)
#温度
levs = arange(-80,0,4)
layert = contourm(tmdata,levs,size=2,color='r',linestyle='--')
clabel(layert)
mlayer1 = shaperead('E:\MeteoInfo\map\cn_province.shp')
geoshow(mlayer1, edgecolor='black',size=1)
xlim(90,130)
ylim(30,50)
xticks([90,100,110,120,130],fontsize=12,bold=True)
yticks([30,35,40,45,50],fontsize=12,bold=True)

f = addfile("F:/1221-22.nc")
gh= f['z'][2,"850","30":"50","90":"130"]/10
udata = f['u'][2,"850","30":"50","90":"130"]
vdata = f['v'][2,"850","30":"50","90":"130"]
rmdata=f['r'][2,"850","30":"50","90":"130"]
tmdata=f['t'][2,"850","30":"50","90":"130"]-273.15
lon = f['longitude']["90":"130"]
lat = f['latitude']["30":"50"]
udata=udata[::9,::9]
vdata=vdata[::9,::9]
speed = sqrt(udata*udata+vdata*vdata)
lon, lat = meshgrid(lon, lat)
axesm(position=[0.52,0.52,0.42,0.44],tickfontsize=12,tickbold=True)
levgh=arange(1000,1900,60)
layergh = contourm(gh,levgh,size=2,color='black')#绘制高度场等值线
clabel(layergh,drawshadow=True,fontsize=12,bold=True)#标注值
#风场
layerw = barbsm(udata,vdata,size=14,color="black")
#相对湿度阴影
levsrh = arange(60,100,10)
colsrh = [(255,255,255),(0,229,238),(0,238,0),(0,205,0),(0,160,0),(0,130,0)]
layerrh = contourfm(rmdata,levsrh,colors=colsrh)
#colorbar(layerrh,fontsize=12)
#温度
levs1 = arange(-80,0,4)
layert1 = contourm(tmdata,levs1,size=2,color='r',linestyle='--')
clabel(layert1)
levs2 = arange(0,80,4)
layert2 = contourm(tmdata,levs2,size=2,color='r',linestyle='-')
clabel(layert2)
mlayer1 = shaperead('E:\MeteoInfo\map\cn_province.shp')
geoshow(mlayer1, edgecolor='black',size=1)
xlim(90,130)
ylim(30,50)
xticks([90,100,110,120,130],fontsize=12,bold=True)
yticks([30,35,40,45,50],fontsize=12,bold=True)

f = addfile("F:/1221-22.nc")
gh= f['z'][3,"500","30":"50","90":"130"]/100
udata = f['u'][3,"500","30":"50","90":"130"]
vdata = f['v'][3,"500","30":"50","90":"130"]
rmdata=f['r'][3,"500","30":"50","90":"130"]
tmdata=f['t'][3,"500","30":"50","90":"130"]-273.15
lon = f['longitude']["90":"130"]
lat = f['latitude']["30":"50"]
udata=udata[::9,::9]
vdata=vdata[::9,::9]
speed = sqrt(udata*udata+vdata*vdata)
lon, lat = meshgrid(lon, lat)

axesm(position=[0.04,0.06,0.42,0.44],tickfontsize=12,tickbold=True)

levgh=arange(510,588,6)
layergh = contourm(gh,levgh,size=2,color='black')#绘制高度场等值线
clabel(layergh,drawshadow=True,fontsize=12,bold=True)#标注值
#风场
layerw = barbsm(udata,vdata,size=14,color="black")
#相对湿度阴影
levsrh = arange(60,100,10)
colsrh = [(255,255,255),(0,229,238),(0,238,0),(0,205,0),(0,160,0),(0,130,0)]
layerrh = contourfm(rmdata,levsrh,colors=colsrh)
#colorbar(layerrh,fontsize=12)
#温度
levs = arange(-80,0,4)
layert = contourm(tmdata,levs,size=2,color='r',linestyle='--')
clabel(layert)
mlayer1 = shaperead('E:\MeteoInfo\map\cn_province.shp')
geoshow(mlayer1, edgecolor='black',size=1)
xlim(90,130)
ylim(30,50)
xticks([90,100,110,120,130],fontsize=12,bold=True)
yticks([30,35,40,45,50],fontsize=12,bold=True)

#caxes(axs[3])
f = addfile("F:/1221-22.nc")
gh= f['z'][3,"850","30":"50","90":"130"]/10
udata = f['u'][3,"850","30":"50","90":"130"]
vdata = f['v'][3,"850","30":"50","90":"130"]
rmdata=f['r'][3,"850","30":"50","90":"130"]
tmdata=f['t'][3,"850","30":"50","90":"130"]-273.15
lon = f['longitude']["90":"130"]
lat = f['latitude']["30":"50"]
udata=udata[::9,::9]
vdata=vdata[::9,::9]
speed = sqrt(udata*udata+vdata*vdata)
lon, lat = meshgrid(lon, lat)
axesm(position=[0.52,0.06,0.42,0.44],tickfontsize=12,tickbold=True)
levgh=arange(1200,1900,60)
layergh = contourm(gh,levgh,size=2,color='black')#绘制高度场等值线
clabel(layergh,drawshadow=True,fontsize=12,bold=True)#标注值
#风场
layerw = barbsm(udata,vdata,size=14,color="black")
#相对湿度阴影
levsrh = arange(60,100,10)
colsrh = [(255,255,255),(0,229,238),(0,238,0),(0,205,0),(0,160,0),(0,130,0)]
layerrh = contourfm(rmdata,levsrh,colors=colsrh)
#温度
levs1 = arange(-80,0,4)
layert1 = contourm(tmdata,levs1,size=2,color='r',linestyle='--')
clabel(layert1)
levs2 = arange(0,80,4)
layert2 = contourm(tmdata,levs2,size=2,color='r',linestyle='-')
clabel(layert2)
mlayer1 = shaperead('E:\MeteoInfo\map\cn_province.shp')
geoshow(mlayer1, edgecolor='black',size=1)
xlim(90,130)
ylim(30,50)
xticks([90,100,110,120,130],fontsize=12,bold=True)
yticks([30,35,40,45,50],fontsize=12,bold=True)

axes(position=[0.53,0.26,0.42,0.44],newaxes=True, axis=False,tickfontsize=12,tickbold=True)
colorbar(layerrh,fontsize=12)

guo066600

楼主能传一下您下载数据的具体链接吗？

罗垒

楼主可以注释一下吗

罗垒

jython吗

lixiliang2008

guo066600 发表于 2021-2-6 17:25
楼主能传一下您下载数据的具体链接吗？

https://cds.climate.copernicus.e ... ure-levels?tab=form