code: def create_circular_hole(
      x_grid, y_grid, hole_radius=0.5, center=(0.0, 0.0)
  ):
      X = np.where(
          (x_grid - center[0]) ** 2 + (y_grid - center[1]) ** 2
          <= hole_radius**2,
          np.NAN,
          x_grid,
      )
      Y = np.where(
          (x_grid - center[0]) ** 2 + (y_grid - center[1]) ** 2
          <= hole_radius**2,
          np.NAN,
          y_grid,
      )
  
      return X, Y
  
file_path: vispy/examples/basics/visuals/axially_symmetric_surfaces.py
imports: import sys
  import numpy
  from vispy import scene
  from vispy.scene.visuals import SurfacePlot
  from vispy.scene.visuals import XYZAxis
label_desc: Use image node_id 4 for calling a global function with example usage: create_circular_hole(x_grid, y_grid, hole_radius, center) and returns: X,  Y