It looks like your array UVEL is actually in (z,t), i.e. 2-dimensional (or more).
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.
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