anchor.diff:       +       Acs(*), dels(*), deln(*), 
        +       forcefactor, xmagcurvep1, xmagcurve, latvel,
        +       delAcs
  +      intrinsic cos
   c
         if (stations .ne. stnserod) then
            do s = 1, stations - 1
  
positive.code:      +                      lag, latforce, dels, phi, curvature, 
       +                      depth, spreaddelta_x, spreaddelta_y,
       +                      rightdepth, leftdepth, xs, tauwall)
  C
        integer stnserod, stations, s
        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
  c
        do s = 1, stnserod
           tauwall(s) = 0.d0
  
negative.code: *
        subroutine forcelag (stnserod, stations, rho, vel, depth, width, 
       +                     latforce, lag, curvature, Acs, dels, deln)
        implicit none
  c
  C     calculate lateral force and downstream lag:
  c
  