mkparfile.c

seismic software,very useful

/* MKPARFILE - convert ascii file to par file format */

#include "par.h"

/*********************** self documentation ******************************/
string sdoc =
" 									\n"
" MKPARFILE - convert ascii to par file format 				\n"
" 									\n"
" mkparfile <stdin >stdout 						\n"
" 									\n"
" Optional parameters:							\n"
" 	string1=\"par1=\"	first par string			\n"
" 	string2=\"par2=\"	second par string			\n"
" 									\n"
" This is a tool to convert values written line by line to parameter 	\n"
" vectors in the form expected by getpar.  For example, if the input	\n"
" file looks like:							\n"
" 	t0 v0								\n"
" 	t1 v1								\n"
"	...								\n"
" then									\n"
"	mkparfile <input >output string1=tnmo string2=vnmo		\n"
" yields:								\n"
"	tnmo=t0,t1,...							\n"
"	vnmo=v0,v1,...							\n"
" 									\n"
;
/**************** end self doc *******************************************/

/* Credits:
 *	CWP: Jack
 */


/* Caveat: A more general tool allowing n1 strings would be desirable. */

main(int argc, char **argv)
{
	int i2, n2 = 0;
	float x1, x2;
	char *string1;
	char *string2;
	char buf[BUFSIZ];
	FILE *datafp;


	/* Hook up getpar */
	initargs(argc, argv);
	askdoc(1);


	/* Get parameters and set up tmpfile */
	if (!getparstring("string1", &string1))	string1 = "par1";
	if (!getparstring("string2", &string2))	string2 = "par2";
	datafp = etmpfile();


	/* Extract x1's from data and save data for later pass over x2's */
	if (!gets(buf))  err("no data found");
	if (2 == sscanf(buf, "%f %f", &x1, &x2)) { /* no comma first time */
		printf("%s=%g", string1, x1);
		efwrite(&x2, FSIZE, 1, datafp);
	} else  err("line #%d: scan failed:\n%s", n2+1, buf);
	++n2;

	while (gets(buf)) {
		if (2 == sscanf(buf, "%f %f", &x1, &x2)) {
			printf(",%g", x1);
			efwrite(&x2, FSIZE, 1, datafp);
		} else  err("line #%d: scan failed:\n%s", n2+1, buf);
		++n2;
	}
	putchar('\n');

	/* Rewind and get the x2's */
	rewind(datafp);

	efread(&x2, FSIZE, 1, datafp);
	printf("%s=%g", string2, x2);

	for (i2 = 1; i2 < n2; ++i2) {
		efread(&x2, FSIZE, 1, datafp);
		printf(",%g", x2);
	}
	putchar('\n');

	return EXIT_SUCCESS;
}