Changeset 22262

Timestamp:

11/15/17 10:45:19 (7 years ago)

Author:

Mathieu Morlighem

Message:

CHG: trying to fix error

File:

• issm/trunk-jpl/src/c/cores/movingfront_core.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/cores/movingfront_core.cpp

@@ -40 +40 @@
         /* determine variables for extrapolation */
         num_extrapol_vars=0;
-        if(isstressbalance) num_extrapol_vars+=3;
+        if(isstressbalance){
+                if(domaintype==Domain3DEnum)
+                 num_extrapol_vars+=3;
+                else
+                 num_extrapol_vars+=2;
+        }
         if(ismasstransport) num_extrapol_vars+=1;
         if(isthermal && domaintype==Domain3DEnum) num_extrapol_vars+=1;
@@ -48 +53 @@
                 extrapol_vars[index]=VxEnum; index++;
                 extrapol_vars[index]=VyEnum; index++;
-                extrapol_vars[index]=VzEnum; index++;
+                if(domaintype==Domain3DEnum){
+                        extrapol_vars[index]=VzEnum; index++;
+                }
         }
-        if(ismasstransport){extrapol_vars[index]=ThicknessEnum; index++;}
+        if(ismasstransport){
+                extrapol_vars[index]=ThicknessEnum; index++;
+        }
         if(isthermal && domaintype==Domain3DEnum){
-                if(isenthalpy){extrapol_vars[index]=EnthalpyEnum;}
-                else{extrapol_vars[index]=TemperatureEnum;}
+                if(isenthalpy){
+                        extrapol_vars[index]=EnthalpyEnum;
+                }
+                else{
+                        extrapol_vars[index]=TemperatureEnum;
+                }
                 index++;
         }