Changeset 27840 for branches/simtest_nebulous_branches/ippScripts/scripts/magic_tree.pl

Timestamp:

May 3, 2010, 8:50:52 AM (16 years ago)

Author:

eugene

Message:

updates from trunk

Location:

branches/simtest_nebulous_branches

Files:

• . (modified) (1 prop)
• ippScripts/scripts/magic_tree.pl (modified) (3 diffs)

branches/simtest_nebulous_branches/ippScripts/scripts/magic_tree.pl

@@ -8 +8 @@
 use Sys::Hostname;
 my $host = hostname();
+my $date = `date`;
 print "\n\n";
-print "Starting script $0 on $host\n\n";
+print "Starting script $0 on $host at $date\n\n";
 
 use vars qw( $VERSION );
@@ -158 +159 @@
 
     # Relative coordinates of centre of the field
-    my $x = $naxis1 - $crpix1;
-    my $y = $naxis2 - $crpix2;
+    my $x = $naxis1/2 - $crpix1;
+    my $y = $naxis2/2 - $crpix2;
 
     # Coordinates on tangent plane
@@ -168 +169 @@
 
     # Coordinates on rotated celestial sphere
-    my $phi = atan2($eta,$xi) + pi/2;
-    my $theta = atan(180 / pi / sqrt($xi**2 + $eta**2));
+    my ($phi, $theta);
+    if ($xi == 0 and $eta == 0) {
+        $phi = 0;
+        $theta = 0;
+    } else {
+        $phi = atan2($eta,$xi) + pi/2;
+        my $denominator = sqrt($xi**2 + $eta**2);
+        &my_die("denominator is zero!!", $magic_id, $PS_EXIT_PROG_ERROR) if $denominator == 0;
+        $theta = atan(180 / pi / $denominator);
+    }
 
     # Coordinates on celestial sphere