Hi -
To make a 2D array with all rows the same you could do
let xb2=xb + 0*yb
And, yes, I do see the problems with contouring here. Try the SHADE plot,
shade/hlimits=-100:400:50 xb2+zb,xb2,zb
which looks fine.
I'm not sure why the CONTOUR command thinks it needs to go draw that bit
over near -100. Possibly there are places in the code that normally handle
modulo-x operations for wrapping around in longitude, that are somehow
coming into play here.
I'll see what I can find out.
Ansley
Patrick McGuire wrote:
thank you, Ansley!!
Yes, I am trying to use these curvilinear grids. (sorry for calling them
as irregular grids).
It looks like I need to make the xb grid be 2D instead of 1D.
I just tried to do this, with this script:
let xb=x[x=0:290:10]
let yb=y[y=0:290:10]
let yb1= if (yb1 GT -1) then 1 else 0
let xb2=xb*yb1 !this makes a 2D array with all the rows the same; is
! there an easier way to do this??
let zb=I[x=0:290:10]*J[y=0:290:10]
contour/hlimits=-100:400:50 xb2+zb,xb2,zb
But even with both xb2 and zb being two-dimensional, I get the same
problems with contour. I will keep studying this problem. But can you
advise further?
PM
---------------------------------------------
http://geosci.uchicago.edu/~mcguire/
---------------------------------------------
Department of the Geophysical Sciences
University of Chicago
On Thu, 17 Sep 2009, Ansley Manke wrote:
hi Patrick -
What Ferret calls an irregular grid, is just an irregularly-spaced
rectangular grid.
We have graphics commands for data on "Curvilinear grids" -
http://ferret.pmel.noaa.gov/Ferret/documentation/users-guide/working-with-special-data-sets/CURVILINEAR-COORDINATE-DATA
In that section is the definition:
A curvilinear grid has longitudes and latitudes defined by
coordinates (lon[i, j],lat[i, j]) in 2D, and the data fields are
also defined on the [i,j] index grid. ...
And the details of graphical commands for these grid types is at
http://ferret.pmel.noaa.gov/Ferret/documentation/users-guide/customizing-plots/MAP-PROJECTIONS-AND-CURVILINEAR-COORDINATES#_VPINDEXENTRY_876
Ansley
Patrick McGuire wrote:
But can't the contour command use irregular grids?
zb is an example of an irregular grid.
PM
---------------------------------------------
http://geosci.uchicago.edu/~mcguire/
---------------------------------------------
Department of the Geophysical Sciences
University of Chicago
On Thu, 17 Sep 2009, Ansley Manke wrote:
Hi Patrick,
I'm confused by your example. The variables zb and (xb+zb) are
2-dimensional arrays, 30x30, but xb is 1-d. The 3-argument contour
command expects all arguments to be 2-D arrays.
Ansley
Patrick McGuire wrote:
Hi Ferret Users
I've been having some problems with ferret recently.
The example below demonstrates the problem. I can give more
examples if you wish. But this problem affects ALL of my FERRET plots
for the methane-hydrate model spongebob, where I plot data from a
netcdf file.
FERRET EXAMPLE:
let xb=x[x=0:290:10]
let zb=I[x=0:290:10]*J[y=0:290:10]
contour/hlimits=-100:400:50 xb+zb,xb,zb
The problem seems to be when I use an irregular zb grid, sometimes
when I do contour or fill, the plot has really weird and large jumps
in the x direction of the contours. We're using ferret 6.2.
Do you have any insights? I have attached a screenshot if you
don't see the problem when you run these commands yourself. If you
don't see it on your end, you might try changing the commands
slightly to see if you can see similar problems. But this has been
reproduced by other ferret 6.2 users.
PM
---------------------------------------------
http://geosci.uchicago.edu/~mcguire/
---------------------------------------------
Department of the Geophysical Sciences
University of Chicago
------------------------------------------------------------------------