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
This archive was generated by hypermail 2.1.8 : Thu Feb 09 2012 - 13:33:26 MST