WIP: Resolve "WENO flux contains bug"

Closes #19 (closed)

WENO.F90

@@ -688,6 +688,7 @@ subroutine consistent_part(nvars, generalized_fluxes, consistent)
     consistent(v) = generalized_fluxes(v,-2) - 8*generalized_fluxes(v,-1) &
       + 37*generalized_fluxes(v,0) + 37*generalized_fluxes(v,1) &
       - 8*generalized_fluxes(v,2) + generalized_fluxes(v,3)
+    consistent(v) = consistent(v)/60.0
   end do
 end subroutine