Re: fatal:Dimension size mismatch, dimension (1) of left hand side reference does not have the same size as the right hand side reference after subscripting. fatal:Execute: Error occurred at or near line 229 in file $NCARG_ROOT/lib/ncarg/nclscripts/c

[1]
The MJO is a zonally propagating wave. The MJO diagnostics are designed
to examine that phenomena.

[2]
You can look at
http://www.ncl.ucar.edu/Document/Functions/Diagnostics/mjo_wavenum_freq_season.shtml
which is contained within

%> less $NCARG_ROOT/lib/ncarg/nclscripts/csm/diagnostics_cam.ncl
<http://www.ncl.ucar.edu/Document/Functions/Diagnostics/diagnostics_cam.shtml>

This does a "poorman's" 2D complex FFT in the longitude direction.
You could copy the function and change the code to apply to the latitude
(meridional) direction.

[3]
PLEASE change the "subject" line. Make the subject
reflect the contents of the question.

[4]
Do *not* send questions to "ncl-talk-request@ucar.edu".
I eliminated that in the response. Send questions
  "ncl-talk@ucar.edu" ... only

On 01/11/2010 12:27 AM, Sabeerali(sebi) wrote:
> I came to know that the below script(mjoclivar_10.ncl) will catch only
> zonally propagating waves. Actually I want to catch meridionally
> propagating waves. Is there any (space time spectra ) script to catch
> the meridionally propagating waves?
>
> Thanks in advance
>
>
>
> 2010/1/10 Sabeerali(sebi) <sabeerl@gmail.com <mailto:sabeerl@gmail.com>>
>
> Dear all,
>
> When I tried to calculate the frequency wavenumber spectra by
> using the below script I am getting the following errors.
>
> fatal:Dimension size mismatch, dimension (1) of left hand side
> reference does not have the same size as the right hand side
> reference after subscripting.
> fatal:Execute: Error occurred at or near line 229 in file
> $NCARG_ROOT/lib/ncarg/nclscripts/csm/diagnostics_cam.ncl
>
> fatal:Execute: Error occurred at or near line 2910 in file
> $NCARG_ROOT/lib/ncarg/nclscripts/csm/diagnostics_cam.ncl
>
> fatal:Execute: Error occurred at or near line 60 in file
> mjoclivar_10.ncl
>
> Any help would be appreciated,
> Thanking you
>
>
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/diagnostics_cam.ncl"
>
> ;++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
> begin
> ; time window
>
> twStrt = 20020102
> twLast = 20341231
>
> diri = "/home/sebi/convert/"
>
> vName = "daily_preci" ; name of variable on the file
> fili = "daily_preci.nc <http://daily_preci.nc>"
> f = addfile(diri+fili, "r")
>
> pltDir = "./"
> pltType = "pdf"
> pltName = "t61_wave_freq"
>
> ;************************************************************
> ; time indices corresponding to the desired time window
> ;************************************************************
>
> date = f->time
> iStrt = ind(date.eq.twStrt) ; desired dates
> iLast = ind(date.eq.twLast)
> delete(date)
>
> ;************************************************************
> ; Read user specified period
> ;************************************************************
>
> if (getfilevartypes(f,vName) .eq. "short") then
> X = short2flt( f->$vName$(iStrt:iLast,{-10:40},{50:150}) ) ;
> all data in window
> else
> X = f->$vName$(iStrt:iLast,{-10:40},{50:150})
> end if
> printVarSummary( X )
> printMinMax(X, True)
>
> time = X&time ; clarity
> ; date = ut_calendar( time , -2 ) ; yyyymmdd
> date=time
> ;************************************************************
> x = dim_avg_Wrap( X(time|:,lon|:,lat|:) ) ; (time,lon)
> ;************************************************************
> ; calculate and plot the spectra
> ;************************************************************
> optWavFrq = False
> optPlot = TrueoptPlot@smth9 = False
>
> nameSeason = (/"winter" , "summer"/)
>
> do ns=0,dimsizes(nameSeason)-1
> wf = mjo_wavenum_freq_season (x, date,nameSeason(ns),
> optWavFrq)
> optPlot@tiMainString ="Wave no:- Freq: spectra: "+
> nameSeason(ns)+": "+twStrt+"-"+twLast
> mjo_wavenum_freq_season_plot (wf, nameSeason(ns), pltDir,
> pltType, pltName, optPlot)
> end do
>
> end
>
>
>
>
>
>
>
>
> --
> **********************************
> Sebi
> Climate and Global Modeling Division
> Indian Institute of Tropical Meteorology
> Pashan, Pune, 411 008
> ****************************************
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>

-- 
======================================================
Dennis J. Shea                  tel: 303-497-1361    |
P.O. Box 3000                   fax: 303-497-1333    |
Climate Analysis Section                             |
Climate&  Global Dynamics Div.                       |
National Center for Atmospheric Research             |
Boulder, CO  80307                                   |
USA                        email: shea 'at' ucar.edu |
======================================================

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Jan 11 09:36:20 2010