Re: [ferret_users]:capricious reshape function

[Ansley Manke <ansley.b.manke@xxxxxxxx>]

Hi Xavier,
If you are just defining a new time axis with the same number of points 
as the existing axis, you can use a regridding operation to put your 
data onto that time axis. The @ASN regridding operation just assigns the 
data to the new axis.

  define AXIS/T/T0="01-jan-0000"/unit="days" taxis=a
  let Var3=vbar[gt=taxis@ASN]

I'm not sure what you are doing with the climatological_axes, but 
generally when you are going to use it for a regridding operation, you 
want to make sure to CANCEL that data set after opening it. All you are 
doing is defining the axes it contains so you can use them later in your 
script. That could possibly be why you are getting a "Variable not 
found" error, if you have not canceled the climatological_axes.cdf  
dataset and it is defined as the default data set. Though here, you seem 
to be specifying [d=1] everywhere so that may not be what's happening.


Xavier Couvelard wrote:
> Hi all
> I'am regriding my model output to obtain climatological axes
> i have a lots of file that leaded me to use a descriptor file .
> and i also use the climatological axes
>
> SET MEM/SIZE=500
>
> set data avg.0.mc
>
> use climatological_axes.cdf
>
> i first define a new time axes and reshape my variable
> Define AXIS/T/T0="01-jan-0000"/unit="days" taxis=a
>
> let Var3=reshape(vbar[d=1],t[gt=taxis]+x[gx=vbar,d=1]+y[gy=vbar,d=1])
> .............
> .............
> So, while it works very well with my two first variable (zeta and 
> ubar) with vbar (2 and the two others variable) it don't
>
> for exemple
> sha var3[l=1]
> ** netCDF error: Variable not found
> while sha var2[l=1] who is ubar make a good plot and vbar[d=1,l=1] 
> also make a good plot.
>
> is someone has an idea?
>
> Xavier

--

~>   ~>   ~>   ~>   ~>   ~>   ~>   ~>   ~>   ~>   ~> 

Ansley Manke, NOAA/PMEL   ansley.b.manke@xxxxxxxx
7600 Sand Point Way NE, Seattle WA
Phone 206-526-6246,  FAX 206-526-6744