Home > custom > colfilt_filters > calvfilt.m

calvfilt

PURPOSE

flag cells with ice and one ice free neighbour 1

SYNOPSIS

function nt=calvfilt(X,calvrate)

DESCRIPTION

 flag cells with ice and one ice free neighbour 1

CROSS-REFERENCE INFORMATION

This function calls:

This function is called by:

SOURCE CODE

0001 function nt=calvfilt(X,calvrate)
0002 % flag cells with ice and one ice free neighbour 1
0003 i=ceil(size(X,1)/2);
0004 nt=X(i,:);
0005 
0006 % find center cell with ice
0007 c = find(X(i,:)>0);
0008 % reduce ice thickness for cells which border on ice free cells
0009 nt(c(any(X(:,c)==0)))=nt(c(any(X(:,c)==0)))*(1-calvrate);