MeteoInfoLab脚本示例：MERRA HDF数据

MERRA是NOAA的一种再分析资料，HDF数据遵循COARDS协议，读取比较简单。

脚本程序：

#Add data file
folder = 'D:/Temp/hdf/'
fns = 'MERRA300.prod.assim.tavg1_2d_slv_Nx.20150101.SUB.hdf'
fn = folder + fns
f = addfile(fn)
vname = 't2m'
t = f[vname]
t2m = t[0,:,:]
#Plot
axesm()
mlayer = shaperead('D:/Temp/map/country1.shp')
geoshow(mlayer)
layer = imshowm(t2m, 20)
colorbar(layer, orientation='horizontal')
title(vname + ' at time=0')
axism()