function phi = upperlip(phi,m,n)
dims = repmat([m],[1 n]);

phi = ceil(phi);

for val = n:-1:1
  xind = find(phi==val);
  %xind = find(ceil(phi)==val);
  for i=xind(:)'
    v = phi(i);
    x = ind2coord(i,dims);
    Y = allrho1(x,m);
    Yind = coord2ind(Y,dims);
    badYi = find(v-phi(Yind) > 1);
    phi(Y(badYi)) = 0*phi(Y(badYi)) + v - 1;
  end
end

return