- 积分
- 10
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2017-4-15
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
各位前辈,我这个读取的再分析资料温度场得出来的温度图明显错误的,但是又不知道应该怎么改正代码。请各位前辈指教!
温度图如下:
代码如下:
fn = 'H:/NOAA/fnl_20050516_00_00.grib1'
f = addfile(fn)
T500 = f['Temperature_isobaric'][0,'50000.0','25:45','100:130'] - 273.16
#Plot
axesm()
lworld = shaperead('H:/meteo/MeteoInfo_Java_1.4.2_Files/MeteoInfo/map/country1.shp')
geoshow(lworld, edgecolor='k')
levs = arange(-30,30,1)
layer = contourm(T500, levs)
clabel(layer)
colorbar(layer)
t = f.gettime(0)
title('Temperature 500hpa (' + t.strftime('%Y-%m-%d %H:00') + ')')
|
|