def METHOD_NAME():
    radar = pyart.testing.make_empty_spectra_radar(5, 5, 2)
    radar.ds["azimuth"] = xr.DataArray(np.array([0, 90, 180, 270, 0]), dims="time")
    radar.ds["elevation"] = xr.DataArray(np.array([0, 0, 0, 0, 10]), dims="time")
    radar.ds["range"] = xr.DataArray(np.array([5, 15, 25, 35, 45]), dims="range")
    radar.init_gate_x_y_z()
    assert radar.gate_x.values.shape == (5, 5)
    assert_allclose(radar.gate_x.values[0], [0, 0, 0, 0, 0], atol=1e-14)
    assert_allclose(radar.gate_x.values[1], [5, 15, 25, 35, 45], atol=1e-14)
    assert_allclose(radar.gate_x.values[2], [0, 0, 0, 0, 0], atol=1e-5)
    assert_allclose(radar.gate_x.values[3], [-5, -15, -25, -35, -45], atol=1e-14)
    assert radar.gate_y.values.shape == (5, 5)
    assert_allclose(radar.gate_y.values[0], [5, 15, 25, 35, 45], atol=1e-14)
    assert_allclose(radar.gate_y.values[1], [0, 0, 0, 0, 0], atol=1e-5)
    assert_allclose(radar.gate_y.values[2], [-5, -15, -25, -35, -45], atol=1e-14)
    assert_allclose(radar.gate_y.values[3], [0, 0, 0, 0, 0], atol=1e-6)
    assert radar.gate_z.values.shape == (5, 5)
    z_sweep0 = np.array([1.47e-6, 1.324e-5, 3.679e-5, 7.210e-5, 1.1919e-4])
    assert_allclose(radar.gate_z.values[0], z_sweep0, atol=1e-3)
    assert_allclose(radar.gate_z.values[1], z_sweep0, atol=1e-3)
    assert_allclose(radar.gate_z.values[2], z_sweep0, atol=1e-3)
    assert_allclose(radar.gate_z.values[3], z_sweep0, atol=1e-3)