Morphological modeling using a fully coupled, total variation diminishing upwind-biased centered scheme