Changeset 4670
- Timestamp:
- 07/20/10 09:27:09 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/modules/ModelProcessorx/Qmu/CreateParametersQmu.cpp
r4491 r4670 203 203 if (strcmp(descriptor,"misfit")==0){ 204 204 205 /*We need the observed velocity: */ 206 IoModelFetchData(&iomodel->vx_obs,NULL,NULL,iomodel_handle,"vx_obs"); 207 IoModelFetchData(&iomodel->vy_obs,NULL,NULL,iomodel_handle,"vy_obs"); 205 /*We need the fit: */ 208 206 IoModelFetchData(&iomodel->fit,NULL,NULL,iomodel_handle,"fit"); 209 210 if(iomodel->vx_obs)for(i=0;i<iomodel->numberofvertices;i++)iomodel->vx_obs[2*i+0]=iomodel->vx_obs[i]/iomodel->yts; 211 if(iomodel->vy_obs)for(i=0;i<iomodel->numberofvertices;i++)iomodel->vy_obs[2*i+1]=iomodel->vy_obs[i]/iomodel->yts; 212 213 parameters->AddObject(new DoubleVecParam(VxObsEnum,iomodel->vx_obs,iomodel->numberofvertices)); 214 parameters->AddObject(new DoubleVecParam(VyObsEnum,iomodel->vy_obs,iomodel->numberofvertices)); 215 parameters->AddObject(new DoubleVecParam(FitEnum,iomodel->fit,iomodel->nsteps)); 216 217 xfree((void**)&iomodel->vx_obs); 218 xfree((void**)&iomodel->vy_obs); 207 parameters->SetParam(iomodel->fit,iomodel->nsteps,FitEnum); 219 208 xfree((void**)&iomodel->fit); 220 209
Note:
See TracChangeset
for help on using the changeset viewer.