NAME

     mkphase - estimate phase correction for seismic data

SYNOPSIS

     mkphase [ -Nntap ] [ -Sofile ] [ -rsirs ] [ -reire ] [ -siws
     ] [ -eiwe ] [ -? ] [ -h ]

DESCRIPTION

     mkphase estimates the constant phase  rotation  required  to
     maximize the kurtosis, or sharpness, of seismic data.  It is
     assumed that  this  phase  is  the  correction  required  to
     correct  the  data  to  zero  phase.   The  theory  and data
     requirements are described in Maximum kurtosis phase correc-
     tion,  R.E.  White,  Geophysical Journal (1988) 95, 371-389.
     The output from the program  consists  of  a  trace-by-trace
     breakdown  of the kurtosis components in the standard system
     output file, and a file for input to XGRAPH, describing  the
     phase  variance of kurtosis and the estimated wavelet phase.
     The phase of the wavelet is reported in the  legend  of  the
     xgraph  file.   The phase correction is the negative of this
     and corresponds to the peak in the graph.

     mkphase gets both its data and its parameters  from  command
     line  arguments.  These arguments specify the input and out-
     put data set names,  the  beginning  and  ending  sequential
     record  numbers for processing, and the beginning and ending
     times for the estimation window.

     -Nntap
          Enter the input data set name or file immediately after
          typing  -N;  entering  nothing  specifies a stdin pipe.
          This input file should include the complete  path  name
          if  the  file resides in a directory other than the one
          from   which   you   start   the   program.    Example:
          -N/work_area/dummy  tells  the program to look for file
          "dummy" in directory "work_area".

     -Sofile
          Enter the name of the output file containing the  phase
          estimation  and  phase variance of the kurtosis immedi-
          ately after typing -S. The output file name is required
          and  may require a fully qualified path name. The phase
          of the wavelet is annotated in the legend of the XGRAPH
          output.  No  output  file  sends  the  output stream to
          stdout and can be piped into xgraph

     -rs irs
          (integer) First (sequential) record  to  process.   All
          records prior to this one are skipped. Default = 1 (the
          first).


     -re ire
          (integer) Last (sequential)  record  to  process.   All
          records after this are ignored. Default = process all.

     -s   Enter the time, in  ms,  for  estimation  window  start
          time. Default is 0 (first sample).

     -e   Enter the time, in ms, for the  estimation  window  end
          time.  Default is maximum time (last sample).

     -?   Enter this (or '-?' if you are running in  c-shell)  to
          get  online  help.   The  program  will terminate after
          printing this.

     -h   Enter this to get online help.  The program  will  ter-
          minate after printing the help.

BUGS

     Unknown

AUTHOR

     R.E. White, ported to USP by Richard Crider, Sue Raikes, and
     Steve Lancaster.

COPYRIGHT

     copyright 2001, Amoco Production Company
               All Rights Reserved
          an affiliate of BP America Inc.