Ignore:
Timestamp:
04/07/10 16:45:45 (15 years ago)
Author:
seroussi
Message:

change CreateParameters with numberofvertices

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/c/ModelProcessorx/Thermal/CreateParametersThermal.cpp

    r3332 r3436  
    4444        IoModelFetchData(&vz,NULL,NULL,iomodel_handle,"vz");
    4545
    46         u_g=(double*)xcalloc(iomodel->numberofnodes*3,sizeof(double));
     46        u_g=(double*)xcalloc(iomodel->numberofvertices*3,sizeof(double));
    4747
    48         if(vx) for(i=0;i<iomodel->numberofnodes;i++)u_g[3*i+0]=vx[i]/iomodel->yts;
    49         if(vy) for(i=0;i<iomodel->numberofnodes;i++)u_g[3*i+1]=vy[i]/iomodel->yts;
    50         if(vz) for(i=0;i<iomodel->numberofnodes;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;
    5151
    5252        count++;
    5353        param= new Param(count,"u_g",DOUBLEVEC);
    54         param->SetDoubleVec(u_g,3*iomodel->numberofnodes,3);
     54        param->SetDoubleVec(u_g,3*iomodel->numberofvertices,3);
    5555        parameters->AddObject(param);
    5656
     
    6565        count++;
    6666        param= new Param(count,"p_g",DOUBLEVEC);
    67         if(pressure) param->SetDoubleVec(pressure,iomodel->numberofnodes,1);
     67        if(pressure) param->SetDoubleVec(pressure,iomodel->numberofvertices,1);
    6868        else param->SetDoubleVec(pressure,0,0);
    6969        parameters->AddObject(param);
     
    8080                count++;
    8181                param= new Param(count,"t_g",DOUBLEVEC);
    82                 if(temperature) param->SetDoubleVec(temperature,iomodel->numberofnodes,1);
     82                if(temperature) param->SetDoubleVec(temperature,iomodel->numberofvertices,1);
    8383                else ISSMERROR(exprintf("Missing initial temperature"));
    8484                parameters->AddObject(param);
Note: See TracChangeset for help on using the changeset viewer.