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

   end;
end; { Calculate_Vorticity_in_Boundary_Points }


    После этого пишутся процедуры для расчета функции вихря на грани-
цах области для различных способов – формулы Тома, Вудса, Йенсена, Та-
рунна или Полежаева и Грязнова.

procedure TomaFormulas;
{
   Назначение: Вычисление функции вихря на границах области реше-
               ния задачи по формулам Тома

}
var
   i, j : byte;
begin
   for i := 0 to Ny do
        begin
            W[i,0] := 2/sqr(Hx)*(Psi[i,0]-Psi[i,1]);
            W[i,Nx] := 2/sqr(Hx)*(Psi[I,Nx]-Psi[i,Nx-1]);
        end;
   for j := 0 to Nx do
        begin
            W[0,j] := 2/sqr(Hy)*(Psi[0,j]-Psi[1,j]);
            W[Ny,j] := -2/sqr(Hy)*(Psi[Ny,j]-Psi[Ny-1,j]);
        end;

end; { TomaFormulas }

procedure WoodFormulas;
{
   Назначение: Вычисление функции вихря на границах области реше-
               ния задачи по формулам Вудса

}
var
   i, j : byte;
begin
   for i := 0 to Ny do
        begin
            W[i,0] := 3/sqr(Hx)*(Psi[i,0]-Psi[i,1])-0.5*W[i,1];
                                        65