MeteoInfo记录篇-MeteoInfoLab

最近使用grads无法实现自己的功能（当然也可能是我没有找到方法），最后使用MeteoInfo解决了问题，下面是使用记录，留个笔记，方便大家一起学习！【气象绘图功能，画等值面和风向杆】

fn='E://RMAPS/20200322/Z_NAFP_C_BEPK_20200322120000_P_EMC-NNC-f03.BIN'
f=addfile(fn)
print(f)
#查看数据中字段后，下面是读取字段，第一个[字段名]，第二个[]里面代表的是4维的，“：”前后代表的是开始结束的index
data=f['VAR0-14-204_FROM_38-0--1_height_above_ground'][:,:,:,:]
axesm()
#设置图例等级
levs=[5,35,75,115,150,250,350,500]
#设置图例颜色
cols=[(250,250,250),(60,250,0),(250,250,0),(250,160,0),(250,100,0),(250,0,0),(168,26,10),(74,4,4),(217,78,143)]
#绘制等值面
layer=contourfm(data,levs,colors=cols,smooth=False)
#显示图例
colorbar(layer)
#获取风的数据
udata=f['u-component_of_wind_height_above_ground'][:,:,:,:]
vdata=f['v-component_of_wind_height_above_ground'][:,:,:,:]
#风数据横纵取的风向杆的间隔数，每隔10个显示
udata=udata[::10,::10]
vdata=vdata[::10,::10]
#绘制风向杆，isuv是代表是否是uv分量，False的时候，说明是风速、风向值
layer = barbsm(udata,vdata,colors='k',isuv=True)
#绘图的经纬度范围
xlim(106,114.5)
ylim(33,39)
#写标题，\r\n表示的是换行，$^{3}$上标，$_{2.5}$下标，只是显示的时候效果比较浅，觉得不太好看
title('RMAPS\r\n2020032220','left',color='b',bold=False,fontsize=18)
#保存图片
savefig("D://a.png",1000,700)
#清除数据
clear()

提示
1、这个里面注释使用的是“#”
2、fn的路径中必须是“/”，如果是""，我遇到的提示utf-8的编码错误。(当时换了短的路径可以，还以为是路径太长了呢)
当初错误截图：

3、print这个可以进行参数的输出，先输出看一下参数名字，后面读取使用
4、最后这个clear不知道是否作用！！！