added derivative of total element flowrate, and added OMP around mathieu funciton initialization. re-factored some io to move if/thens out of loops