From 8770b7470e06f97f162bf9b939991b953cc6e120 Mon Sep 17 00:00:00 2001 From: Jeff Whitaker Date: Sat, 16 Jul 2011 10:50:40 -0600 Subject: [PATCH] fix typo --- examples/makelsmask.py | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/examples/makelsmask.py b/examples/makelsmask.py index 154945f0e..7a3c0a3b8 100644 --- a/examples/makelsmask.py +++ b/examples/makelsmask.py @@ -3,8 +3,8 @@ import numpy as np # make land-sea-lake mask from built-in coastline dataset using is_land method. # this is very slow! -minutes = 5 -resolution = 'i' +minutes = 150 +resolution = 'c' delon = minutes/60. nlons = int(360./delon); nlons = nlons + 1 @@ -18,11 +18,11 @@ m =\ Basemap(llcrnrlon=-180,llcrnrlat=-90,urcrnrlon=180,urcrnrlat=90,resolution=resolution,projection='cyl') for j,lat in enumerate(lats): - print j + #print j for i,lon in enumerate(lons): lsmask[j,i] = m.is_land(lon,lat,lsmask=True) m.drawlsmask(land_color='coral',ocean_color='aqua',lsmask=lsmask,lsmask_lons=lons,lsmask_lats=lats,lakes=True) plt.title('%s minute degree land-sea mask' % minutes) -lsmask.tofile('%sminlsmask_gshhs_%s.dat' % (minutes, resolution)) +#lsmask.tofile('%sminlsmask_gshhs_%s.dat' % (minutes, resolution)) plt.show()