anchor.diff:       +                    curvature, dels, Acs, deln,
        +                    rightdepth, leftdepth, depth, vel, 
        +                    delx, dely, transfactor, transslope)
  +      implicit none
   C
   C  xmagtransslope is the absolute value of the transverse bedslope (transslope)
   C  Acs is cross-sectional area of the half-channel at the inside of the bend
  
positive.code:       subroutine changeposition (stnserod, stations, lerody, rerody,  
       +                     spreaddelta_x, spreaddelta_y, delx, dely, 
       +                     depth, delta_x, delta_y)
  c
  C     change channel position:
  c
  
negative.code: *
        subroutine getcurv (stnserod, stations, delx, dely, dels, 
       +                    curvature)
        implicit none
  c     use law of cosines to find magnitude of angle, use cross 
  c     product to find sign of curvature
  c
        integer stnserod, stations, s
        real*8 delx(*), dely(*), dels(*), curvature(*), mag, sn, a, b, c,
       +       carg
        intrinsic acos
  c
        do s = 2, stnserod - 1
           a = dels(s) 
  