#include <stdio.h>
#define n 9999
FILE *fp2,*fp3;
char a[130],b[130];
main(ac,av)
int ac; char *av[];
{
  int i,j,k;
  char f1[60],f2[60],f3[60];
  int  p[n];
  float x[n],y[n],u[n],v[n],r[n],x1,y1,t,s;
  
  if(ac<2){
    printf("\n\t *****  compare coo & als******\n");
    exit(0);
  }
  strcpy(f1,av[1]);
  k=nindex(f1,"."); if(k>0)f1[k]=0;
  strcpy(f2,f1); strcpy(f3,f1); 
  strcat(f2,".coo"); strcat(f3,".als"); 
  fp2=fopen(f2,"r"); if(fp2==0){ printf("%s not found!\n",f2); exit(0); }
  fgets(a,130,fp2); fgets(a,130,fp2); fgets(a,130,fp2);
  for(i=0;i<n;i++)p[i]=0;
l05:
  fgets(a,130,fp2); if(feof(fp2))goto l06;
  sscanf(a,"%d %f %f",&j,&x1,&y1);
  j--;
  x[j]=x1; y[j]=y1;
  goto l05;
l06:
  fclose(fp2);
  
  fp3=fopen(f3,"r"); if(fp3==0){ printf("%s not found!\n",f3); exit(0); }
  fgets(b,130,fp3); fgets(b,130,fp3); fgets(b,130,fp3);
l10:
  fgets(b,130,fp3);
  if(feof(fp3)!=0)goto l20;
  sscanf(b,"%d %f %f %f",&j,&x1,&y1,&s);
  j--;
  p[j]=1; u[j]=x1; v[j]=y1; r[j]=s;
  goto l10; 
l20:
  fclose(fp3);
  fp2=fopen("coo2als.dat","w");
  for(i=0;i<n;i++){
    if(p[i]){
      t=x[i]-u[i]; s=y[i]-v[i]; t=t*t+s*s;
      if(t>30.){
             printf("%d %f %f_%f %f_%f %f\n",i+1,x[i],y[i],u[i],v[i],r[i],t);
        fprintf(fp2,"%d %f %f %f %f %f %f\n",i+1,x[i],y[i],u[i],v[i],r[i],t);
      }
    }
  }
  fclose(fp2);
}