Go to SIOSEIS introduction. Back to Healy06
#!/bin/csh -f if( $#argv != 1 ) then echo "Usage: mkenv filename (without .sgy suffix)" exit 1 endif set FILE = $1 sioseis << eof procs diskin header t2f f2t gains header2 prout diskoa end diskin ipath /Volumes/Data/Raw/knudsenraw/$FILE.sgy end end header i120 = i58 ! save the original trace length fno 0 lno 9999999 ftr 0 ltr 999 end end header2 i58 = i120 ! restore the original trace length fno 0 lno 9999999 ftr 0 ltr 999 end end prout info 1 ! print plotting info fno 0 lno 999999 noinc 50 end ! print every 50th trace end gains type 7 end ! complex modulus - make envelope from analytic end t2f end ! number of sample is the next power of two larger than the input end f2t type analytic end ! create the complex trace end diskoa opath env-$FILE.sgy end end end eof #chmod 444 /Users/seismic/Data/Chirp/envelope/env-$FILE.sgy #lsh /Users/seismic/Data/Chirp/envelope/env-$FILE.sgy >> info #sort +2 -5 info > info-sortedBack to SIOSEIS Examples. Go to the list of seismic processes. Go to SIOSEIS introduction. Back to Healy06