Введение в технологию программирования. Часть 3. Проект OLYMPUS-D. Коржов Е.Н - 63 стр.

procedure Stepon;
{
   Назначение: вычисление неизвестных функций на одном шаге
}
begin
   Calculate_Vorticity_in_Inner_Points;
   Calculate_StreamFunction_in_Inner_Points;
   Calculate_Vorticity_in_Boundary_Points;
end; { Stepon }

    Можно выполнить разработку процедур, операторы которых образуют
тело процедуры Stepon. Их алгоритмы уже обсуждались в учебной литера-
туре и для случая использования прямоугольных областей решения разо-
браны досконально [22–23]. Тексты этих процедур располагаются в разде-
ле implementation этого же модуля до текста процедуры Stepon.

procedure Calculate_Vorticity_in_Inner_Points;
{
   Назначение: Вычисление функции вихря во внутренних точках области
                   решения задачи
}
var
   R1, R2, R3, R4 : real;
begin
   for i := 1 to Ny-1 do
      for j := 1 to Nx-1 do
          begin
              R1 := Ht_Hx2*(Re_1-Hx_Hy_4*(Psi[i+1,j]-Psi[i-1,j]));
              R2 := Ht_Hx2*(Re_1+Hx_Hy_4*(Psi[i+1,j]-Psi[i-1,j]));
              R3 := Ht_Hx2*(Re_1-Hx_Hy_4*(Psi[i,j+1]-Psi[i,j-1]));
              R4 := Ht_Hx2*(Re_1+Hx_Hy_4*(Psi[i,j+1]-Psl[i,j-1]));
              W[1,i,j] := K_W*W[0,i,j]+R1*W[0,i,j+1]+R2*W[0,i,j-1] +
                                       R3*W[0,i+1,j]+R4*W[0,i-1,j];
          end;
end; { Calculate_Vorticity_in_Inner_Points }

    Здесь использованы коэффициенты K_W, Ht_Hx2, Hx_Hy_4, относя-
щиеся к числу глобальных данных и вычисленные в процедуре Auxval мо-
дуля PROLOG.

                                  63