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