[Thread Prev][Thread Next][Index]
[ferret_users] how to calculate confidence level in fft
- To: ferret_users@xxxxxxxx
- Subject: [ferret_users] how to calculate confidence level in fft
- From: jaya kumar A <jayamet@xxxxxxxxx>
- Date: Thu, 01 Oct 2009 13:50:34 +0530
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:mime-version:received:date:message-id:subject :from:to:content-type; bh=kEUzl14DPRfacHFPhDXXDQLlxYYCG2rZhWXLmSMIiqE=; b=YRFuDY7Bdysf9spKBfDdZ+K0RE+iGdLiaWccFTfjy2cZYFdeg1rq32G9wK+KcvCUuo CuN/FsvN9qwVKda7Pk265kv79NzfYffkGAOJupLTiXGRb3gdk3FzMT1YLEqbWBiTbn9+ sR8lHJabCseWI9WIZlZUeAziGoJjAdK0gRyVg=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:date:message-id:subject:from:to:content-type; b=PP7L5rYZpE5WKIR8q22JE69EZVn8uknlQWhWrywwf0c4brxiW54sIqlBcluCweOkIZ 3avhVJxu5ZILnwLIgyhrx0YGubXOUrea4ayU4xNOqF6att6exg1lwPcApV118vnZPe0H DJmN9cllK3A5clEOPUdT0FKryfrxQLqYX9z8s=
- Sender: owner-ferret_users@xxxxxxxx
dear ferret users,
1. I want to find the peak signals in the JJAS time, how many record in input file to get the significant peak of 90 days in the amplitude spectrum? I use the code below, what it means FFT_nyquist?
2. After evaluating FFT and calculate confidence Level need to check the statistical significance of the peaks, whether any option in ferret
3. If one year spectrum is below confidence how to make the ensemble of the JJAS amplitude spectrum for 10 years.
thanks in advance for the help
with regards
--
JAYAKUMAR A
Senior Research Fellow,
Theoretical Studies Division
IITM, Pune-8
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
USE "2008_mjjaso.nc"
SHOW DATA
LET FFT_etfn = sst[y=-3.6:-5.2@ave,x=60:90@ave]
LET FFT_etfft = FFTA(FFT_etfn[l=1:183])
SET VARIABLE FFT_etfft
! Get the frequency increment used in the FFT.
LET FFT_nf = `FFT_etfft,return=lend`
LET FFT_nyquist = 0.5
LET FFT_freq1 = FFT_nyquist/ FFT_nf
! Define a frequency axis.
DEFINE AXIS/T=`FFT_freq1`:`FFT_nyquist`:`FFT_freq1` faxis
DEFINE GRID/T=faxis gfftfreq
LET a = T[g=gfftfreq]
! Define the period from the frequency axis.
LET per = 1./a
PLOT/VS/nolab/thick=3/LINE=7/HLIMITS=0:90:10/SET_UP/vlimits=0:0.5:0.05 per[l=1:`FFT_nf`], FFT_etfft
!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Privacy Policy | Disclaimer | Accessibility Statement