matlab_class/gcmfaces_calc/calc_zonmedian_T.m

function [FLD]=calc_zonmedian_T(fld);
%object:    compute zonal median
%inputs:    fld is the field of interest
%output:    FLD is the zonal median field
%
%notes:     mygrid.LATS_MASKS is the set of quasi longitudinal lines along which
%               medians will be computed, as computed in gcmfaces_lines_zonal

global mygrid;

%initialize output:
n3=max(size(fld.f1,3),1); n4=max(size(fld.f1,4),1);
FLD=NaN*squeeze(zeros(length(mygrid.LATS_MASKS),n3,n4));

%use array format to speed up computation below:
fld=convert2array(fld.*mygrid.mskC); 
n1=size(fld,1); n2=size(fld,2); 
fld=reshape(fld,n1*n2,n3*n4); 

for iy=1:length(mygrid.LATS_MASKS); 

   %get list ofpoints that form a zonal band:
   mm=convert2array(mygrid.LATS_MASKS(iy).mmt);
   mm=find(~isnan(mm));

   %do the median along this band: 
   tmp1=nanmedian(fld(mm,:),1); 

   %store:
   FLD(iy,:,:)=reshape(tmp1,n3,n4);

end; 


1	gforget	1.3	function [FLD]=calc_zonmedian_T(fld);
2	gforget	1.2	%object: compute zonal median
3			%inputs: fld is the field of interest
4			%output: FLD is the zonal median field
5	gforget	1.3	%
6			%notes: mygrid.LATS_MASKS is the set of quasi longitudinal lines along which
7			% medians will be computed, as computed in gcmfaces_lines_zonal
8	gforget	1.1
9			global mygrid;
10
11			%initialize output:
12			n3=max(size(fld.f1,3),1); n4=max(size(fld.f1,4),1);
13	gforget	1.3	FLD=NaN*squeeze(zeros(length(mygrid.LATS_MASKS),n3,n4));
14	gforget	1.1
15			%use array format to speed up computation below:
16			fld=convert2array(fld.*mygrid.mskC);
17			n1=size(fld,1); n2=size(fld,2);
18			fld=reshape(fld,n1n2,n3n4);
19
20	gforget	1.3	for iy=1:length(mygrid.LATS_MASKS);
21	gforget	1.1
22			%get list ofpoints that form a zonal band:
23	gforget	1.3	mm=convert2array(mygrid.LATS_MASKS(iy).mmt);
24	gforget	1.1	mm=find(~isnan(mm));
25
26			%do the median along this band:
27			tmp1=nanmedian(fld(mm,:),1);
28
29			%store:
30			FLD(iy,:,:)=reshape(tmp1,n3,n4);
31
32			end;
33
34