C-Shell script to plot indicidual moved-out cmp gathers

#! /bin/csh -f
setenv DISPLAY walrus:0.0    # CHANGE THIS FOR YOUR MACHINE!
if( $#argv < 2 ) then
    echo "Usage: chk.nmo filename start_num [ num_inc END_num]"
    exit 1
endif

set FILE = $1
set START_NUM = $2
if( $#argv < 3 ) then
    set NUM_INC = 1
else
    set NUM_INC = $3
endif

if( $#argv < 4 ) then
    set END_NUM = $START_NUM
else
    set END_NUM = $4
endif

echo $START_NUM $END_NUM $NUM_INC
set num = $START_NUM
while ( $num <= $END_NUM )
/home/henkart/src/sioseis/sioseis << eof
procs diskin nmo mute filter agc plot END
diskin
   ipath $FILE fno $num lno $num allno no END
END
mute
  addwb yes xtp 2000 -.02 2500 .6 4000 .7 end
end
nmo
  ! stretc .6
  fno 306701 vtp 1510 3.481 1500 3.504 1520 3.545 1519 3.661 1540 4.020
1538 4.292 1774 5.984 end
  fno 306901 vtp 1510 3.464 1520 3.608 1540 4.031 1543 4.269 1780 5.845
end
  fno 307101 vtp 1510 3.475 1510 3.614 1538 4.020 1521 4.095 1768 5.967
end
  fno 307301 vtp 1509 3.481 1510 3.591 1530 3.852 1539 4.020 1539 4.153
1546 4.431 1747 5.949 end
  fno 307501 vtp 1510 3.493 1519 3.684 1539 4.055 1542 4.136 1774 5.978
end
  fno 307701 vtp 1509 3.504 1510 3.672 1540 4.055 1773 5.961 end
  fno 307901 vtp 1510 3.504 1509 3.585 1509 3.695 1541 4.055 1569 4.472
1758 6.262 end
  fno 308101 vtp 1510 3.533 1510 3.666 1519 3.811 1540 4.084 1569 4.460
1771 5.978 end
  fno 308301 vtp 1510 3.533 1510 3.585 1519 3.811 1540 4.101 1570 4.443
1762 5.978 end
  fno 308501 vtp 1510 3.562 1519 3.805 1530 3.962 1540 4.101 1570 4.443
1759 5.984 end
  ETC.
end
end
filter
   ftype 0 dbdrop 48 pass 20 220 end
end
agc
   winlen .5 end
end
plot
   nsecs 2 nibs 75 def .06 srpath sunfil.pch ann range taginc 2
   trpin 10 vscale 6.6666 END
END
END
eof
xloadimage -r 90 sunfil.pch
@ num = $num + $NUM_INC
end


Return to EW0008 example.               Go to SIOSEIS examples.

Go to the list of seismic processes.    Go to SIOSEIS introduction.