Changeset 6889

Timestamp:

01/03/11 15:06:28 (14 years ago)

Author:

jschierm

Message:

Added calculation of md.thickness_coeff to ThicknessCorrection.m.

File:

• issm/trunk/src/m/model/ThicknessCorrection.m (modified) (5 diffs)

issm/trunk/src/m/model/ThicknessCorrection.m

@@ -5 +5 @@
 %   by forcing the hydrostatic equilibrium.
 %   the thickness is modified as follows:
-%      thickness = coeff * thickness_observation + (1-coeff) * thickness_hydrostatic
+%      thickness = (1-coeff) * thickness_observation + coeff * thickness_hydrostatic
 %   where:
 %      coeff=(d/distance)^2;
@@ -20 +20 @@
 thickness=md.thickness;
 thickness_hydro=md.surface/(1-md.rho_ice/md.rho_water);
+thickness_coeff=zeros(size(md.thickness));
 
 %get nodes on ice sheet and on ice shelf
@@ -51 +52 @@
         if d>distance,
 
-                %if d > 15km, hydrostatique equilibrium
+                %if d > 15km, hydrostatic equilibrium
+                thickness_coeff(pos_shelf(i))=1;
                 thickness(pos_shelf(i))=thickness_hydro(pos_shelf(i));
 
@@ -57 +59 @@
 
                 %else: quadratic combination of hydrostatic equilibrium and observations
-                coeff=(d/distance)^2;
-                thickness(pos_shelf(i))=(1-coeff)*thickness(pos_shelf(i))+coeff*thickness_hydro(pos_shelf(i));
+                thickness_coeff(pos_shelf(i))=(d/distance)^2;
+                thickness(pos_shelf(i))=(1-thickness_coeff(pos_shelf(i)))*thickness(pos_shelf(i))+thickness_coeff(pos_shelf(i))*thickness_hydro(pos_shelf(i));
 
         end
@@ -71 +73 @@
 %change bed to take into account the changes in thickness
 md.thickness=thickness;
+md.thickness_coeff=thickness_coeff;
 md.bed=md.surface-md.thickness;