[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