import netCDF4 as nc
import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap
import numpy as np
filepath = "/run/media/root/*WRF/-5/s/wrfout_d02_2014-09-28_12_00_00"
data = nc.Dataset(filepath,"r")
var = "REFL_10CM"
dbz = data.variables[var][:]
dbzsize = np.shape(data.variables[var])
timelen = dbzsize[0]
heightlen = dbzsize[1]
latlen = dbzsize[2]
lonlen = dbzsize[3]
print(dbzsize)
lats = 30; late = 35;
lons = 114; lone = 120;
m = Basemap(projection="cyl",llcrnrlat=lats,urcrnrlat=late,llcrnrlon=lons,urcrnrlon=lone,resolution="i")
m.drawparallels(np.arange(lats,late+1,1.),labels = [1,0,0,0],fontsize=15,linewidth = 0. )
# set linewidth = 0 to turn off longitude and latitude dash line
<