/*
	wvf_msng.c

	wm1.bat = cl /AL wvf_msng.c getfile

	removes points with missing values

*/
#include <stdio.h>
#include <string.h>
#include "nomouse.h"

#define MAX_PTS 100000
#define MAX_HIT 32000

int huge Asn[MAX_PTS];

main(int argc,char *argv[])

{
	long i,j,k;
	long i1,i2,i3,nummsng=0;
	char string[100],name1[100],name2[100];
	char str1[50],str2[50],str3[50],str4[50],str5[50],cr=0xd;
	int num=0,hit=0;
	double msng,fmsng,vx,s,z0,x,y,z,Z;
	long lz;
	FILE *fpin,*fpout;

	num=get_file_name("*.wvf",name1,num);
	if(num<0)
		exit(0);
	fpin=fopen(name1,"rt");
	if(!fpin)
	{
		printf("Could not open source WVF file '%s'\n\n",name1);
		exit(0);
	}
	printf("Give output file name.\n");
	scanf("%s",string);
	for(k=0;k<strlen(string);k++)
		if(string[k]=='.')
			string[k]='\0';
	sprintf(name2,"%s.wvf",string);
	if(strcmpi(name1,name2)==0)
	{
		printf("Output filename must be different than input name.\n\n");
		exit(0);
	}
	fpout=fopen(name2,"wt");
	if(!fpout)
	{
		printf("Could not open output WVF file '%s'\n",name2);
		exit(0);
	}
	printf("Give missing value for Z im meters.  (Original data value)\n");
	scanf("%lf",&msng);
	printf("Scanning for comments and parameters\n");
	while(fgets(string,100,fpin))
	{
		if(string[0]=='#')
		{
			sscanf(string+1,"%s%s%s%s%s",str1,str2,str3,str4,str5);
			fprintf(fpout,"%s",string);
			if(strcmpi(str1,"Vertical")==0&&strcmpi(str2,"Exaggeration")==0)
			{
				sscanf(str4,"%lf",&vx);
				printf("Vertical Exaggeration = %10.6lf\n",vx);
				hit+=1;
			}
			if(strcmpi(str1,"Scale")==0)
			{
				sscanf(str3,"%lf",&s);
				printf("                Scale = %10.6lf\n",s);
				hit+=2;
			}
			if(strcmpi(str1,"Z0")==0)
			{
				sscanf(str3,"%lf",&z0);
				printf("                   Z0 = %10.6lf\n",z0);
				hit+=4;
			}
		}
	}
	if(hit!=7)
	{
		printf(
			"Could not find all 3 parameters -- Scale, Vert.Exagg., and Z0\n\n");
		exit(0);
	}
	s*=1000.0;
	vx/=s;
	i=j=0;
	printf("\n");
	rewind(fpin);
	printf("\nScanning vertexes.\n");
	while(fgets(string,100,fpin))
	{
		if(string[0]=='v')
		{
			sscanf(string+1,"%lf%lf%lf",&x,&y,&z);
			Z=z*vx+z0;
			lz=Z*1000000.0;
			Z=lz;
			Z/=1000000.0;
			if(Z==msng)
			{
				nummsng+=1;
				Asn[i]=-1;
			}
			else
			{
				Asn[i+1]=j+1;
				fprintf(fpout,"v %0.6lf %0.6lf %0.6lf\n",x,y,z);
				j+=1;
				if(j>=MAX_PTS)
				{
					printf("Too many good points.\n\n");
					exit(0);
				}
			}
			i+=1;
			printf("%c %6ld %6ld %6ld",cr,nummsng,j,i);
		}
	}
	rewind(fpin);
	i=j=0;
	printf("\nScanning tringles\n");
	while(fgets(string,100,fpin))
	{
		if(string[0]=='f')
		{
			sscanf(string+1,"%ld%ld%ld",&i1,&i2,&i3);
			if(Asn[i1]>0&&Asn[i2]>0&&Asn[i3]>0)
			{
				fprintf(fpout,"f %d %d %d\n",Asn[i1],Asn[i2],Asn[i3]);
				j+=1;
			}
			i+=1;
			printf("%c %6ld %6ld",cr,j,i);
		}
	}
	fclose(fpin);
	fclose(fpout);
}


                                                                                                                   