From 23c7de860c5434735705071773dfc6ff3aa82de0 Mon Sep 17 00:00:00 2001
From: jeanbraun <jeanbraun@me.com>
Date: Tue, 30 Jun 2020 13:11:59 +0200
Subject: [PATCH] Fixed bug in estimating erosional flux (contribution from
 Marine.f90 removed)

---
 src/Marine.f90 | 3 +++
 1 file changed, 3 insertions(+)

diff --git a/src/Marine.f90 b/src/Marine.f90
index f9a4f10..6ab7fe2 100644
--- a/src/Marine.f90
+++ b/src/Marine.f90
@@ -148,6 +148,9 @@ subroutine Marine()
   !h=ht+dh
   etot=etot+ht-h
   erate=erate+(ht-h)/dt
+  where (h.lt.sealevel) Sedflux=0.d0
+  where (h.lt.sealevel) etot=0.d0
+  where (h.lt.sealevel) erate=0.d0
 
   ! set the silt fraction in continent
   where (h.ge.sealevel+1.d-3) Fmix=0.d-1