Go to SIOSEIS introduction.          Back to Healy06

Script to make envelopes from Knudsen correlates

#!/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-sorted