function Neff=effectivepressure(md)
%EFFECTIVEPRESSURE - compute effective pressure
%
%   Usage:
%      Neff=effectivepressure(md)
%
%   Example:
%      Neff=effectivepressure(md)

Neff=md.materials.rho_ice*md.constants.g*md.geometry.thickness+md.materials.rho_ice*md.constants.g*md.geometry.bed;
pos=find(Neff<0);
Neff(pos)=0;