Re: [ferret_users] Plotting variable with alternate time axis

[Andrew Wittenberg - NOAA Federal <andrew.wittenberg@xxxxxxxx>]

Hi Paul,

Simple answer is:

define axis/t timeaxnew = t[gt=sos] + 10
let sos_new = sos[gt=timaxnew@asn]

Or here's a fancy version, to keep all the axis and variable attributes:

def ax/t/unit="`sos,r=tunit`"/cal="`sos,r=cal`"/t0="`sos,r=t0`" timeaxnew = 
t[gt=sos] + 10
let/title="`sos,r=title`"/unit="`sos,r=unit`" sos_new = sos[gt=timeaxnew@asn]

Andrew

On 03/16/2013 05:04 PM, Durack, Paul J. wrote:
> Hi there folks..
>
> I have a super noob question which I've been unable to resolve without
> assistance..
>
> I have a variable which has a time coordinate in UTC, I want to plot this
> using ferret with the time indicated automagically to show the local +10
> (not UTC) time, so in effect adding 10hrs to each element of the time
> axis.. I'm sure this is really simple, however I'm really struggling to
> get my syntax correct
>
>
> So I've tried:
>
> [dir]$ ferret -nojnl 	NOAA/PMEL TMAP
> 	FERRET v6.82
> 	Linux 2.6.32-279.1.1.el6.x86_64 64-bit - 08/03/12
> 	16-Jan-13 13:02
>
> yes? use infile.nc
> yes? sh da
>       currently SET data sets:
>      1> ./infile.nc  (default)
> name     title                             I         J         K
> L         M         N
> SOS      Sea Surface Salinity             1:192     1:189     ...
> 1:156     ...       ...
> yes? show grid sos
>      GRID GFG1
> name       axis              # pts   start                end
> LON       LONGITUDE          192mr   0E                   1.875W
> LAT       LATITUDE           189 i   88.195S              87.24N
> normal    Z
> TIME      TIME               156 r   02-JUL-1850 12:00    02-JUL-2005
> 12:00
> normal    E
> normal    F
> yes? let timeax = T[GT=sos]
> yes? define axis/T timeaxnew=timeax+10
> yes? show grid sos
>      GRID GFG1
> name       axis              # pts   start                end
> LON       LONGITUDE          192mr   0E                   1.875W
> LAT       LATITUDE           189 i   88.195S              87.24N
> normal    Z
> TIME      TIME               156 r   02-JUL-1850 12:00    02-JUL-2005
> 12:00
> normal    E
> normal    F
>
>
> Now I have the data (sos variable) and I have my new time axis (timeaxnew,
> which I presume has been incremented by 10 hrs), how do I now plot this
> using the new axis?
>
>
> I've been using:
> yes? sha sos[l=1] ! How do I use a new time axis, or add 10hrs to each
> timestep?
>
>
> Some help cleaning up my syntax here would certainly be appreciated!
>
> Cheers,
>
> P