[Thread Prev][Thread Next][Index]
[ferret_users] cant label t axis
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
[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce /
NOAA /
OAR /
PMEL /
Ferret
Privacy Policy | Disclaimer | Accessibility Statement