Merge pull request #39 from fastscape-lem/FixFluxes

Fixed bug in estimating erosional flux (contribution from Marine.f90 …
fastscape-lem · Jun 30, 2020 · 28c6afc · 28c6afc
2 parents c231a67 + 23c7de8
commit 28c6afc
Showing 1 changed file with 3 additions and 0 deletions.
diff --git 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