Re: [ferret_users] spectrum and FFTA error

On Mon, Nov 11, 2013 at 12:13 PM, Akshay Hegde <akshay.k.hegde@xxxxxxxxx> wrote:

look at the difference here... I suspect there might be some problem with your GMT Environment Settings (path/bash profile etc) if its not working even with default ferret datasets.

Here spectrum1d is working

yes? spawn spectrum1d
spectrum1d 4.5.6 - Compute auto- [and cross- ] spectra from one [or two] timeseries

usage: spectrum1d -S<segment_size> [-C[<xycnpago>]] [-D<dt>] [-H[i][<nrec>]] [-N<name_stem>]
[-V] [-W] [-b[i|o][s|S|d|D[<ncol>]|c[<var1>/...]]] [-f[i|o]<colinfo>]

Here I just tried some test command "spectrum1d_test_command" since it does not exist so getting following error
yes? spawn spectrum1d_test_command
sh: 1: spectrum1d_test_command: not found
Thanks and Regards,
Akshay Hegde 

CSIR-National Institute of Oceanography,
Dona Paula, Panjim, Goa -403004
On Mon, Nov 11, 2013 at 12:02 PM, Akshay Hegde <akshay.k.hegde@xxxxxxxxx> wrote:
Did you try with default dataset which comes with ferret ? whether following commands are working ?

    yes? SET DATA monthly_navy_winds
    yes? GO gmt_power uwnd[x=180,y=0] uwnd x180y0 64
    yes? PLOT x180y0_power
    yes? GO overlay_bars x180y0_power x180y0_error red
Thanks and Regards,
Akshay Hegde 
CSIR-National Institute of Oceanography,
Dona Paula, Panjim, Goa -403004
On Mon, Nov 11, 2013 at 11:55 AM, chinnu sachi <chinnu.sachi7@xxxxxxxxx> wrote:
My time axis is ok,I list it below. even I removed clobber ,still same error is coming.
I can open spectrum1d with man spectrum1d

         VARIABLE : L
                        axis TIME
             FILENAME : v_test.nc
             SUBSET   : 600 points (TIME)
15-JAN-1958 00 /   1:    1.0
15-FEB-1958 00 /   2:    2.0
15-MAR-1958 00 /   3:    3.0
15-APR-1958 00 /   4:    4.0
15-MAY-1958 00 /   5:    5.0
15-JUN-1958 00 /   6:    6.0
15-JUL-1958 00 /   7:    7.0
15-AUG-1958 00 /   8:    8.0
15-SEP-1958 00 /   9:    9.0
15-OCT-1958 00 / 10:   10.0
15-NOV-1958 00 / 11:   11.0
15-DEC-1958 00 / 12:   12.0
15-JAN-1959 00 / 13:   13.0
15-FEB-1959 00 / 14:   14.0
15-MAR-1959 00 / 15:   15.0
15-APR-1959 00 / 16:   16.0
15-MAY-1959 00 / 17:   17.0
15-JUN-1959 00 / 18:   18.0
15-JUL-1959 00 / 19:   19.0
15-AUG-1959 00 / 20:   20.0
15-SEP-1959 00 / 21:   21.0
15-OCT-1959 00 / 22:   22.0
15-NOV-1959 00 / 23:   23.0
15-DEC-1959 00 / 24:   24.0
On Mon, Nov 11, 2013 at 11:48 AM, Akshay Hegde <akshay.k.hegde@xxxxxxxxx> wrote:
Sorry I forgot to mention make sure that you can open man page using following command since you are getting this error "sh: 1: spectrum1d: not found"

man spectrum1d
Thanks and Regards,
Akshay Hegde 
CSIR-National Institute of Oceanography,
Dona Paula, Panjim, Goa -403004
On Mon, Nov 11, 2013 at 11:40 AM, Akshay Hegde <akshay.k.hegde@xxxxxxxxx> wrote:
Try something like this
GO gmt_power variable[x=180,y=0] grid output seg_len
and in your script you used both clobber and append
append ---> If file exists then append, otherwise create and write to it

clobber ---> If the file exists then overwrite it, otherwise create and write to it
in your case you need to use only append
and Hows your Time axis ?
list t[gt=variable]
Thanks and Regards,
Akshay Hegde 
CSIR-National Institute of Oceanography,






Dona Paula, Panjim, Goa -403004
On Mon, Nov 11, 2013 at 10:32 AM, chinnu sachi <chinnu.sachi7@xxxxxxxxx> wrote:

Dear ferret users,

                   I am getting an error error message while trying to run the gmt_power function. The details of the error for power spectrum as follows:

yes? go power.jnl
set mem/size = 650
Cached data cleared from memory
use v_test.nc
repeat/k=1:`20` ( go gmt_power V_NEW GHA1 gmt 256; list/file=test/clobber/append gmt_power;plot gmt_power)
!-> repeat/k=1:20 ( go gmt_power V_NEW GHA1 gmt 256; list/file=test/clobber/append gmt_power;plot gmt_power)
!-> REPEAT: K=1
LISTing to file gmt.tseries
>>> Using GMT command: spectrum1d gmt.tseries -Ngmt -S256 -D30.43583
sh: 1: spectrum1d: not found
**TMAP ERR: non-existent or not on line
             gmt.xpower
SET DATA/EZ/var=freq gmt.xpower
Command file, command group, or REPEAT execution aborted

And also when using FFTA command, error is coming like "Bailing out of external function "ffta":
    Time axis must be a regular axis
**ERROR: error in external function" .But when gone through the mail archives i found a link stating that it happens when there is missing values. But in my data there is no missing value, in that region.

Thanks in advance

--
Regards,
Chinnu Sachidanandan
Physical Oceanography Division
National Institute of Oceanography
csachidanandan@xxxxxxx
--
Regards,
Chinnu Sachidanandan
Physical Oceanography Division
National Institute of Oceanography
csachidanandan@xxxxxxx

--
Regards,
Chinnu Sachidanandan
Physical Oceanography Division
National Institute of Oceanography
csachidanandan@xxxxxxx