Skip to content

Commit

Permalink
fix typo
Browse files Browse the repository at this point in the history
  • Loading branch information
jswhit committed Jul 16, 2011
1 parent b37927c commit 8770b74
Showing 1 changed file with 4 additions and 4 deletions.
8 changes: 4 additions & 4 deletions examples/makelsmask.py
Expand Up @@ -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
Expand All @@ -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()

0 comments on commit 8770b74

Please sign in to comment.