test/em_fire/heatflux/create_input_hfx.m

   1 wd=pwd;                                          % remember work directory
   2 cd ..
   3 startup                                          % set up search path
   4 cd(wd)                                           % return
   5 n1 = 420; n2=420;                                % grid dimensions
   6 width=50;                                        % bump width in grid cells
   7 [i,j]=ndgrid(1:n1,1:n2);                         % create index arrays
   8 hfx=1e6*exp(-((i-n1/2).^2 + 2*((j-n2/3).^2))/width^2);   % elliptic gaussian bump
   9 write_array_2d('input_hfx',hfx)                  % write for ideal.exe
  10 mesh(hfx)