; proj3d file axis algo
;
; axis is F1, F2 or F3
; algo is M or S
;
; equivalent to proj in 3D, but on joined file
;
; see also : PROJ JOIN proc3d

message "Enter file name"
set file = $_

join $file
if ($c_dim != 3) error "available on 3D dataset only"

message "enter axis (F1,F2,F3)"
set axis = $_
set axis = (toupper($axis))

dim 2
if ($axis s= "F1") then
   chsize $c_sizef2 $c_sizef3
   set siz = $c_sizef1
elsif ($axis s= "F2") then
   chsize $c_sizef1 $c_sizef3
   set siz = $c_sizef2
elsif ($axis s= "F3") then
   chsize $c_sizef1 $c_sizef2
   set siz = $c_sizef3
else
   error "Wrong axis"
endif

message "enter algorithm, Skyline or Mean : (S/M)"
set alg = $_
set alg = (toupper($alg))
if ($alg s= 'S') then
   set todo = MAXDATA
elsif ($alg s= 'M') then
   set todo = ADDDATA
else
   error "wrong algorithm"
endif

initinprogress $siz
getc $axis 1 %%  put data
for i = 2 to $siz
    inprogress $i
    getc $axis $i %%
    $todo  put data
endfor

get data   mult (1/$siz)

