- Timestamp:
- 05/08/10 14:37:25 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/ModelProcessorx/Prognostic2/CreateParametersPrognostic2.cpp
r3633 r3673 11 11 #include "../ModelProcessorx.h" 12 12 13 void CreateParametersPrognostic2( DataSet** pparameters,IoModel* iomodel,ConstDataHandle iomodel_handle){13 void CreateParametersPrognostic2(Parameters** pparameters,IoModel* iomodel,ConstDataHandle iomodel_handle){ 14 14 15 15 Param* param = NULL; 16 DataSet* parameters=NULL;16 Parameters* parameters=NULL; 17 17 int count; 18 18 int i; … … 86 86 xfree((void**)&iomodel->bed); 87 87 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; 91 91 92 92 count++; 93 93 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); 96 96 parameters->AddObject(param); 97 97 98 /*Free melting : */99 xfree((void**)&iomodel->melting );98 /*Free melting_rate: */ 99 xfree((void**)&iomodel->melting_rate); 100 100 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; 104 104 105 105 count++; 106 106 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); 109 109 parameters->AddObject(param); 110 110 111 /*Free accumulation : */112 xfree((void**)&iomodel->accumulation );111 /*Free accumulation_rate: */ 112 xfree((void**)&iomodel->accumulation_rate); 113 113 114 114 /*Assign output pointer: */
Note:
See TracChangeset
for help on using the changeset viewer.