[Thread Prev][Thread Next][Index]

Re: [ferret_users] Problem with time series of velocity plot

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

Attachment: uvel on time and z.png
Description: PNG image

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

Privacy Policy | Disclaimer | Accessibility Statement