anchor.diff:       +                      lag, latforce, dels, phi, curvature, 
        +                      depth, spreaddelta_x, spreaddelta_y,
        +                      rightdepth, leftdepth, xs, tauwall)
  -C
  -      integer stnserod, stations, s
  +      implicit none
  +C
  +      integer stnserod, stations, s, sp
         real*8 lambda(*), width(*), lag(*), 
        +       latforce(*), dels(*), phi(*), 
        +       curvature(*), depth(*), spreaddelta_x(*), 
        +       spreaddelta_y(*), rightdepth(*),
        +       leftdepth(*), xs(*),
        +       tauwall(*), gaussfactor, xstrt,
  -     +       xdel, xdepth, gaussian, xdest, tenlambda
  +     +       xdel, xdepth, gaussian, xdest, tenlambda, xtrmnt
  +      intrinsic cos, sin, exp
   c
         do s = 1, stnserod
            tauwall(s) = 0.d0
  
positive.code: *
  *
  *
        real function angle (y, x)
        real*8 y, x
        if (x .ne. 0.0 .or. y .ne. 0.0) then
           angle = atan2(y, x)
        else
           angle = 0.0
  
negative.code:      +       Acs(*), dels(*), deln(*), 
       +       forcefactor, xmagcurvep1, xmagcurve, latvel,
       +       delAcs
        intrinsic cos
  c
        if (stations .ne. stnserod) then
           do s = 1, stations - 1
  