expvel.m

来自「基于matlab的反演程序,用于地球物理勘探中射线追踪及偏移成像程序.」· M 代码 · 共 2 行

function expvel
% export the receiver coordinates, the first and the second layer velocities 
% at the receiver location 
recelev = refdata('get','recelev');
v1rec = 1000 * refdata('get','v1rec');
v2rec = 1000 * refdata('get','v2rec');
reccoord = recelev(1,:);
recnum = 1:length(reccoord);

% Ask for filename
[filename,path] = myuifile(gcf, '*.vel', 'Export Velocity', 'put');
if( filename == 0 )
   return;
end

ind = findstr(filename,'.vel');
if(strcmp(computer,'MAC2'))
%   fullfilename=filename(1:ind(1)-1);
   fullfilename=filename;
else
%   fullfilename=[path filename(1:ind(1)-1)];
   fullfilename=[path filename];
end

fid = fopen(fullfilename, 'w');
if( fid ~= -1 )
   fprintf(fid, '#Rec.   Coord.   V1(m/s)   V2(m/s)\n');
   for i=1:length(reccoord)
      count = fprintf(fid, '%5d %11.5f %11.5f %11.5f\n', recnum(i), ...
                      reccoord(i), v1rec(i), v2rec(i) );
   end
   fclose(fid);
else
   str = sprintf('Could not open file: %s',fullfilename);
   error(str);
end