global_ocean.90x40x15/diags_matlab/mit_getparm.m

function y = mit_getparm(fname,pname);
%function y = mit_getparm(fname,pname);
  
  y = [];
  [fp, msg] = fopen(fname,'r');
  if fp > 0
    str{1} = [];
    notfound = 0;
    while isempty(mystrfind(str{1},pname))
      str{1} = fgetl(fp); 
      if ~ischar(str{1}); % this catches the end of file (not very clean)
        notfound = 1;
        break; 
      end
      % clear again if commented
      if ~isempty(mystrfind(str{1},'#')) 
        str{1} = [];
      end
    end
    if notfound
      disp(['Warning: ' pname ' not found in ' fname])
      y=[];
%      disp(['         setting ' pname ' to zero'])
%      y = 0;
    else
      teststr = [];
      % find the termination of parameter pname 
      % (next line with an '=' sign
      % or with a namelist termination character '&' or '/')
      n = 1;
      while ( isempty(mystrfind(teststr,'=')) & ...
              isempty(mystrfind(teststr,'&')) & ...
              isempty(mystrfind(teststr,'/')) )
        n = n + 1;
        teststr = fgetl(fp);
        str{n}  = teststr;
        if ~ischar(teststr); break; end
      end
      eqind = findstr(str{1},'=');
      y = str{1}(eqind+1:end-1);
      % check whether it is a string in quotes or something else
      quotes = findstr(y,'''');
      if ~isempty(quotes)
        y(quotes(2):end) = [];
        y(1:quotes(1)) = [];
      else
        if ~( strcmpi(y,'.TRUE.') | strcmpi(y,'.FALSE.') )
          y = convert_str2num(str{1}(eqind+1:end-1));
          for k=2:n-1
            y = [y; convert_str2num(str{k}(eqind+1:end-1))];
          end
% $$$     y = str2num(str{1}(eqind+1:end-1))';
% $$$     for k=2:n-1
% $$$       y = [y; str2num(str{k}(eqind+1:end-1))'];
% $$$     end
        end
      end
    end
    fclose(fp);
  else
    error([fname ' could not be opened: ' msg])
  end

  return

function y = convert_str2num(str)
  
% take care of the n*value format
  stars = findstr(str,'*');
  if isempty(stars)
    y = str2num(str)';
  else
    if length(stars)==1
      y=repmat(str2num(str(stars+1:end)),[str2num(str(1:stars-1)) 1]);
    else
      warning('The format n*x1,m*x2, ... cannot be handled')
% $$$       for k=1:length(stars)
% $$$   % find the beginnin and termination of the set
% $$$       end
    end
  end
  
  return

function [ind] = mystrfind(str,pattern)
  
  if length(pattern) > length(str)
    ind = [];
    return
  else
    ind = findstr(str,pattern);
  end
1	function y = mit_getparm(fname,pname);
2	%function y = mit_getparm(fname,pname);
3
4	y = [];
5	[fp, msg] = fopen(fname,'r');
6	if fp > 0
7	str{1} = [];
8	notfound = 0;
9	while isempty(mystrfind(str{1},pname))
10	str{1} = fgetl(fp);
11	if ~ischar(str{1}); % this catches the end of file (not very clean)
12	notfound = 1;
13	break;
14	end
15	% clear again if commented
16	if ~isempty(mystrfind(str{1},'#'))
17	str{1} = [];
18	end
19	end
20	if notfound
21	disp(['Warning: ' pname ' not found in ' fname])
22	y=[];
23	% disp([' setting ' pname ' to zero'])
24	% y = 0;
25	else
26	teststr = [];
27	% find the termination of parameter pname
28	% (next line with an '=' sign
29	% or with a namelist termination character '&' or '/')
30	n = 1;
31	while ( isempty(mystrfind(teststr,'=')) & ...
32	isempty(mystrfind(teststr,'&')) & ...
33	isempty(mystrfind(teststr,'/')) )
34	n = n + 1;
35	teststr = fgetl(fp);
36	str{n} = teststr;
37	if ~ischar(teststr); break; end
38	end
39	eqind = findstr(str{1},'=');
40	y = str{1}(eqind+1:end-1);
41	% check whether it is a string in quotes or something else
42	quotes = findstr(y,'''');
43	if ~isempty(quotes)
44	y(quotes(2):end) = [];
45	y(1:quotes(1)) = [];
46	else
47	if ~( strcmpi(y,'.TRUE.') \| strcmpi(y,'.FALSE.') )
48	y = convert_str2num(str{1}(eqind+1:end-1));
49	for k=2:n-1
50	y = [y; convert_str2num(str{k}(eqind+1:end-1))];
51	end
52	% $$$ y = str2num(str{1}(eqind+1:end-1))';
53	% $$$ for k=2:n-1
54	% $$$ y = [y; str2num(str{k}(eqind+1:end-1))'];
55	% $$$ end
56	end
57	end
58	end
59	fclose(fp);
60	else
61	error([fname ' could not be opened: ' msg])
62	end
63
64	return
65
66	function y = convert_str2num(str)
67
68	% take care of the n*value format
69	stars = findstr(str,'*');
70	if isempty(stars)
71	y = str2num(str)';
72	else
73	if length(stars)==1
74	y=repmat(str2num(str(stars+1:end)),[str2num(str(1:stars-1)) 1]);
75	else
76	warning('The format nx1,mx2, ... cannot be handled')
77	% $$$ for k=1:length(stars)
78	% $$$ % find the beginnin and termination of the set
79	% $$$ end
80	end
81	end
82
83	return
84
85	function [ind] = mystrfind(str,pattern)
86
87	if length(pattern) > length(str)
88	ind = [];
89	return
90	else
91	ind = findstr(str,pattern);
92	end