On Wed, Feb 13, 2008 at 2:16 PM, Jeff Whitaker <[EMAIL PROTECTED]> wrote:
> Evan Mason wrote:
> > Hi Jeff
> >
> > Here are the corners:
> >
> > lon_corners = N.array([-4.09300764,-35.76003475,-43.72330207,
> > -12.05627497])
> > lat_corners = N.array([41.90278813, 49.2136974, 14.7209971, 7.4100878
Evan Mason wrote:
> Hi Jeff
>
> Here are the corners:
>
> lon_corners = N.array([-4.09300764,-35.76003475,-43.72330207,
> -12.05627497])
> lat_corners = N.array([41.90278813, 49.2136974, 14.7209971, 7.41008784])
>
> The reason for the differences is that the matlab script is very
> fiddly, lots o
By losing the memory I mean that the grid is no longer rotated; that the
rotation I introduced through lat1, lon1, lat2, lon2 is lost. If you look
at the latitude of the two bottom corners you see that they are the same,
they should be different - for the matlab script they are different. In
othe
Evan Mason wrote:
> Thanks for the replies. The map you produced, Jeff, looks as it
> should. However, I am trying to make an ocean model grid, and so I
> require two 2d arrays of lon and lat, at my desired grid spacing.
> This is why I try the steps:
>
> dl = 2.
> nx = int((M.xmax - M.xm
Thanks for the replies. The map you produced, Jeff, looks as it should.
However, I am trying to make an ocean model grid, and so I require two 2d
arrays of lon and lat, at my desired grid spacing. This is why I try the
steps:
dl = 2.
nx = int((M.xmax - M.xmin) / dl) + 1
ny = int((M.ymax - M.
On Feb 12, 2008, at 8:39 PM, Evan Mason wrote:
> dl = 2.
> nx = int((M.xmax - M.xmin) / dl) + 1
> ny = int((M.ymax - M.ymin) / dl) + 1
>
> lonr, latr = M.makegrid(nx, ny)
>
> plot(lonr, latr, 'c.')
> show()
I think you might be looking for M.plot() rather than plot()??
plot() will just overw
Evan Mason wrote:
> Hi, I am having some problems using the oblique mercator projection in
> basemap. I want to define a rectangular orthogonal grid, rotated
> clockwise by about 13 degrees. I want to define grid cells of size,
> say, about 20x20 km. The script I have so far is below. The pr
Hi, I am having some problems using the oblique mercator projection in
basemap. I want to define a rectangular orthogonal grid, rotated clockwise
by about 13 degrees. I want to define grid cells of size, say, about 20x20
km. The script I have so far is below. The problem is that at some point
(
