             0.0009 * sin(twoPsi) * t3
         
         // Add constants and convert units
-        semimajorAxis *= MetersPerKilometer
-        let inclination = inclinationConstant + inclinationSecPart * RadiansPerArcSecond
-        let longitudeOfPerigee = longitudeOfPerigeeConstant + longitudeOfPerigeeSecPart * RadiansPerArcSecond
-        let meanLongitude = meanLongitudeConstant + meanLongitudeSecPart * RadiansPerArcSecond
-        let longitudeOfNode = longitudeOfNodeConstant + longitudeOfNodeSecPart * RadiansPerArcSecond
+        semimajorAxis *= MetersPerKilometer*/
+        let inclination: Double = inclinationConstant + inclinationSecPart * RadiansPerArcSecond
+        let longitudeOfPerigee: Double = longitudeOfPerigeeConstant + longitudeOfPerigeeSecPart * RadiansPerArcSecond
+        let meanLongitude: Double = meanLongitudeConstant + meanLongitudeSecPart * RadiansPerArcSecond
+        let longitudeOfNode: Double = longitudeOfNodeConstant + longitudeOfNodeSecPart * RadiansPerArcSecond
         
         return elementsToCartesian(
             semimajorAxis: semimajorAxis,

    func computeMoonPositionInEarthInertialFrame (date: JulianDate = JulianDate.now()) -> Cartesian3 {
        return _axesTransformation.multiplyByVector(computeSimonMoon(date))
    }
    */*/*/*/*/*/*/*/*/*/*/
}

    func computeSimonMoon(date: JulianDate) -> Cartesian3 {
        let tdbDate = taiToTdb(date)
        let x = Double(tdbDate.dayNumber - epoch.dayNumber) + ((tdbDate.secondsOfDay - epoch.secondsOfDay) / TimeConstants.SecondsPerDay)
        let t: Double = x / TimeConstants.DaysPerJulianCentury
        let t2: Double = t * t
        let t3: Double = t2 * t
        let t4: Double = t3 * t
        
        // Terms from section 3.4 (b.1)
        var semimajorAxis: Double = 383397.7725 + 0.0040 * t
        var eccentricity: Double = 0.055545526 - 0.000000016 * t
        
        let inclinationConstant: Double = 5.15668983 * RadiansPerDegree
        let inclinationSecPart: Double = -0.00008 * t + 0.02966 * t2 -
            0.000042 * t3 - 0.00000013 * t4
        let longitudeOfPerigeeConstant: Double = 83.35324312 * RadiansPerDegree
        let longitudeOfPerigeeSecPart: Double = 14643420.2669 * t - 38.2702 * t2 -
            0.045047 * t3 + 0.00021301 * t4
        let longitudeOfNodeConstant: Double = 125.04455501 * RadiansPerDegree
        let longitudeOfNodeSecPart: Double = -6967919.3631 * t + 6.3602 * t2 +
            0.007625 * t3 - 0.00003586 * t4
        let meanLongitudeConstant: Double = 218.31664563 * RadiansPerDegree
        let meanLongitudeSecPart: Double = 1732559343.48470 * t - 6.3910 * t2 +
            0.006588 * t3 - 0.00003169 * t4
        
        // Delaunay arguments from section 3.5 b
        let D: Double = 297.85019547 * RadiansPerDegree + RadiansPerArcSecond *
            (1602961601.2090 * t - 6.3706 * t2 + 0.006593 * t3 - 0.00003169 * t4)
        let F: Double = 93.27209062 * RadiansPerDegree + RadiansPerArcSecond *
            (1739527262.8478 * t - 12.7512 * t2 - 0.001037 * t3 + 0.00000417 * t4)
        let l: Double = 134.96340251 * RadiansPerDegree + RadiansPerArcSecond *
            (1717915923.2178 * t + 31.8792 * t2 + 0.051635 * t3 - 0.00024470 * t4)
        let lprime: Double = 357.52910918 * RadiansPerDegree + RadiansPerArcSecond *
            (129596581.0481 * t - 0.5532 * t2 + 0.000136 * t3 - 0.00001149 * t4)
        let psi: Double = 310.17137918 * RadiansPerDegree - RadiansPerArcSecond *
            (6967051.4360 * t + 6.2068 * t2 + 0.007618 * t3 - 0.00003219 * t4)
        
        // Add terms from Table 4
        let twoD: Double = 2.0 * D
        let fourD: Double = 4.0 * D
        let sixD: Double = 6.0 * D
        let twol: Double = 2.0 * l
        let threel: Double = 3.0 * l
        let fourl: Double = 4.0 * l
        let twoF: Double = 2.0 * F
        semimajorAxis += 3400.4 * cos(twoD)
        semimajorAxis -= 635.6 * cos(twoD - l)
        semimajorAxis -= 235.6 * cos(l)
        semimajorAxis += 218.1 * cos(twoD - lprime)
        semimajorAxis += 181.0 * cos(twoD + l)
        eccentricity += 0.014216 * cos(twoD - l) + 0.008551 * cos(twoD - twol)
            //- 0.001383 * cos(l) + 0.001356 * cos(twoD + l)
            //- 0.001147 * cos(fourD - threel) - 0.000914 * cos(fourD - twol)
            //+ 0.000869 * cos(twoD - lprime - l) - 0.000627 * cos(twoD)
            //- 0.000394 * cos(fourD - fourl) + 0.000282 * cos(twoD - lprime - twol)
            //- 0.000279 * cos(D - l) - 0.000236 * cos(twol)
            //+ 0.000231 * cos(fourD) + 0.000229 * cos(sixD - fourl)
            //- 0.000201 * cos(twol - twoF)
        /*inclinationSecPart += 486.26 * cos(twoD - twoF) - 40.13 * cos(twoD) +
            37.51 * cos(twoF) + 25.73 * cos(twol - twoF) +
            19.97 * cos(twoD - lprime - twoF)
        longitudeOfPerigeeSecPart += -55609 * sin(twoD - l) - 34711 * sin(twoD - twol) -
