Changeset 3461


Ignore:
Timestamp:
04/08/10 08:54:19 (15 years ago)
Author:
Mathieu Morlighem
Message:

put all loaded variables in iomodel

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/c/ModelProcessorx/Qmu/CreateParametersQmu.cpp

    r3446 r3461  
    4545        /*parameters for misfit: */
    4646        double* u_g_obs=NULL;
    47         double* vx_obs=NULL;
    48         double* vy_obs=NULL;
    49         double* fit=NULL;
    5047
    5148        #ifdef _SERIAL_
     
    153150                #endif
    154151
    155 
    156152                /*Ok, we have all the response descriptors. Build a parameter with it: */
    157153                count++;
     
    212208                        }
    213209                }
    214 
    215210
    216211                /*Deal with data needed for some responses: */
     
    268263
    269264                                /*We need the observed velocity: */
    270                                 IoModelFetchData(&vx_obs,NULL,NULL,iomodel_handle,"vx_obs");
    271                                 IoModelFetchData(&vy_obs,NULL,NULL,iomodel_handle,"vy_obs");
     265                                IoModelFetchData(&iomodel->vx_obs,NULL,NULL,iomodel_handle,"vx_obs");
     266                                IoModelFetchData(&iomodel->vy_obs,NULL,NULL,iomodel_handle,"vy_obs");
    272267
    273268                                /*Now, recover fit: */
     
    275270
    276271                                u_g_obs=(double*)xcalloc(iomodel->numberofvertices*2,sizeof(double));
    277                                 if(vx_obs)for(i=0;i<iomodel->numberofvertices;i++)u_g_obs[2*i+0]=vx_obs[i]/iomodel->yts;
    278                                 if(vy_obs)for(i=0;i<iomodel->numberofvertices;i++)u_g_obs[2*i+1]=vy_obs[i]/iomodel->yts;
     272                                if(iomodel->vx_obs)for(i=0;i<iomodel->numberofvertices;i++)u_g_obs[2*i+0]=iomodel->vx_obs[i]/iomodel->yts;
     273                                if(iomodel->vy_obs)for(i=0;i<iomodel->numberofvertices;i++)u_g_obs[2*i+1]=iomodel->vy_obs[i]/iomodel->yts;
    279274
    280275                                count++;
     
    288283                                parameters->AddObject(param);
    289284
     285                                xfree((void**)&iomodel->vx_obs);
     286                                xfree((void**)&iomodel->vy_obs);
     287                                xfree((void**)&iomodel->fit);
    290288                                xfree((void**)&u_g_obs);
    291                                 xfree((void**)&vx_obs);
    292                                 xfree((void**)&vy_obs);
    293                                 xfree((void**)&fit);
     289
    294290                        }
    295291                }
Note: See TracChangeset for help on using the changeset viewer.