[ncl-talk] Plot wavelet with "Months" as X-axis
Dennis Shea
shea at ucar.edu
Fri Jan 13 10:28:43 MST 2017
That should be:
res at tmXBMode = "Explicit"
res at tmXBValues = (/3000,3500,4000,.....,8000/) ; whatever
res at tmXBLabels = (/"Jan","Feb",.....,"Dec"/)
On Fri, Jan 13, 2017 at 9:12 AM, Dennis Shea <shea at ucar.edu> wrote:
> You will have to put in the correct numbers. ncl-talk has no idea to what the 'time' axis refers refers
>
>
> res at tmYLMode = "Explicit"
> res at tmYLValues = (/3000,3500,4000,.....,8000/)
> res at tmYLLabels = (/"Jan","Feb",.....,"Dec"/)
>
>
> On Thu, Jan 12, 2017 at 11:16 PM, Lyndon Mark Olaguera <
> olagueralyndonmark429 at gmail.com> wrote:
>
>> Hi All,
>>
>> I'm plotting a wavelet based from this example https://www.ncl.ucar.e
>> du/Document/Functions/Built-in/wavelet.shtml
>>
>> Is there an easy way to make the x-axis in terms of months (Jan to Dec)?
>>
>> Here's my code:
>>
>> begin
>> ;************************************
>> ; create pointer to file and read in variables
>> ;************************************
>> f = addfile ("hpass_SCS.nc","r")
>> prec = f->PREC(:,0,0) ; one value/year
>> time = prec&time
>> N = dimsizes(time)
>> ;************************************
>> ; Remove mean and linear trend (recommended, not required)
>> ; Missing values are not allowed.
>> ;************************************
>> prec = dtrend(prec,False)
>> ;************************************
>> ; compute wavelet (Missing values not allowed due to use of FFT)
>> ;************************************
>> mother = 0
>> param = 6.0
>> dt = 1.
>> s0 = 0.25
>> dj = 0.25
>> jtot = 44
>> npad = N
>> nadof = new(2,float)
>> noise = 1
>> siglvl = .05
>> isigtest= 1
>>
>> w = wavelet(prec,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
>>
>> power!1 = "time" ; X axis
>> power&time = time
>>
>> power at long_name = "Power Spectrum"
>> power at units = "unit^2"
>>
>>
>> ;************************************
>> ; initial resource settings
>> ;************************************
>> wks = gsn_open_wks("png","wavelet") ; send graphics to PNG
>> file
>> res = True ; plot mods desired
>> 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 ; turn off contour
>> lines
>> res at lbOrientation = "Vertical" ; vertical label bar
>> res at trYReverse = True ; reverse y-axis
>> ;res at tmLabelAutoStride = True
>> res at cnLevelSelectionMode = "ManualLevels" ; set manual contour
>> levels
>> res at cnMinLevelValF = 0.0 ; set min contour level
>> res at gsnLeftString = "Wavelet Power"
>> plot = gsn_csm_contour(wks,power({0:120},:),res)
>> end
>>
>>
>> I'll appreciate any help.
>>
>> Many thanks,
>>
>>
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170113/b23d51ec/attachment.html
More information about the ncl-talk
mailing list