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