anchor.diff:          sparsecol(j+1) = counter + 1
        end if
    enddo ! do j
  -
  +call enter_exit(sub_name,2)
    
   end subroutine Mat_to_CC
   
  
positive.code:                stiff2(i+N_nodes+1,i) = -1.0 * omega
              end if
      enddo
      
      end subroutine Matrix
  
  !
  
negative.code: integer, intent(out) :: NonZeros
  !Local Variables:
  integer :: i, j, counter
  character(len=60) :: sub_name
  
      
  sub_name = 'Mat_to_CC'
  call enter_exit(sub_name,1)
  
  allocate(sparsecol(2*nn+3)) !allocation of sparsecol
  NonZeros = 0 !Initialisation for NonZeros
  