Fix the zero-distance problem also at its root, in geodist.f90.

git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@6651 ab8295b8-cf94-4d9e-aec4-7959e3be5d79
This commit is contained in:
Joe Taylor 2016-04-27 20:26:46 +00:00
parent d9c04c5472
commit b2b53fd23e
1 changed files with 9 additions and 0 deletions

View File

@ -29,6 +29,13 @@ subroutine geodist(Eplat,Eplon,Stlat,Stlon,Az,Baz,Dist)
data D2R/0.01745329251994/ ! degrees to radians conversion factor
data Pi2/6.28318530718/
if(abs(Eplat-Stlat).lt.0.02 .and. abs(Eplon-Stlon).lt.0.02) then
Az=0.
Baz=180.0
Dist=0
go to 999
endif
BOA = BL/AL
F = 1.0 - BOA
! Convert st/end pts to radians
@ -93,4 +100,6 @@ subroutine geodist(Eplat,Eplon,Stlat,Stlon,Az,Baz,Dist)
az = 360.0 - az
baz = 360.0 - baz
999 return
end subroutine geodist