Фиктивный аргумент INTENT(IN) в контексте определения переменной

Question

Я получаю сообщение об ошибке в моей подпрограмме, когда я запускаю свой код.

Этот код взят из упражнения из книги Кинкейда и Чейни о методах Гаусса Зейделя для решения эллиптических уравнений в частных производных.

Сообщение об ошибке:

фиктивный аргумент 'u' с INTENT(IN) в контексте определения переменной (присваивания) в (1).

Я ссылаюсь на (1) в коде ниже. Как я могу исправить подпрограмму, чтобы сообщение об ошибке не появлялось?

subroutine seidel(ax,ay,nx,ny,h,itmax,u)     
  real, dimension(0:nx,0:ny), intent(in) :: u        
  real, intent(in) :: ax,ay, h                 
  integer, intent(in) :: nx, ny, itmax         
  integer:: i,j,k

  do  k = 1,itmax      
    do  j = 1,ny-1     
      y = ay + real(j)*h
      do  i = 1,nx-1   
        x = ax + real(i)*h      
        v = u(i+1,j) + u(i-1,j) + u(i,j+1) + u(i,j-1) 
        u(i,j) = (v - h*h*g(x,y))/(4.0 - h*h*f(x,y))           (1)
      end do
    end do
  end do
end subroutine seidel

Answer 1

intent(in) это обещание компилятору, что подпрограмма не будет пытаться обновить значение аргумента. Код, который вы показываете разрывы, которые обещают в строке:

u(i,j) = (v - h*h*g(x,y))/(4.0 - h*h*f(x,y))

Исправьте это, не нарушая свои обещания или, возможно, более уместно в этом случае, сделав намерение inout, как это

real, dimension(0:nx,0:ny), intent(inout) :: u        

inout сообщает компилятору, что подпрограмме будет передан аргумент, и он может внести в него изменения.

(Я думаю, что это, вероятно, дубликат, но пока не могу его найти.)