#define	TEMP_FILE	"\\temp.ray"
#define THROWAWAY	"\\throwa.way"

/* OPRA Data Acquisition Program */
#define ZOFFSET  20.0
#define	ZLOW	-20.0
#define ZHIGH 	 60.0
#define ZSCALE	16./(ZHIGH - ZLOW)	

#define NOTAPE

#include        <math.h>
#include        <conio.h>       /* for getch, kbhit */
#include        <stdio.h>       /* for malloc */
#include        <dos.h>
#include        <signal.h>
#include	<time.h>
#include	<unistd.h>
#include	<allegro.h>

#define EXTERN
#include        "..\include\globals.h"


#define         PPI     0
#define         RHI     1

int nodisplay = 1, radiate = 0;
double DACVOLTAGE = -4.2;
DISPLAY *remdisp;

main(argc,argv)
char    *argv[];   
int     argc;
   {
   char    c, fname[20];
   DWELL   *dwell;   
   RADAR   *radar;
   CONFIG  *config;
   PIRAQ   *piraq;
   DISPLAY *disp;
   PACKET  *packet;
   float   *prods, *prptr;
   double  zsave;
   int     rayv, gate,oldbutton=0, ascptr;
   int     type=PPI,dispnum='1',dchange=0,pmode,singlefilter;   
   int 	   togcnt = 0, counter = 0, runcnt = 0;


   int       gates_noclut, gates_clut;
   int       noafc = 0;
   double    dopp = 0.0;
   FILE	*tmpfp, *thrwfp;


int i, j, test=1,sync,oldtime,bidelay,leftover,xmit,event;
unsigned int    itime;
unsigned short  subsec;
double crap,vconst;
struct tm *ascz;

   
ascptr = 0;

   if(argc > 1)
   	{
	if(!strncmp("-G", argv[1]))nodisplay = 0;
	if(!strncmp("-g", argv[1]))nodisplay = 0;
	if(!strncmp("-A", argv[1]))nodisplay = 2;
	if(!strncmp("-a", argv[1]))nodisplay = 2;
   	}
		
   gates_noclut = gates_clut = 200;
   
   signal(SIGFPE,SIG_IGN);

   if((prods = (float *)malloc(sizeof(float) * 1200 * 8)) == NULL)
      {puts("not enough memory"); exit(1);}

   if((dwell = (DWELL *)malloc(sizeof(DWELL))) == NULL)
      {puts("not enough memory"); exit(1);}

   if((disp = (DISPLAY *)malloc(sizeof(DISPLAY))) == NULL)
      {puts("not enough memory"); exit(1);}

   if((radar = (RADAR *)malloc(sizeof(RADAR))) == NULL)
      {puts("not enough memory"); exit(1);}

   if((config = (CONFIG *)malloc(sizeof(CONFIG))) == NULL)
      {puts("not enough memory"); exit(1);}

   if((piraq = (PIRAQ *)malloc(sizeof(PIRAQ))) == NULL)
      {puts("not enough memory"); exit(1);}

   if((packet = (PACKET *)malloc(sizeof(PACKET))) == NULL)
      {puts("not enough memory"); exit(1);}

   /* read in the dsp configuration file */
   readconfig("",config);
   singlefilter = config->clutterfilter;

   /* set the pointers in the piraq structure based on the configuration */
   setpiraqptr(config,piraq);

   /* load the dsp code into the dual port ram */
   initsystem(config);
   cleardpram(config,piraq);
   load(config,piraq);

   /* put all configuration parameters into the piraq's dual port ram */
   setdpram(config,piraq);
   
   /* read in the radar configuration file */
   readradar("",radar);
   
   /* read in the display configuration file */
   readdisplay("",disp);

   if(!timerset(config)) {printf("EXEC: Error setting timer\n"); exit(1);}      /* set up the range gate timer */
   if(!start(config))  {printf("EXEC: Could not start PIRAQ\n"); exit(1);}        /* start the range gate timer and DSP */

   /* initialize functions such as gps,serial,display,tape,disk after this */

   /* get a dwell structure header filled so the display can initialize */
   /* by using the header of a bogus dwell */
   fillheader(dwell,piraq,config);  /* fill the dwell header */
   
   /* initialize SCSI tape drive */
   
#ifndef NOTAPE
   init_tape();
#endif
   
   /* Force Vert Pointing mode */
   disp->type |=2; /* Put on the grid bit */
   disp->type |=4; /* Put on the VP bit */
   
   if(!nodisplay)
   	{
   	if(!init_display(radar,dwell,disp))        
      		{
		printf("could not open display lookup file\n"); 
		closeall(); 
		exit(1);
		}
	}
   if(nodisplay == 1)
   	printf("EXEC running, hit Q to exit\n");
   
   /* calculate this constant just once */
   vconst = 128.0 * (2.0 * radar->frequency * 2.0 * M_PI * dwell->header.prt) / C / M_PI;
   
   /* Ensure that wrtollg.trg does NOT exist */
   if(exists("\\wrtoggle.trg"))
    	{
	unlink("\\wrtoggle.trg");
	}
    
   
   pmode = *piraq->dataformat;
   zsave = disp->rgpp/(3*1.414);
/*  GRG 4/19/98 disp->rgpp =zsave; */


   do
    {
    /* Test for disk write toggle flag */
    if(exists("\\wrtoggle.trg"))  /* Test for existence of the trigger file */
    	{
	unlink("\\wrtoggle.trg"); /* Delete the toggle trigger file */
	
	if(togcnt == 0) /* Don't let it trigger more than once */
		{
		disp->recording = (disp->recording == 2) ? 0 : 2;        
		togcnt++;
		}
	}
    else	/* If trigger file doesn't exist, zero the toggle count */
    	{
	togcnt = 0;
	}
    
      /* if c = 0 and there is already is another character, then */ 
      c = kbhit() ? toupper(getch()) : 1;
      
      /* a function key was hit */
      if(c == 0) 
	 {
	 c = getch() - 0x3B;    /* parse the function key */
	 if(c < 6  || (c < 8 && dwell->header.dataformat == DATA_POL1))
	    {
	    disp->displayparm = c;  
	    ascptr = (int)c;
	    if(!nodisplay)selections(disp);
	    }
	    /* This line eliminates crash if tape selected with F10 */
	    if(c==9)c=8;
	 switch(c) /* recording to tape (3) or disk (2) */
	    {   case 8:     
			disp->recording = (disp->recording == 2) ? 0 : 2;        
			write_disk(1,dwell,radar,config->outfilename);  /* write RADAR to disk */
			break;
		case 9:     
			disp->recording = (disp->recording == 3) ? 0 : 3;        
			write_tape(1,dwell,radar);  /* write a RADAR struct to tape */
			break;
	    }
	 c = 1;        /* reset c to benign value */
	 }
      
      if(c == 'G')   savescr(dwell->header.time);
      
      if(c == '1')   disp->rgpp *= 2;
      if(c == '2')   disp->rgpp /= 2;
      
/* Remove these options for VP */
/*      if(c == '3')   disp->type ^= 2;*/   /* flip the grid bit */
/*      if(c == '4')   disp->type ^= 4;*/   /* vertical profile bit */


      /* dac output voltage */
      if(c == '>')     dac(DACVOLTAGE+=0.05);
      if(c == '<')     dac(DACVOLTAGE-=0.05);
      
      /* change gate spacing on the fly */
      if (c == '[' || c == ']')
	 {
	 switch(disp->recording)  /* write a closing radar structure*/
	   {
	   case 2:      write_disk(1, dwell, radar, config->outfilename);  break;
	   case 3:      write_tape(1, dwell, radar);  break;
	   default:     break;
	   }
	 if (c == ']')   config->pulsewidth +=2;
	 if (c == '[')   config->pulsewidth -=2;
	 if (config->pulsewidth < 2)   config->pulsewidth = 2;
	 /* the following code will set up a scan for full amb range */
	 if (config->pulsewidth > 92 && config->clutterfilter == 0) 
		{
		config->pulsewidth = 94; /* set up xxxm gates */
		config->gates = gates_noclut;
		}
	 if (config->pulsewidth > 62 && config->clutterfilter == 1) 
		{
		config->pulsewidth = 64; /* set up xxxm gates */
		config->gates = gates_noclut;
		}
	 switch (disp->recording)  /* write a closing radar structure*/
	   {
	   case 2:      write_disk(1, dwell, radar, config->outfilename);  break;
	   case 3:      write_tape(1, dwell, radar);  break;
	   default:     break;
	   }
	 stop(config);
	 /* put all configuration parameters into the piraq's dual port ram */
	 setdpram(config, piraq);
	 timerset(config);   /* set up the range gate timer */
	 if (!start(config))  exit(1); /* start the range gate timer */
					/* get a dwell structure header */
					/* filled so the display can */
					/* initialize by using the header */
					/* of a bogus dwell */
	 fillheader(dwell, piraq, config);     /* fill the dwell header */
//         init_disp_scales(radar, dwell, disp);
	 }
      if (c == 'S' || c == 'D')
	 {
	 switch (disp->recording)
	   {
	   case 2:      write_disk(1, dwell, radar, config->outfilename);  break;
	   case 3:      write_tape(1, dwell, radar);  break;
	   default:     break;
	   }
	 if (c == 'D')   config->hits *=2;
	 if (c == 'S')   config->hits /=2;
	 if (config->hits < 32)     config->hits = 32;
	 if (config->hits > 1024)   config->hits = 1024;
	 setdpram(config, piraq);
	 }
      
      if (c == ',' || c == '.')
	 {
	 switch (disp->recording)
	   {
	   case 2:      write_disk(1, dwell, radar, config->outfilename);break;
	   case 3:      write_tape(1, dwell, radar); break;
	   default:     break;
	   }
	 if (c == ',')   config->gates -=10;
	 if (c == '.')   config->gates +=10;
	 setdpram(config, piraq);
	 }


      if (c == 'M' || c == 'N')
	 {
	 if (c == 'M') config->tsgate +=1;
	 if (c == 'N') config->tsgate -=1;
	 if (config->tsgate < 1)   config->tsgate = 1;
	 if (config->tsgate > config->gates -1)  config->tsgate = config->gates -1;
	 config->timeseries = 1;
	 setdpram(config, piraq);
	 }
      if (c == 'H' || c == 'J')
	 {
	 if (c == 'J')   config->tsgate +=10;
	 if (c == 'H')   config->tsgate -=10;
	 if (config->tsgate < 1)   config->tsgate = 1;
	 if (config->tsgate > config->gates -1)  config->tsgate = config->gates -1;
	 config->timeseries = 1;
	 setdpram(config, piraq);
	 }

      
      if(c == 'F')   disp->fakeangles = 1 - disp->fakeangles;
      if(c == 'T')   {disp->threshold += 0.05; if(disp->threshold > .99) disp->threshold=.99;}
      if(c == 'U')   disp->threshold = 0.0;
      
      /*chang*/
      if (c == 'C')   
	{
	switch (disp->recording)  /* write a closing radar structure*/
	   {
	   case 2:      write_disk(1, dwell, radar, config->outfilename);break;
	   case 3:      write_tape(1, dwell, radar);break;
	   default:     break;
	   }
	config->clutterfilter ^= 1;
	if (config->clutterfilter == 1)   config->gates = gates_clut;
	if (config->clutterfilter == 0)   config->gates = gates_noclut;
	switch (disp->recording)  /* write a closing radar structure*/
	   {
	   case 2:      write_disk(1, dwell, radar, config->outfilename); break;
	   case 3:      write_tape(1, dwell, radar); break;
	   default:     break;
	   }
	stop(config);
	/* put all configuration parameters into the piraq's dual port ram */
	setdpram(config, piraq);
	timerset(config);   /* set up the range gate timer */
	if (!start(config))  exit(1); /* start the range gate timer */
					/* get a dwell structure header */
					/* filled so the display can */
					/* initialize by using the header */
					/* of a bogus dwell */
	fillheader(dwell, piraq, config);     /* fill the dwell header */
//         init_disp_scales(radar, dwell, disp);
	}

      if (c == 'A' && config->pcorrect)        
	 {
	 if(!nodisplay)close_display();
	 afclock(config, piraq);
	 
	 if (!nodisplay && !init_display(radar, dwell, disp))        
	    {
	    closeall();
	    printf("could not open display lookup file\n");
	    exit(1);
	    }
	 }

/* some of these parameters are no longer used */      
/* Remove option for VP radars
      if(c == 'P' && disp->type != PPI)              {disp->type = PPI; dchange = 1; c = dispnum;}
      if(c == 'R' && disp->type != RHI)              {disp->type = RHI; dchange = 1; c = dispnum;}
*/
/* DEBUG This permanently kills the TX...can't restart it 
   The reason is: there is no trigger to the DSP and so the wait cycle
   stalls it out.  Need to set flags to bypass wait when it's in stop
   mode 
   
	if(c == '!')
		{
		if(radiate)
			{
			radiate = 0;
			Prf(8.0E+06/(float)config->prt, 1.0, 0);
			}
		else
			{
			radiate = 1;
			Prf(8.0E+06/(float)config->prt, 1.0, 1);
			}
		}
*/
      if(c >= '0' && c <= '9' && c != dispnum) dchange = 1;
      
      if (c == 'O')
	 {
	 config->pcorrect = 1 - config->pcorrect;
	 setdpram(config, piraq);
	 }
      if(c == 'K') /* Kill OPRA */
      	{
	closeall();
	exit(2);
	}


# ifdef notdef
      /* if left mouse button was just clicked */
      gate = disp->mouserange / (C * 0.5 * dwell->header.rcvr_pulsewidth);
      if(disp->mousebutton & ~oldbutton & LEFTBUTTON)        
	if(gate < config->gates)
	   {
	   config->tsgate = gate; 
	   setdpram(config,piraq);
	   }
      if(disp->mousebutton & ~oldbutton & RIGHTBUTTON)        
	 {
	 config->timeseries ^= 1;
	 setdpram(config,piraq);
	 }
      /* update the mouse button memory */
      oldbutton = disp->mousebutton;      
# endif
      
      if(!wait(config,piraq))  break;
	
      if(config->pcorrect && config->gate0mode)    afc(config,piraq);
     
      filldwell(dwell,piraq,config); /* create a dwell structure in full */

      products(dwell,radar,prods);  /* compute scientific radar parameters */
/* 
This section writes a data ray for remote and a throwaway file
to let remote know that it is done writing.   Remote will then read the
data file and delete the throwaway.
*/
      
      tmpfp = fopen(TEMP_FILE, "wb");
      fwrite(radar,radar->recordlen,1,tmpfp);
      fwrite(dwell,dwell->header.recordlen,1,tmpfp);
      fwrite(disp,sizeof(DISPLAY),1,tmpfp);
      fclose(tmpfp);
      
      thrwfp = fopen(THROWAWAY, "wb");
      fwrite(radar,radar->recordlen,1,thrwfp);
      fclose(thrwfp);
/* End of remote communication section */

      if(disp->recording == 2)     write_disk(2,dwell,radar,config->outfilename);  /* write to disk */
      if(disp->recording == 3)     write_tape(2,dwell,radar);  /* write to tape */
      
	if(nodisplay == 0)
      		display(dwell,disp,prods);  /* write one dwell to display */
	if(nodisplay == 2) 
      		zprint(prods, dwell, disp);	
		
    /* Test for kill flag (existence of "killopra.fil") */
    if(exists("\\killopra.fil"))  /* Test for existence of the trigger file */
    	{
	unlink("\\killopra.fil"); /* Delete the kill trigger file */
	closeall();	
	exit(2);
	}
    /* Test for quit flag (existence of "qexec.tog") */
    if(exists("\\qexec.tog"))  /* Test for existence of the trigger file */
    	{
	unlink("\\qexec.tog"); /* Delete the toggle trigger file */
	closeall();	
	exit(0);
	}
    }while(c != 'Q');
   closeall();
   exit(0);
   }

