Changeset 16475

Timestamp:

10/21/13 09:37:01 (11 years ago)

Author:

Mathieu Morlighem

Message:

CHG: reverting to previous version

File:

• issm/trunk-jpl/src/m/geometry/NowickiProfile.m (modified) (2 diffs)

issm/trunk-jpl/src/m/geometry/NowickiProfile.m

@@ -12 +12 @@
 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 +36 @@
 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