Add bridge.f90 which calculates the bridging probability