Changeset 1310 for issm/trunk/src/c/ModelProcessorx/Thermal/CreateElementsNodesAndMaterialsThermal.cpp

Timestamp:

07/14/09 14:56:05 (16 years ago)

Author:

Eric.Larour

Message:

Updated model processors to take elementonwater into account

File:

• issm/trunk/src/c/ModelProcessorx/Thermal/CreateElementsNodesAndMaterialsThermal.cpp (modified) (5 diffs)

issm/trunk/src/c/ModelProcessorx/Thermal/CreateElementsNodesAndMaterialsThermal.cpp

@@ -82 +82 @@
         double penta_viscosity_overshoot;
         double penta_stokesreconditioning;
+        bool   penta_onwater;
 
         /*matpar constructor input: */
@@ -197 +198 @@
         ModelFetchData((void**)&model->B,NULL,NULL,model_handle,"B","Matrix","Mat");
         ModelFetchData((void**)&model->n,NULL,NULL,model_handle,"n","Matrix","Mat");
+        ModelFetchData((void**)&model->elementonwater,NULL,NULL,model_handle,"elementonwater","Matrix","Mat");
         
         for (i=0;i<model->numberofelements;i++){
@@ -232 +234 @@
                 penta_meanvel=model->meanvel;
                 penta_epsvel=model->epsvel;
+                penta_onwater=(bool)*(model->elementonwater+i);
 
                 /*We need the field collapse for transient, so that we can use compute B with the average temperature*/
@@ -246 +249 @@
                                 penta_p,penta_q,penta_shelf,penta_onbed,penta_onsurface,penta_meanvel,penta_epsvel,
                                 penta_collapse,penta_melting,penta_accumulation,penta_geothermalflux,penta_artdiff,
-                                penta_thermal_steadystate,penta_viscosity_overshoot,penta_stokesreconditioning,0); 
+                                penta_thermal_steadystate,penta_viscosity_overshoot,penta_stokesreconditioning,penta_onwater); 
 
                 /*Add penta element to elements dataset: */
@@ -303 +306 @@
         xfree((void**)&model->n);
         xfree((void**)&model->B);
+        xfree((void**)&model->elementonwater);
 
         /*Add one constant material property to materials: */