#include <stdio.h>
char f1[50],a[72][80],b[72][80];
int i,j,k,ipos,jpos;
FILE *fp,*fp1;
main(ac,av)
int ac; char *av[];
{
  if(ac<3){
    printf("\n\tUsage: pchead refence_fits batch_file\n");
    printf("\n\t MUST BE CAREFUL, same sky_area MUST!!\n");
    printf("\n\tmodify fits_head RA. Dec, according a refence fits in batch\n");
    printf("\n\t no matter of file size\n");
    printf("\n\t\t\t jiang, 61120\n\n");
    exit(0);
  }
  fp=fopen(av[1],"rb");
  if(fp==0){
    printf("\n\trefence_fits_file not found!\n");
    exit(0);
  }
  fread(a,72,80,fp); 
  ipos=indexpos(a,"RA      ",72);
  fclose(fp);
  
  fp=fopen(av[2],"r");
  if(fp==0){
    printf("\n\tbatch_file not found!\n");
    exit(0);
  }
  k=0;
l10:
  k++;
  fgets(f1,50,fp);
  if(feof(fp)!=0){ printf("\nok !\n"); exit(0); }
  i=strlen(f1)-1; f1[i]=0;
  fp1=fopen(f1,"r+b");
  if(fp1==0){
    printf("\n\tfits_file not found!\n");
    goto l10;
  }
  fread(b,72,80,fp1);
  jpos=indexpos(b,"RA      ",72);
  for(i=10;i<30;i++)b[jpos][i]=a[ipos][i];
  for(i=10;i<30;i++)b[jpos+1][i]=a[ipos+1][i];
  fseek(fp1,0,0);
  fwrite(b,72,80,fp1);
  fclose(fp1);
  printf("%5d: %s\n",k,f1);
  goto l10; 
}