gen_dune_params.m

function fn = gen_dune_params(params)
%function fn = gen_dune_params(params)
%
% function generating new params file from default file for dune-detailed
% simulations. The new file is written to $RBMATLABTEMP/burgers_fem
% the default settings are taken from params.
%
% required fields of params:
%         dune_params_fn: the filename of the params-file with default
%         settings. Further the simulation parameters are required:
%         udir_height, k, cx, cy, udir_amplitude, udir_xscale, T, nt,
%         max_newton_steps, newton_eps, dune_grid_fn

% Bernard Haasdonk 18.1.2007
  
  
  if ~exist(fullfile(rbmatlabresult,'burgers_fem'))
    mkdir(rbmatlabresult,'burgers_fem');
  end
  
  [p,f,e] = fileparts(params.dune_params_fn);
  new_params_fn = fullfile(rbmatlabtemp,'burgers_fem',[f,'_new.dat']);
  pnames = {'mu4','mu1','mu2','mu3','mu5','mu6',...
            'EndTime','deltaT','TimeStepNumber',...
            'maxNewtonStep','IteratedNewtonSolverEpsilon',...
            'Task','GrapeDisplay','VisualizeRBSpace','VisualizeEachNewtonStep',...
            'VisualizeEachTimeStep',...
            'OutputPath','OutputPrefix','Grid'};
  
  % generate outputprefix for snapshot-files, e.g.
  % mu1_0.1,mu2_0,mu3_0,mu4_0.2,mu5_1,mu6_1
  mu = get_mu(params);
  outputprefix = ...
      ['mu1_',num2str(mu(1)),...
       ',mu2_',num2str(mu(2)),...
       ',mu3_',num2str(mu(3)),...
       ',mu4_',num2str(mu(4)),...
       ',mu5_',num2str(mu(5)),...
       ',mu6_',num2str(mu(6)) ...
      ];
  
  pvalues = {[num2str(params.udir_height),' % = mu_4'],...
             [num2str(params.k),' % = mu_1'],...
             [num2str(params.cx),' % = mu_2'],...
             [num2str(params.cy),' % = mu_3'],...
             [num2str(params.udir_amplitude),' % = mu_5'],...
             [num2str(params.udir_xscale),' % = mu_6'],...
             num2str(params.T), ...
             num2str(params.T/params.nt), ...
             num2str(params.nt), ...
             num2str(params.max_newton_steps),...
             num2str(params.newton_eps),...
             '1  % == Task',...
             '0  % == GrapeDisplay',...
             '0  % == VisualizeRBSpace',...
             '0  % == Visualize each newton step',...
             '0  % == Visualize each time step',...
             [fullfile(rbmatlabtemp,'burgers_fem'),' % = outputpat for snapshots'],...
             outputprefix,...
             params.dune_grid_fn
            };
  
  disp(['generating new params file']);
  disp(new_params_fn);
  disp(['from old ',params.dune_params_fn]);
  change_params_file(params.dune_params_fn, new_params_fn, ...
                     pnames, pvalues);
  
  fn = new_params_fn;
% TO BE ADJUSTED TO NEW SYNTAX
%| \docupdate