[Thread Prev][Thread Next][Index]

Re: [ferret_users] add (degenerate) axis to variable




Satyesh,

My axis shows up as a single point:

yes? sho ax t1
 name       axis              # pts   start                end
 T1        TIME                 1 r   15-MAR-2000 00:00    15-MAR-2000 00:00
T0 = 1-JAN-1970
   Axis span (to cell edges) = 1

That seems right.

-Andy

On 3/31/21 12:43 PM, Satyesh Ghetiya wrote:
You can enter:
show grid thist;

Also in axis definition,  try t=1-jan-1970-1-jan-1970
Also try same as yours but add how many days are there, some syntax is there, I m forgetting it. Just add that you have 1 day in your time axis  .
At last check by " show grid T1" command


On Wed, 31 Mar 2021, 23:45 Andy Jacobson, <andy.jacobson@xxxxxxxx <mailto:andy.jacobson@xxxxxxxx>> wrote:

    Satyesh,


    Thanks, that's interesting. However, I get the same behavior using whatever syntax. No complaint from ferret, but time axis is also not present in the output:

     > let thist = days1900(${yr},${mon},15) - days1900(1970,1,1)   ! ${yr} and ${mon} are shell variables
     > define axis/T/T0=1-jan-1970/UNIT=days time=thist

     > let/title="GOSIF v2"/units="W m−2  microm−1 sr−1" sif=sif0[gx=regions[d=2],gy=regions[d=2],gt=time@asn]*0.0001

    or

     > let/title="GOSIF v2"/units="W m−2  microm−1 sr−1" sif=sif0[gx=regions[d=2],gy=regions[d=2],gt=time]*0.0001

    or

     > let/title="GOSIF v2"/units="W m−2  microm−1 sr−1" sif=sif0[gx=regions[d=2],gy=regions[d=2],t=time]*0.0001

     > save/clobber/file="$outfile" sif      ! $outfile is a shell variable

    -Andy

    On 3/31/21 11:59 AM, Satyesh Ghetiya wrote:
     > No, not =@asn. Just =T1. (@asn works on some already defined time axis of same L steps)
     >
     > On Wed, 31 Mar 2021, 22:31 Andy Jacobson, <andy@xxxxxxxx <mailto:andy@xxxxxxxx> <mailto:andy@xxxxxxxx <mailto:andy@xxxxxxxx>>> wrote:
     >
     >     Yes, I tried [...gt=t1@asn] syntax. Ferret didn't complain, but it didn't output that time axis either.
     >
     >     Thanks,
     >
     >     Andy
     >
     >     On 3/31/21 10:44 AM, Satyesh Ghetiya wrote:
     >      > Have you tried giving that time axis to variable by GT ? -->
     >      > let your variable is Temp and time axis name is T1, then -->
     >      > let Temp_new=Temp[gt=T1]
     >      >
     >      > ^ This should work.
     >      >
     >      > On Wed, Mar 31, 2021 at 10:04 PM Andy Jacobson <andy.jacobson@xxxxxxxx <mailto:andy.jacobson@xxxxxxxx> <mailto:andy.jacobson@xxxxxxxx <mailto:andy.jacobson@xxxxxxxx>> <mailto:andy.jacobson@xxxxxxxx <mailto:andy.jacobson@xxxxxxxx> <mailto:andy.jacobson@xxxxxxxx <mailto:andy.jacobson@xxxxxxxx>>>> wrote:
     >      >
     >      >     Howdy,
     >      >
     >      >     I have a 2-d variable (lat, lon) that I read from a well-formatted netCDF file. This represents one month out of a time series of many such files. I'd like to add a time axis to this variable. I've defined my single-point time axis, but I can't figure out how to convince ferret to add this time axis to my variable. Normally I'd do this in NCO but I'd like to stay in ferret this time. Is there some reshape trickery involved?
     >      >
     >      >     Thanks,
     >      >
     >      >     Andy
     >      >
     >      >     --
     >      >     Andy Jacobson
     >      > andy.jacobson@xxxxxxxx <mailto:andy.jacobson@xxxxxxxx> <mailto:andy.jacobson@xxxxxxxx <mailto:andy.jacobson@xxxxxxxx>> <mailto:andy.jacobson@xxxxxxxx <mailto:andy.jacobson@xxxxxxxx> <mailto:andy.jacobson@xxxxxxxx <mailto:andy.jacobson@xxxxxxxx>>>
     >      >
     >      >     NOAA Global Monitoring Laboratory
     >      >     325 Broadway R/GMD1
     >      >     Boulder, Colorado 80305
     >      >
     >      >     303/497-4916
     >      >
     >      >
     >      >
     >      > G Satyesh,
     >      > India
     >      >
     >
     >     --
     >     Andy Jacobson
     > andy@xxxxxxxx <mailto:andy@xxxxxxxx> <mailto:andy@xxxxxxxx <mailto:andy@xxxxxxxx>>
     >

-- Andy Jacobson
    andy.jacobson@xxxxxxxx <mailto:andy.jacobson@xxxxxxxx>

    NOAA Earth System Research Lab
    Global Monitoring Division
    325 Broadway R/GMD1
    Boulder, Colorado 80305

    303/497-4916


--
Andy Jacobson
andy.jacobson@xxxxxxxx

NOAA Earth System Research Lab
Global Monitoring Division
325 Broadway R/GMD1
Boulder, Colorado 80305

303/497-4916
	

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

Privacy Policy | Disclaimer | Accessibility Statement