Re: Locate good data point

[Mick Spillane <spillane@xxxxxxxxxxxxx>]

Hi Jens-Olaf,

Try this demo

! define a 3-d variable
def axis/x=1:10:1 xax ; def axis/y=1:10:1 yax ; def axis/depth/z=1:5:1 zax
def grid/x=xax/y=yax/z=zax grd ; let v=x[g=grd]+y[g=grd]-z[g=grd]

! and suppose only values greater than 10 are good
let vg=if(v gt 10)then v

! look at the first layer or two to see what its like
shade/k=1 vg
...

! define a mask function for "good" values ...
let mask=if(vg)then 1

! ... and an x-edge function that locates the first good point
let xedge=mask[x=@loc:1]

! see if it is what you want
shade/k=1 vg ; list/k=1 xedge

! you could also define a y-edge function
let yedge=mask[x=@loc:1]
list/k=1 yedge

Hope this helps,
Mick


|--****--****-*---*---***--***--|____spillane@pmel.noaa.gov____|
|-*__---*-----*--*-*--*--*-*--*-|_SCIENCE APPLICATIONS SUPPORT_|
|--***--*-----*-*---*-***--***--|____EPIC/Ferret/PlotPlus______|
|-----*-*-----*-*****-*----*----|__Room 2070 Bldg#3 NOAA/PMEL__|
|-****---****-*-*---*-*----*----|____Phone_:_(206)526-6780_____|