Changeset 879

Timestamp:

06/10/09 10:22:26 (16 years ago)

Author:

seroussi

Message:

missing temperature for 3d transient

File:

• issm/trunk/src/c/ModelProcessorx/Prognostic/CreateParametersPrognostic.cpp (modified) (2 diffs)

issm/trunk/src/c/ModelProcessorx/Prognostic/CreateParametersPrognostic.cpp

@@ -26 +26 @@
         double* u_g=NULL;
         double* pressure=NULL;
+        double* temperature=NULL;
         double* thickness=NULL;
         double* surface=NULL;
@@ -72 +73 @@
                 /*Free pressure: */
                 xfree((void**)&pressure);
+        }
+
+        /*Get temperature if 3d model: */
+        parameters->FindParam((void*)&dim,"dim");
+        if (dim==3){ 
+                ModelFetchData((void**)&temperature,NULL,NULL,model_handle,"temperature","Matrix","Mat");
+                
+                count++;
+                param= new Param(count,"t_g",DOUBLEVEC);
+                if(temperature) param->SetDoubleVec(temperature,model->numberofnodes,1);
+                else param->SetDoubleVec(temperature,0,0);
+                parameters->AddObject(param);
+
+                /*Free temperature: */
+                xfree((void**)&temperature);
         }