[Thread Prev][Thread Next][Index]
Re: [ferret_users] Problem with time series of velocity plot
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
[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce /
NOAA /
OAR /
PMEL /
Ferret
Privacy Policy | Disclaimer | Accessibility Statement