void closeall(void)
   {
   Prf(1000,1,0);
   if(!nodisplay)close_display();
   printf("EXEC terminated\n");
   close_disk();
   close_tape();
   savev();
   }

zprint(float *prods, DWELL *dwell, DISPLAY *disp)
{
int i, rayv, ascptr;
struct tm *ascz;
float *prptr;
char curchar, colchars[16] = {'a','b','c','d','0','1','2','3','4','5','6','7',
			'8','9','A','B'};

ascz = gmtime((time_t *)&dwell->header.time);
prptr = prods;
ascptr = 4;
printf("%02d%02d%02d:", ascz->tm_hour, ascz->tm_min, ascz->tm_sec);
for(i=0; i < 70; i++)
	{
	rayv = ((int)((*(prptr + ascptr) + ZOFFSET)*ZSCALE))&0xF;
	curchar = colchars[rayv];
	if(*(prptr + ascptr) > ZHIGH)curchar = '^';
	if(*(prptr + ascptr) < ZLOW)curchar = '.';
	if(*(prptr + 2) < disp->threshold)curchar = '-';
	printf("%1c", curchar);
	prptr += PRODS_ELEMENTS;
	} 
printf("\n");
if(ascz->tm_sec%10 == 0)
	printf("  km-> 0----^---1.5---^---3.0---^---4.5---^---6.0---^---7.5---^---9.0 %s\n", 
		disp->recording != 2? "Rec off": "To disk");
}