def METHOD_NAME():
    def spherical_limits(x, y, z):
        """
        Checks that the field is inside a sphere of radius 2
        """
        return np.sqrt(x**2 + y**2 + z**2) <= 2
    mips_sl = MercuryiPS(
        "mips_sl",
        address="GPIB::1::INSTR",
        pyvisa_sim_file="MercuryiPS.yaml",
        field_limits=spherical_limits,
    )
    yield mips_sl
    mips_sl.close()