Re: Spectral Analysis

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Mon Nov 19 2012 - 13:03:38 MST

The following could be used as a template.

   sdof = specx_anal(soi,d,sm,pct)

   printVarSummary(sdof)
   print(sdof@frq+" "+(1.0/sdof@frq)+" "+sdof@spcx)
   nfrq = dimsizes(sdof@frq)

----------------
Use info above to select the index values
of the frequencies to be plotted

    period = 1.0/sdof@frq

    r@tmXBMode = "Explicit"
    r@tmXBValues = (/0.1, 0.2, 0.3, 0.4, 0.5/)
   ;r@tmXBLabels = (/0.1+"~C~"+sprintf("%3.1f",period(167)) \
   ; ,0.2+"~C~"+sprintf("%3.1f",period(333)) \
   ; ,0.3+"~C~"+sprintf("%3.1f",period(501)) \
   ; ,0.4+"~C~"+sprintf("%3.1f",period(668)) \
   ; ,0.5+"~C~"+sprintf("%3.1f",period(nfrq-1)) /)

    r@tmXBLabels = (/sprintf("%3.1f",period(167))+"~C~"+0.1 \
                             ,sprintf("%3.1f",period(333))+"~C~"+0.2 \
                             ,sprintf("%3.1f",period(501))+"~C~"+0.3 \
                             ,sprintf("%3.1f",period(668))+"~C~"+0.4 \
                             ,sprintf("%3.1f",period(nfrq-1))+"~C~"+0.5 /)

    plot = gsn_csm_xy(wks,sdof@frq, splt,r)

On 11/19/2012 10:31 AM, Renata Tedeschi wrote:
> Hi everyone.
>
> I use the spec_4.ncl (available on the website of ncl), but I need to do
> two x-axis: one (the inferior) is the same that have in spec_4.ncl
> (frequency), and other (the superior) its values should be the period
> (inverse of frequency). How can I do this?
>
> Thanks
> Renata Tedeschi
>
>
> _______________________________________________
> 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 Mon Nov 19 13:03:47 2012

This archive was generated by hypermail 2.1.8 : Wed Nov 21 2012 - 11:16:05 MST