sac.h

解吸SEED格式的源代码

/*===========================================================================*/
/* SEED reader     |               sac.h                   |     header file */
/*===========================================================================*/
/*
	Name:		sac.h
	Purpose:	structure for header of a SAC (Seismic Analysis Code) data file
	Usage:		#include "sac.h"
	Input:		not applicable
	Output:		not applicable
	Warnings:	not applicable
	Errors:		not applicable
	Called by:	output_data
	Calls to:	none
	Algorithm:	not applicable
	Notes:		Key to comment flags describing each field:
				Column 1:
					R = required by SAC
					  = (blank) optional
				Column 2:
					A = settable from a priori knowledge
					D = available in data
					F = available in or derivable from SEED fixed data header
					T = available in SEED header tables
					  = (blank) not directly available from SEED data, header
					    tables, or elsewhere
	Problems:	none known
	References:	O'Neill, D. (1987).  IRIS Interim Data Distribution Format
		(SAC ASCII), Version 1.0 (12 November 1987).  Incorporated
		Research Institutions for Seismology, 1616 North Fort Myer
		Drive, Suite 1440, Arlington, Virginia 22209.  11 pp.
		Tull, J. (1987).  SAC User's Manual, Version 10.2, October 7,
			1987.  Lawrence Livermore National Laboratory, L-205,
			Livermore, California 94550.  ??? pp.
	Language:	C, hopefully ANSI standard
	Author:	Dennis O'Neill
	Revisions:07/15/88  Dennis O'Neill  Initial preliminary release 0.9
		11/21/88  Dennis O'Neill  Production release 1.0
		09/19/89  Dennis O'Neill  corrected length of char strings
*/

struct sac
{
	float	delta;			/* RF time increment, sec    */
	float	depmin;			/*    minimum amplitude      */
	float	depmax;			/*    maximum amplitude      */
	float	scale;			/*    amplitude scale factor */
	float	odelta;			/*    observed time inc      */
	float	b;			/* RD initial value, ampl.   */
	float	e;			/* RD final value, amplitude */
	float	o;			/*    event start, sec > 0   */
	float	a;			/*    1st arrival time       */
	float	internal1;		/*    internal use           */
	float	t0;			/*    user-defined time pick */
	float	t1;			/*    user-defined time pick */
	float	t2;			/*    user-defined time pick */
	float	t3;			/*    user-defined time pick */
	float	t4;			/*    user-defined time pick */
	float	t5;			/*    user-defined time pick */
	float	t6;			/*    user-defined time pick */
	float	t7;			/*    user-defined time pick */
	float	t8;			/*    user-defined time pick */
	float	t9;			/*    user-defined time pick */
	float	f;			/*    event end, sec > 0     */
	float	resp0;			/*    instrument respnse parm*/
	float	resp1;			/*    instrument respnse parm*/
	float	resp2;			/*    instrument respnse parm*/
	float	resp3;			/*    instrument respnse parm*/
	float	resp4;			/*    instrument respnse parm*/
	float	resp5;			/*    instrument respnse parm*/
	float	resp6;			/*    instrument respnse parm*/
	float	resp7;			/*    instrument respnse parm*/
	float	resp8;			/*    instrument respnse parm*/
	float	resp9;			/*    instrument respnse parm*/
	float	stla;			/*  T station latititude     */
	float	stlo;			/*  T station longitude      */
	float	stel;			/*  T station elevation, m   */
	float	stdp;			/*  T station depth, m      */
	float	evla;			/*    event latitude         */
	float	evlo;			/*    event longitude        */
	float	evel;			/*    event elevation        */
	float	evdp;			/*    event depth            */
	float	unused1;		/*    reserved for future use*/
	float	user0;			/*    available to user      */
	float	user1;			/*    available to user      */
	float	user2;			/*    available to user      */
	float	user3;			/*    available to user      */
	float	user4;			/*    available to user      */
	float	user5;			/*    available to user      */
	float	user6;			/*    available to user      */
	float	user7;			/*    available to user      */
	float	user8;			/*    available to user      */
	float	user9;			/*    available to user      */
	float	dist;			/*    stn-event distance, km */
	float	az;			/*    event-stn azimuth      */
	float	baz;			/*    stn-event azimuth      */
	float	gcarc;		/*    stn-event dist, degrees*/
	float	internal2;	/*    internal use           */
	float	internal3;	/*    internal use           */
	float	depmen;		/*    mean value, amplitude  */
	float	cmpaz;		/*  T component azimuth     */
	float	cmpinc;		/*  T component inclination */
	float	unused2;	/*    reserved for future use*/
	float	unused3;	/*    reserved for future use*/
	float	unused4;	/*    reserved for future use*/
	float	unused5;	/*    reserved for future use*/
	float	unused6;	/*    reserved for future use*/
	float	unused7;	/*    reserved for future use*/
	float	unused8;	/*    reserved for future use*/
	float	unused9;	/*    reserved for future use*/
	float	unused10;	/*    reserved for future use*/
	float	unused11;	/*    reserved for future use*/
	float	unused12;	/*    reserved for future use*/
	long	nzyear;		/*  F zero time of file, yr  */
	long	nzjday;		/*  F zero time of file, day */
	long	nzhour;		/*  F zero time of file, hr  */
	long	nzmin;		/*  F zero time of file, min */
	long	nzsec;		/*  F zero time of file, sec */
	long	nzmsec;		/*  F zero time of file, msec*/
	long	internal4;	/*    internal use           */
	long	internal5;	/*    internal use           */
	long	internal6;	/*    internal use           */
	long	npts;		/* RF number of samples      */
	long	internal7;	/*    internal use           */
	long	internal8;	/*    internal use           */
	long	unused13;	/*    reserved for future use*/
	long	unused14;	/*    reserved for future use*/
	long	unused15;	/*    reserved for future use*/
	long	iftype;		/* RA type of file          */
	long	idep;		/*    type of amplitude      */
	long	iztype;		/*    zero time equivalence  */
	long	unused16;	/*    reserved for future use*/
	long	iinst;		/*    recording instrument   */
	long	istreg;		/*    stn geographic region  */
	long	ievreg;		/*    event geographic region*/
	long	ievtyp;		/*    event type             */
	long	iqual;		/*    quality of data        */
	long	isynth;		/*    synthetic data flag    */
	long	unused17;	/*    reserved for future use*/
	long	unused18;	/*    reserved for future use*/
	long	unused19;	/*    reserved for future use*/
	long	unused20;	/*    reserved for future use*/
	long	unused21;	/*    reserved for future use*/
	long	unused22;	/*    reserved for future use*/
	long	unused23;	/*    reserved for future use*/
	long	unused24;	/*    reserved for future use*/
	long	unused25;	/*    reserved for future use*/
	long	unused26;	/*    reserved for future use*/
	long	leven;		/* RA data-evenly-spaced flag*/
	long	lpspol;		/*    station polarity flag  */
	long	lovrok;		/*    overwrite permission   */
	long	lcalda;		/*    calc distance, azimuth */
	long	unused27;	/*    reserved for future use*/
	char	kstnm[8];	/*  F station name           */
	char	kevnm[16];	/*    event name             */
	char	khole[8];	/*    man-made event name    */
	char	ko[8];		/*    event origin time id   */
	char	ka[8];		/*    1st arrival time ident */
	char	kt0[8];		/*    time pick 0 ident      */
	char	kt1[8];		/*    time pick 1 ident      */
	char	kt2[8];		/*    time pick 2 ident      */
	char	kt3[8];		/*    time pick 3 ident      */
	char	kt4[8];		/*    time pick 4 ident      */
	char	kt5[8];		/*    time pick 5 ident      */
	char	kt6[8];		/*    time pick 6 ident      */
	char	kt7[8];		/*    time pick 7 ident      */
	char	kt8[8];		/*    time pick 8 ident      */
	char	kt9[8];		/*    time pick 9 ident      */
	char	kf[8];		/*    end of event ident     */
	char	kuser0[8];	/*    available to user      */
	char	kuser1[8];	/*    available to user      */
	char	kuser2[8];	/*    available to user      */
	char	kcmpnm[8];	/*  F component name         */
	char	knetwk[8];	/*    network name           */
	char	kdatrd[8];	/*    date data read         */
	char	kinst[8];	/*    instrument name        */
};

