Changeset 3654 in Sophya for trunk/SophyaExt/XephemAstroLib/dbfmt.c

Timestamp:

Jul 16, 2009, 12:34:39 PM (16 years ago)

Author:

cmv

Message:

mise a niveau Xephem 3.7.4, cmv 16/07/2009

File:

• trunk/SophyaExt/XephemAstroLib/dbfmt.c (modified) (5 diffs)

trunk/SophyaExt/XephemAstroLib/dbfmt.c

@@ -18 +18 @@
 #define SUBFLD          '|'     /* subfield separator */
 #define MAXFLDS 20              /* must be more than on any expected line */
+#define MAXESGOOD       100     /* max earth satellite good, days */
 
 static char *enm (char *flds[MAXFLDS]);
@@ -183 +184 @@
  * we detect nonconformance as efficiently as possible.
  * name ends at first '\0', '\r' or '\n'.
+ * set startok/endok.
  * if ok return 0 else return -1
  */
@@ -244 +246 @@
         op->es_M = (float)tle_fld (l2, 44, 51);
         op->es_orbit = (int)tle_fld (l2, 64, 68);
+
+        /* limit date range to decay period that changes period by 1% but
+         * never more than MAXESGOOD.
+         * es_n is rev/day, es_decay is (rev/day)/day
+         */
+        if (fabs(op->es_decay) > 0) {
+            double dt = 0.01*op->es_n/fabs(op->es_decay);
+            if (dt > MAXESGOOD)
+                dt = MAXESGOOD;
+            op->es_startok = op->es_epoch - dt;
+            op->es_endok = op->es_epoch + dt;
+        }
 
         /* yes! */
@@ -576 +590 @@
         if (nf == 12)
             op->es_drag = (float) atod (flds[11]);
+
+        /* if not already specified, limit date range to decay period that
+         * changes period by 1% but never longer than MAXESGOOD.
+         * es_n is rev/day, es_decay is (rev/day)/day
+         */
+        if (op->es_startok == 0 && op->es_endok == 0 && fabs(op->es_decay) > 0){
+            double dt = 0.01*op->es_n/fabs(op->es_decay);
+            if (dt > MAXESGOOD)
+                dt = MAXESGOOD;
+            op->es_startok = op->es_epoch - dt;
+            op->es_endok = op->es_epoch + dt;
+        }
 
         return (0);
@@ -1000 +1026 @@
 
 /* For RCS Only -- Do Not Edit */
-static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: dbfmt.c,v $ $Date: 2008-03-25 17:45:12 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
+static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: dbfmt.c,v $ $Date: 2009-07-16 10:34:36 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};