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