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)",