How about if you send the function enough data around the location you want, so that it can compute the result. Not knowing what it needs, say,
> let radial = my_function(SST[j=190:210])> plot/j=100 radial
Ferret just loads the data requested in the argument to the function, whether it is all of SST, or SST[j=1:200]. The argument subscripts are computed from the context of the arguments that are sent in.
Possibly you could set up the function to work on only some of the data by passing more arguments,
let radial = my_function(SST, xout, yout)
where xout and yout are the desired output range in two directions (or single values). The function would have access to the full data in SST, but it could use xout and yout to determine the amount of computation it needs to do.
On 1/9/2015 8:03 AM, Serena Illig wrote: