Changeset 11314

Timestamp:

02/03/12 10:55:53 (13 years ago)

Author:

seroussi

Message:

added error message if unphysical water friction value

File:

• issm/trunk-jpl/src/c/objects/Elements/Penta.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/objects/Elements/Penta.cpp ¶

@@ -4271 +4271 @@
         if(converged){
                 /*Convert enthalpy into temperature and water fraction*/
-                for(i=0;i<numdof;i++) matpar->EnthalpyToThermal(&temperatures[i],&waterfraction[i],values[i],pressure[i]);
+                for(i=0;i<numdof;i++){
+                        matpar->EnthalpyToThermal(&temperatures[i],&waterfraction[i],values[i],pressure[i]);
+                        if(waterfraction[i]<0) _error_("Negative water fraction found in solution vector");
+                        if(waterfraction[i]>1) _error_("Water fraction >1 found in solution vector");
+                }
                         
                 this->inputs->AddInput(new PentaP1Input(EnthalpyEnum,values));