Hi Patrick,
This was a bug, and we have a fix for it, which will be in the upcoming
Ferret v6.3 release.
let xb = x[x=0:290:10]
let yb = y[y=0:290:10]
let xb2 = xb + 0*yb
let zb = I[x=0:290:10] * J[y=0:290:10]
contour/hlimits=-100:400:50 xb2+zb,xb2,zb
Ansley Manke wrote:
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
------------------------------------------------------------------------
|