femwind/params_defaults.m

   1 function params=params_defaults
   2     params.run_fortran=0;
   3     params.run_matlab=1;
   4     params.femwind_fortran_test=0;
   5     params.test_fortran=0;
   6     params.level=1;  % starting from 1
   7     params.graphics=2;  % -1: none 1=basic, 2=all
   8     params.expand=1.2;  % exponential grid expansion in the vertical
   9     params.mesh_top=1000; % if given, ignore params_expand
  10     params.sc_all=[1,2]; % mesh refinements for tests at multiple scales
  11     params.sc2_all=[1,2,4];  % additional factors for horizonal mesh extent
  12     params.nelem3=[50,50,10]; % base size in elements, horizontal=2*odd
  13     params.h=[30,30,2]; % base mesh spacing before scaling
  14     params.a=[1 1 1]; % penalty factors in x y z directions
  15     params.initial_wind='log'; % or uniform
  16     params.roughness_height=0.5;
  17     params.terrain_shape='hill'; % terrain for add_terrain_to_mesh
  18     params.terrain_top='shift'; % mesh top treatment for add_terrain_to_mesh
  19     params.terrain_height=0.1; % terrain height as part of domain height
  20     params.solver='2-level' ; % see sparse_solve.m
  21     params.maxit=50; % max iterations
  22     params.maxit_coarse=8; % 2 smoothing, coarse, 2 smoothing, coarse, 2 smoothing
  23     params.coarsest_iter=100; % 0 = direct solver n>0 number of iterations
  24     params.nsmooth=3; % smoothing iterations before correcton
  25     params.coarsening='2 linear';
  26     % params.coarse_P='variational';
  27     params.coarse_K='assembly';
  28     params.P_by_x=1;  % prolongation by geometrically linear interpolation
  29     params.smoothing='vertical sweeps';
  30     % params.smoothing='3D red-black';
  31     params.restol=1e-6;
  32     params.exact=0; % compare with exact solution to compute error
  33     params.slice=0.5; % vertical y slice of error to display, 0 to 1
  34     params.err_slice_fig=12; % figure number for residual slice
  35     params.res_slice_fig=13; % figure number for error slice
  36     params.iterations_fig=14; % figure number for iterations progress
  37     params.maxaspect=3;  % do not coarsen vertically if vertical layer is too thick
  38     params.minaspect=1/3; % do not coarsen horizontally if the layer is too thin
  39     params.levels=15;
  40     params.debug_level=0;
  41     params.apply_coarse_boundary_conditions=1;
  42     params.nsmooth_coarse=2;
  43     params.save_files=0; % save progress levels=3, workspace=2 params only=1
  44     params.save_file_prefix='femwind';
  45     %Define Streamline Starting Points: Defined in terms of scale*nelem
  46     params.in_height_stream = [175];
  47     params.time_stream  = 0;
  48     params.st_contour = 1; %Produce contour planes in streamlines plot 0 off, 1 on
  49     params.st_quiver = 1; %Produce vectors along streamlines 0 off, 1 on
  50 end