/*
	wvf2ec.c

	reads a WVF file and makes an earth coordinate file of the 
	verticies of the triangles in 9 floats each

	lat0 lon0 z0 lat1 lon1 z1 lat2 lon2 z2

	we1.bat = cl /AL wvf2ec.c proj1

*/
#include <stdio.h>
#include <string.h>

struct coef
{
	double X0,Y0,Z0;
	double scale,zscale;
};


struct control
{
	char proj;
	int pnum;
	double lon_0,lon_1,lon_2,lon_p;
	double lat_0,lat_1,lat_2,lat_p;
	double pixsize,scale;
	double x0,y0;
	double c,C,F,H,k,k0,k1,k2,n,P,rho_0;
	double omega; /* elevation (tilt of surface plane) */
	double gamma; /* azimuth (east of north) */
	double horiz; /* radius of horizon arc */
}Map;


main(int argc,char *argv[])

{
	int i,j,k;
	char infile[100],string[100];
	FILE *fpwvf;

	if(argc<2)
	{
		printf("\n\nwvf2ec [infile]\n\n");
		printf("    where:  infile = name of WVF file\n\n");
		exit(0);
	}
	strcpy(infile,argv[1]);
	for(i=0;i<strlen(infile);i++)
		if(infile[i]=='.')
			infile[i]='\0';
	sprintf(string,"%s.wvf",infile);
	fpwvf=fopen(string,"rt");
	if(!fpwvf)
	{
		printf("Could not open '%s' to read.\n\n",string);
		exit(0);
	}
	read_map_data(fpwvf);
}

/**********************************************************************
**
**
**
**********************************************************************/

int read_map_data(FILE *fp)

{
	char string[100];
	char str[7][50];
	int num;

	rewind(fp);
	do
	{
		fgets(string,100,fp);
		if(string[0]=='#')
		{
			num=sscanf(string+2,"%s%s%s%s%s%s%s",
				str[0],str[1],str[2],str[3],str[4],str[5],str[6]);
			if(strcmpi(str[0],"Projection")==0)
			{
				Map.proj=str[num-1][2];
printf("Projection = '%c'\n",Map.proj);
			}
			if(strcmpi("longitude_of_central_meridian_of_map",str[0])==0)
			{
				sscanf(str[2],"%lf",&Map.lon_0);
printf("longitude_of_central_meridian_of_map = %lf\n",Map.lon_0);
			}
			if(strcmpi("central_latitude_of_map",str[0])==0)
			{
				sscanf(str[2],"%lf",&Map.lat_0);
printf("central_latitude_of_map = %lf\n",Map.lat_0);
			}
			if(strcmpi("longitude_(degrees)_of_central_meridian_of_map",str[0])==0)
			{
				sscanf(str[2],"%lf",&Map.lon_0);
printf("longitude_(degrees)_of_central_meridian_of_map = %lf\n",Map.lon_0);
			}
			if(strcmpi("latitude_and_longitude_of_first_point",str[0])==0)
			{
				sscanf(str[2],"%lf",&Map.lat_1);
				sscanf(str[3],"%lf",&Map.lon_1);
printf("latitude_and_longitude_of_first_point = %lf %lf\n",Map.lat_1,Map.lon_1);
			}
			if(strcmpi("latitude_and_longitude_of_second_point",str[0])==0)
			{
				sscanf(str[2],"%lf",&Map.lat_2);
				sscanf(str[3],"%lf",&Map.lon_2);
printf("latitude_and_longitude_of_second_point = %lf\n",Map.lat_2,Map.lon_2);
			}
			if(strcmpi("latitude_and_longitude_of_origin",str[0])==0)
			{
				sscanf(str[2],"%lf",&Map.lat_0);
				sscanf(str[3],"%lf",&Map.lon_0);
printf("latitude_and_longitude_of_origin = %lf %lf\n",Map.lat_0,Map.lon_0);
			}
			if(strcmpi("latitude_of_2_standard_parallels",str[0])==0)
			{
				sscanf(str[2],"%lf",&Map.lat_1);
				sscanf(str[3],"%lf",&Map.lat_2);
printf("latitude_of_2_standard_parallels = %lf %lf\n",Map.lat_1,Map.lat_2);
			}
			if(strcmpi("latitude_of_standard_parallel",str[0])==0)
			{
				sscanf(str[2],"%lf",&Map.lat_1);
printf("latitude_of_standard_parallel = %lf\n",Map.lat_1);
			}
			if(strcmpi("longitude_of_central_meridian",str[0])==0)
			{
				sscanf(str[2],"%lf",&Map.lon_0);
printf("longitude_of_central_meridian = %lf\n",Map.lon_0);
			}
			if(strcmpi("scale_factor",str[0])==0)
			{
				sscanf(str[2],"%lf",&Map.k0);
printf("scale_factor = %lf\n",Map.k0);
			}
			if(strcmpi("altitude_of_view_in_kilometers",str[0])==0)
			{
				sscanf(str[2],"%lf",&Map.H);
printf("altitude_of_view_in_kilometers = %lf\n",Map.H);
			}
			if(strcmpi("viewing_azimuth_(N=0)",str[0])==0)
			{
				sscanf(str[2],"%lf",&Map.gamma);
printf("viewing_azimuth_(N=0) = %lf\n",Map.gamma);
			}
			if(strcmpi("viewing_angle_up_from_straight_down",str[0])==0)
			{
				sscanf(str[2],"%lf",&Map.omega);
printf("viewing_angle_up_from_straight_down = %lf\n",Map.omega);
			}
			if(strcmpi("latitude_of_origin",str[0])==0)
			{
				sscanf(str[2],"%lf",&Map.lat_0);
printf("latitude_of_origin = %lf\n",Map.lat_0);
			}
		}
	}while(string[0]=='#');
}



                     