Go to SIOSEIS introduction.          Back to Healy06

Script to make a plot from Knudsen envelopes

#!/bin/csh -f
rm sunfil.ras
sioseis << eof
procs diskin prout geom stack wbt avenor filter mix gains plot end
diskin
   fgmt 0300 lgmt 320
    ipath env-2006_226_0158_LF_028.sgy
   end
end
geom
   type 17 dbrps 3 end
end
gains
    subwb yes type 5 alpha 5 end
end
avenor
   sets 0 .1 addwb yes end
end
wbt
   vel 1475 end
end
mix
   weight 1  1 end
end
filter
    ftype 0 pass 2 500 dbdrop 48 end
end
prout
    fno 0 lno 9999999 noinc 100 end
end
plot
   stime 1 dir ltr nsecs .2 dptr 1 tlines .05 .1 
    colors gray opath siofil wiggle 0  ann gmtint anninc 5 ann2 range
    trpin 300 def .01 tlines .05 nibs 2859  vscale 5 end
end
end
eof
sio2sun siofil sunfil.ras
#convert -rotate 90 sunfil.ras $FILE.png
display -rotate 90 sunfil.ras &