construct_grid.m

function grid = construct_grid(params)
%function grid = construct_grid(params)
%
% function giving a simple common constructor syntax to different
% grid types. See the constructors of the grid types with a single
% params option for details on specifications
%
% required fields of params:
%    gridtype : 'rectgrid' or 'triagrid'

% Bernard Haasdonk 18.5.2007

if isfield(params,'gridtype')
  switch params.gridtype
   case 'rectgrid'
    grid = rectgrid(params);
   case 'onedgrid'
    grid = onedgrid(params);
   case 'triagrid'
    grid = triagrid(params);
    grid = set_boundary_types(grid, params);
   case 'none'
    % nothing to do
    grid = [];
   otherwise
    error('gridtype not known')
  end;
else
  grid = [];
end;