Hi Sara,
Instead of RESHAPE, use a regridding transformation.
Regridding in T
will keep the grid of the variable the same in all the
other dimensions.
RESHAPE is for changing from one entire grid to another.
It will just
take the first 7969 points of temp and put them on the
new grid.
Because you are going from one time axis to another with
the same number
of points, but not from one calendar to another, you
want to use the
@ASN regridding transformation.
Let temp_on_time=temp[gt=
days_axis@ASN]
The axis you have defined will be labeled with just
T=1:365. If you are
asking how to get calendar time labels, then you need a
/T0 qualifier
on your DEFINE AXIS command.
If all of this doesn't answer the question then please
write back with
more information. I'm not sure which previous message
you're talking
about. Instead of referring to your last email, please
show the grid of
your variable.
Ansley
On 12/16/2014 9:25 PM, 'Sara Sari' via _OAR PMEL
Ferret Users wrote:
> Hi everyone,
> I forgot to mention in my last email that I
actually want to reset the time axis from 7969 steps
to 365 days while plotting for a specific depth. I use
this code but it ignores the depth specifications.
> Def axis/t=1:365/npoint=7969/units=days Days_axis
>
> Yes? Let
temp_on_time=RESHAPE(temp,[gt=temp_on_time])
>
> Yes? Plot temp_on_time
>
>
>
> ...now i want to plot temp graph for water layer
20. (Z=20), but don't know how to include it! Normally
if I write
> plot temp[z=20] it plots for this layer but taxis
is in time step. I want to keep the same command but
just on taxis in days unit.
> Could any one please help me.
> Thank you,
> Sara