remove q_c and direction from functions and make them as a initial parameter
There is no benefit of having q_c and direction being set in the functions since the bse calculations are not remade when new values of q_c and direction are given to the functions. Therefore we have moved them to the initialization.