anchor.diff:      integer :: MatrixSize,NonZeros,submatrixsize,ngen,i
       real(dp) :: area,Q01_mthrees,sheet_number
       character(len=60) :: sub_name
  -    real(dp) :: ha,hv,omega
  +
       
       sub_name = 'cap_flow_ladder'
       call enter_exit(sub_name,1)
   
  -    ha = 4.0
  -    hv = 1.0
  -    omega = 10.0
  -    call calc_cap_imped(ha,hv,omega)
   !     Number of non-zero entries in solution matrix. 
         NonZeros=3
         do i=2,cap_param%num_symm_gen
  
positive.code:   write ( *, '(a)' ) ' '
  
  write (*,*) 'C_3 = ',RHS2(N_Nodes+1),'+ i(',RHS2(2*N_Nodes+2),')'
  
  !
  !  Free memory.
  
negative.code:   write ( *, '(a)' ) '  Computed solution:'
    write ( *, '(a)' ) ' '
  
  write (*,*) 'C_1 = ',RHS1(2),'+ i(',RHS1(N_Nodes+2),')'
  deriv1(1) = (h(2)*RHS1(3) - h(1)*RHS1(1))/(2*dx)
  !c3 = (h(3) * H1R(3) - h(1) * H1R(1))/(2*dx);
  !c4 = (h(3) * H1I(3) - h(1) * H1I(1))/(2*dx);
  deriv1(2) = (h(2)*RHS1(N_nodes+4) - h(1)*RHS1(N_Nodes+2))/(2*dx)
  write (*,*) 'C_2 = ',deriv1(1),'+ i(',deriv1(2),')'
  
  !
  !  Free memory.
  