% Script

clear all,close all,cd bin,path(path,pwd),cd ..

% Rør her
Rmax = 50;
M = 200; % Antal radier
N = 200; % Antal vinkler
gamma = 0; % Brug gamma=0 ved beregning af flowfelter

Reynolds = [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 15, 20]; % 1.58, 3.5, 4.88, 4.93, 7.52, 9.77];
% Alternativer: [ 0,1,3,6,10 ]   

% Rør ikke ved mere

psi = NumStreamFct;

cd Userdata/Flows

for i = 1:length(Reynolds),
   if length(Rmax)>1,
      Rm = Rmax(i);
   else
      Rm = Rmax;
   end
   
   if length(M)>1,
      Mi = M(i);
   else
      Mi = M;
   end
   
   if length(N)>1,
      Ni = N(i);
   else
      Ni = N;
   end
   
   Re = Reynolds(i);
   
   lg = loggrid(Rm,Mi,Ni,'Logarithmic grid',gamma);
	Set_grid(psi,lg);
   Set_Re(psi,Re);
   Solve(psi);
   ff = flowfield(psi);
   save(strcat('FlowRe',num2str(Re),'R',int2str(Rm), ...
      'M',int2str(Mi),'N',int2str(Ni)),'psi','Re','ff')
end

cd ..
cd ..

return