Refactor trimbd routine

Refactor trimbd routine to fortran 95 and also reduce used array to be
integer(kind=1) since it is simply a flagging array.  Also uses SUM
intrinsic for initial all set check.
Refactor infil routine

Includes some enhancements such as vectorization of some statements.
Also adds new routines to test examples.
Convert flaguse to use integer(kind=1) instead
Remove 2d_edges plot type due to error
Adjust comments
Refactor to f90 of filpatch (filrecur) routine

Initial refactor of basic filrecur routine including minor changes to
names.  Moved statement functions to pure functions defined in the
subroutine.
Add new filpatch source to tests
Change variable names
Variable updates

Continue to put more indices and coordinate rectangles into small arrays
with the corners.  This commit changes a number of additional routines
so that the function signatures are shortened and the rectangles are
passed around rathter than the corners or indices.
Final name changes for filrecur
Refactor of filval routine
Add new routine to tests
Remove cpp directives
Remove pre-processing directives
Use thread-safe Riemann solvers, output times
Rename timing source to new source name
Update timing routines
Remove extraneous timing files
Correct check for Richardson error estimation
Change file input name