Changeset 3436 for issm/trunk/src/c/ModelProcessorx/Thermal/CreateParametersThermal.cpp

Timestamp:

04/07/10 16:45:45 (15 years ago)

Author:

seroussi

Message:

change CreateParameters with numberofvertices

File:

• issm/trunk/src/c/ModelProcessorx/Thermal/CreateParametersThermal.cpp (modified) (3 diffs)

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

@@ -44 +44 @@
         IoModelFetchData(&vz,NULL,NULL,iomodel_handle,"vz");
 
-        u_g=(double*)xcalloc(iomodel->numberofnodes*3,sizeof(double));
+        u_g=(double*)xcalloc(iomodel->numberofvertices*3,sizeof(double));
 
-        if(vx) for(i=0;i<iomodel->numberofnodes;i++)u_g[3*i+0]=vx[i]/iomodel->yts;
-        if(vy) for(i=0;i<iomodel->numberofnodes;i++)u_g[3*i+1]=vy[i]/iomodel->yts;
-        if(vz) for(i=0;i<iomodel->numberofnodes;i++)u_g[3*i+2]=vz[i]/iomodel->yts;
+        if(vx) for(i=0;i<iomodel->numberofvertices;i++)u_g[3*i+0]=vx[i]/iomodel->yts;
+        if(vy) for(i=0;i<iomodel->numberofvertices;i++)u_g[3*i+1]=vy[i]/iomodel->yts;
+        if(vz) for(i=0;i<iomodel->numberofvertices;i++)u_g[3*i+2]=vz[i]/iomodel->yts;
 
         count++;
         param= new Param(count,"u_g",DOUBLEVEC);
-        param->SetDoubleVec(u_g,3*iomodel->numberofnodes,3);
+        param->SetDoubleVec(u_g,3*iomodel->numberofvertices,3);
         parameters->AddObject(param);
 
@@ -65 +65 @@
         count++;
         param= new Param(count,"p_g",DOUBLEVEC);
-        if(pressure) param->SetDoubleVec(pressure,iomodel->numberofnodes,1);
+        if(pressure) param->SetDoubleVec(pressure,iomodel->numberofvertices,1);
         else param->SetDoubleVec(pressure,0,0);
         parameters->AddObject(param);
@@ -80 +80 @@
                 count++;
                 param= new Param(count,"t_g",DOUBLEVEC);
-                if(temperature) param->SetDoubleVec(temperature,iomodel->numberofnodes,1);
+                if(temperature) param->SetDoubleVec(temperature,iomodel->numberofvertices,1);
                 else ISSMERROR(exprintf("Missing initial temperature"));
                 parameters->AddObject(param);