GPM swath data

This example code illustrates how to access and visualize a GPM L1C swath data.

#Add data file
folder = 'D:/Temp/hdf/'
fns = '1C.F19.SSMIS.XCAL2015-P.20160105-S214106-E232259.009078.V03A.HDF5'
fn = folder + fns
f = addfile(fn)
lon = f['Longitude'][:,:]
lat = f['Latitude'][:,:]
vname = 'Tc'
v_data = f[vname]
data = v_data[:,:,0]
data[data<=-9999.9] = nan
long_name = v_data.attrvalue('LongName')[0]
units = v_data.attrvalue('Units')[0]
#Plot
axesm()
geoshow('country', edgecolor='k')
levs = arange(40, 90, 1)
layer = pcolor(lon, lat, data, 20)
colorbar(layer, orientation='horizontal', aspect=40, label=units)
title('{0}{1}'.format(fns, long_name+' (nchannel1=0)'))
axism()