Bug fix from Jean

src/CASCADE/diffusion_erosion.f

@@ -28,7 +28,7 @@ c                     Voronoi cells
 c OUTPUT: a range of arrays/variables are updated:
 c        h, dh, dhmin, dhmax, outflux
 
-c subroutines calles:
+c subroutines called:
 c - none
 
       common /vocal/ ivocal
@@ -40,7 +40,7 @@ c - none
       integer nb(nnode),nn(nbmax,nnode)
       real    sides(nbmax,nnode)
       real    outflux
-      double precision  dhmin,dhmax, dhh
+      double precision dhmin,dhmax,dhh
 
 c finds height change
 
@@ -68,9 +68,14 @@ c finds height change
             outflux=outflux+dhh*surf
             dhh=0.
             endif
-            if (h(i)+dhh.lt.sea_level) then
-            dhh=sea_level-h(i)
-            endif
+!            if (h(i)+dhh.lt.sea_level) then
+!            dhh=sea_level-h(i)
+!            endif
+           if ((h(i)-dhh.lt.sea_level).and.(h(i).ge.sea_level)) then
+           dhh=sea_level-h(i)            
+           elseif ((h(i)-dhh.lt.sea_level).and.(h(i).lt.sea_level)) then
+           dhh=0.0        
+           endif
           dh(i)=-dhh
           h(i)=h(i)-dhh
           dhmin=min(dhmin,dhh)
@@ -80,3 +85,4 @@ c finds height change
 
       return
       end
+