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	cnh	1.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);