code: def make_lat_lons(cvects):
      
      lats = np.degrees(np.arcsin(cvects[2]))
      lons = np.degrees(np.arctan2(cvects[0], cvects[1]))
      return np.hstack([lats, lons])
docs: Convert from directional cosines to latitidue and longitude
  
      Parameters
      ----------
      cvects : directional cosine (i.e., x,y,z component) values
  
      returns (np.ndarray(2,nsrc)) with the directional cosine (i.e., x,y,z component) values