- 积分
- 9
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2017-5-9
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
fn = 'H:/meteo/GRIB/
fnl_20121005_12_00.grib2'
f = addfile(fn)
T850 =
f['Temperature_isobaric'][0,'85000.0','15:55','40:70'] - 273.16
#Plot
axesm()
lworld = shaperead('H:/meteo/MeteoInfo_Java_1.4.2_Files/MeteoInfo/map/country1.shp')
geoshow(lworld, edgecolor='k')
levs = arange(-40,40,4)
layer = contourm(T850, levs)
clabel(layer)
colorbar(layer)
t = f.gettime(0)
title('Temperature 850hpa (' + t.strftime('%Y-%m-%d %H:00') + ')')
可以读出来中国的图,但是想要一个点的温度和相对湿度,需要改哪里,求大神指教
|
|