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 {