Re: space-time spectra over a specified domain

This and other space-time functions in this package
assume the data are periodic in longitude. I suggest
altering the data that are input to just the region of interest.

   latS = -10. ; latitude band
   latN = 10. ; symmetric about equator

   lonL = 30. ; longitude segment
   lonR = 120.

Then later,

   x = dim_avg_n_Wrap( X, 1 ) ; (time,lon)
   printVarSummary(x)

;************************************************************
; Create necessary regional data
; (a) Remove zonal (longitudinal) means at each time
; (b) Set all data outside chosen longitudes to 0.0
; (c) Taper the data in the region
;************************************************************
   xZon = dim_avg_n_Wrap( x, 1 )
   printVarSummary(xZon) ; xZon(time)

   x = x - conform(x, xZon, 0) ; remove zonal mean at each time step
   printVarSummary(x) ; [time | 9497] x [lon | 144]

   LON = conform(x, x&lon, 1)
   x = where(LON.lt.lonL .or. LON.gt.lonR, 0.0, x)

   x(:,{lonL:lonR}) = taper(x(:,{lonL:lonR}), 0.1, 0)

On 2/7/12 10:56 PM, Bidyut Bikash Goswami wrote:
> Hi,
>
> I want to plot the space-time spectra over a specified domain(say indian region) using
> http://www.ncl.ucar.edu/Document/Functions/Diagnostics/mjo_wavenum_freq_season.shtml (using the script mjoclivar_10.ncl).
>
> In order to do this, please suggest me regarding the required modifications needed to be done in the file http://www.ncl.ucar.edu/Document/Functions/Diagnostics/diagnostics_cam.shtml. Please specify the places where I need to put the longitude boundaries, in diagnostics_cam file.
>
> Regards,
>
> Bidyut Bikash Goswami
> Research Fellow
> Indian Institute of Tropical Meteorology
> Dr. Homi Bhabha Road
> Pashan
> Pune-411008
> Telephone:+91 20 25904226 (O)
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Feb 8 10:40:45 2012