make_sources_file.m

来自「Spectral Element Method for wave propaga」· M 代码 · 共 13 行

ngll=5;
nex=20;
zmin=-1e3;
xmin=-1e3; xmax=1e3;

ns=nex*(ngll-1)+1;
clear table
table(:,1) = xmin+(xmax-xmin)/(ns-1)*(0:ns-1)';
table(:,2) = zmin;
table(:,3) = (table(:,1)-xmin)/4e3;  % = x / horizontal_speed

save('input_sources.tab', 'table', '-ascii')