[no subject]

[Marcelo Barreiro <marcelo@xxxxxxxxxxxxxx>]

Hi Emilio,

Thanks for your suggestion, it worked. I did

? set data etopo120
? set region/x=59W:50W/y=35S:27S
? let Prec_filt = if rose[d=2] ge -500 then Prec[d=1]
? contour/fill/l=5/I=141:145/J=28:32 Prec_filt           !to plot

I had to put "I=141:145/J=28:32" because otherwise ferret sends the error 
message

 **ERROR: inconsistent sizes of data regions: X axis
           PREC[D=1] has 4 points (I=108:111)
           expression has 5 points (I=141:145)

and an analog for J.

Then I calculated the mean with

? let/quiet Prec_ave=Prec_filt[d=1,x=@AVE,y=@AVE]

and try to plot it using 

? plot Prec_ave

and I got the same error as above. I can not understand why do I get this 
message because Prec_ave is one-dimensional (only in time). I tried to 
use 

plot/I=141:145/J=28:32 Prec_ave

but ferret displays again the same error message.

Do you have any idea on how to plot the time series of the average (of 
Prec_ave)???

Thanks. (Muchas gracias!)

Marcelo



================================================
Marcelo Barreiro
Dept. of Oceanography, Physical Section
College of Geosciences and Maritime Studies
Texas A&M University
College Station, TX 77843-3146
Phone: (409) 845-3238
Email: marcelo@ocean.tamu.edu
URL: http://www-ocean.tamu.edu/~marcelo/
================================================


On Wed, 25 Aug 1999, Emilio Mayorga wrote:

> 
> Marcelo,
> 
> Ferret doesn't support non-rectangular regions per se. Depending on what
> you want to do, "masks" may work for you. Let's say you have a variable
> that has positive values in your continental region, and negative values
> (or nodata) outside of it. Set your rectangular region:
> 
> SET REGION ... 
> 
> then create a masked copy of your variable of interest, MYVAR:
> 
> LET MYVARMASK = IF continentvar ge 0 THEN MYVAR
> 
> Now use the new variable MYVARMASK instead of the original MYVAR. Areas
> outside of the continent have been set to nodata, and will not be included
> in calculations.
> 
> Hope that helps. Saludos,
> 
> Emilio Mayorga
> 
> ________________________________________________________________
> University of Washington	
> School of Oceanography			emiliom@u.washington.edu
> Box 357940				ph. (206) 543-5334
> Seattle, WA  98195-7940	  USA		FAX (206) 685-3351
> http://boto.ocean.washington.edu/staff_pages/emilio
> 
> On Wed, 25 Aug 1999, Marcelo Barreiro wrote:
> 
> > 
> > Hi all, 
> > 
> > I would like to work in a continental region that is not rectangular. Is 
> > there any way to set an XY-region which is not rectangular?
> > 
> > Any help will be greatly appreciated.
> > 
> > Marcelo
> > 
> > 
> > ================================================
> > Marcelo Barreiro
> > Dept. of Oceanography, Physical Section
> > College of Geosciences and Maritime Studies
> > Texas A&M University
> > College Station, TX 77843-3146
> > Phone: (409) 845-3238
> > Email: marcelo@ocean.tamu.edu
> > URL: http://www-ocean.tamu.edu/~marcelo/
> > ================================================
> > 
> > 
> > 
>