Changeset 3462


Ignore:
Timestamp:
04/08/10 09:02:06 (15 years ago)
Author:
Mathieu Morlighem
Message:

minor bug fix

Location:
issm/trunk/src/c/ModelProcessorx
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/c/ModelProcessorx/Balancedthickness/CreateParametersBalancedthickness.cpp

    r3460 r3462  
    1  * \brief driver for creating parameters dataset, for prognostic analysis.
     1 /* \brief driver for creating parameters dataset, for prognostic analysis.
    22 */
    33
     
    5050        count++;
    5151        param= new Param(count,"m_g",DOUBLEVEC);
    52         if(iomodel->melting) param->SetDoubleVec(melting,iomodel->numberofvertices,1);
    53         else param->SetDoubleVec(melting,0,1);
     52        if(iomodel->melting) param->SetDoubleVec(iomodel->melting,iomodel->numberofvertices,1);
     53        else param->SetDoubleVec(iomodel->melting,0,1);
    5454        parameters->AddObject(param);
    5555
  • issm/trunk/src/c/ModelProcessorx/Prognostic/CreateParametersPrognostic.cpp

    r3460 r3462  
    5050        parameters->FindParam(&dim,"dim");
    5151        if (dim==3){
    52                 IoModelFetchData(&pressure,NULL,NULL,iomodel_handle,"pressure");
     52                IoModelFetchData(&iomodel->pressure,NULL,NULL,iomodel_handle,"pressure");
    5353               
    5454                count++;
    5555                param= new Param(count,"p_g",DOUBLEVEC);
    56                 if(pressure) param->SetDoubleVec(pressure,iomodel->numberofvertices,1);
    57                 else param->SetDoubleVec(pressure,0,0);
     56                if(iomodel->pressure) param->SetDoubleVec(iomodel->pressure,iomodel->numberofvertices,1);
     57                else param->SetDoubleVec(iomodel->pressure,0,0);
    5858                parameters->AddObject(param);
    5959
    6060                /*Free pressure: */
    61                 xfree((void**)&pressure);
     61                xfree((void**)&iomodel->pressure);
    6262        }
    6363
     
    6565        parameters->FindParam(&dim,"dim");
    6666        if (dim==3){
    67                 IoModelFetchData(&temperature,NULL,NULL,iomodel_handle,"temperature");
     67                IoModelFetchData(&iomodel->temperature,NULL,NULL,iomodel_handle,"temperature");
    6868               
    6969                count++;
    7070                param= new Param(count,"t_g",DOUBLEVEC);
    71                 if(temperature) param->SetDoubleVec(temperature,iomodel->numberofvertices,1);
    72                 else param->SetDoubleVec(temperature,0,0);
     71                if(iomodel->temperature) param->SetDoubleVec(iomodel->temperature,iomodel->numberofvertices,1);
     72                else param->SetDoubleVec(iomodel->temperature,0,0);
    7373                parameters->AddObject(param);
    7474
    7575                /*Free temperature: */
    76                 xfree((void**)&temperature);
     76                xfree((void**)&iomodel->temperature);
    7777        }
    7878
     
    114114
    115115        /*Get melting: */
    116         IoModelFetchData(&iiomodel->melting,NULL,NULL,iomodel_handle,"melting");
     116        IoModelFetchData(&iomodel->melting,NULL,NULL,iomodel_handle,"melting");
    117117        if(iomodel->melting) for(i=0;i<iomodel->numberofvertices;i++)iomodel->melting[i]=iomodel->melting[i]/iomodel->yts;
    118118       
Note: See TracChangeset for help on using the changeset viewer.