!---(1)      make a sweep trace ...
        !---(2)      ... and get its autocorrelation
        !---(3)      list-plot them in the time     -domain
        !---(4)      list-plot them in the frequency-domain
        !
        !---(1)
        call = sweep                            output   =  sweep
               f1          =  5
               f2          =  20
               t_max       =  4.0
               ramp_on     =  0.6
               ramp_off    =  1.2
               ramp_type   =  hanning
        !---(2)
        call = ccorf2                           input(1) =  sweep
               lags        = -25  25            input(2) =  sweep
                                                output   =  klauder
        !---(3)
        call = list2                            input(1) =  sweep
                                                input(2) =  klauder
        !---(4)
        call = listspec                         input    =  sweep
        call = listspec                         input    =  klauder

        !---(1)      read some segd data
        !---(2)      create/modify some trace headers
        !---(3)      write as segy format
        !
        !---(1)
        call = read_sgd                         input    =  e:\rmt0\segd.tap
               accept if (chan_set)  =  2
        !---(2)
        call = headit

               if   (channel) = 1    120
               then (offset)  = 2480 100

               if   (file)    = 101, 200  202, 273
               then (file)    = 101, 200  201, 272
        !---(3)
        call = write_sy                         output   =  segy.dsk

               ebcdic(1)  =  C_Line_1234
               ebcdic(2)  =  C
               ebcdic(3)  =  C_Re-format_from_Segd
               ebcdic(4)  =  C_blah_blah_blah

               header( 9,12) = file
               header(13,16) = channel
               header(37,40) = offset

Home