Update get_reflection_spherical_fujimura.m

Fix bug on computing grazing angle at 90 degrees
ufrgs-gnss-lab · Nov 28, 2022 · 59e46fb · 59e46fb
1 parent d6e34a1
commit 59e46fb
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/geo-alhazen-aux/get_reflection_spherical_fujimura.m b/geo-alhazen-aux/get_reflection_spherical_fujimura.m
@@ -84,7 +84,11 @@
 geo_ang = rad2deg(angle(u));
 
 %% Grazing Angle
-graz_ang = atand(cotd(geo_ang)-(Rs./Ra)./sind(geo_ang));
+if e==90
+    graz_ang = get_grazing_angle_vector (pos_ant_loc,pos_spec_loc,pos_trans_loc); 
+    return
+end
 
+graz_ang = atand(cotd(geo_ang)-(Rs./Ra)./sind(geo_ang));
 end