created python driving script to compute fractures solutions for cartesian product of input parameters
added OMP statements to parallelize boundary flux integral calcs, tweaked makefile
modified input to allow 0 calc locations (only compute element flowrates), and consolidated element flowrate calc stuff
added comment pointing out mostly repeated code, after re-factor
working on plotting comparison of WR vs multiporosity solutions