Create fluxcalculation.f