Fixed a concurrency bug in flux function calculation
Added a parameter Bmin to control minimum pair-creating B field