home *** CD-ROM | disk | FTP | other *** search
/ PC World 1999 April / PCWorld_1999-04_cd.bin / Software / Vyzkuste / HomePlan / ASTRELEM.C < prev    next >
C/C++ Source or Header  |  1994-11-06  |  6KB  |  221 lines
  1. /*
  2.  
  3.     Convert asteroid elements in the 8-line format used
  4.     by the Minor Planet Center (and in the Minor Planet
  5.     Electronic Circulars) to our CSV format.
  6.     
  7.     Designed and implemented by John Walker in October of 1994.
  8.     
  9.      (4646) Kwee
  10.     Epoch 1993 Jan. 13.0 TT = JDT 2449000.5  (M-P)          Oishi
  11.     M 220.20030              (2000.0)            P               Q
  12.     n   0.26003234     Peri.  350.28914     +0.79481320     +0.60665475
  13.     a   2.4309973      Node   332.34432     -0.55450494     +0.71559595
  14.     e   0.1921938      Incl.    1.92064     -0.24656897     +0.34625491
  15.     P   3.79           H   14.0           G   0.15
  16.     From observations at 4 oppositions, 1960-1990.  Ref. MPC 17195.
  17.               1         2         3         4         5         6         7
  18.     01234567890123456789012345678901234567890123456789012345678901234567890
  19.     
  20.     Name,Magnitude H,Magnitude G,Mean anomaly,Arg. perihelion,Long. node,Inclination,Eccentricity,Semimajor axis,Epoch (MJD)
  21.     AALTJE 677,9.70,0.15,47.68385,272.29565,272.71469,8.48713,0.04530949,2.95528649,48600
  22.  
  23.     
  24. */
  25.  
  26. #include <stdio.h>
  27. #include <ctype.h>
  28. #include <string.h>
  29. #include <math.h>
  30.  
  31. static char s[8][132];                    /* Buffer to hold complete entry */
  32.  
  33. static int rs(char *s, FILE *fp)
  34. {
  35.     if (fgets(s, 132, fp) == NULL) {
  36.         return 0;
  37.     }
  38.     while (strlen(s) > 0 && s[strlen(s) - 1] < ' ') {
  39.         s[strlen(s) - 1] = 0;
  40.     }
  41.     return 1;
  42. }
  43.  
  44. static void sc(char *b, int line, int start, int end) {
  45.     char *c = s[line];
  46.     int i, n;
  47.     
  48.     for (i = 0, n = start; n <= end; n++) {
  49.         if (c[n] == 0) {
  50.             break;
  51.         }
  52.         if (!isspace(c[n])) {
  53.             b[i++] = c[n];
  54.         }
  55.     }
  56.     b[i] = 0;
  57. }
  58.  
  59. int main(int argc, char *argv[])
  60. {
  61.     int i, p, running = 1, recov, written = 0;
  62. #define IBN    4    
  63.     static char ib[IBN][132];
  64.     FILE *fi = stdin, *fo = stdout;            /* Input and output file names */
  65. #ifdef EDITMODE
  66.     FILE *fn = NULL;
  67. #endif    
  68.     static char anum[80], aname[80], epoch[80], meanan[80], semimaj[80],
  69.         eccen[80], argperi[80], longnode[80], inclin[80], magh[80], magg[80];
  70.     
  71.     if (argc > 1) {
  72.         fi = fopen(argv[1], "r");            /* Input file */
  73.         if (fi == NULL) {
  74.             fprintf(stderr, "Cannot open input file %s\n", argv[1]);
  75.             return 2;
  76.         }
  77.         if (argc > 2) {
  78.             fo = fopen(argv[2], "w");        /* Output file */
  79.             if (fo == NULL) {
  80.                 fprintf(stderr, "Cannot open output file %s\n", argv[2]);
  81.                 return 2;
  82.             }
  83. #ifdef EDITMODE            
  84.             if (argc > 3) {
  85.                 fn = fopen(argv[3], "w");    /* Edit mode file */
  86.                 if (fn == NULL) {
  87.                     fprintf(stderr, "Cannot open numeric format output file %s\n", argv[3]);
  88.                     return 2;
  89.                 }
  90.             }
  91. #endif            
  92.         }
  93.     }
  94.  
  95. #ifdef EDITMODE    
  96.     fprintf(fo, "Number,Name,Number Name,Name Number,Magnitude H,Magnitude G,Mean anomaly,Arg. perihelion,Long. node,Inclination,Eccentricity,Semimajor axis,Epoch (MJD)\n");
  97.     if (fn != NULL) {
  98.         fprintf(fn, "Number,Name,Magnitude H,Magnitude G,Mean anomaly,Arg. perihelion,Long. node,Inclination,Eccentricity,Semimajor axis,Epoch (MJD)\n");
  99.     }
  100. #else
  101.     fprintf(fo, "Name,Magnitude H,Magnitude G,Mean anomaly,Arg. perihelion,Long. node,Inclination,Eccentricity,Semimajor axis,Epoch (MJD)\n");
  102. #endif
  103.     
  104.     while (1) {
  105.         char *u, *v;
  106.  
  107.         p = 0;
  108.         while (1) {
  109. #define ibm(x)    ib[(p + (IBN - (x))) % IBN]
  110.             if (!rs(ib[p], fi)) {
  111.                 running = 0;
  112.                 break;
  113.             }
  114.             if (strncmp(ib[p], "Epoch ", 6) == 0 &&
  115.                 strstr(ib[p], "TT") != NULL &&
  116.                 strchr(ib[p], '=') != NULL &&
  117.                 strstr(ib[p], "JDT") != NULL) {
  118.         
  119.                 /* MPECs for recovered objects contain an "Id." line after the name
  120.                    and before the Epoch, instead of a "From" line at the end of the
  121.                    elements.  If this is the case with these elements, move the ID
  122.                    to the end of the elements and shift everything else up one line. */
  123.                    
  124.                 recov = strncmp(ibm(1), "Id", 2) == 0;
  125.  
  126.                 strcpy(s[0], ibm(recov ? 2 : 1));
  127.                 strcpy(s[1], ib[p]);
  128.                 for (i = 0; i < 5; i++) {
  129.                     if (!rs(s[i + 2], fi)) {
  130.                         running = 0;
  131.                         break;
  132.                     }
  133.                 }
  134.                 p = (p + 1) % IBN; 
  135.                 break;
  136.             }
  137.             p = (p + 1) % IBN; 
  138.         }
  139.         if (!running) {
  140.             break;
  141.         }        
  142.         
  143.         /*    Scan asteroid number, if one has been assigned. */
  144.         
  145.         if ((u = strchr(s[0], '(')) != NULL && (v = strchr(u + 1, ')')) != NULL) {
  146.             u += 1;
  147.             *v++ = 0;
  148.             strcpy(anum, u);
  149.         } else {
  150.             v = s[0];
  151.             strcpy(anum, "(Unnamed)");
  152.         }
  153.         
  154.         /*    Scan name or other designation */
  155.         
  156.         while (isspace(*v)) {
  157.             v++;
  158.         }
  159.         while (strlen(v) > 0 && isspace(v[strlen(v) - 1])) {
  160.             v[strlen(v) - 1] = 0;    
  161.         }
  162.         strcpy(aname, v);
  163.         
  164.         /* Epoch isn't always precisely aligned in MPC postings.
  165.            Allow for some slack. */
  166.            
  167.         v = strstr(s[1], "JDT");
  168.         if (v == NULL) {
  169.             strcpy(epoch, "0");            /* Zero should get somebody's attention */
  170.         } else {
  171.             sprintf(epoch, "%g", atof(v + 3) - 2400000.5);
  172.         }
  173.         
  174.         sc(meanan, 2, 2, 15);
  175.         sc(semimaj, 4, 2, 15);
  176.         sc(eccen, 5, 2, 15);
  177.         sc(argperi, 3, 26, 35);
  178.         sc(longnode, 4, 26, 35);
  179.         sc(inclin, 5, 26, 35);
  180.         sc(magg, 6, 39, 48);
  181.         sc(magh, 6, 20, 30);
  182.         
  183.         /* Mean anomaly is blank in some MPECs, meaning zero.
  184.            Home Planet correctly interprets the blank field as zero,
  185.            but it might confuse other programs.  Replace it with 0.0 if
  186.            blank. */
  187.            
  188.         if (strlen(meanan) == 0) {
  189.             strcpy(meanan, "0.0");
  190.         }
  191.         
  192.         /* Output information in CSV format */
  193.  
  194. #ifdef EDITMODE        
  195.         fprintf(fo, "%s,%s,%s %s,%s %s,%s,%s,%s,%s,%s,%s,%s,%s,%s\n",
  196.                     anum, aname, anum, aname, aname, anum, magh, magg, meanan, argperi,
  197.                     longnode, inclin, eccen, semimaj, epoch);
  198.         if (fn != NULL) {
  199.             fprintf(fn, "%s,%s %s,%s,%s,%s,%s,%s,%s,%s,%s,%s\n",
  200.                         anum, anum, aname, magh, magg, meanan, argperi,
  201.                         longnode, inclin, eccen, semimaj, epoch);
  202.         }
  203. #else
  204.         fprintf(fo, "%s %s,%s,%s,%s,%s,%s,%s,%s,%s,%s\n",
  205.                     aname, anum, magh, magg, meanan, argperi,
  206.                     longnode, inclin, eccen, semimaj, epoch);
  207. #endif
  208.         written++;                    
  209.     }
  210.     
  211.     fclose(fi);
  212.     fclose(fo);
  213. #ifdef EDITMODE    
  214.     if (fn != NULL) {
  215.         fclose(fn);
  216.     }
  217. #endif
  218.     fprintf(stderr, "%d record%s written.\n", written, written == 1 ? "" : "s");        
  219.     
  220.     return 0;
  221. }