X-Git-Url: http://dxcluster.org/gitweb/gitweb.cgi?a=blobdiff_plain;f=perl%2FSun.pm;fp=perl%2FSun.pm;h=09ebdbac32248e2969ea4335ade79749c628416b;hb=8dfb46f118199934c742e1fcaee0ceb8084a6998;hp=09af995baf0a0dc9a981bd936841d04109b3f52b;hpb=02ab82e5a26f6055cc6379ef2ff1a16fc93be2e9;p=spider.git diff --git a/perl/Sun.pm b/perl/Sun.pm index 09af995b..09ebdbac 100644 --- a/perl/Sun.pm +++ b/perl/Sun.pm @@ -608,9 +608,9 @@ sub get_satellite_pos my $epoch_day=$epoch-int(1000*$epoch_year); #printf("epoch_year = %10.2f\n",$epoch_year); #printf("epoch_day = %17.12f\n",$epoch_day); - $epoch_year=$epoch_year+2000 if ($epoch_year < 57); - $epoch_year=$epoch_year+1900 if ($epoch_year >= 57); - my $jt_epoch=Julian_Date_of_Year($epoch_year); + my $ep_year=$epoch_year+2000 if ($epoch_year < 57); + $ep_year=$epoch_year+1900 if ($epoch_year >= 57); + my $jt_epoch=Julian_Date_of_Year($ep_year); $jt_epoch=$jt_epoch+$epoch_day; #printf("JT for epoch = %17.12f\n",$jt_epoch); my $tsince=($jtime-$jt_epoch)*24*60; @@ -673,16 +673,15 @@ sub get_satellite_pos my $xl=mod2p($xls-$c5/$p*$axnsl); my $u=mod2p($xl-$xnodes); - my $item3=0; + my $item3; my $eo1=$u; my $tem5=1; my $coseo1=0; my $sineo1=0; - while ( abs($tem5) >= 1e-6 && $item3 < 10 ) + for ($item3=0; abs($tem5) >= 1e-6 && $item3 < 10; $item3++ ) { $sineo1=sin($eo1); $coseo1=cos($eo1); - $item3 = $item3+1; $tem5=1-$coseo1*$axnsl-$sineo1*$aynsl; $tem5=($u-$aynsl*$coseo1+$axnsl*$sineo1-$eo1)/$tem5; my $tem2=abs($tem5);