matlab-grid-generator/bin/pert.m

function [Tp] = pert(T);
% Calculates the perturbation about the mean of the first direction.

N=size(T);
T=reshape(T,[N(1) prod(N)/N(1)]);

q=T;q(isnan(T))=0;
qm=0*q+1;qm(find(q==0))=0;
n=sum(qm,1); n(find(n==0))=1;
tx=sum(q,1)./n;
n=sum(qm,2); n(find(n==0))=1;
ty=sum(q,2)./n;
[Ty,Tx]=ndgrid(ty,tx);

Tp=reshape( T-Tx.*qm ,N);
1	function [Tp] = pert(T);
2	% Calculates the perturbation about the mean of the first direction.
3
4	N=size(T);
5	T=reshape(T,[N(1) prod(N)/N(1)]);
6
7	q=T;q(isnan(T))=0;
8	qm=0*q+1;qm(find(q==0))=0;
9	n=sum(qm,1); n(find(n==0))=1;
10	tx=sum(q,1)./n;
11	n=sum(qm,2); n(find(n==0))=1;
12	ty=sum(q,2)./n;
13	[Ty,Tx]=ndgrid(ty,tx);
14
15	Tp=reshape( T-Tx.*qm ,N);