<div dir="ltr">Hi All,<div><br></div><div>I&#39;m plotting a wavelet based from this example <a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/wavelet.shtml">https://www.ncl.ucar.edu/Document/Functions/Built-in/wavelet.shtml</a></div><div><br></div><div>Is there an easy way to make the x-axis in terms of months (Jan to Dec)?</div><div><br></div><div>Here&#39;s my code:</div><div><br></div><div><div>begin</div><div>;************************************</div><div>; create pointer to file and read in variables</div><div>;************************************</div><div>  f    = addfile (&quot;hpass_SCS.nc&quot;,&quot;r&quot;)</div><div>  prec  = f-&gt;PREC(:,0,0)                     ; one value/year</div><div>  time = prec&amp;time</div><div>  N    = dimsizes(time)</div><div>;************************************</div><div>; Remove mean and linear trend (recommended, not required)</div><div>; Missing values are not allowed.</div><div>;************************************</div><div>  prec  = dtrend(prec,False)</div><div>;************************************</div><div>; compute wavelet (Missing values not allowed due to use of FFT)</div><div>;************************************</div><div>  mother  = 0</div><div>  param   = 6.0</div><div>  dt      = 1.</div><div>  s0      = 0.25</div><div>  dj      = 0.25</div><div>  jtot    = 44</div><div>  npad    = N</div><div>  nadof   = new(2,float)</div><div>  noise   = 1</div><div>  siglvl  = .05</div><div>  isigtest= 1</div><div><br></div><div>  w = wavelet(prec,mother,dt,param,s0,dj,jtot,npad,noise,isigtest,siglvl,nadof)</div><div>;************************************</div><div>; create coodinate arrays for plot</div><div>;************************************</div><div>  power            = onedtond(w@power,(/jtot,N/))</div><div>  power!0          = &quot;period&quot;                        ; Y axis</div><div>  power&amp;period     = w@period</div><div><br></div><div>  power!1          = &quot;time&quot;                          ; X axis</div><div>  power&amp;time       = time</div><div><br></div><div>  power@long_name  = &quot;Power Spectrum&quot;</div><div>  power@units      = &quot;unit^2&quot;</div></div><div><br></div><div><br></div><div><div>;************************************</div><div>; initial resource settings</div><div>;************************************</div><div>  wks = gsn_open_wks(&quot;png&quot;,&quot;wavelet&quot;)             ; send graphics to PNG file</div><div>  res                     = True                  ; plot mods desired</div><div>  res@cnFillOn            = True                  ; turn on color</div><div>  res@cnFillPalette       = &quot;BlAqGrYeOrReVi200&quot;   ; set color map</div><div>  res@cnFillMode          = &quot;RasterFill&quot;          ; turn on raster mode</div><div>  res@cnRasterSmoothingOn = True                  ; turn on raster smoothing</div><div>  res@cnLinesOn           = False                 ; turn off contour lines</div><div>  res@lbOrientation       = &quot;Vertical&quot;            ; vertical label bar</div><div>  res@trYReverse          = True                  ; reverse y-axis</div><div>  ;res@tmLabelAutoStride   = True</div><div>  res@cnLevelSelectionMode = &quot;ManualLevels&quot;       ; set manual contour levels</div><div>  res@cnMinLevelValF       = 0.0                  ; set min contour level</div><div>  res@gsnLeftString       = &quot;Wavelet Power&quot;</div><div>  plot = gsn_csm_contour(wks,power({0:120},:),res)</div><div>end</div></div><div><br></div><div><br></div><div>I&#39;ll appreciate any help.</div><div><br></div><div>Many thanks,</div><div><br></div><div><br></div></div>