Changeset 5513

Timestamp:

Nov 14, 2005, 1:44:35 PM (21 years ago)

Author:

drobbin

Message:

Updated Precession

File:

• trunk/psLib/src/astro/psEarthOrientation.c (modified) (8 diffs)

trunk/psLib/src/astro/psEarthOrientation.c

@@ -8 +8 @@
 *  @author Robert Daniel DeSonia, MHPCC
 *
-*  @version $Revision: 1.12 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2005-11-12 03:37:34 $
+*  @version $Revision: 1.13 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2005-11-14 23:44:35 $
 *
 *  Copyright 2005 Maui High Performance Computing Center, University of Hawaii
@@ -349 +349 @@
     // Check for null parameter
     PS_ASSERT_PTR_NON_NULL(time, NULL);
+    if (time->type == PS_TIME_UT1) {
+        psError(PS_ERR_BAD_PARAMETER_VALUE, true,
+                "Invalid time input.  Time cannot be of type UT1.\n");
+        return NULL;
+    }
 
     // Convert psTime to MJD
@@ -529 +534 @@
             // Date too early for tables. Get default time delta value from metadata, and issue warning.
             psLogMsg(__func__,PS_LOG_WARN,PS_ERRORTEXT_psTime_TIME_PREDATES_TABLES,mjd,"xp");
-     
+
             // Lookup value from time metadata loaded from pslib.config
             tableMetadataItem = psMetadataLookup(timeMetadata, "psLib.time.before.xp");
@@ -538 +543 @@
             }
             x = tableMetadataItem->data.F64;
-     
+
         } else if(status == PS_LOOKUP_PAST_BOTTOM) {
             tableMetadataItem = psMetadataLookup(timeMetadata, "psLib.time.predict.mjd");
@@ -547 +552 @@
             }
             mjdPred = tableMetadataItem->data.F64;
-     
+
             // Get xp
             tableMetadataItem = psMetadataLookup(timeMetadata, "psLib.time.predict.xp");
@@ -556 +561 @@
             xp = (psVector*)tableMetadataItem->data.V;
             PS_ASSERT_PTR_NON_NULL(xp,NULL);
-     
+
             // Get yp
             tableMetadataItem = psMetadataLookup(timeMetadata, "psLib.time.predict.yp");
@@ -565 +570 @@
             yp = (psVector*)tableMetadataItem->data.V;
             PS_ASSERT_PTR_NON_NULL(yp,NULL);
-     
+
             // Calculate "a" and "c" constants
             a = TWOPI*(mjd - mjdPred)/365.25;
             c = TWOPI*(mjd - mjdPred)/435.0;
-     
+
             // Calculate x and y polar coordinates
             x = xp->data.F64[0] +
@@ -576 +581 @@
                 xp->data.F64[3]*cos(c) +
                 xp->data.F64[4]*sin(c);
-     
-     
-     
-     
-     
+
+
+
+
+
         }
     */