Changeset 3043

Timestamp:

02/16/10 14:22:28 (15 years ago)

Author:

Eric.Larour

Message:

Added parameters necessary for misfit computation in qmu_analysis=1.

File:

• issm/trunk/src/c/ModelProcessorx/Qmu/CreateParametersQmu.cpp (modified) (2 diffs)

issm/trunk/src/c/ModelProcessorx/Qmu/CreateParametersQmu.cpp

@@ -47 +47 @@
         int num_qmu_mass_flux_segments=0;
         int my_num_qmu_mass_flux_segments=0;
+        
+        /*parameters for misfit: */
+        double* u_g_obs=NULL;
+        double* vx_obs=NULL;
+        double* vy_obs=NULL;
+        double* fit=NULL;
 
         #ifdef _SERIAL_
@@ -272 +278 @@
                                 xfree((void**)&my_qmu_mass_flux_segments);
                         }
+                        if (strcmp(descriptor,"misfit")==0){
+
+                                /*We need the observed velocity: */
+                                IoModelFetchData(&vx_obs,NULL,NULL,iomodel_handle,"vx_obs");
+                                IoModelFetchData(&vy_obs,NULL,NULL,iomodel_handle,"vy_obs");
+
+                                /*Now, recover fit: */
+                                IoModelFetchData(&iomodel->fit,NULL,NULL,iomodel_handle,"fit");
+
+                                u_g_obs=(double*)xcalloc(iomodel->numberofnodes*2,sizeof(double));
+                                if(vx_obs)for(i=0;i<iomodel->numberofnodes;i++)u_g_obs[2*i+0]=vx_obs[i]/iomodel->yts;
+                                if(vy_obs)for(i=0;i<iomodel->numberofnodes;i++)u_g_obs[2*i+1]=vy_obs[i]/iomodel->yts;
+
+                                count++;
+                                param= new Param(count,"velocity_obs",DOUBLEVEC);
+                                param->SetDoubleVec(u_g_obs,2*iomodel->numberofnodes,2);
+                                parameters->AddObject(param);
+
+                                count++;
+                                param= new Param(count,"fit",DOUBLEVEC);
+                                param->SetDoubleVec(iomodel->fit,iomodel->nsteps);
+                                parameters->AddObject(param);
+
+                                xfree((void**)&u_g_obs);
+                                xfree((void**)&vx_obs);
+                                xfree((void**)&vy_obs);
+                                xfree((void**)&fit);
+                        }
                 }