meteoinfoLab 绘制grads 格式PM2.5 遇到问题！

bearwoo

格点pm2.5的ctl用grads打开，里面每个时次的变量都可以画出来，但是我用MeteoinfoLab 把ctl 导入后，有的时次可以画出来而有的直接报错，不懂这是什么原因，还请大神帮忙看一下！    以下贴出我的ctl 以及 Lab中的脚本：
ctl:

1. dset ^aero_norad2016121300_11d
   
2. options sequential big_endian
   
3. title model output from grapes
   
4. undef -9.99E+33
   
5. xdef  300 linear  100.0000    0.1500
   
6. ydef  200 linear   20.0000    0.1500
   
7. zdef   33 linear 1 1
   
8. tdef   792 linear 00z13DEC2016    60mn
   
9. vars   8
   
10. rh    3 99 relative humidity
    
11. aero_norad_ext   30 99 ext
    
12. pm25  3 99 pm2.5
    
13. vis_vert   0 99 vertical visibility
    
14. vis  0 99 old visibility
    
15. aod  0 0 aod
    
16. pblh 0 0 pblh
    
17. fktm 15 99 diffusion coefficient
    
18. endvars

复制代码

meteoinfo脚本：

1. fn1 = addfile('H:/model/output/mrf/case1/aero_hour/aero.ctl')
   
2. pm = fn1['pm25'][168,0,'30:45','110:125']
   
3. #Set data folders
   
4. basedir = 'D:/py/APP/MeteoInfo_Java_1.4R8_Files/MeteoInfo'
   
5. mapdir = os.path.join(basedir, 'map')
   
6. #Read shape files
   
7. bou2_layer = shaperead(os.path.join(mapdir, 'bou2_4p.shp'))
   
8. bou1_layer = shaperead(os.path.join(mapdir, 'bou1_4l.shp'))
   
9. city_layer = shaperead(os.path.join(mapdir, 'res1_4m.shp'))
   
10. 
    
11. #plot
    
12. axesm()

复制代码

MeteoInfo

是一个bug，你可以加入MeteoInfo QQ群下载最新版本。

bearwoo

MeteoInfo 发表于 2017-11-12 14:36
是一个bug，你可以加入MeteoInfo QQ群下载最新版本。

原来是这样，谢谢王老师！