[ncl-talk] wavelet period
Shaimah ali
shaimah.ali at gmail.com
Wed Sep 4 00:40:56 MDT 2019
Hi,
I am plotting a wavelet based on this example
https://www.ncl.ucar.edu/Applications/Scripts/wavelet_2.ncl.
I need to convert the Y-axis periods in to units of months(like
12,20,40..months),instead of units of years.. So what are the changes i
have to be done in my code.
I have monthly mean data of AOD for 19.7 years.(235 values).
here is my code:
aod = asciiread("aod.dat",-1,"float")
aod!0 = "time"
ntime = dimsizes(aod)
timeo = fspan(2000.19,2019.,ntime)
aod&time = timeo
time = timeo
N = dimsizes(time)
;WAVELET
mother = 0
param = 6.0
dt = 0.25 ;timesteps in units of years
s0 = dt
dj = 0.25
jtot = 1+floattointeger(((log10(N*dt/s0))/dj)/log10(2.))
npad = N
nadof = 0
noise = 1
siglvl = .05
isigtest= 0
w =
wavelet(aod,mother,dt,param,s0,dj,jtot,npad,noise,isigtest,siglvl,nadof)
; create coodinate arrays for plot
power = onedtond(w at power,(/jtot,N/))
power!0 = "period" ; Y axis
power&period = w at period ; convert period to
units of years
power!1 = "time" ; X axis
power&time = time
power at long_name = "Power Spectrum"
power at units = "1/unit-freq"
; compute significance ( >= 1 is significant)
SIG = power ; transfer meta data
SIG = power/conform (power,w at signif,0)
SIG at long_name = "Significance"
SIG at units = " "
; initial resource settings
wks = gsn_open_wks("png","01") ; send graphics to PNG file
res = True ; plot mods desired
res at gsnDraw = True ;False ; Do not draw plot
res at gsnFrame = True ;False ; Do not advance
frome
res at gsnMaximize = True
res at cnFillOn = True ; turn on color
res at cnFillPalette = "BlAqGrYeOrReVi200" ; set color map
res at cnFillMode = "RasterFill" ; turn on raster mode
res at cnRasterSmoothingOn = True ; turn on raster smoothing
res at cnLinesOn = False ;True ; turn off contour
lines
res at cnLineLabelsOn = False
res at cnInfoLabelOn = False
res at trYReverse = False ;True ; reverse y-axis
res at tmLabelAutoStride = False ;True
res at tmYROn = False
res at tmXTOn = False
res at tmYLMode = "Explicit"
res at tmYLValues = (/1,2,4,8,16,32,64,128/)
res at tmYLLabels = (/"1","2","4","8","16","32","64","128"/)
res at tmXBMode = "Explicit"
res at tmXBValues = (/20,40,60,80,100,120,140,160,180,200,220,240/)
res at tmXBLabels =
(/"20","40","60","80","100","120","140","160","180","200","220","240"/)
res at cnLevelSelectionMode = "ExplicitLevels" ; set manual contour
levels
res at cnLevels =
(/0.01,0.02,0.04,0.08,0.1,0.12,0.14,0.16,0.20,0.24,0.28,0.32/)
res at tiMainString = ""
res at gsnRightString = ""
res at gsnCenterString = "AOD"
res at gsnLeftString = ""
res at tiYAxisString = "Periodicities (years)"
res at tiXAxisString = "Months(2000-2019)"
plot = new(2,graphic)
plot(0) = gsn_csm_contour(wks,power,res)
THANK YOU
Shaima.N
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190904/3d03e70f/attachment.html>
More information about the ncl-talk
mailing list