Fix a problem where the dot product might possibly be greater than 1.0 causing a...

diff --git a/src/Time/sunpos.cxx b/src/Time/sunpos.cxx
index 81f2d87b121744abef61d839ee236cdea385e38d..0b1e03f595c496fb517912640dbc63b74bd5f25e 100644 (file)
--- a/src/Time/sunpos.cxx
+++ b/src/Time/sunpos.cxx
@@ -341,6 +341,12 @@ void fgUpdateSunPos( void ) {
     dot = sgScalarProductVec3( surface_to_sun, v->get_surface_south() );
     // cout << "  Dot product = " << dot << endl;
 
+    if (dot > 1) {
+        SG_LOG( SG_ASTRO, SG_WARN,
+                "Dot product  = " << dot << " is greater than 1.0" );
+        dot = 1.0;
+    }
+
     if ( east_dot >= 0 ) {
        l->set_sun_rotation( acos(dot) );
     } else {