code: def get_xyz_from_ll(projected, reference):
      
  
      azims = geod.azimuth(reference.longitude, reference.latitude,
                           projected.longitude, projected.latitude)
      depths = np.subtract(reference.depth, projected.depth)
      dists = geod.geodetic_distance(reference.longitude,
                                     reference.latitude,
                                     projected.longitude,
                                     projected.latitude)
      return (dists * math.sin(math.radians(azims)),
              dists * math.cos(math.radians(azims)),
              depths)
docs: This method computes the x, y and z coordinates of a set of points
      provided a reference point
  
      :param projected:
          :class:`~openquake.hazardlib.geo.point.Point` object
          representing the coordinates of target point to be projected
      :param reference:
          :class:`~openquake.hazardlib.geo.point.Point` object
          representing the coordinates of the reference point.
  
      :returns:
              x
              y
              z