#include <stdio.h>
char head[72][80];
FILE *fp,*fp1;
main(ac,av)
int ac; char **av;
{
  int i,j,k;
  float ep;
  char a[40];
  if(ac<2){ printf("\n\tUsage: addep1 t0001\n\n"); exit(0); }
  sprintf(a,"ls %s/*/*.fit* >1.1",av[1]);
  system(a);
  fp1=fopen("1.1","r");
l10:
  fgets(a,40,fp1); if(feof(fp1))goto l20;
  k=strlen(a)-1; a[k]=0;
  printf("%s  ",a);
  fp=fopen(a,"rb"); if(fp==0){ printf("file not found !\n"); goto l20; }
  fread(head,72,80,fp); fclose(fp);
  k=indexpos(head,"EPOCH   ",72);
  if(k!=72){ printf("already done !\n"); goto l10; }
  k=indexpos(head,"RA      ",72);
  sscanf(&head[k][50],"%f",&ep);
  if(ep<1993. || ep >1995.){ printf("undo !\n"); goto l10; }
  k=indexpos(head,"PZERO   ",72);
  if(k==72)k=indexpos(head,"RECORDS ",72);
  if(k==72)k=indexpos(head,"        ",72);
  if(k==72)k=indexpos(head,"COMMENT ",72);
  if(k==72)goto l10;
  sprintf(&head[k][0],"EPOCH   =   %18.1f /",ep); 
  for(i=32;i<70;i++)head[k][i]=32;
  fp=fopen(a,"rb+"); fwrite(head,72,80,fp); fclose(fp);
  printf("\n");
  goto l10;
l20:
  fclose(fp1);
  system("rm 1.1");
}