nan bug

StrayLight.cpp

@@ -588,6 +588,12 @@ void EarthShine(double ju, double sat[3], double ob[3], double py[3], double E[7
                     XYZ ves = XYZ(xs[w][l], ys[w][l] - ye[i][j], zs[w][l] - ze[i][j]);//镜头分块到盖板小块
                     double tyj = acos(ves.X / ves.Length);
                     double fwj = acos(ves.Y / (sqrt(ves.Y*ves.Y + ves.Z*ves.Z)));
+                    if (isfinite(tyj) == false){
+                        tyj = 0;
+                    }
+                    if (isfinite(fwj) == false){
+                        fwj = 0;
+                    }
                     double ppp = pow(10, PST(fwj, tyj));
                     Econe[i][j] += Eshield[w][l] / 2 / PI / ves.Square()*0.03*pow(10, PST(fwj, tyj));
                 }