matlab_class/sample_processing/example_smooth.m

function []=example_smooth(fld);

gcmfaces_global;

%use the field produced in sample_proccessing/example_bin_average.m
if myenv.verbose>0;
    gcmfaces_msg('* apply land mask (1/NaN) to gridded before smoothing');
end;
fld=fld.*mygrid.mskC(:,:,1);

%%%%%%%% isotropic diffusion %%%%%%%%

%choose smoothing scale: here 3 X grid spacing
if myenv.verbose>0;
    gcmfaces_msg('* set the smoothing scale to 3 grid points');
end;
distXC=3*mygrid.DXC; distYC=3*mygrid.DYC;

%do the smoothing:
if myenv.verbose>0;
    gcmfaces_msg(['* call diffsmoth2D : apply smoothing operator ' ...
        'that consists in time stepping a diffusion equation ' ...
        'with accordingly chosen diffusivity and duration. In particular ' ...
        'diffsmoth2D illustrates gradient computations (using calc_T_grad) ' ...
        'and convergence computations (using calc_UV_conv).']);
end;
fld_smooth=diffsmooth2D(fld,distXC,distYC);

%display results:
if myenv.verbose>0;
    gcmfaces_msg('* crude display of results in array format');
end;

figureL; gcmfaces_sphere(fld_smooth); caxis([-1 1]*0.4); colorbar;
title('smoothed data')

%%%%%%%% rotated diffusion %%%%%%%%

%choose smoothing scale: here 3 X grid spacing
if myenv.verbose>0;
    gcmfaces_msg('* set anisotropic and rotated smoothing');
end;
distLarge=3*sqrt(mygrid.RAC); distSmall=1*sqrt(mygrid.RAC);
fldRef=mygrid.YC;

%do the smoothing:
if myenv.verbose>0;
    gcmfaces_msg(['* call diffsmoth2Drotated : apply anisotropic smoothing operator ' ...
        'that acts preferentially along contours of a reference field (here latitude).']);
end;
fld_smooth2=diffsmooth2Drotated(fld,distLarge,distSmall,mygrid.YC);

%display results:
if myenv.verbose>0;
    gcmfaces_msg('* crude display of results in array format');
end;

figureL; gcmfaces_sphere(fld_smooth2); caxis([-1 1]*0.4); colorbar;
title('zonally smoothed data')


1	gforget	1.6	function []=example_smooth(fld);
2
3			gcmfaces_global;
4	gforget	1.1
5			%use the field produced in sample_proccessing/example_bin_average.m
6	gforget	1.2	if myenv.verbose>0;
7			gcmfaces_msg('* apply land mask (1/NaN) to gridded before smoothing');
8			end;
9	gforget	1.6	fld=fld.*mygrid.mskC(:,:,1);
10	gforget	1.1
11	gforget	1.3	%%%%%%%% isotropic diffusion %%%%%%%%
12
13	gforget	1.1	%choose smoothing scale: here 3 X grid spacing
14	gforget	1.2	if myenv.verbose>0;
15			gcmfaces_msg('* set the smoothing scale to 3 grid points');
16			end;
17	gforget	1.1	distXC=3mygrid.DXC; distYC=3mygrid.DYC;
18
19			%do the smoothing:
20	gforget	1.2	if myenv.verbose>0;
21			gcmfaces_msg(['* call diffsmoth2D : apply smoothing operator ' ...
22			'that consists in time stepping a diffusion equation ' ...
23			'with accordingly chosen diffusivity and duration. In particular ' ...
24			'diffsmoth2D illustrates gradient computations (using calc_T_grad) ' ...
25	gforget	1.4	'and convergence computations (using calc_UV_conv).']);
26	gforget	1.2	end;
27	gforget	1.6	fld_smooth=diffsmooth2D(fld,distXC,distYC);
28	gforget	1.1
29			%display results:
30	gforget	1.2	if myenv.verbose>0;
31			gcmfaces_msg('* crude display of results in array format');
32			end;
33	gforget	1.1
34	gforget	1.6	figureL; gcmfaces_sphere(fld_smooth); caxis([-1 1]*0.4); colorbar;
35	gforget	1.2	title('smoothed data')
36	gforget	1.3
37			%%%%%%%% rotated diffusion %%%%%%%%
38
39			%choose smoothing scale: here 3 X grid spacing
40			if myenv.verbose>0;
41			gcmfaces_msg('* set anisotropic and rotated smoothing');
42			end;
43			distLarge=3sqrt(mygrid.RAC); distSmall=1sqrt(mygrid.RAC);
44			fldRef=mygrid.YC;
45
46			%do the smoothing:
47			if myenv.verbose>0;
48			gcmfaces_msg(['* call diffsmoth2Drotated : apply anisotropic smoothing operator ' ...
49			'that acts preferentially along contours of a reference field (here latitude).']);
50			end;
51	gforget	1.6	fld_smooth2=diffsmooth2Drotated(fld,distLarge,distSmall,mygrid.YC);
52	gforget	1.3
53			%display results:
54			if myenv.verbose>0;
55			gcmfaces_msg('* crude display of results in array format');
56			end;
57
58	gforget	1.6	figureL; gcmfaces_sphere(fld_smooth2); caxis([-1 1]*0.4); colorbar;
59	gforget	1.3	title('zonally smoothed data')
60
61