Index: /issm/trunk-jpl/src/m/geometry/NowickiProfile.m
===================================================================
--- /issm/trunk-jpl/src/m/geometry/NowickiProfile.m	(revision 16474)
+++ /issm/trunk-jpl/src/m/geometry/NowickiProfile.m	(revision 16475)
@@ -12,4 +12,5 @@
 delta = 0.1;          % ratio of water density and ice density -1
 hg    = 1;            % ice thickness at grounding line
+sea   = hg/(1+delta); % sea level
 lamda = 0.1;          % ration of deviatoric stress and water pressure
 beta  = 5;            % friction coefficient
@@ -35,4 +36,4 @@
 for i = ceil(numel(x)/2):numel(x)
 	h(i) = (x(i)/(4*(delta+1)*q)+hg^(-2))^(-0.5); % ice thickness for ice shelf from (3.1)
-	b(i) = -h(i)*(1/(1+delta));
+	b(i) = sea-h(i)*(1/(1+delta));
 end