/*
	makecol.c

	read '*.mcl' file as follows:

d:\slar\es115a.slr   		:	image file name
0.013     						:	minimum pixel size 
9195   4224   					:	rows and columns
512                        :  file header size
16                         :  row header
4                          :  row trailer
5   								:	number of latitude coefficients
	60.198 						:	b[0] -- intercept
	 7.1776e-006 				:	b[1] -- y
	-1.8923e-010 				:	b[2] -- y*y
	-1.1064e-004 				:	b[3] -- x
	-1.2695e-010 				:	b[4] -- x*x
5 									:	number of longitude coefficients
	-133.18150 					:	b[0] -- intercept
	  -1.0837e-005 			:	b[1] -- x	
	  -2.3141e-020 			:	b[2] -- x*x
	  -2.2208e-004 			:	b[3] -- y
	  -2.3462e-011 			:	b[4] -- y*y

*/
#include <stdio.h>
#include <math.h>
#include <string.h>
#include "getfile3.c"

#define STR_LEN 100
#define MIN_DIM 25

int NumImage=5;
struct dir
{
	char id;
	float scale;
	int iscale;
}D[5]=
{
	{'\0', 1.0,1},
	{'w', 3.0, 3},
	{'x', 9.0, 9},
	{'y',27.0,27},
	{'z',81.0,81}
};

main(argc,argv)

int argc;
char *argv[];

{
	int i,j,k;
	char string[100],path[100],name[100],type[100],outfile[100];
	FILE *fpin,*fpout;
	float minpix,a[8],b[8];
	int row,col,hdr,row_hdr,row_tlr;
	int nlat,nlon;

	if(get_file_name("*.mcl",string)<0)
	{
		printf("No '.mcl' files\n\n");
		exit(0);
	}
	fpin=fopen(string,"rb");
	if(!fpin)
	{
		printf("Could not open '%s'\n");
		exit(0);
	}

	fpout=fopen("mcl.col","wt");
	fprintf(fpout,"begin_vector\n");
	fprintf(fpout,"begin_collage\n");
	fprintf(fpout,"color_lut = d:\gloria\grey.lut\n");
	fprintf(fpout,"assignment = d:\gloria\grey.asn\n");
	fprintf(fpout,"missing = 0\n");
	fprintf(fpout,"minimum_pixel = %f",minpix);
	fprintf(fpout,"begin_image\n");
	fprintf(fpout,
"---------------------------------------------------------------------\n");
	fprintf(fpout,
"name         maxpix minpix  row  col hdr r_hdr r_tlr N   a0        a1\n");
	fprintf(fpout,
"------------ ------ ------ ---- ---- --- ----- ----- - ------ ----------*\n");
	for(;;)
	{
		if(fscanf(fpin,"%s",string)<1)
			break;
		split(string,path,name,type);
		printf("path = '%s'  name = '%s'  type = '%s'\n",path,name,type);
		fscanf(fpin,"%f",&minpix);
		printf("     minpix = %f\n",minpix);
		fscanf(fpin,"%d",&row);
		printf("        row = %d\n",row);
		fscanf(fpin,"%d",&col);
		printf("        col = %d\n",col);
		fscanf(fpin,"%d",&hdr);
		printf("file header = %d\n",hdr);
		fscanf(fpin,"%d",&row_hdr);
		printf(" row header = %d\n",row_hdr);
		fscanf(fpin,"%d",&row_tlr);
		printf("row trailer = %d\n",row_tlr);
		fscanf(fpin,"%d",&nlat);
		printf("nlat = %d\n",nlat);
		if(nlat<0&&nlat>7)
			exit(0);
		for(i=0;i<nlat;i++)
		{
			fscanf(fpin,"%f",a+i);
			printf("a[%d] = %20.10f\n",i,a[i]);
		}
		fscanf(fpin,"%d",&nlon);
		printf("nlon = %d\n",nlon);
		if(nlon<0&&nlon>7)
			exit(0);
		for(i=0;i<nlon;i++)
		{
			fscanf(fpin,"%f",b+i);
			printf("b[%d] = %20.10f\n",i,b[i]);
		}
		for(i=0;i<NumImage;i++)
		{
			if(row/(int)D[i].iscale<MIN_DIM&&col/(int)D[i].iscale<MIN_DIM)
				break;
			if(i==0)
				fprintf(fpout,"%s%s.%s  ",path,name,type);
			else
				fprintf(fpout,"%s%s%c.%s  ",path,name,D[i].id,type);
			if(i<NumImage-1)
				fprintf(fpout,"%f %f ",minpix*3.0*D[i].scale,minpix*D[i].scale);
			else
				fprintf(fpout,"100.000 %f ",minpix*D[i].scale);
			fprintf(fpout,"%5d %5d",row/(int)D[i].iscale,col/(int)D[i].iscale);
			if(i==0)
			{
				col-=row_hdr+row_tlr;
				fprintf(fpout," %4d %3d %3d\n",hdr,row_hdr,row_tlr);
			}
			else
				fprintf(fpout," 0 0 0\n");
			fprintf(fpout,"%d\n",nlat);
			for(k=0;k<nlat;k++)
			{
				if(k==0)
					fprintf(fpout,"%15.7e\n",a[k]);
				else
					fprintf(fpout,"%15.7e\n",a[k]*D[i].scale);
			}
			fprintf(fpout,"%d\n",nlon);
			for(k=0;k<nlon;k++)
			{
				if(k==0)
					fprintf(fpout,"%15.7e\n",b[k]);
				else
					fprintf(fpout,"%15.7e\n",b[k]*D[i].scale);
			}
		}
	}
	fprintf(fpout,"end_image\n");
	fprintf(fpout,"end_collage\n");
	fclose(fpout);
}

/**************************************************************************
**
**
**
************************************************************************* */

int split(string,path,name,type)

char *string,*path,*name,*type;

{
	int i,j,k;
	int dot=-1,last_slash=-1,len=strlen(string),num=0;

	if(strlen(string)>STR_LEN)
	{
		printf("In split() -- string 1 too long!\n\n");
		exit(0);
	}

	for(i=0;i<len;i++)
	{
		if(string[i]=='.')
			dot=i;
		if(string[i]=='\\')
			last_slash=i;
	}
	last_slash+=1;
	if(dot==-1)
		dot=len-1;
	if(last_slash>0)
	{
		for(i=0;i<last_slash;i++)
			path[i]=string[i];
	}
	path[last_slash]='\0';
	for(i=last_slash;i<dot;i++)
		name[i-last_slash]=string[i];
	name[dot-last_slash]='\0';
	dot+=1;
	for(i=dot;i<len;i++)
		type[i-dot]=string[i];
	type[len-dot]='\0';
	if(strlen(name)>0)
		num+=1;
	if(strlen(type)>0)
		num+=1;

	return(num);
}

                                                                                                         