RDI_Coords.pl
changeset 14 8c79b38a7086
parent 13 b176da8559b3
child 16 68a9fc5e7d45
--- a/RDI_Coords.pl
+++ b/RDI_Coords.pl
@@ -1,9 +1,9 @@
 #======================================================================
 #                    R D I _ C O O R D S . P L 
 #                    doc: Sun Jan 19 17:57:53 2003
-#                    dlm: Wed Aug  7 11:18:51 2013
+#                    dlm: Wed Nov 27 11:21:49 2013
 #                    (c) 2003 A.M. Thurnherr
-#                    uE-Info: 199 23 NIL 0 0 72 0 2 4 NIL ofnI
+#                    uE-Info: 287 70 NIL 0 0 72 2 2 4 NIL ofnI
 #======================================================================
 
 # RDI Workhorse Coordinate Transformations
@@ -33,6 +33,7 @@
 #	Jan 15, 2012: - replaced defined(@...) by (@...) to get rid of warning
 #	Aug  7, 2013: - BUG: &velBeamToBPInstrument did not return any val unless
 #						 all beam velocities are defined
+#	Nov 27, 2013: - added &RDI_pitch(), &tilt_azimuth()
 
 use strict;
 use POSIX;
@@ -274,6 +275,18 @@
 	return deg(atan(tan(rad($RDI_pitch)) * cos(rad($RDI_roll))));
 }
 
+sub RDI_pitch($$)
+{
+	my($gimbal_pitch,$roll) = @_;
+	return deg(atan(tan(rad($gimbal_pitch))/cos(rad($roll))));
+}
+
+sub tilt_azimuth($$)
+{
+	my($gimbal_pitch,$roll) = @_;
+	return angle(deg(atan2(sin(rad($gimbal_pitch)),sin(rad($roll)))));
+}
+
 # - angle from vertical is home grown and should be treated with caution
 # - angle between two unit vectors given by acos(v1 dot v2)
 # - vertical unit vector v1 = (0 0 1) => dot product = z-component of v2