#Add data file

fn = 'AMSR_E_L3_DailyLand_V06_20091231.hdf'

f = addfile(os.path.join('D:/Temp/hdf', fn))

#vname = 'D_Soil_Moisture'

vname = 'A_TB36.5H_(Res_1)'

data = f[vname][::-1,:]

data[data==0] = -9999

data.fill_value = -9999

data = data * 0.1

yn = data.dimlen(0)

xn = data.dimlen(1)

#Project data

toproj = projinfo()    #longlat projection

x = linspace(-180, 180, xn)

y = linspace(-90, 90, yn)

lon, lat = meshgrid(x, y)

data = data.project(lon, lat, toproj=toproj)

#Plot

axesm(tickfontsize=12)

mlayer = shaperead('D:/Temp/map/country1.shp')

geoshow(mlayer, edgecolor='k')

levs = arange(140, 320, 10)

layer = imshowm(x, y, data, levs)

colorbar(layer, fontsize=12)

xticks(arange(-180, 181, 30))

yticks(arange(-90, 91, 30))

title([fn, vname])