[Thread Prev][Thread Next][Index]

[ferret_users] Randomly distributed time series plots



Dear all,

Sometimes I need to do a 'quick check' across a large area and visualize the time series of a particular variable at random locations to get an idea of the general trend through time.

At the moment I am using a not-so-flexible approach, such as

let indice_a = {1, 1, 1, 2, 2, 2, 3, 3, 3, 4, 4, 4}
let indice_b = {1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3}
let indice_x = {73,75,77,79}
let indice_y = {11,13,15}
go multi_view 4,3

repeat/RANGE=1:10/NAME=m (let a = indice_a[i=`m`] ;\
let b = indice_b[i=`m`] ;\
let x = indice_x[i=`m`] ;\
let y = indice_y[i=`m`] ;\
set view V`a``b`;\
plot my_variable[x=`x`,y=`y`,k=2]

This nested loop is not working, but it's a concept I tried to develop based on a single loop. As you can see it's rather tedious to write these indices and not easily adaptable to any region when lat and long can vary. Also I need two or three nested loops, which for now I can't make them work.

Additionally, I want to use the maximum allowed by multi_view, which I think is 15 windows, and do it multiple times in a loop if I need to sample a large region.

Is there anyway to make this more automatic? Such as sampling randomly lat and long from the available range? For example In R I would use the sampling function to create indices for x and y, and then use them in a loop.

I am sure other people have a need to check spatial change in a variable across time. Any suggestion or idea is well appreciated.

Thanks

[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce / NOAA / OAR / PMEL / Ferret

Privacy Policy | Disclaimer | Accessibility Statement