uk2iod.c

#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <math.h>
#include <ctype.h>

#define LIM 128

int fgetline(FILE *file,char *s,int lim);
int find_satno(char *desig0)
{
  FILE *file;
  int satno=99999,status;
  char desig[16];
  char *env,filename[LIM];

  env=getenv("ST_DATADIR");
  sprintf(filename,"%s/data/desig.txt",env);
  file=fopen(filename,"r");
  if (file==NULL) {
    fprintf(stderr,"Designation file not found!\n");
    exit(0);
  }
  while (!feof(file)) {
    status=fscanf(file,"%d %s",&satno,desig);
    if (strcmp(desig,desig0)==0) 
      break;
  }
  fclose(file);

  return satno;
}


int main(int argc,char *argv[])
{
  FILE *file;
  char line[LIM];
  int intidy,intido,piece,site,year,month,day,hour,min,sec,fsec,satno;
  char desig[16],pdesig[16];
  int format,epoch;
  float csec,cang,x;
  int rah,ram,rafm,ded,dem,defm;
  float tm,tx,am,ax;
  char sign;

  file=fopen(argv[1],"r");
  while (fgetline(file,line,LIM)>0) {
    // Skip wrong lines
    if (!isdigit(line[0])) 
      continue;
    // Skip short lines
    if (strlen(line)<55)
      continue;

    // Scan line
    sscanf(line,"%02d%03d%02d%04d%02d%02d%02d%02d%02d%02d%03d",&intidy,&intido,
	   &piece,
	   &site,
	   &year,
	   &month,
	   &day,
	   &hour,
	   &min,
	   &sec,
	   &fsec);
    sscanf(line+27,"%f",&csec);
    sscanf(line+33,"%1d",&format);
    if (format==2) {
      sscanf(line+34,"%02d%02d%d",&rah,&ram,&rafm);
      sscanf(line+42,"%c%02d%02d%d",&sign,&ded,&dem,&defm);
    } else if (format==3) {
      sscanf(line+34,"%02d%02d%d",&rah,&ram,&rafm);
      sscanf(line+42,"%c%02d%02d",&sign,&ded,&dem);
    }
    sscanf(line+50,"%f",&cang);
    sscanf(line+54,"%d",&epoch);

    // Year switch
    if (year>50)
      year+=1900;
    else
      year+=2000;

    // Format designation
    if (piece<26) {
      sprintf(desig,"%02d %03d%c",intidy,intido,piece+'A'-1);
      sprintf(pdesig,"%02d%03d%c",intidy,intido,piece+'A'-1);
    } else {
      fprintf(stderr,"Failed to understand designation!\n");
      fprintf(stderr,"%s\n",line);
      continue;
    }

    // Test data format
    if (format==3) {
      x=dem*0.6;
      dem=(int) floor(x);
      defm=(int) (100.0*(x-dem));
    } else if (format!=2) {
      fprintf(stderr,"Angle format not implemented!\n");
      fprintf(stderr,"%s\n",line);
      continue;
    }

    // Fractional seconds
    if (fsec<10)
      fsec*=100;
    else if (fsec<100)
      fsec*=10;

    // Time accuracy
    if (csec<10)
      csec*=0.1;
    else if (csec<100)
      csec*=0.01;
    tx=floor(log10(csec))+8;
    tm=floor(csec/pow(10.0,tx-8));

    // angle accuracy
    if (cang<10)
      cang*=1;
    else if (cang<100)
      cang*=0.1;
    ax=floor(log10(cang))+8;
    am=floor(cang/pow(10.0,ax-8));

    // Fractional RA
    if (rafm<10)
      rafm*=100;
    else if (rafm<100)
      rafm*=10;

    // Fractional DE
    if (defm<10)
      defm*=10;
    else if (defm<100)
      defm*=1;

    // Get satellite number
    satno=find_satno(pdesig);

    // Format IOD line
    printf("%05d %s   %04d G %04d%02d%02d%02d%02d%02d%03d %1.0f%1.0f %d%d ",satno,desig,site,year,month,day,hour,min,sec,fsec,tm,tx,format,epoch);
    printf("%02d%02d%03d%c%02d%02d%02d %1.0f%1.0f\n",rah,ram,rafm,sign,ded,dem,defm,am,ax);
  }
  fclose(file);

  return 0;
}

// Read a line of maximum length int lim from file FILE into string s
int fgetline(FILE *file,char *s,int lim)
{
  int c,i=0;
 
  while (--lim > 0 && (c=fgetc(file)) != EOF && c != '\n')
    s[i++] = c;
  if (c == '\n')
    s[i++] = c;
  s[i] = '\0';
  return i;
}