anchor.diff:      ! 2*pi in this calculation
   #else
       ! use the external pressure solver
  -    call create_Hypre_A_vector(MPI_COMM_WORLD, hypre_grid, hypre_A, &
  +    call create_Hypre_A_matrix(MPI_COMM_WORLD, hypre_grid, hypre_A, &
             a, nx, ny, freesurfFac, dt, ierr)
   #endif
   
  
positive.code: 
  ! ---------------------------------------------------------------------------
  
  subroutine create_Hypre_A_vector(MPI_COMM_WORLD, hypre_grid, hypre_A, &
            a, nx, ny, freesurfFac, dt, ierr)
    implicit none
  
  
negative.code: #endif
  
    return
  end subroutine create_Hypre_A_matrix
  
  ! ---------------------------------------------------------------------------
  
  