anchor.diff:  integer :: nflt, nslip, ndsp, ndsp_dof, nlos, nsts, nnbr, nfixed, nfree
   integer :: i, j, ierr, icomp, iflt, inbr
   double precision, allocatable :: A(:,:), b(:), x(:), atmp(:,:), btmp(:)
  -double precision :: damping_squared, smoothing_squared
  +double precision :: damping_squared, smoothing_squared, factor
   logical, allocatable :: isSlipFixed(:)
  +double precision, parameter :: spy = 60.0d0*60.0d0*24.0d0*365.25d0
   
   
   ! Initialize solution
  
positive.code:     enddo
  endif
  
  ! Load solution into euler pole vector
  do i = 1,npoles
      euler_pole(i,1) = x(j+i-1)
      euler_pole(i,2) = x(j+i  )
      euler_pole(i,3) = x(j+i+1)
  enddo
  
  return
  
negative.code:                   gf_euler, &
                    displacement, &
                    disp_components, &
                    input_disp_unit, &
                    los, &
                    prestress, &
                    cov_matrix, &
  