[Thread Prev][Thread Next][Index]

Re: [ferret_users] Re: inconsistent sizes of data regions: Y axis of X position array



Hi Simon,
When  you say you found it, do you mean you found the answer to your first question?
 
I'm not sure I understand what you are asking here. SET GRID is not relevant when you have data on curvilinear grids like this. If vel is a variable whose curvilinear grid is described by LON from one dataset and LAT, then you should be able to

SHADE/palete=land_sea/title=""/levels=16 vel[d=1], lon[d=2], lat[d=3]

with the values of d=  set coorectly, of course.

Ansley

On 9/29/2010 10:16 AM, Szymon Roziewski wrote:
Hi,
I found it,
but I have another question,
How can I make a grid from two matrices:um_c3_lon_448x615.tx - longitude, um_c3_lat_448x615.tx - latitude - 448x615 both.
This matrix I would use to make a plot on this grid like that:
set grid gridlonlat !(the grid made from two matrices which I discussed above)
SHADE/palete=land_sea/title=""/levels=16 vel ! wind velocity in geographical coordinates

Best regards,
Simon

2010/9/29 Szymon Roziewski <szymon.roziewski@xxxxxxxxx>
hello there,

I have such a problem, I want to draw vectors (wind) in geographical coordinates which are in files (448 columns and 615 rows - longitude, latitude)
The problem is that when I try to draw by this commend:
vector/over/title="WIND"/length=25/xskip=10/yskip=10 wlon[d=1], wlat[d=2],lon[d=3],lat[d=4]
I get the message:
     currently SET data sets:
    1> /home/szymon/WAM/CHECK/03225_2010092100+03000000c3hs000000000000000
 name     title                             I         J         K         L
 WLON     WLON                             1:448     1:615     ...       ...
 
    2> /home/szymon/WAM/CHECK/03226_2010092100+03000000c3hs000000000000000
 name     title                             I         J         K         L
 WLAT     WLAT                             1:448     1:615     ...       ...
 
    3> /home/szymon/WAM/CHECK/um_c3_lon_448x615.txt
 name     title                             I         J         K         L
 LON      LON                              1:448     1:615     ...       ...
 
    4> /home/szymon/WAM/CHECK/um_c3_lat_448x615.txt  (default)
 name     title                             I         J         K         L
 LAT      LAT                              1:448     1:615     ...       ...
 
 **ERROR: inconsistent sizes of data regions: Y axis of X position array
vector/over/title="WIND"/length=25/xskip=10/yskip=10 wlon[d=1], wlat[d=2],lon[d=3],lat[d=4]
Command file, command group, or REPEAT execution aborted
STOP -script mode, ERROR RUNNING SCRIPT

but I dont know why... I think dimensions are correct.

Below I post my script. And I attach appropriate files to this letter.
I appreciate any solutions!

CANCEL MODE LOGO
DEFINE AXIS/x=1:448:1/unit=degree xaxis
DEFINE AXIS/y=1:615:1/unit=degree yaxis

DEFINE GRID/x=xaxis/y=yaxis gridlonlat

!set grid gridlonlat

file/var=wlon/grid=gridlonlat/format=stream/type=r4/swap "/home/szymon/WAM/CHECK/03225_2010092100+03000000c3hs000000000000000"
file/var=wlat/grid=gridlonlat/format=stream/type=r4/swap "/home/szymon/WAM/CHECK/03226_2010092100+03000000c3hs000000000000000"

let vel=(wlon[d=1]^2+wlat[d=2]^2)^0.5

FILE/VARIABLES=lon/COLUMNS=448/GRID=gridlonlat "/home/szymon/WAM/CHECK/um_c3_lon_448x615.txt"
FILE/VARIABLES=lat/COLUMNS=448/GRID=gridlonlat "/home/szymon/WAM/CHECK/um_c3_lat_448x615.txt"

!SHADE/palete=land_sea/title=""/levels=16 vel

show data

!list lat[d=4]

vector/over/title="WIND"/length=25/xskip=10/yskip=10 wlon[d=1], wlat[d=2],lon[d=3],lat[d=4]


frame/file=umwind.gif


Kind regards,
Szymon Roziewski



--
Z wyrazami szacunku,
Szymon Roziewski

[Thread Prev][Thread Next][Index]

Contact Us
Dept of Commerce / NOAA / OAR / PMEL / TMAP

Privacy Policy | Disclaimer | Accessibility Statement