/* a SAC structure containing all null values */
static struct sac sac_null = {
-12345., -12345., -12345., -12345., -12345.,
-12345., -12345., -12345., -12345., -12345.,
-12345., -12345., -12345., -12345., -12345.,
-12345., -12345., -12345., -12345., -12345.,
-12345., -12345., -12345., -12345., -12345.,
-12345., -12345., -12345., -12345., -12345.,
-12345., -12345., -12345., -12345., -12345.,
-12345., -12345., -12345., -12345., -12345.,
-12345., -12345., -12345., -12345., -12345.,
-12345., -12345., -12345., -12345., -12345.,
-12345., -12345., -12345., -12345., -12345.,
-12345., -12345., -12345., -12345., -12345.,
-12345., -12345., -12345., -12345., -12345.,
-12345., -12345., -12345., -12345., -12345.,
-12345, -12345, -12345, -12345, -12345,
-12345, -12345, -12345, -12345, -12345,
-12345, -12345, -12345, -12345, -12345,
-12345, -12345, -12345, -12345, -12345,
-12345, -12345, -12345, -12345, -12345,
-12345, -12345, -12345, -12345, -12345,
-12345, -12345, -12345, -12345, -12345,
-12345, -12345, -12345, -12345, -12345,
{ '-','1','2','3','4','5',' ',' ' },
{ '-','1','2','3','4','5',' ',' ',' ',' ',' ',' ',' ',' ',' ',' ' },
{ '-','1','2','3','4','5',' ',' ' }, { '-','1','2','3','4','5',' ',' ' },
{ '-','1','2','3','4','5',' ',' ' }, { '-','1','2','3','4','5',' ',' ' },
{ '-','1','2','3','4','5',' ',' ' }, { '-','1','2','3','4','5',' ',' ' },
{ '-','1','2','3','4','5',' ',' ' }, { '-','1','2','3','4','5',' ',' ' },
{ '-','1','2','3','4','5',' ',' ' }, { '-','1','2','3','4','5',' ',' ' },
{ '-','1','2','3','4','5',' ',' ' }, { '-','1','2','3','4','5',' ',' ' },
{ '-','1','2','3','4','5',' ',' ' }, { '-','1','2','3','4','5',' ',' ' },
{ '-','1','2','3','4','5',' ',' ' }, { '-','1','2','3','4','5',' ',' ' },
{ '-','1','2','3','4','5',' ',' ' }, { '-','1','2','3','4','5',' ',' ' },
{ '-','1','2','3','4','5',' ',' ' }, { '-','1','2','3','4','5',' ',' ' },
{ '-','1','2','3','4','5',' ',' ' }
};

