[Thread Prev][Thread Next][Index]

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



Hi Júnior, Ansley
Yes it works - I attach a picture for anybody who would like to know how it looks like.
So it means as You said Júnior that I should do a bilinear interpolation to velocity data for properly representation.


many thanks!
Simon

2010/9/30 Júnior <juniorphy@xxxxxxxxx>
Hi Simon and Ansley,

I'm also not sure if I understood your issue, but you can use the Ansley's suggestion.
Note that you have to have lat and lon dependent each other because wind also does, i.e. lat(448,615) , lon(448,615)..and you seem have already, so Ferret will plot exactly you give to him, like Ashley said :


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

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


However, if you want to change the your grid (velocity data set) to lat and lon curvilinear you could do a bilinear interpolation, for example.

Best,

Júnior

 

On Thu, Sep 30, 2010 at 7:35 AM, Szymon Roziewski <szymon.roziewski@xxxxxxxxx> wrote:
Hi,
thank You for your response.
Yes I found the answer to the first question, i.e. I can draw vectors in geographical coordinates which are curvilinear and the wind velocity.
I'm thinking of whether ferret properly changes velocity values respectively to the curvilinear grid, the input data was on uniform grid or maybe I need to change these values myself by some calculations.
Kind regards,
Simon

2010/9/29 Ansley Manke <Ansley.B.Manke@xxxxxxxx>
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



--
Z wyrazami szacunku,
Szymon Roziewski



--
----------------------------------------------------------------
Francisco das Chagas Vasconcelos Júnior
Grupo de Modelagem Atmosférica
Departamento de Meteorologia e Oceanografia
Fundação Cearense de Meteorologia e Recursos Hídricos
Av. Rui Barbosa 1246 - CEP 60115-221
Fortaleza, Ceará
Fone: (85) 3101-1106 / 3101-1126
Email: francisco.vasconcelos@xxxxxxxxxx



--
Z wyrazami szacunku,
Szymon Roziewski

Attachment: wind_2010080201_gc.gif
Description: GIF image


[Thread Prev][Thread Next][Index]

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

Privacy Policy | Disclaimer | Accessibility Statement