anchor.diff:  c
         if (stations .ne. stnserod) then
            do s = 1, stations - 1
  -            if (width(s) .ne. 0.0 .and. 
  +            if (width(s) .ne. 0.d0 .and. 
        +          width(s) .gt. 2.d0 * depth(s)) then
                  forcefactor = rho * vel(s) ** 2.d0 / depth(s) 
        +              * (1.d0 - 2.d0 * depth(s) / width(s))
  
positive.code:       end do
        do s = 1, stnserod
           xdepth = rightdepth(s)
           if (curvature(s) .lt. 0.0) xdepth = leftdepth(s)
           if (xdepth .ne. 0.0) then
              tauwall(s) = tauwall(s) / xdepth
              spreaddelta_x(s) = tauwall(s) * (-1.d0) * sin(phi(s))
              spreaddelta_y(s) = tauwall(s) * cos(phi(s))
  
negative.code:                   leftdepth(s) = 0.d0
                 end if
              end if
              if (curvature(s) .lt. 0.d0) then
                 deln(s) = 0.5d0 * width(s) + Acs(s) * 2.d0 
       +                   / (depth(s) + rightdepth(s))
              else
  