Changeset 870

Timestamp:

06/09/09 15:32:17 (16 years ago)

Author:

seroussi

Message:

no pressure in 2d

File:

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

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

@@ -18 +18 @@
         DataSet* parameters=NULL;
         int      count;
-        int i;
+        int      i;
+        int      dim;
 
         double* vx=NULL;
@@ -57 +58 @@
         xfree((void**)&vz);
         xfree((void**)&u_g);
-        
-        /*Get pressure: */
-        ModelFetchData((void**)&pressure,NULL,NULL,model_handle,"pressure","Matrix","Mat");
-        
-        count++;
-        param= new Param(count,"p_g",DOUBLEVEC);
-        if(pressure) param->SetDoubleVec(pressure,model->numberofnodes,1);
-        else param->SetDoubleVec(pressure,0,0);
-        parameters->AddObject(param);
 
-        /*Free pressure: */
-        xfree((void**)&pressure);
+        /*Get pressure if 3d model: */
+        parameters->FindParam((void*)&dim,"dim");
+        if (dim==3){ 
+                ModelFetchData((void**)&pressure,NULL,NULL,model_handle,"pressure","Matrix","Mat");
+                
+                count++;
+                param= new Param(count,"p_g",DOUBLEVEC);
+                if(pressure) param->SetDoubleVec(pressure,model->numberofnodes,1);
+                else param->SetDoubleVec(pressure,0,0);
+                parameters->AddObject(param);
+
+                /*Free pressure: */
+                xfree((void**)&pressure);
+        }
 
         /*Get thickness: */