[Thread Prev][Thread Next][Index]

[ferret_users] Problems using external function SCAT2GRIDGAUSS_ZT



Hi all,

I'm getting an error, trying to use the external function SCAT2GRIDGAUSS_ZT in ferret v6.62.

I have some chlorophyll profiles from HOT as ASCII and want to put them on a depth-time grid for further processing and plotting. My ferret script so far:

SET DATA/SKIP=2/FORMAT=DELIMITED/TYPE=" numeric , text, text, numeric, numeric"/VAR=nr,date,time,dbar,chla hplc_chla_noinvalid.txt

let mm=STRFLOAT(SUBSTRING(date[d=1],0,2))
let dd=STRFLOAT(SUBSTRING(date[d=1],3,2))
let yy=STRFLOAT(SUBSTRING(date,5,2))

let year=IF yy GT 50 THEN yy + 1900 ELSE yy + 2000

let date_new= DAYS1900(year, mm, dd) - DAYS1900(1988, 01, 01)



Overview:

yes? sho dat
     currently SET data sets:
    1> ./hplc_chla_noinvalid.txt  (default)
 name     title                             I         J         K         L
 NR       NR                               1:20480   ...       ...       ...
 DATE     DATE                             1:20480   ...       ...       ...
 TIME     TIME                             1:20480   ...       ...       ...
 DBAR     DBAR                             1:20480   ...       ...       ...
 CHLA     CHLA                             1:20480   ...       ...       ...

yes? list date_new[i=1:5]
             VARIABLE : DAYS1900(YEAR, MM, DD) - DAYS1900(1988, 01, 01)
             FILENAME : hplc_chla_noinvalid.txt
             SUBSET   : 5 points (X)
 1   / 1:  337.0
 2   / 2:  337.0
 3   / 3:  337.0
 4   / 4:  337.0
 5   / 5:  337.0

yes? list dbar[i=1:5]
             VARIABLE : DBAR
             FILENAME : hplc_chla_noinvalid.txt
             SUBSET   : 5 points (X)
 1   / 1:   2.10
 2   / 2:  22.80
 3   / 3:  48.60
 4   / 4:  73.00
 5   / 5:  98.20

yes? list chla[i=1:5]
             VARIABLE : CHLA
             FILENAME : hplc_chla_noinvalid.txt
             SUBSET   : 5 points (X)
 1   / 1:   23.0
 2   / 2:   34.0
 3   / 3:   55.0
 4   / 4:   56.0
 5   / 5:  135.0


Trying to use SCAT2GRIDGAUSS_XY works:

define axis/x=0:300:1/units=meters xax
define axis/y=300:400:1 yax
let chla_reg = scat2gridgauss_xy(dbar[i=1:5],date_new[i=1:5],chla[i=1:5],x[gx=xax],y[gy=yax],1.,1.,1.,0)
shade chla_reg

SCAT2GRIDGAUSS_ZT doesn't:

define axis/t=300:400:1 tax
define axis/z=0:300:1 zax
let chla_reg = scat2gridgauss_zt(dbar[i=1:5],date_new[i=1:5],chla[i=1:5],z[gz=zax],t[gt=tax],1.,1.,1.,0)
shade chla_reg

Bailing out of external function "scat2gridgauss_zt":
         output axes must be REGULARLY spaced
 **ERROR: error in external function

But I don't know why. The output axes are regularly...

sho ax/all
 name       axis              # pts   start                end
 PSXT      LONGITUDE          160mr   130.5E               70.5W
 PSXU      LONGITUDE          160mr   131E                 70W
 PSYT      LATITUDE           100 i   28.836S              48.567N
 PSYU      LATITUDE           100 i   27.672S              50N
 PSZT      DEPTH (m)           27 i-  5                    3824
 PSZW      DEPTH (m)           27 i-  10                   4149
 ABSTRACT                 9999999 r   1                    10000000
 EZ                         20480 r   1                    20480
 ZPRESS    DEPTH (m)         2284 i-  1.9                  257.4
 YAX       Y                  101 r   300                  400
 XAX       X (METERS)         301 r   0                    300
 TAX       T                  101 r   300                  400
 ZAX       Z                  301 r   0                    300

It doesn't matter how I define the t-axis (using t="startdate":"enddate":step ; using to="01-JAN-1988" ; using different calendar=... also doesn't help).

Any suggestions appreciated.

Thank you,
Bastian Kern

[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce / NOAA / OAR / PMEL / Ferret

Privacy Policy | Disclaimer | Accessibility Statement