ma.drawcoastlines()
ma.drawparallels(np.arange(-90.0, 90.0, 30.0))
ma.drawmeridians(np.arange(-180.0, 180.0, 30.0))
# revhaxby=cmap_map(lambda x: x[::-1], cm.GMT_haxby)
# tlvs=np.arange(-20.,22.,1.)
# tlvs=(-5000,-4500,-4000,-3500,-3000,-2500,-2000,-1500,-1000,-500,0,500,1000,1500,2000,2500,3000,3500,4000,4500,5000)
# tlvs=(0,500,1000,1500,2000,2500,3000,3500,4000,4500,5000)
# tlvs=np.arange(0,5000,100)
tlvs = np.arange(0, 2000, 100)
# cs1 = m1.contourf(x1, y1, obsnaug, lvs, cmap=bcm.GMT_haxby_r)
# csa = ma.contourf(xa, ya, etopo, tlvs, cmap=bcm.GMT_no_green)
# csa = ma.contourf(xa, ya, etopo, tlvs, cmap=bcm.GMT_globe)
# csa = ma.contourf(xa, ya, etopo, tlvs, cmap=bcm.GMT_haxby_r)
csa = ma.contourf(xa, ya, etopo, tlvs, cmap=cm.RdYlGn_r)
ma.drawlsmask(land_color="lightgrey", ocean_color="steelblue", lakes=False)
# draw colorbar.
axa = plt.gca()
# pos = axa.get_position()
# l, b, w, h = pos.bounds
# caxa = plt.axes([l+w+0.01+baroffset, b, 0.02, h]) # setup colorbar axes
# plt.colorbar(csa, caxa,ticks=csa.levels[::2], format='%g', drawedges=False)
# plt.axes(axa) # make the original axes current again
# plt.axes(rect1)
plt.title("ETOPO 10 min (meters)", horizontalalignment="left", position=(0.0, 1.0))
llcrnrlon=-180,urcrnrlon=180,resolution='c')
xa, ya = ma(*np.meshgrid(lon, lat))
ma.drawcoastlines()
ma.drawparallels(np.arange(-90.,90.,30.))
ma.drawmeridians(np.arange(-180.,180.,30.))
#revhaxby=cmap_map(lambda x: x[::-1], cm.GMT_haxby)
#tlvs=np.arange(0.,160.,10.)
tlvs=np.arange(-160.,0.,10.)
#tlvs=(-30,-25,-20,-15,-10,-5,0,5,10,15,20,25,30)
#cs1 = m1.contourf(x1, y1, obsnaug, lvs, cmap=cm.GMT_haxby_r)
csa = ma.contourf(xa, ya, hf, tlvs, cmap=cm.GMT_no_green_r)
#csa = ma.contourf(xa, ya, obstaugavg, tlvs, cmap=cm.GMT_globe)
ma.drawlsmask(land_color='lightgrey',ocean_color='steelblue',lakes=False)
# draw colorbar.
axa = plt.gca()
pos = axa.get_position()
l, b, w, h = pos.bounds
#caxa = plt.axes([l+w+0.01, b+0.18, 0.02, h-0.35]) # setup colorbar axes
caxa = plt.axes([l+w+0.01, b, 0.02, h]) # setup colorbar axes
plt.colorbar(csa, caxa,ticks=csa.levels[::2], format='%g', drawedges=False)
plt.axes(axa) # make the original axes current again
#plt.axes(rect1)
plt.title("Heat Flux (mW/m**2)", horizontalalignment='left', position=(0.0,1.0))
ma.drawcoastlines()
ma.drawparallels(np.arange(-90., 90., 30.))
ma.drawmeridians(np.arange(-180., 180., 30.))
#revhaxby=cmap_map(lambda x: x[::-1], cm.GMT_haxby)
#tlvs=np.arange(-20.,22.,1.)
#tlvs=(-5000,-4500,-4000,-3500,-3000,-2500,-2000,-1500,-1000,-500,0,500,1000,1500,2000,2500,3000,3500,4000,4500,5000)
#tlvs=(0,500,1000,1500,2000,2500,3000,3500,4000,4500,5000)
#tlvs=np.arange(0,5000,100)
tlvs = np.arange(0, 2000, 100)
#cs1 = m1.contourf(x1, y1, obsnaug, lvs, cmap=bcm.GMT_haxby_r)
#csa = ma.contourf(xa, ya, etopo, tlvs, cmap=bcm.GMT_no_green)
#csa = ma.contourf(xa, ya, etopo, tlvs, cmap=bcm.GMT_globe)
#csa = ma.contourf(xa, ya, etopo, tlvs, cmap=bcm.GMT_haxby_r)
csa = ma.contourf(xa, ya, etopo, tlvs, cmap=cm.RdYlGn_r)
ma.drawlsmask(land_color='lightgrey', ocean_color='steelblue', lakes=False)
# draw colorbar.
axa = plt.gca()
#pos = axa.get_position()
#l, b, w, h = pos.bounds
#caxa = plt.axes([l+w+0.01+baroffset, b, 0.02, h]) # setup colorbar axes
#plt.colorbar(csa, caxa,ticks=csa.levels[::2], format='%g', drawedges=False)
#plt.axes(axa) # make the original axes current again
#plt.axes(rect1)
plt.title("ETOPO 10 min (meters)",
horizontalalignment='left',
position=(0.0, 1.0))
llcrnrlon=-180,urcrnrlon=180,resolution='c')
xa, ya = ma(*np.meshgrid(lon, lat))
ma.drawcoastlines()
ma.drawparallels(np.arange(-90., 90., 30.))
ma.drawmeridians(np.arange(-180., 180., 30.))
#revhaxby=cmap_map(lambda x: x[::-1], cm.GMT_haxby)
#tlvs=np.arange(0.,160.,10.)
tlvs = np.arange(-160., 0., 10.)
#tlvs=(-30,-25,-20,-15,-10,-5,0,5,10,15,20,25,30)
#cs1 = m1.contourf(x1, y1, obsnaug, lvs, cmap=cm.GMT_haxby_r)
csa = ma.contourf(xa, ya, hf, tlvs, cmap=cm.GMT_no_green_r)
#csa = ma.contourf(xa, ya, obstaugavg, tlvs, cmap=cm.GMT_globe)
ma.drawlsmask(land_color='lightgrey', ocean_color='steelblue', lakes=False)
# draw colorbar.
axa = plt.gca()
pos = axa.get_position()
l, b, w, h = pos.bounds
#caxa = plt.axes([l+w+0.01, b+0.18, 0.02, h-0.35]) # setup colorbar axes
caxa = plt.axes([l + w + 0.01, b, 0.02, h]) # setup colorbar axes
plt.colorbar(csa, caxa, ticks=csa.levels[::2], format='%g', drawedges=False)
plt.axes(axa) # make the original axes current again
#plt.axes(rect1)
plt.title("Heat Flux (mW/m**2)",