matlab_class/sample_processing/example_bin_average.m

function [fld]=example_bin_average();
% EXAMPLE_BIN_AVERAGE generates a sample of random
%    numbers distributed over the globe and 
%    bin averages to the default grid.
 
%%%%%%%%%%%%%%%%%
%load grid:
%%%%%%%%%%%%%%%%%

gcmfaces_global;

if isempty(mygrid);
  grid_load;
end;

if myenv.verbose>0;
    gcmfaces_msg('===============================================');
    gcmfaces_msg(['*** entering example_bin_average: generate ' ... 
        'randomly distributed data samples (position and value), ' ... 
        'and bin average this data set to a gcmfaces grid'],'');
end;

warning('off','MATLAB:dsearch:DeprecatedFunction');
warning('off','MATLAB:delaunay:DuplicateDataPoints');

%%%%%%%%%%%%%%%%%%%%%%%
%generate random data

if myenv.verbose>0;
    gcmfaces_msg('* generate random data');
end;
nn=1e6;
lat=(rand(nn,1)-0.5)*2*90;
lon=(rand(nn,1)-0.5)*2*180; 
%needed for 0-360 longitude convention
if mygrid.nFaces==1;
    xx=find(lon<0);lon(xx)=lon(xx)+360;
end;
sample=(rand(nn,1)-0.5)*2;

%%%%%%%%%%%%%%%%%%%%%%%
%generate delaunay triangulation

if myenv.verbose>0;
    gcmfaces_msg('* call gcmfaces_bindata : generate delaunay triangulation');
end;
gcmfaces_bindata;

%%%%%%%%%%%%%%%%%%%%%%%%
%bin average random data

if myenv.verbose>0;
    gcmfaces_msg('* call gcmfaces_bindata : bin average data');
end;
fld=gcmfaces_bindata(lon,lat,sample);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%quick display
if myenv.verbose>0;
    gcmfaces_msg('* display results on sphere');
end;
figureL; set(gca,'FontSize',16);
gcmfaces_sphere(fld); axis xy; caxis([-1 1]*0.4); colorbar;
title('bin averaged random data');

if myenv.verbose>0;
    gcmfaces_msg('*** leaving example_bin_average');
    gcmfaces_msg('===============================================');
end;


1	function [fld]=example_bin_average();
2	% EXAMPLE_BIN_AVERAGE generates a sample of random
3	% numbers distributed over the globe and
4	% bin averages to the default grid.
5
6	%%%%%%%%%%%%%%%%%
7	%load grid:
8	%%%%%%%%%%%%%%%%%
9
10	gcmfaces_global;
11
12	if isempty(mygrid);
13	grid_load;
14	end;
15
16	if myenv.verbose>0;
17	gcmfaces_msg('===============================================');
18	gcmfaces_msg(['*** entering example_bin_average: generate ' ...
19	'randomly distributed data samples (position and value), ' ...
20	'and bin average this data set to a gcmfaces grid'],'');
21	end;
22
23	warning('off','MATLAB:dsearch:DeprecatedFunction');
24	warning('off','MATLAB:delaunay:DuplicateDataPoints');
25
26	%%%%%%%%%%%%%%%%%%%%%%%
27	%generate random data
28
29	if myenv.verbose>0;
30	gcmfaces_msg('* generate random data');
31	end;
32	nn=1e6;
33	lat=(rand(nn,1)-0.5)290;
34	lon=(rand(nn,1)-0.5)2180;
35	%needed for 0-360 longitude convention
36	if mygrid.nFaces==1;
37	xx=find(lon<0);lon(xx)=lon(xx)+360;
38	end;
39	sample=(rand(nn,1)-0.5)*2;
40
41	%%%%%%%%%%%%%%%%%%%%%%%
42	%generate delaunay triangulation
43
44	if myenv.verbose>0;
45	gcmfaces_msg('* call gcmfaces_bindata : generate delaunay triangulation');
46	end;
47	gcmfaces_bindata;
48
49	%%%%%%%%%%%%%%%%%%%%%%%%
50	%bin average random data
51
52	if myenv.verbose>0;
53	gcmfaces_msg('* call gcmfaces_bindata : bin average data');
54	end;
55	fld=gcmfaces_bindata(lon,lat,sample);
56
57	%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
58	%quick display
59	if myenv.verbose>0;
60	gcmfaces_msg('* display results on sphere');
61	end;
62	figureL; set(gca,'FontSize',16);
63	gcmfaces_sphere(fld); axis xy; caxis([-1 1]*0.4); colorbar;
64	title('bin averaged random data');
65
66	if myenv.verbose>0;
67	gcmfaces_msg('*** leaving example_bin_average');
68	gcmfaces_msg('===============================================');
69	end;
70
71
72
73