[WIP] Add example for a subduction model

Add example for subduction - still in progress
replace endl with std::endl