Changeset 5466 for trunk/psLib/src/astro

Timestamp:

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

Author:

drobbin

Message:

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

Location:

trunk/psLib/src/astro

Files:

• psEarthOrientation.c (modified) (3 diffs)
• psEarthOrientation.h (modified) (2 diffs)
• psTime.h (modified) (2 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;
 }
 

trunk/psLib/src/astro/psEarthOrientation.h

@@ -9 +9 @@
 *  @author Robert Daniel DeSonia, MHPCC
 *
-*  @version $Revision: 1.7 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2005-11-02 01:07:25 $
+*  @version $Revision: 1.8 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2005-11-03 04:31:22 $
 *
 *  Copyright 2005 Maui High Performance Computing Center, University of Hawaii
@@ -77 +77 @@
 );
 
-/** Calculate the parallax factor for the given position and time.
+
+/* Calculate the parallax factor for the given position and time.
  *
  *  @return double:     the calculated parallax factor.
  */
+/*
 double psEOC_ParallaxFactor(
     const psSphere *coords,            ///< specified position
     const psTime *time                 ///< specified time
 );
+*/
 
 /** Calculates the components of the rotation between the CEO and GCRS frames, X, Y,

trunk/psLib/src/astro/psTime.h

@@ -11 +11 @@
  *  @author Ross Harman, MHPCC
  *
- *  @version $Revision: 1.41 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-09-21 21:37:20 $
+ *  @version $Revision: 1.42 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-11-03 04:31:22 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -139 +139 @@
 );
 
-/** Convert psTime to UTC or TAI time.
- *
- *  Converts psTime to UTC or TAI time based on the psTimeType argument.
+/** Convert psTime to UTC, TAI, UT1, or TT time.
+ *
+ *  Converts psTime to UTC, TAI, UT1, or TT time based on the psTimeType argument.
  *
  *  @return  psTime*: Pointer to psTime.