We present a method for computing first order asymptotics of semiclassical spectra for 1-D Bogoliubov-de Gennes (BdG) Hamiltonian from Supraconductivity, which models the electron/hole scattering through two SNS junctions. This involves: 1) reducing the system to Weber equation near the branching point at the junctions; 2) constructing local sections of the fibre bundle of microlocal solutions; 3) normalizing these solutions for the “flux norm” associated to the microlocal Wronskians; 4) finding the relative monodromy matrices in the gauge group that leaves invariant the flux norm; 5) from this we deduce Bohr-Sommerfeld (BS) quantization rules that hold precisely when the fibre bundle of microlocal solutions (depending on the energy parameter E) has trivial holonomy. Such a semi-classical treatement reveals interesting continuous symetries related to monodromy.