--- 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: Sat Feb 22 09:43:27 2014
# (c) 2003 A.M. Thurnherr
-# uE-Info: 199 23 NIL 0 0 72 0 2 4 NIL ofnI
+# uE-Info: 37 0 NIL 0 0 72 0 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