Re: [ferret_users] cant label t axis

[Russell Fiedler <russell.fiedler@xxxxxxxx>]

Hi,

One further thing. Are you really wanting to restrict your time region from 1-jan-1979 to 1-dec-2006, i.e. the end of November?
If you don't have data for Dec 2006 then fine but you may really want to define your variables as 

let uu=uwnd[d=1,t=1-jan-1979:1-jan-2007]

etc.

Russ

On Thursday 21 July 2011 03:21, Ansley Manke wrote:
> Hi -
> Are the monthly time axes of the original three datasets the same?  If 
> not then when you mix the time axes from 2 datasets in
> 
>     let py=coeff1*(dy2+dy3)
> 
> there will be an inconsistency. You compute climatologies, then put the 
> climatology back on the monthly time axes.  Can you instead put them 
> onto a common monthly axis?
> 
> Ansley
> 
> On 7/14/2011 4:02 AM, nuncio@xxxxxxxxx wrote:
> > Hi,
> >    I am trying to run the following script but don't understand why the script stops with this error when I try to plot the x component px.
> >
> > " **ERROR: inconsistent sizes of data regions: can't label T axis of XXU"
> >
> > But I am able to plot xxu alone.  Error arises when it is multiplied with uclim_new, which is the result of transformation of climatological u velocity using @ASN.
> >
> > what is going wrong?
> >
> > regards
> > nuncio
> >
> >
> >
> >
> >
> > !----------------------------------------------------------------------------------------------
> > ! Horizontal Rossby wave-activity flux
> > !(Takaya and Nakamura (1997, 2001))
> >
> > ! NCEP/NCAR winds and geopotential heights
> >
> > ! level any
> >
> > let g=9.8
> >
> >
> > use uwnd.mon.mean_001.nc
> > use vwnd.mon.mean_001.nc
> > use hgt.mon.mean.nc
> >
> > let uu=uwnd[d=1,t=1-jan-1979:1-dec-2006]
> > let vv=vwnd[d=2,t=1-jan-1979:1-dec-2006]
> > let hh=hgt[d=3,t=1-jan-1979:1-dec-2006]
> >
> >
> > use climatological_axes
> > cancel data climatological_axes
> > let uclim=uu[gt=month_irreg@mod]
> > let vclim=vv[gt=month_irreg@mod]
> > let hclim=hh[gt=month_irreg@mod]
> >
> > let uclim_new=uclim[gt=uwnd[d=1]@asn]
> > let vclim_new=vclim[gt=vwnd[d=2]@asn]
> > let hclim_new=hclim[gt=hgt[d=3]@asn]
> >
> >
> >
> > let hhano=hh[t=1-jan-1979:1-dec-2006]-hclim_new[t=1-jan-1979:1-dec-2006] ! geopotential height anomalies
> >
> >
> > let sinlat= sin(y[g=hh]*(3.14/180.0))
> >
> > let ff=2.0*(7.24/100000)*sinlat
> >
> > let psiaa=(g/ff)*hhano                 ! stream function
> >
> >
> > let W=(uu^2+vv^2)^0.5
> >
> > let dpsi_dx=psiaa[x=@ddc]
> > let dpsi_dx_dx=dpsi_dx[x=@ddc]
> > let dpsi_dy=psiaa[y=@ddc]
> > let dpsi_dy_dy=dpsi_dy[y=@ddc]
> > let dpsi_dx_dy=dpsi_dx[y=@ddc]
> >
> > let xxu=(dpsi_dx*dpsi_dx)-(psiaa*dpsi_dx_dx)
> > let xxv=(dpsi_dx*dpsi_dy)-(psiaa*dpsi_dx_dy)
> > let yyv=(dpsi_dy*dpsi_dy)-(psiaa*dpsi_dy_dy)
> >
> > let coeff1 = z[g=hh]/(1000*2*W)
> >
> >
> >
> > let dx2=uclim_new*xxu
> > let dx3=vclim_new*xxv
> >
> > let px=coeff1*(dx2+dx3)         ! x component
> >
> >
> >
> >
> > let dy2=uclim_new*xxv
> > let dy3=vclim_new*yyv
> > let py=coeff1*(dy2+dy3)         ! y component
> >
> > !-------------------------------------------------------------------------------------------------
> >
> > National Centre for Antarctic and Ocean Research
> > Headland Sada, Vasco Da Gama
> > Goa, 403804
> > India.
> > Ph: Cell-91 9890357423
> >      off- 91 832 2525428
>