#include <stdio.h>
char   b[146],c1[14],c2[14],c3[14];
float a[4],d[4],d1,d2,x,y;
double r1,r2,r3;
FILE *fp;
int g[180];
//  0  2
//  1  3
main(int ac, char **av)
{
  int i,j,k;
  if(ac<2){ printf("\n\t**** status fits.pos9 ******\n");
            printf("\n\tstatpos9 !  ( status DEC)\n");
            printf("\n\tstatpos9 10 -10 ( status RA1 RA2)\n");
            exit(0);
          }
  fp=fopen("/line3/uband-data/cat/fits.pos9","r");
  for(i=0;i<180;i++)g[i]=0;  r1=99; r2=-99;
  if(ac>2){ sscanf(av[1],"%f",&d1); sscanf(av[2],"%f",&d2); goto l51; }
l10:
  fgets(b,146,fp); if(feof(fp))goto l20;
  sscanf(&b[46],"%f %f %f %f %f %f %f %f",
   &a[0],&d[0],&a[1],&d[1],&a[3],&d[3],&a[2],&d[2]); 
  y=0; for(i=0;i<4;i++)y+=d[i]; y*=.25; j=y+90.5; g[j]++;
  goto l10;
l20:
  fclose(fp);  j=0;
  for(i=0;i<180;i++)if(g[i]){ 
    j++; printf("%3d:%4d  ",i-90,g[i]); if(j%5==0)printf("\n"); }
  exit(0);
l51:
  k=0;
  printf("%f %f\n",d1,d2);
l50:
  fgets(b,146,fp); if(feof(fp))goto l60;
  sscanf(&b[46],"%f %f %f %f %f %f %f %f",
   &a[0],&d[0],&a[1],&d[1],&a[3],&d[3],&a[2],&d[2]); 
  y=0; for(i=0;i<4;i++)y+=d[i]; y*=.25; 
  if(y>d1 || y<d2)goto l50;
  k++;
  for(i=0;i<4;i++)if(a[i]<1)goto l50;
  x=0; for(i=0;i<4;i++)x+=a[i]; x*=.25;
  if(x<8.)x+=24.;
//   printf("%8.2f",x);
  if(r1>x)r1=x; if(r2<x)r2=x;
  goto l50;
l60:
  r3=(r1+r2)/2.;  if(r3>24.)r3-=24.;
  printf("%d:  %lf %lf\n",k,r1,r2);  
  if(r2>24.)r2-=24.;
  toms2(r1,c1,1); toms2(r2,c2,1);  toms2(r3,c3,1);
  printf("%s %s   %s\n",c2,c1,c3);
  
  exit(0);
}   
  

toms2(aa,c,k)
double aa;           // input
char c[14];          // output
int k;               // input
{
  int  i,j;
  double a,b,x;
  a=aa;
  if(k==1){ if(a>24.)a-=24.; if(a<0.)a+=24.; }
  c[0]=32; if(a<0.){ a=-a; c[0]='-'; }
  i=a; b=(a-i)*60.;
  j=b; x=(b-j)*60.;
  if(j==60){ j=0; i++; }
  if(k==1)sprintf(&c[1],"%2.2d:%2.2d:%6.3f",i,j,x);
  if(k==0)sprintf(&c[1],"%2.2d:%2.2d:%5.2f",i,j,x);
  if(c[7]==32)c[7]='0';   if(c[8]==32)c[8]='0';
  if(c[7]=='6'){
    c[7]='0'; j++; sprintf(&c[4],"%2.2d",j);
    if(c[4]=='6'){
      c[4]='0'; i++; sprintf(&c[1],"%2.2d",i);
    }
    c[3]=':'; c[6]=':';
  }
}
