Changeset 5466 for trunk/psLib/src/astro/psEarthOrientation.c

Timestamp:

Nov 2, 2005, 6:31:22 PM (21 years ago)

Author:

drobbin

Message:

Updated EOC files(fxns) and tests, specifically TEOtoCEO sphere rotations

File:

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

trunk/psLib/src/astro/psEarthOrientation.c

@@ -9 +9 @@
 *  @author Robert Daniel DeSonia, MHPCC
 *
-*  @version $Revision: 1.8 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2005-11-02 01:07:25 $
+*  @version $Revision: 1.9 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2005-11-03 04:31:22 $
 *
 *  Copyright 2005 Maui High Performance Computing Center, University of Hawaii
@@ -288 +288 @@
 }
 
-
+/*
 double psEOC_ParallaxFactor(const psSphere *coords,
                             const psTime *time)
 {
-
-
+ 
+ 
     return NAN;
 }
+*/
 
 psEarthPole *psEOC_PrecessionModel(const psTime *time)
@@ -455 +456 @@
 psSphereRot* psSphereRot_TEOtoCEO(const psTime *time)
 {
-    return NULL;
+    PS_ASSERT_PTR_NON_NULL(time,NULL);
+    psTime *in = psTimeAlloc(time->type);
+    *in = *time;
+    if (in->type != PS_TIME_UT1) {
+        in = psTimeConvert(in, PS_TIME_UT1);
+    }
+    double T = (double)(in->sec) + (double)(in->nsec / 1e9);
+    printf("\nThe Value of T is = %.13g\n", T);
+    T += -2451545.0;
+    printf("\nThe Value of T is = %.13g\n", T);
+    double theta = 2.0 * M_PI * (0.7790572732640 + 1.00273781191135448 * T);
+    psSphereRot *out = psSphereRotInvert(theta, 0.0, 0.0);
+
+    psFree(in);
+    return out;
 }