Added a conversion factor for contour elevations (contour_conversion) --> defaults to unity
<<NB Prior to merging this pull request, we need to verify that I implemented this correctly. Not yet sure>>

Implemented the compute_zonal flag so it actually determined if compute zonal will be run

Handled the special case where user has already clipped NHD flowlines to the model domain (thus no edge intersection issues which was throwing exceptions)