[cig-commits] r18187 - in seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/seis_process: . sac2asc

[danielpeter at geodynamics.org]

Author: danielpeter
Date: 2011-04-06 12:45:34 -0700 (Wed, 06 Apr 2011)
New Revision: 18187

Added:
   seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/seis_process/sac2asc/
   seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/seis_process/sac2asc/Makefile
   seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/seis_process/sac2asc/sac2asc.c
Log:
adds sac2asc processing utility in UTILS/seis_process/sac2asc/ directory

Added: seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/seis_process/sac2asc/Makefile
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/seis_process/sac2asc/Makefile	                        (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/seis_process/sac2asc/Makefile	2011-04-06 19:45:34 UTC (rev 18187)
@@ -0,0 +1,26 @@
+#########################################################
+# Makefile
+#
+# needs SAC libraries installed 
+########################################################
+
+##
+## modify to include your SAC library path
+##
+INCLUDE_SACPATH=-I/opt/seismo-util/lib/Sacio/single
+LIBRAY_SACPATH=-L/opt/seismo-util/lib
+
+CFLAGS = -g -Wall $(INCLUDE_SACPATH)
+
+PROG=sac2asc
+
+OBJS=$(PROG).o
+
+LIBS=$(LIBRARY_SACPATH) -lSacio -lSacTools
+
+
+$(PROG): $(OBJS)
+	$(CC) $(CFLAGS) -o $(PROG) $(OBJS) $(LIBS)
+
+clean:
+	rm -f $(PROG) $(OBJS) core *~ 

Added: seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/seis_process/sac2asc/sac2asc.c
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/seis_process/sac2asc/sac2asc.c	                        (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/seis_process/sac2asc/sac2asc.c	2011-04-06 19:45:34 UTC (rev 18187)
@@ -0,0 +1,31 @@
+
+#include <stdlib.h>
+#include <stdio.h>
+#include <sacio.h>
+
+int
+main(int argc, char *argv[]) {
+    SACHEAD sachead;
+    float *data;
+    char* filename;    
+    int i;
+    float time;
+
+    if(argc < 2) {
+	fprintf(stderr, "sac2asc: sacfile\n");
+        exit(-1);
+    }
+    filename = argv[1];
+    if((data = read_sac(filename, &sachead) ) == 0) {
+	fprintf(stderr, "Error reading sacfile: %s\n", filename);
+        exit(-1);
+    }
+    time = sachead.b;
+    for(i = 0; i < sachead.npts; i++) {
+        fprintf(stdout, "%f %e\n", time, data[i]);
+        time = time + sachead.delta;
+    }
+    free(data);
+
+    return(0);
+}