      program latlon_to_xyz
c
c     ------------------------------------------------------------------
c *   Conversion lat, longitude to x y z                                   
c     ------------------------------------------------------------------

c
      print*, 'lat, long, dist ? (deg. and any units)' 
      read *, rlat, rlon, rdis
      print*, rlat, rlon, rdis
c
      rlat=rlat*3.1415927/180.
      rlon=rlon*3.1415927/180.
c
      call t_sph_to_car(rdis,rlat,rlon,x,y,z)
c
      print*, x,y,z
c       
      stop 'latlon_to_xyz: normal termination'
      end
c
c     ------------------------------------------------------------------
c
