[Thread Prev][Thread Next][Index]

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



Hi,
Here are the info.

use pycnoB_2.tsout3N
yes? sho grid uvel
    GRID GSF4
 name       axis              # pts   start                end
 XDIM      X                    1 r   1                    1
 YDIM      Y                    1 r   1                    1
 ZDIM      Z                  200 r   1                    200
 TDIM      T                    1 r   1                    1
 normal    E
 normal    F
yes? stat uvel
 
             X-component_of_current
             X: 1
             Y: 1
             Z: 0.5 to 200.5
             T: 1
             E:  N/A
             F:  N/A
             DATA SET: ./pycnoB_2.tsout3N
 
 Total # of data points: 200 (1*1*200*1*1*1)
 # flagged as bad  data: 0
 Minimum value: -0.0077385
 Maximum value: 0.12014
 Mean    value: 0.0047735 (unweighted average)
 Standard deviation: 0.016147
yes? shade uvel

Actually, I'm trying to plot the time series of uvel in a column of water with 200 layers, and there is not any x or y.( I only have z). Also, I attached the uvel time series that I plotted at first and I had problem with in the second file. and normal uvel on depth on the 3rd enclosed.

Thank you,


On Tuesday, July 8, 2014 6:55 PM, Billy Kessler <william.s.kessler@xxxxxxxx> wrote:


What is the shape of uvel?
sho grid uvel
stat uvel
shade uvel
Are there a lot of blanked values?
Does the shape of uvel make sense with the reshaping you're doing? Remember, RESHAPE is just wrapping one grid onto another. 

Similarly, what is the corresponding info for uvel_on_time?

I don't have enough information to diagnose this, but I'm sure ferret is doing something logical. By plotting and listing and using sho grid and stat, you should be able to figure it out. 

Billy K

Sent from my phone .... Please excuse my brevity .... and my typing. 

On Jul 8, 2014, at 4:28 PM, "'Sara Sari' via _OAR PMEL Ferret Users" <ferret_users@xxxxxxxx> wrote:

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


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

Privacy Policy | Disclaimer | Accessibility Statement