#include <stdio.h>
FILE *fp;
char a[80],b1[5],b2[5];
main(int ac, char **av)
{
  int i,j,k;
  if(ac<2){ 
    printf("\n\t******* process dome,twilight FLAT ********");
    printf("\n\t                                  110930");
    printf("\n\tUsage: flat !\n");
    printf("\n\tcp /vega2/rhbin/bias*.fit .");
    printf("\n\td4 *.fits");
    printf("\n\tjob_e4b");
    printf("\n\tls e*_?.fit > ?.f");
    printf("\n\tf4 ?.f fdateobje_?.fit");
    printf("\n\trm 1.1 *.tmp d*.fit\n\n");
    exit(0);
  }
  strcpy(a,"cp /vega2/rhbin/bias*.fit .");
  printf("%s\n",a); system(a);
  strcpy(a,"d4 *.fits");
  printf("%s\n",a); system(a);
  strcpy(a,"job_e4b");
  printf("%s\n",a); system(a);
  strcpy(a,"ls e*_1.fit >1.f");
  printf("%s\n",a); system(a);
  strcpy(a,"ls e*_2.fit >2.f");
  printf("%s\n",a); system(a);
  strcpy(a,"ls e*_3.fit >3.f");
  printf("%s\n",a); system(a);
  strcpy(a,"ls e*_4.fit >4.f");
  printf("%s\n",a); system(a);
  
  system("dirss >1.1");
  fp=fopen("1.1","r");
  fgets(a,80,fp); fgets(a,80,fp); fgets(a,80,fp); 
  for(i=0;i<4;i++)b1[i]=a[6+i]; b1[i]=0;
  for(i=0;i<4;i++)b2[i]=a[45+i]+32; b2[i]=0;

  sprintf(a,"f4 1.f f%s%s_1.fit",b1,b2);
  printf("%s\n",a); system(a);
  sprintf(a,"f4 2.f f%s%s_2.fit",b1,b2);
  printf("%s\n",a); system(a);
  sprintf(a,"f4 3.f f%s%s_3.fit",b1,b2);
  printf("%s\n",a); system(a);
  sprintf(a,"f4 4.f f%s%s_4.fit",b1,b2);
  printf("%s\n",a); system(a);

  system("rm d*.fit 1.1");
}