sattools/src/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;
}