Changeset 18843

Timestamp:

11/25/14 02:26:42 (10 years ago)

Author:

jbondzio

Message:

CHG: Use element boundaries as ice front instead of ZeroLevelset. Change this back when consistent stress balance treatment of front elements is given.

File:

• issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp (modified) (4 diffs)

issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp

@@ -1625 +1625 @@
         IssmDouble gravity     = element->GetMaterialParameter(ConstantsGEnum);
         element->GetVerticesCoordinates(&xyz_list);
-        //element->ZeroLevelsetCoordinates(&xyz_list_front,xyz_list,MaskIceLevelsetEnum);
         element->GetIcefrontCoordinates(&xyz_list_front,xyz_list,MaskIceLevelsetEnum);
         element->NormalSection(&normal[0],xyz_list_front);
@@ -2077 +2076 @@
         IssmDouble gravity     = element->GetMaterialParameter(ConstantsGEnum);
         element->GetVerticesCoordinates(&xyz_list);
-        element->ZeroLevelsetCoordinates(&xyz_list_front,xyz_list,MaskIceLevelsetEnum);
+        element->GetIcefrontCoordinates(&xyz_list_front,xyz_list,MaskIceLevelsetEnum);
         element->NormalSection(&normal[0],xyz_list_front);
 
@@ -2595 +2594 @@
         IssmDouble gravity   = element->GetMaterialParameter(ConstantsGEnum);
         element->GetVerticesCoordinates(&xyz_list);
-        element->ZeroLevelsetCoordinates(&xyz_list_front,xyz_list,MaskIceLevelsetEnum);
+        element->GetIcefrontCoordinates(&xyz_list_front,xyz_list,MaskIceLevelsetEnum);
         element->NormalSection(&normal[0],xyz_list_front);
 
@@ -4167 +4166 @@
         /*Retrieve all inputs and parameters*/
         element->GetVerticesCoordinates(&xyz_list);
-        element->ZeroLevelsetCoordinates(&xyz_list_front,xyz_list,MaskIceLevelsetEnum);
+        element->GetIcefrontCoordinates(&xyz_list_front,xyz_list,MaskIceLevelsetEnum);
         element->NormalSection(&normal[0],xyz_list_front);
         Input* surface_input  = element->GetInput(SurfaceEnum); _assert_(surface_input);