/* This program is a Sun version of the FITS format converter to
   take NICMOS2 images extracted from ramp runs. 
   The files are assumed to contain floating point nos.  All files in a
   directory will be converted and will have .fits appended to indicate
   which are the FITS format images.
   usage: fitsarch  directory
   created: 08/11/89  mjr  { drawing from T. Greene's b2f and Irene's
   							 FITS }
   modified: 11/03/97 mjr revised for FSPEC
   modified: 11/03/97 mjr revised for newer 256x256 format
*/

#include <stdio.h>
#include <string.h>
#include <sys/types.h>
#include <dirent.h>

#define TRUE 	1
#define FALSE 	0
#define MAXSTR 180
#define FOREVER TRUE
#define NBYTES 2880		/* Fits block size */

#define NICMOS3  TRUE

/* Definitions to load for 256x256 version		*/

#if NICMOS3 

#define NUM_PIX 65536			/* 256x256 w/ no overscans */
#define B_SIZE 131072 
#define R_WFRAME TRUE
#define MAXCOLS 256 
#define MAXROWS 256
#define NBLOCKS 46		/* No. of 2880-byte blocks for biggest image */
#endif

short int libuf[MAXCOLS*MAXROWS];
short int sibuf[MAXCOLS*MAXROWS];
char buf[NBYTES],buf2[NBYTES];
char ext[] =  ".fits";
char ext2[] = ".dat";
char *name[3000];
char dummy[] = {"000000000"};

