Re: [ferret_users] problems reading a GMT-style grid

[Hein Zelle <hein.zelle@xxxxxxxxxxxxx>]

Hi Billy,

> let zz=reshape('Z',x[gx=xlon]+y[gy=ylat])
> - - - - - - -
> This works fine except for one crucial thing: the latitudes are
> reversed!


Would it work if you define your y axis in the reverse direction?
Ferret is typically clever enough to understand reversed axes.  You
could also do it in 2 steps, assigning the first attempt to a reversed
y axis with @asn.

A very ugly method could be to create matching arrays of x and y
coordinates, e.g. xcoords = xsequence(x[gx=xlon]+0*y[gy=ylat]).  You
would probably also need a reversed y axis for the ycoords variable.
You could then use the scat2grid functions to put bathymetry onto your
target grid.  Or do it in 2 steps the way Ansley said - that's
definitely cleaner.

Kind regards,
     Hein


-- 

Dr. Hein Zelle
Advisor Meteorology & Oceanography

Tel:    +31 (0)527-242299
Fax:    +31 (0)527-242016
Email:  hein.zelle@xxxxxxxxxxxxx
Web:    www.bmtargoss.com

BMT ARGOSS
P.O. Box 61, 8325 ZH Vollenhove
Voorsterweg 28, 8316 PT Marknesse
The Netherlands

----Confidentiality Notice & Disclaimer---- 

The contents of this e-mail and any attachments are intended for the
use of the mail addressee(s) shown. If you are not that person, you
are not allowed to read it, to take any action based upon it or to
copy it, forward, distribute or disclose the contents of it and you
should please delete it from your system. BMT ARGOSS does not accept
liability for any errors or omissions in the context of this e-mail or
its attachments which arise as a result of internet transmission, nor
accept liability for statements which are those of the author and
clearly not made on behalf of BMT ARGOSS.