建议你用MeteoInfoLab (MeteoInfoLab脚本汇总贴
http://bbs.06climate.com/forum.php?mod=viewthread&tid=36151&fromuid=106
(出处: 气象家园)
),示例脚本:
- fn = 'D:/Temp/grib/fnl_20110416_00_00'
- f = addfile(fn)
- v = f['Temperature_isobaric']
- zn = v.dimlen(1)
- data = v[0,zn-1,:,:] #First level from ground
- data = data - 273.16
- #Plot
- axesm()
- lworld = shaperead('D:/Temp/Map/country1.shp')
- lchina = shaperead('D:/Temp/Map/bou2_4p.shp')
- geoshow(lworld, edgecolor='k')
- geoshow(lchina, edgecolor='k')
- layer = contourfm(data, 20)
- colorbar(layer)
- title('Temperature')
|