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

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

 

脚本程序：

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

 

转载于:https://www.cnblogs.com/yaqiang/p/4586156.html