def METHOD_NAME():
    radar = pyart.testing.make_empty_spectra_radar(5, 4, 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]), dims="range")
    assert radar.gate_longitude.data.shape == (5, 4)
    assert radar.gate_latitude.data.shape == (5, 4)
    assert_almost_equal(radar.gate_longitude.data[0, 0], -97.5, 1)
    assert_almost_equal(radar.gate_latitude.data[0, 0], 36.5, 1)
    # reset and try again with a non-default lat_0/lon_0
    radar.ds.projection.pop("_include_lon_0_lat_0")
    radar.ds.projection["lat_0"] = 20.0
    radar.ds.projection["lon_0"] = 60.0
    radar.init_gate_longitude_latitude()
    assert_almost_equal(radar.gate_longitude.data[0, 0], 60.0, 1)
    assert_almost_equal(radar.gate_latitude.data[0, 0], 20.0, 1)