MITgcm_contrib/jmc_script/c_levs.m

function [ctl]=c_levs(mn,Mx,dc,dBug);
% [ctl]=c_levs(mn,Mx,dc,[option]);
%- return contour levels
%  dc > 0 : select contour-levels as multiple of "dc" within [mn,Mx]
%  dc < 0 : select aproximatively |dc| contour levels, within [mn,Mx]

% $Header: /u/gcmpack/MITgcm_contrib/jmc_script/c_levs.m,v 1.1 2008/10/09 01:08:27 jmc Exp $
% $Name:  $

if nargin < 3,
  error('Not enough argument (need at least 3)');
  return
end
if nargin < 4, dBug=0; end

if dc > 0,
  dd=dc;
else
  d0=(mn-Mx)/dc;
  if d0 <= 0,
    ctl=0; return
  end
  k0=log(d0)/log(10);
  kd=round(k0);
  pp=10.^kd;
  dd=round(d0/pp);
  if dd == 0,
    kd=floor(k0);
    pp=10.^kd;
    dd=round(d0/pp);
  end
  dd=pp*dd;
  if dBug > 0, fprintf(' -dc= %i ; d0 = %e ; kd,k0= %i,%f ; pp=%e ; dd=%e\n', ...
               -dc,d0,kd,k0,pp,dd); end
end

ctl=[ceil(mn/dd):floor(Mx/dd)]*dd;

return
1	function [ctl]=c_levs(mn,Mx,dc,dBug);
2	% [ctl]=c_levs(mn,Mx,dc,[option]);
3	%- return contour levels
4	% dc > 0 : select contour-levels as multiple of "dc" within [mn,Mx]
5	% dc < 0 : select aproximatively \|dc\| contour levels, within [mn,Mx]
6
7	% $Header: /u/gcmpack/MITgcm_contrib/jmc_script/c_levs.m,v 1.1 2008/10/09 01:08:27 jmc Exp $
8	% $Name: $
9
10	if nargin < 3,
11	error('Not enough argument (need at least 3)');
12	return
13	end
14	if nargin < 4, dBug=0; end
15
16	if dc > 0,
17	dd=dc;
18	else
19	d0=(mn-Mx)/dc;
20	if d0 <= 0,
21	ctl=0; return
22	end
23	k0=log(d0)/log(10);
24	kd=round(k0);
25	pp=10.^kd;
26	dd=round(d0/pp);
27	if dd == 0,
28	kd=floor(k0);
29	pp=10.^kd;
30	dd=round(d0/pp);
31	end
32	dd=pp*dd;
33	if dBug > 0, fprintf(' -dc= %i ; d0 = %e ; kd,k0= %i,%f ; pp=%e ; dd=%e\n', ...
34	-dc,d0,kd,k0,pp,dd); end
35	end
36
37	ctl=[ceil(mn/dd):floor(Mx/dd)]*dd;
38
39	return