rdseed.h

解吸SEED格式的源代码

/*                 +=======================================+                 */
/*=================|   Volume and station header tables    |=================*/
/*                 +=======================================+                 */

/* volume header blockettes */
EXTERN struct type10	type10;		/* type 010 blockette */
EXTERN struct type11	type11;		/* type 011 blockette */
EXTERN struct type12	type12;		/* type 012 blockette */

/* pointers to starts of abbreviation dictionary tables (linked lists) */
EXTERN struct type30	*type30_head EQNULL;	/* type 030 blockette */
EXTERN struct type31	*type31_head EQNULL;	/* type 031 blockette */
EXTERN struct type32	*type32_head EQNULL;	/* type 032 blockette */
EXTERN struct type33	*type33_head EQNULL;	/* type 033 blockette */
EXTERN struct type34	*type34_head EQNULL;	/* type 034 blockette */
EXTERN struct type35	*type35_head EQNULL;	/* type 035 blockette */

/* pointers to starts of response dictionaries */
EXTERN struct type41	*type41_head EQNULL;	/* type 041 blockette */
EXTERN struct type42	*type42_head EQNULL;	/* type 042 blockette */
EXTERN struct type43	*type43_head EQNULL;	/* type 043 blockette */
EXTERN struct type44	*type44_head EQNULL;	/* type 044 blockette */
EXTERN struct type45	*type45_head EQNULL;	/* type 045 blockette */
EXTERN struct type46	*type46_head EQNULL;	/* type 046 blockette */
EXTERN struct type47	*type47_head EQNULL;	/* type 047 blockette */
EXTERN struct type48	*type48_head EQNULL;	/* type 048 blockette */

/* pointers to start of station header tables (linked list) */
EXTERN struct type50	*type50_head EQNULL;	/* type 050 blockette */
												/* other type 05x blockettes */
												/* are subsumed under this */

/* pointers to starts of time span control header structures (linked lists) */
EXTERN struct type70	*type70_head EQNULL;	/* type 070 blockette */
EXTERN struct type71	*type71_head EQNULL;	/* type 071 blockette */
EXTERN struct type72	*type72_head EQNULL;	/* type 072 blockette */
EXTERN struct type73	*type73_head EQNULL;	/* type 073 blockette */
EXTERN struct type74	*type74_head EQNULL;	/* type 074 blockette */

/*                 +=======================================+                 */
/*=================|  Seismic data and related variables   |=================*/
/*                 +=======================================+                 */

EXTERN struct type50	*current_station EQNULL;/* ptr to current station */
EXTERN struct type52	*current_channel EQNULL;/* ptr to current channel */
EXTERN struct type71	*current_origin  EQNULL;/* ptr to current origin record */

EXTERN struct data_hdr	*data_hdr;				/* data header */
EXTERN struct time		last_time;				/* time of previous record */
EXTERN double			last_sample_rate;		/* sample rate, previous rcd */
EXTERN int				last_nsamples;			/* # samples, previous rcd */
EXTERN int		seis_buffer_length;

EXTERN float		*seismic_data;				/* seismic data */
EXTERN float			*seismic_data_ptr;		/* ptr to above */
EXTERN char				decode[_DECODEKEY];		/* data decoder key */

EXTERN int				numitems;				/* number of sgms wanted */
EXTERN char				*item[_MAXSGMS];		/* list of sgms wanted */
EXTERN long int			start_record;			/* starting rec # of sgm */
EXTERN double sample_rate_accum;				/* average sample rate accum */
EXTERN int sample_rate_count;					/* count for average sample rate calc */
EXTERN double adj_tolerance;					/* tolerance adjust, read from enviornment */

/*                 +=======================================+                 */
/*=================|          process definitions          |=================*/
/*                 +=======================================+                 */

/* Utilities */
void read_logical_record ();					/* process a logical rec */
int memncpy ();									/* cp, null-term. n chars */
void read_blockette ();							/* read a blockette */
void pagetop ();								/* top-of-page printer */
void allocation_error ();						/* memory allocation error */
int parse_int ();								/* parse an integer field */
long int parse_long ();							/* parse a long integer field */
double parse_double ();							/* parse a double field */
char *parse_nchar ();							/* parse n-char wide ASCII */
char *parse_varlstr ();							/* parse vbl-length ASCII */
struct time timeadd ();							/* add seconds to a time */
struct time timesub();

struct time timeadd_double ();						/* add double seconds to a time */
struct time timeaddphase ();							/* add seconds to a time */
int timecmp ();									/* compare two times */
double timedif ();								/* difference of two times */
void snoop ();									/* list block contents */
void snoop_data ();								/* sub for snoop for data */
unsigned long int swap_4byte ();				/* swap bytes in 4-byte words */
unsigned short int swap_2byte ();				/* swap bytes in 2-byte words */
int find_wordorder ();							/* determine word order */
struct optstruct *cmdlineproc ();				/* command line processor */

