MeteoInfoLab脚本示例：获取一维数据并绘图

MeteoInfo

气象数据基本为多维数据（通常是4维，空间3维加时间维），只让数据中一维可变，其它维均固定即可提取一维数据。比如此例中固定了时间维、高度维、纬度维，只保留经度维可变：hgt = f['hgt'][0,'500','4','180:360']。固定某一维可以用序号（从0开始），比如此例中的时间维为0，既第一个时次。也可以用该维中真实的数据，比如高度维希望是500 hPa，此时需要用引号：'500'，纬度维的固定也是类似。维的范围可用序号来限定，比如：5:20，也可用该维的真实值限定，比如此例中的经度维限定为180-360度，需要用引号：'180:360'。一维线条图绘制用plot命令。


脚本程序：

1. f = addfile('D:/Temp/test/model.ctl')
   
2. hgt = f['hgt'][0, '500','40','180:360']
   
3. plot(hgt, 'r-o', linewidth=2)
   
4. title('HGT')
   
5. grid(True)
   
6. ylabel('HGT')
   
7. xlabel('Longitude')
   
8. show()

f = addfile('D:/Temp/test/model.ctl')<br /> hgt = f['hgt'][0, '500','40','180:360']<br /> plot(hgt, 'r-o', linewidth=2)<br /> title('HGT')<br /> grid(True)<br /> ylabel('HGT')<br /> xlabel('Longitude')<br /> show()

兰溪之水

王老师的语法既像NCL又像MATLAB

MeteoInfo

兰溪之水 发表于 2015-6-18 16:28
王老师的语法既像NCL又像MATLAB

嘿嘿，各取所长