Added flux function to the initial conditions