parse_53.c

解吸SEED格式的源代码

/*===========================================================================*/
/* SEED reader     |               parse_type53              |  station header */
/*===========================================================================*/
/*
	Name:		parse_type53
	Purpose:	parse a response (poles and zeroes) (type 53) blockette
				into a globally-available table
	Usage:		void parse_type53 ();
				char *blockette;
				parse_type53 (blockette);
	Input:		blockette = pointer to beginning of blockette
	Output:		none
	Externals:	temp_char
				type53_head - 1st entry in table; defined in structures.h,
					allocated in globals.h
	Warnings:	none
	Errors:		none
	Called by:	process_stationh
	Calls to:	parse_int - get an integer from a blockette 
				parse_double - get a double from a blockette 
				parse_nchar - get a fixed-length string from a blockette 
	Algorithm:	parse the blockette into a globally-available table
	Notes:		none
	Problems:	none known
	References:	Halbert et al, 1988; see main routine
	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
				08/19/89  Dennis O'Neill  streamlined linked list construction
				05/11/94  CL			  took out the malloc if num==0
*/

#include "rdseed.h"

void parse_type53 (blockette)
char *blockette;								/* ptr to start of blockette */
{
	char *blockette_ptr;						/* ptr to inside blockette */
	struct type53 *type53; 							/* genl member of table */
	int i;										/* counter */
	struct response *response;					/* channel response ptr */

	/* point to beginning of information, past type and length */
	blockette_ptr = blockette;
	blockette_ptr += 7;

/*                 +=======================================+                 */
/*=================| allocate space and initialize a type53  |=================*/
/*                 +=======================================+                 */

	type53 = (struct type53 *) alloc_linklist_element (sizeof (struct type53),
		"parse_type53");	

/*                 +=======================================+                 */
/*=================|   recover information from the type53   |=================*/
/*                 +=======================================+                 */

	/* recover transfer function type */
	type53->transfer = parse_nchar (&blockette_ptr, 1);

	/* recover cascade sequence number */
	type53->stage = parse_int (&blockette_ptr, 2);

	/* recover response in units (0 if none) */
	type53->input_units_code = parse_int (&blockette_ptr, 3);

	/* recover response out units */
	type53->output_units_code = parse_int (&blockette_ptr, 3);

	/* recover AO normalization factor (1.0 if none) */
	type53->ao_norm = parse_double (&blockette_ptr, 12);

	/* recover normalization frequency */
	type53->norm_freq = parse_double (&blockette_ptr, 12);

	/* recover number of complex zeroes */
	type53->number_zeroes = parse_int (&blockette_ptr, 3);

	/* allocate space for complex zeroes */
	if (type53->number_zeroes != 0)
		type53->zero = (struct type53sub *) alloc_linklist_element 
			(sizeof (struct type53sub) * type53->number_zeroes, "parse_type53->zero");

	/* recover complex zeroes */
	for (i = 0; i < type53->number_zeroes; i++)
	{
		type53->zero[i].real = parse_double (&blockette_ptr, 12);
		type53->zero[i].imag = parse_double (&blockette_ptr, 12);
		type53->zero[i].real_error = parse_double (&blockette_ptr, 12);
		type53->zero[i].imag_error = parse_double (&blockette_ptr, 12);
	}

	/* recover number of complex poles */
	type53->number_poles = parse_int (&blockette_ptr, 3);

	/* allocate space for complex poles */
	if (type53->number_poles != 0)
		type53->pole = (struct type53sub *) alloc_linklist_element
			(sizeof (struct type53sub) * type53->number_poles, "parse_type53->pole");

	/* recover complex poles */
	for (i = 0; i < type53->number_poles; i++)
	{
		type53->pole[i].real = parse_double (&blockette_ptr, 12);

		type53->pole[i].imag = parse_double (&blockette_ptr, 12);
		type53->pole[i].real_error = parse_double (&blockette_ptr, 12);
		type53->pole[i].imag_error = parse_double (&blockette_ptr, 12);

		if ((!strcmp(type53->transfer,"A") || !strcmp(type53->transfer,"B")) &&
		     type53->pole[i].real > 0.0) 
			fprintf(stderr, "Warning! Station/channel: %s/%s has a pole with a real value > 0.0: real value=%6.2f\n", 
				current_station->station, 
				current_channel->channel,
				type53->pole[i].real);
#if 0 
=========================

		if ((!strcmp(type53->transfer,"A") || !strcmp(type53->transfer,"B")) &&
		     type53->pole[i].real_error > 0.0)  
			fprintf(stderr, "Warning! Station/channel: %s/%s has a pole with a error value > 0.0\n", 
				current_station->station, 
                current_channel->channel,   
                type53->pole[i].real_error);
#endif


	}

/*                 +=======================================+                 */
/*=================|   link the type53 into the type53 table   |=================*/
/*                 +=======================================+                 */

	/* initialize the pointer to the next record as end-of-list marker */
	type53->next = NULL;

	response = get_response ('P');
	response->ptr.type53 = type53;	
}