/* Functions to read and process the volume identifier block */
void process_volh ();							/* process volume header */
void parse_type10 ();							/* parse volume id blkette */
void parse_type11 ();							/* parse vol stn hdr index */
void parse_type12 ();							/* parse vol time span index */
void print_volh ();								/* print volume header */
void print_type10 ();							/* print volume id blkette */
void print_type11 ();							/* print vol stn hdr index */
void print_type12 ();							/* print vol time span index */

/* Function to process a blank or unknown type block */
void process_blank ();

/* Functions to read and process the abbreviation dictionary block */
void process_abrvd ();							/* process abbreviation hdr */
void parse_type30 ();							/* parse type 30 blockette */
void parse_type31 ();							/* parse type 31 blockette */
void parse_type32 ();							/* parse type 32 blockette */
void parse_type33 ();							/* parse type 33 blockette */
void parse_type34 ();							/* parse type 34 blockette */
void parse_type35 ();							/* parse type 35 blockette */
void print_abrvd ();							/* print abbrev dictionaries */
void print_type30 ();							/* print type30 table */
void print_type31 ();							/* print type31 table */
void print_type32 ();							/* print type32 table */
void print_type33 ();							/* print type33 table */
void print_type34 ();							/* print type34 table */
void print_type35 ();							/* print type35 table */

/* Functions to read and process the response dictionary blocks */
void parse_type41 ();							/* parse sym coefficients blkt */
void parse_type42 ();							/* parse polynomial blkt */
void parse_type43 ();							/* parse poles & zeroes blkt */
void parse_type44 ();							/* parse coefficients blkt */
void parse_type45 ();							/* parse rspns list blkt */
void parse_type46 ();							/* parse generic rspns blkt */
void parse_type47 ();							/* parse generic rspns blkt */
void parse_type48 ();							/* parse ch sensitivity */

/* Functions to read and process the station header blocks */
int process_stnh ();							/* process station headers */
int parse_type50 ();		/* parse station id blkt */
void parse_type51 ();		/* parse station comment */
int parse_type52 ();		/* parse channel id */
struct response *get_response ();	/* get a response for a chl */
void parse_type53 ();			/* parse poles & zeroes blkt */
void parse_type54 ();			/* parse coefficients blkt */
void parse_type55 ();			/* parse rspns list blkt */
void parse_type56 ();			/* parse generic rspns blkt */
void parse_type57 ();			/* parse decimation blkt */
void parse_type58 ();			/* parse ch sensitivity */
void parse_type59 ();			/* parse ch comment */
void parse_type60 ();			/* parse ch comment */
void parse_type61 ();			/* parse sym coefficients blkt */
void parse_type62 ();			/* parse polynomial blkt */

void print_stnh ();			/* print station headers */
void print_type50 ();			/* print station id blkt */
void print_type51 ();			/* print station comment */
void print_channel ();			/* print all channel info */
void print_type52 ();			/* print channel id */
void print_response ();			/* print all response info */
void print_type53 ();			/* print poles & zeroes blkt */
void print_type54 ();			/* print coefficients blkt */
void print_type55 ();			/* print rspns list blkt */
void print_type56 ();			/* print generic rspns blkt */
void print_type57 ();			/* print decimation blkt */
void print_type58 ();			/* print ch sensitivity */
void print_type59 ();			/* print ch comment */
void print_type60 ();			/* print ch comment */
void print_type61 ();			/* print sym coefficients blkt */
void print_type62 ();			/* print polynomial blkt */

/* Functions to read and process the time span header blocks */
void process_timeh ();							/* process time span headers */
void parse_type70 ();							/* parse time span id blkt */
void parse_type71 ();							/* parse hypocenter id blkt */
void parse_type72 ();							/* parse event phases blkt */
void parse_type73 ();							/* parse t.s. data st index */
void print_timeh ();							/* print time span headers */
void print_type70 ();							/* print time span id blkt */
void print_type71 ();							/* print hypocenter id blkt */
void print_type72 ();							/* print event phases blkt */
void print_type73 ();							/* print t.s. data st index */

/* Functions to read and process the data */
void process_data ();							/* process data blocks */
void init_data_hdr ();							/* initialize a data header */
void convert_seedhdr ();						/* parse fixed size data hdr */
void decode_asro ();							/* decode asro data */
void decode_cdsn ();							/* decode cdsn data */
void decode_dwwssn ();							/* decode dwwssn data */
void decode_steim ();							/* decode steim data */
void decode_sro ();								/* decode sro data */
void decode_32bit ();							/* decode 32-bit int data */
void output_data (); 							/* write output data */

/*                 +=======================================+                 */
/*=================|          end of header file           |=================*/
/*                 +=======================================+                 */

/* the next line should be the end of the header file */
#endif