anchor.diff:  real(dp), intent(in) :: ha,hv,omega
   
   ! Local Variables:
  -integer, parameter :: N_nodes = 2000
  +integer, parameter :: N_nodes = 5000
   integer :: n,ldb
   integer :: iopt, info
   real(dp),  allocatable :: sparseval(:)
  
positive.code:   write ( *, '(a)' ) '  Computed solution:'
    write ( *, '(a)' ) ' '
  
  write (*,*) 'C_1 = ',RHS1(1),'+ i(',RHS1(N_Nodes+2),')'
  
  !
  !  Free memory.
  
negative.code:     integer :: MatrixSize,NonZeros,submatrixsize,ngen,i
      real(dp) :: area,Q01_mthrees,sheet_number
      character(len=60) :: sub_name
  
      
      sub_name = 'cap_flow_ladder'
      call enter_exit(sub_name,1)
  
  !     Number of non-zero entries in solution matrix. 
        NonZeros=3
        do i=2,cap_param%num_symm_gen
  