Changeset 25289

Timestamp:

07/18/20 23:17:43 (5 years ago)

Author:

adhikari

Message:

CHG: skip computation for non-icy elements

File:

• issm/trunk-jpl/src/c/classes/Elements/Tria.cpp (modified) (3 diffs)

issm/trunk-jpl/src/c/classes/Elements/Tria.cpp

@@ -4957 +4957 @@
 
         /*gia solution parameters:*/
-        IssmDouble lithosphere_thickness,mantle_viscosity;
+        IssmDouble lithosphere_thickness,mantle_viscosity,ice_mask;
 
         /*output: */
@@ -4971 +4971 @@
         this->parameters->FindParam(&gsize,MeshNumberofverticesEnum);
         this->parameters->FindParam(&yts,ConstantsYtsEnum);
-
+        
         /*recover gia solution parameters: */
         int cross_section_shape;
@@ -4986 +4986 @@
         IssmDouble mantle_density            = FindParam(MaterialsMantleDensityEnum);
         IssmDouble rho_ice                   = FindParam(MaterialsRhoIceEnum);
+
+        /*recover ice mask*/ 
+        Input2* ice_mask_input=this->GetInput2(MaskIceLevelsetEnum); 
+        ice_mask_input->GetInputAverage(&ice_mask);
+
+        /*early return if we are NOT on an icy element:*/
+        if (ice_mask==1) return; 
 
         /*pull thickness averages! */