program Pr2;

{$APPTYPE CONSOLE}

uses
  Fil_XYZ_1u;

const
  {SaveFileName={'gamma.prn''u_matr3.prn''Uq.prn';  }
  SaveFileName1='u_zz1.prn';
  SaveFileName2='u_yy1.prn';
  SaveFileName3='u_xx1.prn';
 { SaveFileName4='u_Rj.prn';
  El_Potential=5;
  FileName1={'u_ii3.prn''U_iq.prn';
  FileName2={'u_jj3.prn''U_jq.prn';
  File1={'u_ii3.prn''J_iq.prn';
  File2={'u_jj3.prn''J_jq.prn'; }

procedure SetInitial;

var
  i, j, k : word;
  uu, t : real;

begin
 for k:=50{60} to Mk do
  begin
    t:=(k-1)/(Mk-1);
    uu:=1*(1-t)-1*t;
   for i:=1 to Mr do
    for j:=1 to Mz do
      u[i,j,k]:=uu;
  end;

  for k:=1 to 50{60} do
  begin
     t:=(k-1)/(50{60}-1);
    uu:=0*(1-t)+1*t;
    for i:=1 to Mr do
    for j:=1 to Mz do
      u[i,j,k]:=uu;
  end;


 {for k:=1 to Mk do
  begin
    t:=(k-1)/(Mk-1);
    uu:=-1*(1-t)+0*t;
   for i:=1 to Mr do
    for j:=1 to Mz do
      u[i,j,k]:=uu;
  end; }

  {for j:=1 to Mz-1 do
    for i:=1 to Mr-1 do
    begin
      t:=(j-1)/(Mz-1);
      gamma[i,j]:=t+0.1;
    end; }
end;{--------------------SetInitial--------------------}


procedure SaveMatrixToFile;

const
  Tab=chr(9);

var
  i, j, k : word;
  f : text;

begin
  assign(f,SaveFileName1);
  rewrite(f);

  for i:=1 to Mr  do
  begin
    for k:=1 to Mz  do
    begin
      write(f,u[i,k,137]:16 {gamma[i*2,(Mz+1-j)*2]*sqrt(sqr(u[(i-1)*2,(Mz+1-j)*2]-u[i*2,(Mz+1-j)*2])+
                        sqr(u[i*2,(Mz-j)*2]-u[i*2,(Mz+1-j)*2]))/0.02
                        u[i,Mz+1-j]:16});
      if i<Mr  then write(f,Tab);
    end;
    writeln(f,'');
  end;

  close(f);


  assign(f,SaveFileName2);
  rewrite(f);

  for i:=1 to Mr  do
  begin
    for j:=1 to Mz  do
    begin
      write(f,u[i,j,120]:16 {gamma[i*2,(Mz+1-j)*2]*sqrt(sqr(u[(i-1)*2,(Mz+1-j)*2]-u[i*2,(Mz+1-j)*2])+
                        sqr(u[i*2,(Mz-j)*2]-u[i*2,(Mz+1-j)*2]))/0.02
                        u[i,Mz+1-j]:16});
      if i<Mr  then write(f,Tab);
    end;
    writeln(f,'');
  end;

  close(f);


  assign(f,SaveFileName3);
  rewrite(f);

  for i:=1 to Mr  do
  begin
    for k:=1 to Mk  do
    begin
      write(f,u[i,21,k]:16 {gamma[i*2,(Mz+1-j)*2]*sqrt(sqr(u[(i-1)*2,(Mz+1-j)*2]-u[i*2,(Mz+1-j)*2])+
                        sqr(u[i*2,(Mz-j)*2]-u[i*2,(Mz+1-j)*2]))/0.02
                        u[i,Mz+1-j]:16});
      if i<Mr  then write(f,Tab);
    end;
    writeln(f,'');
  end;

  close(f);

 { assign(f,FileName1);
  rewrite(f);
  for i:=0 to 39 do
  begin
    for j:=30 downto 0 do
    begin
      write(f,(-u[1+i*2,j*2+1]+u[1+(i+1)*2,j*2+1]):16, tab);
{      if i<40 then write(f,Tab); for i:=2 to 40 do
  begin
    for j:=1 to 30 do
    begin
      write(f,(u[(i-1)*2,(Mz+1-j)*2]-u[i*2,(Mz+1-j)*2]):16);
      if i<40 then write(f,Tab);
    end;
    writeln(f,'');
  end;
  close(f);

  assign(f,FileName2);
  rewrite(f);
  for i:=0 to 40 do
  begin
    for j:=29 downto 0 do
    begin
      write(f,(u[1+i*2,j*2+1]-u[1+i*2,(j+1)*2+1]):16, tab);
{      if i<40 then write(f,Tab); for i:=1 to 40 do
  begin
    for j:=1 to 30 do
    begin
      write(f,(-u[i*2,(Mz+1-(j-1))*2]+u[i*2,(Mz+1-j)*2]):16);
      if i<40 then write(f,Tab);
    end;
    writeln(f,'');
  end;
  close(f);

 {----------Current-----------------

   assign(f,File1);
  rewrite(f);
  for i:=1 to Mr-1 do
  begin
    for j:=1 to Mz-1 do
    begin
      write(f,(gamma[i,(Mz+1-j)]*(-u[i,(Mz+1-j)]+u[i+1,(Mz+1-j)])/0.005):16);
      if i<Mr-1 then write(f,Tab);
    end;
    writeln(f,'');
  end;
  close(f);

  assign(f,File2);
  rewrite(f);
  for i:=1 to Mr-1 do
  begin
    for j:=1 to Mz-1 do
    begin
      write(f,(gamma[i,(Mz+1-j)]*(u[i,(Mz+1-j)]-u[i,(Mz+1-(j+1))])/0.005):16);
      if i<Mr-1 then write(f,Tab);
    end;
    writeln(f,'');
  end;
  close(f);  }

end;{--------------------SaveMatrixToFile--------------------}


{procedure SaveBoundaryToFile;

const
  Tab=chr(9);

var
  i, j : word;
  f : text;

begin
  assign(f,SaveFileName1); {U
  rewrite(f);
  for i:=32 to 53 do
    writeln(f,u[i,1]);
  close(f);

  assign(f,SaveFileName2);  {B
  rewrite(f);
  for i:=32 to 53 do
    writeln(f,u[i,16]);
  close(f);

  assign(f,SaveFileName3);  {L
  rewrite(f);
  for i:=1 to 16 do
    writeln(f,u[32,i]);
  close(f);

  assign(f,SaveFileName4);  {R
  rewrite(f);
  for i:=1 to 16 do
    writeln(f,u[53,i]);
  close(f);
end;{--------------------SaveBoundaryToFile--------------------}


Begin

   SetParameters;
   SetInitial;
   Compute;
  Current;
 { SaveBoundarytoFile;
  SaveMatrixToFile;    }
End.