Ignore:
Timestamp:
05/08/10 14:37:25 (15 years ago)
Author:
Eric.Larour
Message:

Redid Params and Inputs

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/c/ModelProcessorx/Prognostic2/CreateParametersPrognostic2.cpp

    r3633 r3673  
    1111#include "../ModelProcessorx.h"
    1212
    13 void CreateParametersPrognostic2(DataSet** pparameters,IoModel* iomodel,ConstDataHandle iomodel_handle){
     13void CreateParametersPrognostic2(Parameters** pparameters,IoModel* iomodel,ConstDataHandle iomodel_handle){
    1414       
    1515        Param*   param = NULL;
    16         DataSet* parameters=NULL;
     16        Parameters* parameters=NULL;
    1717        int      count;
    1818        int      i;
     
    8686        xfree((void**)&iomodel->bed);
    8787
    88         /*Get melting: */
    89         IoModelFetchData(&iomodel->melting,NULL,NULL,iomodel_handle,"melting");
    90         if(iomodel->melting) for(i=0;i<iomodel->numberofvertices;i++)iomodel->melting[i]=iomodel->melting[i]/iomodel->yts;
     88        /*Get melting_rate: */
     89        IoModelFetchData(&iomodel->melting_rate,NULL,NULL,iomodel_handle,"melting_rate");
     90        if(iomodel->melting_rate) for(i=0;i<iomodel->numberofvertices;i++)iomodel->melting_rate[i]=iomodel->melting_rate[i]/iomodel->yts;
    9191       
    9292        count++;
    9393        param= new Param(count,"m_g",DOUBLEVEC);
    94         if(iomodel->melting) param->SetDoubleVec(iomodel->melting,iomodel->numberofvertices,1);
    95         else param->SetDoubleVec(iomodel->melting,0,1);
     94        if(iomodel->melting_rate) param->SetDoubleVec(iomodel->melting_rate,iomodel->numberofvertices,1);
     95        else param->SetDoubleVec(iomodel->melting_rate,0,1);
    9696        parameters->AddObject(param);
    9797
    98         /*Free melting: */
    99         xfree((void**)&iomodel->melting);
     98        /*Free melting_rate: */
     99        xfree((void**)&iomodel->melting_rate);
    100100
    101         /*Get accumulation: */
    102         IoModelFetchData(&iomodel->accumulation,NULL,NULL,iomodel_handle,"accumulation");
    103         if(iomodel->accumulation) for(i=0;i<iomodel->numberofvertices;i++)iomodel->accumulation[i]=iomodel->accumulation[i]/iomodel->yts;
     101        /*Get accumulation_rate: */
     102        IoModelFetchData(&iomodel->accumulation_rate,NULL,NULL,iomodel_handle,"accumulation_rate");
     103        if(iomodel->accumulation_rate) for(i=0;i<iomodel->numberofvertices;i++)iomodel->accumulation_rate[i]=iomodel->accumulation_rate[i]/iomodel->yts;
    104104       
    105105        count++;
    106106        param= new Param(count,"a_g",DOUBLEVEC);
    107         if(iomodel->accumulation) param->SetDoubleVec(iomodel->accumulation,iomodel->numberofvertices,1);
    108         else param->SetDoubleVec(iomodel->accumulation,0,0);
     107        if(iomodel->accumulation_rate) param->SetDoubleVec(iomodel->accumulation_rate,iomodel->numberofvertices,1);
     108        else param->SetDoubleVec(iomodel->accumulation_rate,0,0);
    109109        parameters->AddObject(param);
    110110
    111         /*Free accumulation: */
    112         xfree((void**)&iomodel->accumulation);
     111        /*Free accumulation_rate: */
     112        xfree((void**)&iomodel->accumulation_rate);
    113113
    114114        /*Assign output pointer: */
Note: See TracChangeset for help on using the changeset viewer.