WSJT-X/geocentric.f

      subroutine geocentric(alat,elev,hlt,erad)

      implicit real*8 (a-h,o-z)

C  IAU 1976 flattening f, equatorial radius a
      f = 1.d0/298.257d0
      a = 6378140.d0
      c = 1.d0/sqrt(1.d0 + (-2.d0 + f)*f*sin(alat)*sin(alat))
      arcf = (a*c + elev)*cos(alat)
      arsf = (a*(1.d0 - f)*(1.d0 - f)*c + elev)*sin(alat)
      hlt = datan2(arsf,arcf)
      erad = sqrt(arcf*arcf + arsf*arsf)
      erad = 0.001d0*erad

      return
      end
First good compile with map65a running inside map65. git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/map65@334 ab8295b8-cf94-4d9e-aec4-7959e3be5d79 2007-01-11 16:25:52 -05:00			`subroutine geocentric(alat,elev,hlt,erad)`

			`implicit real*8 (a-h,o-z)`

			`C IAU 1976 flattening f, equatorial radius a`
			`f = 1.d0/298.257d0`
			`a = 6378140.d0`
			`c = 1.d0/sqrt(1.d0 + (-2.d0 + f)fsin(alat)*sin(alat))`
			`arcf = (ac + elev)cos(alat)`
			`arsf = (a(1.d0 - f)(1.d0 - f)c + elev)sin(alat)`
			`hlt = datan2(arsf,arcf)`
			`erad = sqrt(arcfarcf + arsfarsf)`
			`erad = 0.001d0*erad`

			`return`
			`end`