def METHOD_NAME(Nphi, Ntheta, Nr, dtype, dealias):
    """Tests ghost broadcasting in colatitude."""
    c, d, b, phi, theta, r, x, y, z = build_ball(Nphi, Ntheta, Nr, dtype, dealias, mesh=(2,2))
    b_S2 = b.S2_basis()
    f = field.Field(dist=d, dtype=dtype)
    g = field.Field(dist=d, bases=(b_S2,), dtype=dtype)
    h = field.Field(dist=d, bases=(b_S2,), dtype=dtype)
    f['g'] = 6
    g['g'] = (5*np.cos(theta)**2 - 1) * np.sin(theta) * np.cos(phi)
    h['g'] = 6 * (5*np.cos(theta)**2 - 1) * np.sin(theta) * np.cos(phi)
    h_op = (f * g).evaluate()
    h_op.change_scales(1)
    assert np.allclose(h_op['g'], h['g'])