We study operation of a new device, the superconducting differential double contour interferometer (DDCI), in application for the ultra sensitive detection of magnetic flux and for digital read out of the state of the superconducting flux qubit. DDCI consists of two superconducting contours weakly coupled by Josephson Junctions. In such a device a change of the critical current and the voltage happens in a step-like manner when the angular momentum quantum number changes in one of the two contours. The DDCI may outperform traditional Superconducting Quantum Interference Devices when the change of the quantum number occurs in a narrow magnetic field region near the half of the flux quantum due to thermal fluctuations, quantum fluctuations, or the switching a loop segment in the normal state for a while by short pulse of an external current. Higher sensitivity of DDCI compared to the conventional SQUID is provided by strong discreteness of the energy spectrum of the continuous superconducting loop [6]. According to the conventional theory [7] the total energy of the persistent currentin a loop with small cross section s ≪ λ 2 L (T ) is determined mainly by the kinetic energy [8]:L /s)L is the kinetic inductance of the loop of side l; L ≈ µ 0 l is the magnetic inductance; s is the cross section of superconducting wires; n s is the density of the Cooper pairs;is the London penetration depth. Two permitted states n and n + 1 have equal energy at Φ = (n + 0.5)Φ 0 according to (2) and thus equal probability P ∝ exp(−E k /k B T ). The probability of other permitted states is negligible and P (n+1) ≈ 1−P (n) at Φ ≈ (n+0.5)Φ 0 when Φ 2 0 /2L k = I p,A Φ 0 ≫ k B T . The probability of the n state may be described with the relationat the magnetic flux Φ = (n + 0.5)Φ 0 + δΦ, when δΦ ≪ Φ 0 . Here I p,A = Φ 0 /2L k is the persistent current (1) at |n − Φ/Φ 0 | = 1/2 and ǫ = I p,A Φ 0 /k B T . The probability (3) changes from P (n) ≈ 1 to P (n) ≈ 0 in a region of the magnetic flux from δΦ ≈ −Φ 0 /2ǫ to δΦ ≈ Φ 0 /2ǫ. This region may be very narrow due to a big value ǫ = I p,A Φ 0 /k B T ≫ 1 equal, for example ǫ ≈ 1500 at the temperature of measurement T ≈ 1 K and a typical value I p,A = 10 µA measured, for example in [9]. Measurements [10] of flux qubit (superconducting loop with three