#include <stdio.h>
#define n 9999
#define m  999         // readin coo2als.dat
int cn[m],nc;
float cx[m],cy[m];
FILE *fp0,*fp2,*fp3,*fp4;
char a[130];
main(ac,av)
int ac; char *av[];
{
  int i,j,k,kk;
  float x,y,z,v;
  float c1,c2;
  char f1[60],f2[60],f3[60],f4[60];
  if(ac<2){ printf("\n\t ********get cosmos from als, reorder coo ******\n");
            printf("\n\t cosmos=miss.dat+coo2als.dat+uals(<10)\n");
            printf("\n\tUsage: getcosmos uedo\n");
            exit(0);
          }
  strcpy(f1,av[1]);
  k=nindex(f1,"."); if(k>0)f1[k]=0;
  strcpy(f2,f1); strcpy(f3,f1); strcpy(f4,"cosmos.dat");
  strcat(f2,".coo");  strcat(f3,".als"); 
  fp3=fopen(f3,"r"); if(fp3==0){ printf("%s not found!\n",f3); exit(0); }
  fp2=fopen(f2,"w");
  fp4=fopen(f4,"w");
  fp0=fopen("sort.par","r");
  fgets(a,130,fp0); sscanf(a,"%f",&c1); 
  fgets(a,130,fp0); sscanf(a,"%f",&c2); 
  fclose(fp0);
  k=0;
  fp0=fopen("miss.dat","r");
  if(fp0){
l03:
    fgets(a,130,fp0); if(feof(fp0))goto l04;
    fputs(a,fp4); k++;
    goto l03;
l04:
    fclose(fp0);    
  }
  fp0=fopen("coo2als.dat","r"); i=0;
  if(fp0){
l05:
    fgets(a,130,fp0); if(feof(fp0))goto l06;
    sscanf(a,"%d %f %f",&cn[i],&cy[i],&cx[i]);
    fputs(a,fp4);
    i++; k++;
    goto l05;
l06:    
    fclose(fp0);
  }
  nc=i; kk=0; 
  for(i=0;i<3;i++){ fgets(a,130,fp3);  fputs(a,fp2); }
l10:
  fgets(a,130,fp3); if(feof(fp3))goto l20;
  sscanf(a,"%d %f %f %f %f %f %f %f %f",&j,&y,&x,&z,&z,&z,&z,&z,&v);
  if(v<c1){ fprintf(fp4,"%5d %9.3f %9.3f\n",j,y,x); k++; }
  else {
    for(i=0;i<nc;i++)if(j==cn[i])goto l10;
    fputs(a,fp2); kk++;
  }
  goto l10; 
l20:
  fclose(fp3);
  printf("prodced  cosmos.dat  %d\n",k);
  printf("prodced  %s %d\n",f2,kk);
}