matlab_class/gcmfaces_smooth/diffsmooth2D.m

function [FLD]=diffsmooth2D(fld,dxCsm,dyCsm);

%object: implementation (gforget@mit.edu) of a diffusive smoother (Weaver and Courtier 2001)

%input: fld     field to be smoothed (masked with NaN)
%       dxCsm,dyCsm     scale in first/second direction
%output:FLD     smoothed field

%asumption: dxCsm/dyCsm are assumed to be given at the positions of U/V points

global mygrid;

dxC=mygrid.DXC; dyC=mygrid.DYC; 
dxG=mygrid.DXG; dyG=mygrid.DYG; 
rA=mygrid.RAC;

%scale the diffusive operator:
tmp0=dxCsm./dxC; tmp0(isnan(fld))=NaN; tmp00=nanmax(tmp0);
tmp0=dyCsm./dyC; tmp0(isnan(fld))=NaN; tmp00=max([tmp00 nanmax(tmp0)]);
smooth2D_nbt=tmp00;
smooth2D_nbt=ceil(1.1*2*smooth2D_nbt^2);

smooth2D_dt=1;
smooth2D_T=smooth2D_nbt*smooth2D_dt;
smooth2D_Kux=dxCsm.*dxCsm/smooth2D_T/2;
smooth2D_Kvy=dyCsm.*dyCsm/smooth2D_T/2;

%time-stepping loop:
FLD=fld;

for it=1:smooth2D_nbt;
%    if mod(it,ceil(smooth2D_nbt/50))==0; fprintf([num2str(it) '/' num2str(smooth2D_nbt) ' done\n']); end;

    [dTdxAtU,dTdyAtV]=calc_T_grad(FLD,0);
    tmpU=dTdxAtU.*smooth2D_Kux;
    tmpV=dTdyAtV.*smooth2D_Kvy;
    [fldDIV]=calc_UV_div(tmpU,tmpV);
    dFLDdt=smooth2D_dt*fldDIV./rA;
    FLD=FLD-dFLDdt;

end;


1	function [FLD]=diffsmooth2D(fld,dxCsm,dyCsm);
2
3	%object: implementation (gforget@mit.edu) of a diffusive smoother (Weaver and Courtier 2001)
4
5	%input: fld field to be smoothed (masked with NaN)
6	% dxCsm,dyCsm scale in first/second direction
7	%output:FLD smoothed field
8
9	%asumption: dxCsm/dyCsm are assumed to be given at the positions of U/V points
10
11	global mygrid;
12
13	dxC=mygrid.DXC; dyC=mygrid.DYC;
14	dxG=mygrid.DXG; dyG=mygrid.DYG;
15	rA=mygrid.RAC;
16
17	%scale the diffusive operator:
18	tmp0=dxCsm./dxC; tmp0(isnan(fld))=NaN; tmp00=nanmax(tmp0);
19	tmp0=dyCsm./dyC; tmp0(isnan(fld))=NaN; tmp00=max([tmp00 nanmax(tmp0)]);
20	smooth2D_nbt=tmp00;
21	smooth2D_nbt=ceil(1.12smooth2D_nbt^2);
22
23	smooth2D_dt=1;
24	smooth2D_T=smooth2D_nbt*smooth2D_dt;
25	smooth2D_Kux=dxCsm.*dxCsm/smooth2D_T/2;
26	smooth2D_Kvy=dyCsm.*dyCsm/smooth2D_T/2;
27
28	%time-stepping loop:
29	FLD=fld;
30
31	for it=1:smooth2D_nbt;
32	% if mod(it,ceil(smooth2D_nbt/50))==0; fprintf([num2str(it) '/' num2str(smooth2D_nbt) ' done\n']); end;
33
34	[dTdxAtU,dTdyAtV]=calc_T_grad(FLD,0);
35	tmpU=dTdxAtU.*smooth2D_Kux;
36	tmpV=dTdyAtV.*smooth2D_Kvy;
37	[fldDIV]=calc_UV_div(tmpU,tmpV);
38	dFLDdt=smooth2D_dt*fldDIV./rA;
39	FLD=FLD-dFLDdt;
40
41	end;
42
43