Re: [ferret_users] spectrum and FFTA error

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