void main(argc,argv)
int argc;
char *argv[];
{
	int t_ok,t_stat,handle,n_samp,n_line,err_no,n,obs_no;
	int nframes,hdr_len,num_read,i,l,j,wframe,deg,min;
	long int nbyte;
	float dmax,dmin,itime,am,fsec;
	char cur_path[MAXSTR],*p,*hdstr,time[30],date[30],temp[80];
	char c, title[80], data[MAXSTR], filename[MAXSTR], filter[13];
	char swap_bytes = TRUE, t_append = TRUE, answer[5];
	char l_int = TRUE, *fname[3000],outname[105],gain[MAXSTR],sign[2];
	char revs[MAXSTR], lst[30], tele_ra[30],tele_dec[30];
	char telescope[4],tele[10];
	FILE *fpin,*fopen(),*fptest,*fplog;
   	char *malloc();
	void float_swap();	
	void strip(),append(),fixname();

	if ( argc != 3 ) {
	   fprintf(stderr,"Error in fspecfit. Usage: fspecfit filelist telescope.\n");
	   exit(0);
	   }
	sscanf(argv[2],"%s",telescope);
	sscanf(argv[2],"%s",tele);
  	append(tele);
	
   /* Pad data buffers with zeros to end of NBLOCKS*2880 */
	
	for (i=0; i<NBLOCKS*NBYTES/2; i++) {
		sibuf[i] = 0;
		}
	
 	swap_bytes = TRUE;


/*--- Get list of files to process -------------------*/

	fpin = fopen(argv[1],"r");
	printf(" File  %s being read.\n",argv[1]);
	
	if ( fpin == NULL ) {
		fprintf(stderr,"Error in fspecfit. File %s not found.\n",
			argv[1]);
		exit(0);
		}
	nframes = 0;
	fplog = fopen("datalog.lst","w");

	fname[0] = malloc(400000); /* Room for 1000 200 char names */

	while (! feof(fpin)) {
	    fscanf(fpin,"%s",filename);
		if ( nframes>0 ) fname[nframes]=fname[nframes-1]+l+1; /* Inc name ptr */
		l = strlen(filename);
		strcpy(fname[nframes],filename);
		nframes++;
		}
	fclose(fpin);
    nframes--;


	for (j=0; j<nframes; j++) {
		strcpy(outname,fname[j]);
		strcat(outname,ext);
		fptest = fopen(outname,"wb");
		fpin = fopen(fname[j],"r");
		fprintf(stderr,"File %s is being moved.\n",fname[j]);		
		if ( fpin == NULL ) printf(" Null pointer on %s.\n",fname[j]);

		/* Get data values read in. */

		num_read = fread((char *)sibuf,sizeof(short int),
			MAXCOLS*MAXROWS,fpin);
		nbyte = num_read * sizeof(short int);
		printf("num_read=%d nbyte=%d.\n",num_read,nbyte);

		if ( num_read != MAXCOLS*MAXROWS ) {
	       	   fprintf(stderr,"Unexpected end-of-file on %s %ld.\n",
	        	fname[j],num_read);
	           exit(0);
		   }
				
		/* Get trailing info to place in the FITS header. */

		fgets(temp,80,fpin);
		sscanf(temp,"%f %d",&itime,&wframe); 
		fgets(time,30,fpin);
		l = strlen(time);
		time[l-2] = '\0';
		fgets(date,30,fpin);
		l = strlen(date);
		date[l-2] = '\0';
		fgets(title,82,fpin);
		l = strlen(title);
		title[l-2] = '\0';
		fgets(filter,12,fpin);
        	l = strlen(filter);
		filter[l-2] = '\0';
		fgets(temp, 80, fpin);
		sscanf(temp, "%d", &obs_no);
		fgets(temp, 80, fpin);
		sscanf(temp, "%2d%2d%f", &deg,&min,&fsec);
	    sprintf(tele_ra,"%02d:%02d:%05.2f",deg,min,fsec);
		fgets(temp, 80, fpin);
		sscanf(temp, "%1s%2d%2d%f", sign,&deg,&min,&fsec);
		sprintf(tele_dec,"%s%02d:%02d:%04.1f",sign,deg,min,fsec);
		fgets(temp, 80, fpin);
		sscanf(temp,"%s",lst);
		fgets(temp, 80, fpin);
		sscanf(temp,"%f", &am);
		fclose(fpin);

		if ( strlen(title) == 0 ) {
			title[0] = 'X';
			title[1] = '\0';
			}	
		if ( strlen(filter) == 0 ) {
			filter[0] = 'X';
			filter[1] = '\0';
			}	
		fprintf(fplog,"%s 256 %s %s %s %f %d %s %s %s %f %s  \"%s\"\n",	fname[j],date,time,filter,itime,wframe,telescope,tele_ra,tele_dec,am,lst,title);

		
		append(date);
		append(time);
		append(title);
		append(tele_ra);
		append(tele_dec);
		append(lst);
		
		strcpy(filename,fname[j]);
		append(filename);

		dmax = -32767.;
		dmin = -dmax;
		
		for (i=0; i<MAXCOLS*MAXROWS; i++) {
			if ( (float)sibuf[i] > dmax ) dmax = (float)sibuf[i];
			if ( (float)sibuf[i] < dmin ) dmin = (float)sibuf[i];
			}

		
		if ( swap_bytes ) {
		    swab(sibuf, libuf, MAXCOLS*MAXROWS*2); 
			}

         for (i=0; i<2880; buf[i++] = ' ');

	     hdstr = "SIMPLE  =                    T";
	     p = buf;
	     strncpy(p, hdstr, 30) ;


		hdstr = "BITPIX  =                   16";
		p = buf+80;
	        strncpy(p, hdstr, 30) ;

	     	hdstr = "NAXIS   =                    2";
		p = buf+160;
		strncpy(p, hdstr, 30) ;

	     	hdstr = "NAXIS1  =                  256";
	        p = buf+240;
	        strncpy(p, hdstr, 30) ;

	        hdstr = "NAXIS2  =                  256";
	        p = buf+320;
	        strncpy(p, hdstr, 30) ;
	     
		hdstr = "BSCALE  =               1.0E00";
    	 	p = buf+400;
     	 	strncpy(p, hdstr, 30) ;
		
	     	hdstr = "BZERO   =               0.0E00";
    	 	p = buf+480;
     	 	strncpy(p, hdstr, 30) ;
		
	     	hdstr = "BUNIT   = 'ADU/PIX '          ";
     	 	p = buf+560;
     	 	strncpy(p, hdstr, 30) ;

		hdstr = "BLANK   =               -32767";
		p = buf+640;
            	strncpy(p, hdstr, 30) ;

	 	hdstr = "COMMENT = '                                        ";	
	 	p = buf+720 ;
	 	strncpy(p, hdstr, 30);
	 	p = buf + 731;
		sprintf(data,"Coadds = %d'",wframe);
		l = strlen(data);
		strncpy(p, data, l);

	 	hdstr = "DATAMAX =                                          ";
    	 	p = buf+800;
     	 	strncpy(p, hdstr, 50) ;
     	 	sprintf(&(data[0]), "%9.1f", (float)dmax);
         	p = buf+821;
         	strncpy(p, data, 9) ;  

	 	hdstr = "DATAMIN =                                          ";
    	 	p = buf+880;
     	 	strncpy(p, hdstr, 50) ;
     	 	sprintf(&(data[0]), "%9.1f", (float)dmin);
         	p = buf+901;
         	strncpy(p, data, 9) ;  

         	hdstr = "DATE-OBS= '                   ";
         	p = buf + 960;
         	strncpy(p, hdstr, 30) ;
     	 	p = buf + 971;
	 	i = strlen(date);
	 	strncpy(p, date, i);

	 	hdstr = "OBJECT  = '                   ";
    	 	p = buf+1040;
     	 	strncpy(p, hdstr, 30) ;
	 	p = buf + 1051;
	 	l = strlen(title);
	 	if ( l > 70 ) l = 70;
	 	strncpy(p, title, l);

		hdstr = "INSTRUME= '256x256 '                              ";
         	p = buf+1120;
         	strncpy(p, hdstr, 50) ;

		hdstr = "FILTER  = '      '                                ";
         	p = buf+1200;
         	strncpy(p, hdstr, 30) ;
     	 	p = buf + 1211;
	 	i = strlen(filter);
	 	strncpy(p, filter, i);
		
	 	hdstr = "ITIME   =                     ";
    	 	p = buf+1280;
     	 	strncpy(p, hdstr, 30) ;
	 	sprintf(&(data[0]),"%9.1f", itime);
	 	p = buf + 1301;
	 	strncpy(p, data, 9);

	       hdstr = "IRAFNAME= '                                       ";
    	 	p = buf+1360;
     	 	strncpy(p, hdstr, 50) ;
     	 	p = buf+1371;
	 	i = strlen(filename);
     	 	strncpy(p, filename, i) ;  

         	hdstr = "TIME-OBS= '                   ";
	 	p = buf + 1440;
	 	strncpy(p, hdstr, 30);
	 	p = buf + 1451;
	 	l = strlen(time);
	 	strncpy(p, time, l);
         	
		hdstr = "TELE RA = '                   ";
	 	p = buf + 1520;
	 	strncpy(p, hdstr, 30);
	 	p = buf + 1531;
	 	l = strlen(tele_ra);
	 	strncpy(p, tele_ra, l);

         	hdstr = "TELE DEC= '                   ";
	 	p = buf + 1600;
	 	strncpy(p, hdstr, 30);
	 	p = buf + 1611;
	 	l = strlen(tele_dec);
	 	strncpy(p, tele_dec, l);

         	hdstr = "LST     = '                   ";
	 	p = buf + 1680;
	 	strncpy(p, hdstr, 30);
	 	p = buf + 1691;
	 	l = strlen(lst);
	 	strncpy(p, lst, l);
	 	
         	hdstr = "AIRMASS =                     ";
    	 	p = buf+1760;
     	 	strncpy(p, hdstr, 30) ;
	 	sprintf(&(data[0]),"%9.2f", am);
	 	p = buf + 1771;
	 	strncpy(p, data, 9);

         	hdstr = "TELESCOP= '                   ";
	 	p = buf + 1840;
	 	strncpy(p, hdstr, 30);
	 	p = buf + 1851;
	 	l = strlen(tele);
	 	strncpy(p,tele , l);

	 	hdstr = "END                                               ";
         	p = buf+1920;
     	 	strncpy(p, hdstr, 40) ;
       

         	hdr_len = 2880;
		
	 	fwrite(buf,sizeof(char),NBYTES,fptest);

	 /* Now work on transferring data to disk. */

	 	n_samp = NBLOCKS;
		 		 	
	 	for (i=0; i<n_samp; i++) {
	      		memcpy(buf,((char *)libuf)+i*NBYTES,NBYTES);
	      		fwrite(buf,sizeof(char),NBYTES,fptest);
	     		}
	  	fclose(fptest);		 
	  	} 
	fclose(fplog);	
	fprintf(stderr,"%d files written to disk.\n", j);
	fprintf(stderr,"FITS archiving is complete.\n");
}
			
