Hi Sara,
What are you trying to accomplish? In your original question you
had these commands:
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
Your data is on a grid with 200 levels in Z and one T. The above
commands will define a grid which consists only of a time series,
and treat the data as a time series of 200 points. So yes, as far as
the plot is concerned you are just relabeling the axis on which the
data varies, and plotting it as if it varies in T not in Z.
If this data is intended to be a time series at 200 layers, then
perhaps are there other files which contain other timesteps?
Ansley
On 7/9/2014 3:13 AM, 'Sara Sari' via
_OAR PMEL Ferret Users wrote:
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,
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.
Thank you William, but still I
don't get the time series. please see
the result. it shows that T is 0 (hour)!
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>
|