- Timestamp:
- 04/07/10 16:45:45 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/ModelProcessorx/Thermal/CreateParametersThermal.cpp
r3332 r3436 44 44 IoModelFetchData(&vz,NULL,NULL,iomodel_handle,"vz"); 45 45 46 u_g=(double*)xcalloc(iomodel->numberof nodes*3,sizeof(double));46 u_g=(double*)xcalloc(iomodel->numberofvertices*3,sizeof(double)); 47 47 48 if(vx) for(i=0;i<iomodel->numberof nodes;i++)u_g[3*i+0]=vx[i]/iomodel->yts;49 if(vy) for(i=0;i<iomodel->numberof nodes;i++)u_g[3*i+1]=vy[i]/iomodel->yts;50 if(vz) for(i=0;i<iomodel->numberof nodes;i++)u_g[3*i+2]=vz[i]/iomodel->yts;48 if(vx) for(i=0;i<iomodel->numberofvertices;i++)u_g[3*i+0]=vx[i]/iomodel->yts; 49 if(vy) for(i=0;i<iomodel->numberofvertices;i++)u_g[3*i+1]=vy[i]/iomodel->yts; 50 if(vz) for(i=0;i<iomodel->numberofvertices;i++)u_g[3*i+2]=vz[i]/iomodel->yts; 51 51 52 52 count++; 53 53 param= new Param(count,"u_g",DOUBLEVEC); 54 param->SetDoubleVec(u_g,3*iomodel->numberof nodes,3);54 param->SetDoubleVec(u_g,3*iomodel->numberofvertices,3); 55 55 parameters->AddObject(param); 56 56 … … 65 65 count++; 66 66 param= new Param(count,"p_g",DOUBLEVEC); 67 if(pressure) param->SetDoubleVec(pressure,iomodel->numberof nodes,1);67 if(pressure) param->SetDoubleVec(pressure,iomodel->numberofvertices,1); 68 68 else param->SetDoubleVec(pressure,0,0); 69 69 parameters->AddObject(param); … … 80 80 count++; 81 81 param= new Param(count,"t_g",DOUBLEVEC); 82 if(temperature) param->SetDoubleVec(temperature,iomodel->numberof nodes,1);82 if(temperature) param->SetDoubleVec(temperature,iomodel->numberofvertices,1); 83 83 else ISSMERROR(exprintf("Missing initial temperature")); 84 84 parameters->AddObject(param);
Note:
See TracChangeset
for help on using the changeset viewer.