#include <stdio.h>
#define n1 4096
#define n2 4032
FILE *fp2,*fp3,*fp4;
float u[n2][n1];
char head[72][80];
char a[130];
main(ac,av)
int ac; char *av[];
{
  int i,j,k,ip;
  int ii,jj,kk;
  float x,y;
  char f0[60],f1[60],f2[60],f3[60],f4[60];
  if(ac<2){ printf("\n\t ********putcoo ******\n");
            printf("\n\tUsage: putcoo uedo cosmos.dat or uedo.als,coo\n");
            exit(0);
          }
  strcpy(f1,av[1]);
  k=nindex(f1,"."); if(k>0)f1[k]=0;
  strcpy(f2,f1);  strcpy(f4,f1); 
  strcat(f2,".fit"); strcpy(f3,av[2]); strcat(f4,f3);
  k=nindex(f4,".");
  ip=0; if(f4[k+1]=='d' && f4[k+2]=='a' && f4[k+3]=='t')ip=1;
  if(ip==1)f4[k]=0;
  else { strcpy(f0,&f4[k+1]); strcpy(f4,f1); strcat(f4,f0); }
  strcat(f4,".fit");
  fp2=fopen(f2,"rb"); if(fp2==0){ printf("%s not found!\n",f2); exit(0); }
  fp3=fopen(f3,"r"); if(fp3==0){ printf("%s not found!\n",f3); exit(0); }

  fread(head,72,80,fp2);
  fread(u,n1*n2,4,fp2);
  fclose(fp2);
  swap4(u,n1*n2*4);
  if(ip==0)for(i=0;i<3;i++)fgets(a,130,fp3);
l10:
  fgets(a,130,fp3); if(feof(fp3))goto l20;
  sscanf(a,"%d %f %f",&j,&x,&y);
  i=y-0.5; 
  j=x-0.5;
  for(ii=0;ii<=9;ii++)for(jj=0;jj<=9;jj++){
    x=ii*ii+jj*jj; x=sqrt(x); kk=x+0.5;
    if(kk==9)u[i+ii][j+jj]=u[i-ii][j+jj]=u[i+ii][j-jj]=u[i-ii][j-jj]=50000;
  } 
  goto l10; 
l20:
  fclose(fp3);
  swap4(u,n1*n2*4);
  fp4=fopen(f4,"wb"); fwrite(head,72,80,fp4); 
  fwrite(u,n1*n2,4,fp4); fclose(fp4);
  printf("%s produced!\n",f4);
}