anchor.diff:        end do
         do s = 1, stnserod
            xdepth = rightdepth(s)
  -         if (curvature(s) .lt. 0.0) xdepth = leftdepth(s)
  -         if (xdepth .ne. 0.0) then
  +         if (curvature(s) .lt. 0.d0) xdepth = leftdepth(s)
  +         if (xdepth .ne. 0.d0) then
               tauwall(s) = tauwall(s) / xdepth
               spreaddelta_x(s) = tauwall(s) * (-1.d0) * sin(phi(s))
               spreaddelta_y(s) = tauwall(s) * cos(phi(s))
  
positive.code:             do while (xs(sp) .le. xtrmnt .and. sp .le. stnserod)
                 if (xs(sp) .ge. xstrt) then
                    xdel = sign(xdest - xs(sp), 1.d0)
                    if (lambda(s) .ne. 0.0) then
                       gaussian = exp(-1.d0 * xdel ** 2.d0 / 2.d0 
       +                    / lambda(s) ** 2.d0) / sqrt(2.d0 * 3.1416d0) 
       +                    / lambda(s) 
  
negative.code: c     stress calculation
  c                 1.7  8/11: debugged version SL
  c
  c     $Id: meander.f,v 1.4 2002-08-14 10:34:26 arnaud Exp $
  c
        subroutine meander (stations, stnserod, x, y, xs, dels, flow,  
       +                    rerody, lerody, slope, width, depth, 
  