sucommand.c

su 的源代码库

/* Copyright (c) Colorado School of Mines, 2006.*/
/* All rights reserved.                       */

/* Copyright (c) Virginia Tech */
/* All rights reserved.		       */

/* SUCOMMAND: $Revision: 1.3 $ ; $Tue Aug 31 11:16:11 EDT 1999 $	*/

#include "su.h"
#include "segy.h"

/*********************** self documentation **********************/
char *sdoc[] = {
" 									",
" SUCOMMAND - pipe traces having the same key header word to command	",
" 									",
"     sucommand <stdin >stdout [Optional parameters]			",
" 									",
" Required parameters:							",
" 	none								",
" 									",
" Optional parameters: 							",
" 	verbose=0		wordy output				",
" 	key=cdp			header key word to pipe on		",
" 	command=\"suop nop\"    command piped into			",
"	dir=0		0:  change of header key			",
"			-1: break in monotonic decrease of header key	",
"			+1: break in monotonic increase of header key	",
" 									",
" 									",
"Notes:									",
" This program permits limited parallel processing capability by opening",
" pipes for processes, signaled by a change in key header word value.	",
" 									",
" 									",
NULL};

/* Credits:
 *	VT: Matthias Imhof
 *
 * Note:
 *	The "valxxx" subroutines are in su/lib/valpkge.c.  In particular,
 *      "valcmp" shares the annoying attribute of "strcmp" that
 *		if (valcmp(type, val, valnew) {
 *			...
 *		}
 *	will be performed when val and valnew are different.
 *
 */
/**************** end self doc ***********************************/


segy trace;

int
main(int argc, char **argv)
{
	cwp_String key;	/* header key word from segy.h		*/
	cwp_String command;	/* pipe command			*/
	cwp_String type;/* ... its type				*/
	int indx;	/* ... its index			*/
	int nsegy;	/* number of bytes in the segy		*/
	Value val;	/* value of key in current gather	*/
	Value valnew;	/* value of key in trace being treated	*/
	int fold;	/* number of traces per gather		*/
	int verbose;	/* verbose flag				*/
	int dir;	/* direction flag			*/
	FILE *pipe;

	/* Initialize */
	initargs(argc, argv);
	requestdoc(1);

	/* Set parameters */
	if (!getparint   ("verbose", &verbose))	 verbose = 0;
	if (!getparstring("command", &command))	 command = "suop nop";
	if (!getparstring("key", &key))		 key = "cdp";
	if (!getparint   ("dir", &dir))	 	 dir = 0;
	dir=-dir;

	type = hdtype(key);
	indx = getindex(key);

	/* Set up for first trace (must compare new key field each time) */
	pipe = popen(command, "w");
	fold = 1;
	
	nsegy = gettr(&trace);
	gethval(&trace, indx, &val);
	fputtr(pipe, &trace);

	/* Loop over traces */
	for(;;)
	{	int	cmp;
		nsegy = gettr(&trace);
		if (! nsegy) break;
		
		gethval(&trace, indx, &valnew);
		cmp = valcmp(type, val, valnew);
		val = valnew;
		
		if (((dir==0) && cmp) || ((dir != 0) && (dir!=cmp))) {	
			/* val and valnew differ, indicating a  */
			/* new gather ? 				*/
			if (verbose) {
				fprintf(stderr, "val=");
				fprintfval(stderr, type, val);
				fprintf(stderr, "\tgather size=%d\n", fold);
			}

			pclose(pipe);

			/* Set up for next gather */
			
			pipe = popen(command, "w");
			fold = 1;
			fputtr(pipe, &trace);
		}
		
		else {	/* still in same gather */
			fputtr(pipe, &trace);
			++fold;
		}
	}
	
	if (verbose) {  fprintf(stderr, "val=");
			fprintfval(stderr, type, val);
			fprintf(stderr, "\tgather size=%d\n", fold);
	}
	pclose(pipe);
	
	return(CWP_Exit());
}