/* definitions of constants for SAC enumerated data values */
/* undocumented == couldn't find a definition for it (denio, 07/15/88) */
#define IREAL   0								/* undocumented              */
#define ITIME   1								/* file: time series data    */
#define IRLIM   2								/* file: real&imag spectrum  */
#define IAMPH   3								/* file: ampl&phas spectrum  */
#define IXY     4								/* file: gen'l x vs y data   */
#define IUNKN   5								/* x data: unknown type      */
												/* zero time: unknown        */
												/* event type: unknown       */
#define IDISP   6								/* x data: displacement (nm) */
#define IVEL    7								/* x data: velocity (nm/sec) */
#define IACC    8								/* x data: accel (cm/sec/sec)*/
#define IB      9								/* zero time: start of file  */
#define IDAY   10								/* zero time: 0000 of GMT day*/
#define IO     11								/* zero time: event origin   */
#define IA     12								/* zero time: 1st arrival    */
#define IT0    13								/* zero time: user timepick 0*/
#define IT1    14								/* zero time: user timepick 1*/
#define IT2    15								/* zero time: user timepick 2*/
#define IT3    16								/* zero time: user timepick 3*/
#define IT4    17								/* zero time: user timepick 4*/
#define IT5    18								/* zero time: user timepick 5*/
#define IT6    19								/* zero time: user timepick 6*/
#define IT7    20								/* zero time: user timepick 7*/
#define IT8    21								/* zero time: user timepick 8*/
#define IT9    22								/* zero time: user timepick 9*/
#define IRADNV 23								/* undocumented              */
#define ITANNV 24								/* undocumented              */
#define IRADEV 25								/* undocumented              */
#define ITANEV 26								/* undocumented              */
#define INORTH 27								/* undocumented              */
#define IEAST  28								/* undocumented              */
#define IHORZA 29								/* undocumented              */
#define IDOWN  30								/* undocumented              */
#define IUP    31								/* undocumented              */
#define ILLLBB 32								/* undocumented              */
#define IWWSN1 33								/* undocumented              */
#define IWWSN2 34								/* undocumented              */
#define IHGLP  35								/* undocumented              */
#define ISRO   36								/* undocumented              */
#define INUCL  37								/* event type: nuclear shot  */
#define IPREN  38								/* event type: nuke pre-shot */
#define IPOSTN 39								/* event type: nuke post-shot*/
#define IQUAKE 40								/* event type: earthquake    */
#define IPREQ  41								/* event type: foreshock     */
#define IPOSTQ 42								/* event type: aftershock    */
#define ICHEM  43								/* event type: chemical expl */
#define IOTHER 44								/* event type: other source  */
												/* data quality: other problm*/
#define IGOOD  45								/* data quality: good        */
#define IGLCH  46								/* data quality: has glitches*/
#define IDROP  47								/* data quality: has dropouts*/
#define ILOWSN 48								/* data quality: low s/n     */
#define IRLDTA 49								/* data is real data         */
#define IVOLTS 50								/* file: velocity (volts)    */
#define INIV51 51								/* undocumented              */
#define INIV52 52								/* undocumented              */
#define INIV53 53								/* undocumented              */
#define INIV54 54								/* undocumented              */
#define INIV55 55								/* undocumented              */
#define INIV56 56								/* undocumented              */
#define INIV57 57								/* undocumented              */
#define INIV58 58								/* undocumented              */
#define INIV59 59								/* undocumented              */
#define INIV60 60								/* undocumented              */

/* True/false definitions */
#ifndef TRUE
#define FALSE 0
#define TRUE !FALSE
#endif

/* Format strings for writing headers for SAC ASCII files */
#define FCS "%15.7g%15.7g%15.7g%15.7g%15.7g\n"	/* for floats */
#define ICS "%10d%10d%10d%10d%10d\n"			/* for integers */
#define CCS1 "%-8.8s%-8.8s%-8.8s\n"				/* for strings */
#define CCS2 "%-8.8s%-16.16s\n"					/* for strings */

/* Construct for writing data to SAC ASCII files */
/*
	for (i = 0; i < data_hdr->nsamples; i++)
	{
		fprintf (outfile, "%15.7f", *seismic_data_ptr++);
		if ( ( ((i+1)%5) == 0) && (i > 0) ) fprintf (outfile, "\n");
	}
*/