Thank you William, but still I don't get the time series. please see the result. it shows that T is 0 (hour)!On Tuesday, July 8, 2014 2:08 PM, William S. Kessler <william.s.kessler@xxxxxxxx> wrote:
It looks like your array UVEL is actually in (z,t), i.e. 2-dimensional (or more).
In that case you must tell RESHAPE the complete structure of the new array. That is, if the z-axis is called DEPTH, then:
> let uvel_on_time = RESHAPE(uvel,t[gt=hours_axis]+z[gz=depth]) ! plus any other axes
Remember the 2nd argument to RESHAPE does nothing more than specify the _shape_ of the output array, but that shape must be complete.
Billy K
On Jul 8, 2014, at 12:46 PM, 'Sara Sari' via _OAR PMEL Ferret Users wrote:
> Dear Ferret Users,
> I use the following command to plot time series of x-component of velocity, UVEL, (in 200 layers). But, it plots velocity on y-axis and layers on x-axis, only changing the x-axis lable to time while it actually is layers instead of time (seems only the label has changed). Could anyone let me know what the problem is and how I can fix it to have a real time series plot of velocity for 24 hours?
>
> def axis/t=0:24/npoints=200/units=hours hours_axis
> let uvel_on_time = RESHAPE(uvel,t[gt=hours_axis])
> plot uvel_on_time
>
> Regards,
> Sara
<uvel on time and z.png>
Attachment:
shade uvel.png
Description: PNG image
Attachment:
uvel time series.png
Description: PNG image
Attachment:
uvel.png
Description: PNG image