sac-101.4 for linux

SAC(Seismic Analysis Code)是美国加州大学Lawrence Livermore国家实验室(LLNL)开发研制的，用于研究连续信号，特别是时间序列数据的通用程序

#include 
#include 
#include 

#include 
#include 

#define MAX        4000
#define ERROR_MAX  256

int 
main（int argc char *argv[]） {

    /* Local variables */
    int i j;
    int nlen nlen1 nlen2 nerr max;

    float beg delta end;
    char *kname;

    float yarray1[MAX] yarray2[MAX] ytmp[MAX] xarray[1];
    float *out;

    int nwin wlen nfft leven;

    char error[ERROR_MAX];

    max = MAX;

    for（i = 0; i < MAX; i++） {
      yarray1[i] = 0.0;
      yarray2[i] = 0.0;
      ytmp[i] = 0.0;
    }
    /* Read in the first file  */        
    kname = strdup（“convolvec_in1.sac“）;
    rsac1（kname ytmp &nlen1 &beg &delta &max &nerr SAC_STRING_LENGTH）;

    if （nerr != 0） {
      fprintf（stderr “Error reading in file（%d）: %s\n“ nerr kname）;
      exit（-1）;
    }


    /* Read in the second file  */
    kname = strdup（“convolvec_in2.sac“）;
    rsac1（kname yarray2 &nlen2 &beg &delta &max &nerr SAC_STRING_LENGTH）;

    if （nerr != 0） {
      fprintf（stderr “Error reading in file: %s\n“ kname）;
      exit（-1）;
    }

    /* Reverse the First Signal */
    j = 0;
    for（i = nlen1 - 1; i >= 0; i--） {
      yarray1[j] = ytmp[i];
      j++;
    }
    
    nlen = nlen1;
    if（nlen2 > nlen） {
      nlen = nlen2;
    }
    /* Allocate space for the correlation of yarray1 and yarray2 */
    max = next2（（2 * nlen） - 1） * 2;
    out = （float *） malloc（sizeof（float） * max）;
    if（out == NULL） {
      fprintf（stderr “Error allocating memory for correlation\n“）;
      exit（-1）;
    }

    /* Set up values for the cross correlation */
    nwin = 1;
    wlen = nlen;
    nfft = 0;
    
    /*     Call crscor （ Cross Correlation no wait uh Convolution ）
     *        - yarray1 - First  Input array to correlate
     *        - yarray2 - Second Input array to correlate
     *        - nlen    - Number of points in yarray and yarray2
     *        - nwin    - Windows to use in the correlation
     *        - wlen    - Length of the windows
     *        - type    - Type of Window （SAC_RECTANGLE）
     *        - out     - output sequence 
     *        - nfft    - Length of the output sequence
     *        - error   - Error Message
     *        - err_len - Length of Error Message （on input）
     */
    crscor（yarray1 yarray2 nlen 
           nwin wlen SAC_RECTANGLE
           out &nfft error ERROR_MAX）;

    /* Zero out the tmp signal */
    for（i = 0; i < MAX; i++） {
      ytmp[i] = 0.0;
    }
    
    /* Reconstruct the signal from the “cross correlation“ back to front
     *  
     *  ytmp[0         : nlen1 - 2         ] <- out[nfft-nlen1+1 : nfft  - 1 ] 
     *  ytmp[nlen1 - 1 : nlen1 + nlen2  -2 ] <- out[0            : nlen2 - 1 ]
     *
     *  nfft-1 is the last point of the output sequence
     */
    for（i = 0; i <= nlen1 - 2; i++） {  
      ytmp[i] = out[nfft - nlen1 + i + 1];
    }
    for（i = 0; i <= nlen2 - 1; i++） {
      ytmp[nlen1 + i - 1] = out[i];
    }