Add new field to the mesh object to indicate whether the model has rivers. Also update the river functions to use that field and also populate the existing field with the number of rivrs in it.