s_gu.m

基于Matlab的地震数据处理显示和测井数据显示于处理的小程序

function units=s_gu(seismic,mnem)
% Get units of measurement of the header with mnemonic "mnem" from seismic data set "seismic".
% If S4M.case_sensitive is set to false, the case of the header mnemonic is disregarded.
%
% Written by: E. R., December 30, 2000
% Last updated: December 30, 2000; compact header representation
%
%            units=s_gu(seismic,mnem)
% INPUT
% seismic    seismic data set
% mnem       header mnemonic
% OUTPUT
% units      string with units of measurements of header with mnemonic "mnem"

global S4M

if strcmpi(mnem,'trace_no')       % Implied header "trace_no"
   units='n/a';
   return
end

mnems=seismic.header_info(:,1);

if S4M.case_sensitive
   idx=find(ismember(mnems,mnem));
else 
   idx=find(ismember(lower(mnems),lower(mnem)));
end

if length(idx) == 1
   units=seismic.header_info{idx,2};
   return
end

% Handle error condition
if isempty(idx)
   disp([' Header "',mnem,'" not found. Available headers (in addition to the pseudo-header) are:'])
   disp(mnems')
else
   disp([' More than one header found: ',cell2str(mnems(idx),', ')])
end

error(' Abnormal termination.')