Hi
You are using "SET_UP per" as if it's a variable name. The error message is telling you that variable names in Ferret cannot have a space in them.
Do you mean to say this?
!Save fft in period unitlet A = per[l=1:`FFT_nf`], fft_uwnd
On 8/2/2012 10:58 AM, antonio avio wrote:
Hallo,I wonder why I got this error when i wanted to save the results of FFT analysis (period unit not frequency) in different nc files:
!Save fft in period unitlet A = SET_UP per[l=1:`FFT_nf`], fft_uwnd!-> DEFINE VARIABLE A = SET_UP per[l=1:378], fft_uwnd**ERROR: command syntax: SET_UP per[l=1:378], fft_uwndSET_UP before per[l=1:378] is illegalDEFINE VARIABLE A = SET_UP per[l=1:378], fft_uwndCommand file, command group, or REPEAT execution aborted
Here is my command:
[SNIP]!converting to Period! Get the frequency increment used in the FFT.LET FFT_nf = `fft_uwnd,return=lend`LET FFT_nyquist = 0.5LET FFT_freq1 = FFT_nyquist/ FFT_nf
! Define a frequency axis.DEFINE AXIS/T=`FFT_freq1`:`FFT_nyquist`:`FFT_freq1` FAXISDEFINE GRID/T=FAXIS gfftfreqLET a = T[g=gfftfreq]
! Define the period from the frequency axis.LET per = 1./a
! Plot period vs FFT Amplitudes showing the first 24 months where most of the energy is.! The PPL commands clean up the appearance of the plot.
SET VIEW urPLOT/VS/LINE/HLIMITS=0:32:2/TITLE="Amplitude Spectrum"/SET_UP per[l=1:`FFT_nf`], fft_uwndPPL XFOR (I2)PPL XLAB Period, months/cyclePPL YLABPPL PLOT!****************************************************************************************************! Compute and plot the phase using fftp.
LET FFT_uwndfftp = fftp(fft_u[l=1:757])SET VARIABLE/TITLE="FFT Phase"/UNITS="deg" FFT_uwndfftp
SET VIEW llPLOT FFT_uwndfftp
SET VIEW lrPLOT/VS/LINE/HLIMITS=0:32:2/TITLE="FFT Phase"/SET_UP per[l=1:`FFT_nf`],FFT_uwndfftpPPL XFOR (I2)PPL XLAB Period, months/cyclePPL YLAB DegPPL PLOT
!Save fft period unitlet A = SET_UP per[l=1:`FFT_nf`], fft_uwndsave/file=fft_period.nc A
!Save phase in period unitlet B = SET_UP per[l=1:`FFT_nf`], FFT_uwndfftpsave/file=phase_period.nc B
!**********************************************************************************Thanks for your help.
Antonio.
Attachment:
test.jnl
Description: Binary data