Re: [ferret_users] how to do this with "repeat"?

[Ansley Manke <ansley.b.manke@xxxxxxxx>]

Hi Antonio,
If I understand correctly what you want to do, it would be this, 
averaging 1:10 in time

let sst_fav_tav =sst13[X=10w:4w@FAV:7,Y=35n:38n@FAV:7,L=1:10@AVE]

Or, if you want the x and y averaging at each time, you would just 
define the fill in x and y and then use the desired time steps to save 
or in whatever other calculations you want to do.

  let sst13_fav = sst13[X=10W:4w@FAV:7,Y=35n:38n@FAV:7]
  save/file=myfile.nc/L=1:10 sst13_fav

In any case this shows how to apply the @FAV transformation to a region 
in x and y.

Most often, the REPEAT command is not needed. Ferret operates over 
entire grids or axes, or the subsets you specify with qualifiers /X= 
/Y=  or SET REGION/X=/Y=  or in square brackets [T=]. 

Ansley


arv@xxxxxxx wrote:
> Hi,
>
> I have a parameter called sst13 where I want to perform a @FAV 
> for the X and Y axes averaging along the time axis so I get I new 
> parameter, i.e.
> "sst13_fav" with the new values averaged. But don't know how to put
> this into words. I've tried something like this (and around this), but
> doesn't work:
>
> repeat/l=1:10 let s=sst13[X=10W:4w,Y=35n:38n,[X=@FAV:7,Y=@FAV:7]
>
> yes? show data/variables
> currently SET data sets:
>     1> ./sst13.nc  (default)
> name     title                             I         J        
> K         L
> SST13    PHSSTA[X=12W:2W,Y=34N:49N,L=120  1:201     1:301     1:
> 1       1:99
>                on grid GHJ1 with -1.000000E+34 for missing data
>              X=12W:2W(358)  Y=34N:49N  Z=-0.5:0.5
>
>   time range: 17-SEP-2004 00:00 to 21-NOV-2006 00:00
>
> Best regards,
> Antonio
>   

--

~>   ~>   ~>   ~>   ~>   ~>   ~>   ~>   ~>   ~>   ~> 

Ansley Manke, NOAA/PMEL   ansley.b.manke@xxxxxxxx
7600 Sand Point Way NE, Seattle WA
Phone 206-526-6246,  FAX 206-526-6744