output:
Description: This function calculates upwinded normal flux normal including the mhd terms and put in Uf[n][0]. The interior flux is also put in Us[n]
