%navier_testproblem    sets up reference Examples 7.1 to 7.4
%   IFISS scriptfile: DJS; 30 March 2005. 
% Copyright (c) 2005 D.J. Silvester, H.C. Elman, A. Ramage 
gohome
clear 
fprintf('\nspecification of reference Navier-Stokes problem.\n')
fprintf('\nchoose specific example (default is cavity)');
fprintf('\n     1  Channel domain')
fprintf('\n     2  Flow over a backward facing step')
fprintf('\n     3  Lid driven cavity')
fprintf('\n     4  Flow over a plate\n');
sn = default('',3);
if sn==1
   !copy .\stokes_flow\test_problems\poiseuille_flow.m .\stokes_flow\specific_flow.m
   !copy .\stokes_flow\test_problems\poiseuille_bc.m .\diffusion\specific_bc.m
   pipe_stokes, solve_navier
elseif sn==2
   !copy .\stokes_flow\test_problems\backwardstep_flow.m .\stokes_flow\specific_flow.m
   !copy .\stokes_flow\test_problems\backwardstep_bc.m .\diffusion\specific_bc.m
   longstep_stokes, solve_step_navier
elseif sn==3
   lid_model=default('cavity type leaky/tight/regularised 1/2/3 (regularised)',3); 
   if lid_model ==3    
      !copy .\stokes_flow\test_problems\regcavity_flow.m .\stokes_flow\specific_flow.m
   elseif lid_model ==2    
      !copy .\stokes_flow\test_problems\tightcavity_flow.m .\stokes_flow\specific_flow.m
   else
      !copy .\stokes_flow\test_problems\leakycavity_flow.m .\stokes_flow\specific_flow.m
   end
   !copy .\stokes_flow\test_problems\zero_bc.m .\diffusion\specific_bc.m
   square_stokes, solve_navier
elseif sn==4
   !copy .\stokes_flow\test_problems\plate_flow.m .\stokes_flow\specific_flow.m
   !copy .\stokes_flow\test_problems\plate_bc.m .\diffusion\specific_bc.m
   plate_stokes, solve_plate_navier
else
   error('reference problem datafile not found!')
end