<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Hi,<br></div>I am plotting a wavelet based on this example<br><a href="https://www.ncl.ucar.edu/Applications/Scripts/wavelet_2.ncl">https://www.ncl.ucar.edu/Applications/Scripts/wavelet_2.ncl</a>.<br><br></div><div>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.<br></div><div>I have monthly mean data of AOD for 19.7 years.(235 values).<br>here is my code:<br><br>aod = asciiread("aod.dat",-1,"float")<br>  aod!0 = "time"<br>  ntime = dimsizes(aod)<br>  timeo = fspan(2000.19,2019.,ntime)<br>  aod&time = timeo<br><br>  time  = timeo<br>  N     = dimsizes(time)<br>  <br>;WAVELET<br><br>  mother  = 0<br>  param   = 6.0<br>  dt      = 0.25   ;timesteps in units of years<br>  s0      = dt<br>  dj      = 0.25 <br>  jtot    = 1+floattointeger(((log10(N*dt/s0))/dj)/log10(2.))<br>  npad    = N<br>  nadof   = 0<br>  noise   = 1<br>  siglvl  = .05<br>  isigtest= 0<br><br>  w = wavelet(aod,mother,dt,param,s0,dj,jtot,npad,noise,isigtest,siglvl,nadof)<br><br>; create coodinate arrays for plot<br><br>  power            = onedtond(w@power,(/jtot,N/))<br>  power!0          = "period"                        ; Y axis<br>  power&period     = w@period                      ; convert period to units of years<br><br>  power!1          = "time"                          ; X axis<br>  power&time       = time<br><br>  power@long_name  = "Power Spectrum"<br>  power@units      = "1/unit-freq"<br><br>  <br>; compute significance ( >= 1 is significant)<br>  SIG              = power                            ; transfer meta data<br>  SIG              = power/conform (power,w@signif,0)<br>  SIG@long_name    = "Significance"<br>  SIG@units        = " "<br><br>; initial resource settings<br><br>  wks = gsn_open_wks("png","01")             ; send graphics to PNG file<br><br>  res                     = True                  ; plot mods desired<br>  res@gsnDraw             = True ;False                 ; Do not draw plot<br>  res@gsnFrame            = True ;False                 ; Do not advance frome<br>  res@gsnMaximize         = True <br>  res@cnFillOn            = True                  ; turn on color<br>  res@cnFillPalette       = "BlAqGrYeOrReVi200"   ; set color map<br>  res@cnFillMode          = "RasterFill"          ; turn on raster mode<br>  res@cnRasterSmoothingOn = True                  ; turn on raster smoothing<br>  res@cnLinesOn           = False ;True                 ; turn off contour lines<br>  res@cnLineLabelsOn      = False<br>  res@cnInfoLabelOn       = False<br>  res@trYReverse          = False  ;True                  ; reverse y-axis <br>  res@tmLabelAutoStride   = False   ;True<br>  res@tmYROn              = False <br>  res@tmXTOn              = False<br><br>  res@tmYLMode = "Explicit"<br><br>  res@tmYLValues = (/1,2,4,8,16,32,64,128/)<br>  res@tmYLLabels = (/"1","2","4","8","16","32","64","128"/)<br>  <br>  res@tmXBMode = "Explicit"<br>  res@tmXBValues = (/20,40,60,80,100,120,140,160,180,200,220,240/)<br>  res@tmXBLabels = (/"20","40","60","80","100","120","140","160","180","200","220","240"/)<br>  <br>  <br>  res@cnLevelSelectionMode = "ExplicitLevels"       ; set manual contour levels<br>  res@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/)<br><br>  res@tiMainString       = ""<br>  res@gsnRightString      = ""<br>  res@gsnCenterString     = "AOD"<br>  res@gsnLeftString       = ""<br>  res@tiYAxisString       = "Periodicities (years)"<br>  res@tiXAxisString       = "Months(2000-2019)"<br><br>  plot = new(2,graphic)<br>  plot(0) = gsn_csm_contour(wks,power,res)<br><br><br></div><div>THANK YOU<br></div><div>Shaima.N<br></div></div></div></div>