void strip(str)		/* ctrip CR LFs from trailing strings */
char *str;
{
	int i;

	i = strlen(str);
	str[i-1] = '\0';
	str[i-2] = '\0';
}

void append(str)	/* append single quote to string */
char *str;
{
	int i;

	i = strlen(str);
	while ( i<8 ) 
		str[i++] = ' ';
	str[i] = '\'';
	str[i+1] = '\0';
}

void fixname(fname,n,uname)
char fname[],uname[];
int n;
{
      int i,nlength,next;
      char ext[4];
      extern void itoa();

      nlength = strlen(fname);

      strncpy(uname,fname,nlength);
      uname[nlength] = '.';

      itoa(n,ext,10);
      next = strlen(ext);

      for (i=0; i<3-next; i++) {
	 uname[nlength+i+1] = '0';
	 }
      for (i=nlength+4-next; i<nlength+4; i++) {
	 uname[i] = ext[i-nlength-4+next];
	 }  
      uname[nlength+4] = '\0';

}

void itoa(n,s)	/* convert n to chars in s. K&R p.60. */
char s[];
int n;
{
	int i, sign;
	extern void reverse();

	if (( sign = n) < 0)  	/* record sign */
		n = -n; 	/* make n positive */
	i = 0;
	do {			/* generate digits in reverse order. */
		s[i++] = n % 10 + '0';		/* get next digit */
	} while ( ( n /= 10 ) > 0);		/* delete it */
	if ( sign < 0 )
		s[i++] = '-';
	s[i] = '\0';
	reverse(s);
}

void reverse(s)	/* reverse string s in place. K&R p.59 */
char s[];
{
	int c, i, j;

	for (i = 0, j =strlen(s) - 1; i < j; i++, j--) {
		c = s[i];
		s[i] = s[j];
		s[j] = c;
		}
}


void float_swap(array)
float *array;
{
	int j;
	char temp[4];
	union {			/* This hack is used for byte-swapping. */
 	char t_num[4];
	    float num;
		} nums;

	for (j=0; j<MAXROWS*MAXCOLS; j++) {
		nums.num = *(array+j);
		temp[0] = nums.t_num[0];
		temp[1] = nums.t_num[1];
		temp[2] = nums.t_num[2];
		temp[3] = nums.t_num[3];
    	nums.t_num[0] = temp[3];
		nums.t_num[1] = temp[2];
		nums.t_num[2] = temp[1];
		nums.t_num[3] = temp[0];
	 	*(array+j)= (float)nums.num;
		}
}