Changeset 4430 for issm/trunk/src/m/dakota/qmuresponse.m

Timestamp:

07/07/10 11:31:55 (15 years ago)

Author:

jschierm

Message:

Qmu: Calculate mass fluxes on multiple profiles.

File:

• issm/trunk/src/m/dakota/qmuresponse.m (modified) (1 diff)

issm/trunk/src/m/dakota/qmuresponse.m

@@ -8 +8 @@
 elseif strcmpi(descriptor,'max_vx'),
         response=max(processedresults.vx);
+elseif strcmpi(descriptor,'max_abs_vx'),
+        response=max(abs(processedresults.vx));
 elseif strcmpi(descriptor,'min_vx'),
         response=min(processedresults.vx);
 elseif strcmpi(descriptor,'max_vy'),
         response=max(processedresults.vy);
+elseif strcmpi(descriptor,'max_abs_vy'),
+        response=max(abs(processedresults.vy));
 elseif strcmpi(descriptor,'min_vy'),
         response=min(processedresults.vy);
-elseif strcmpi(descriptor,'mass_flux'),
+elseif strncmpi(descriptor,'mass_flux',9),
+        indx=str2int(descriptor(10:end));
+        if isempty(indx) || ~indx
+                indx=1;
+        end
+
         %call mass flux module.
         m_dh=models.dh;
         m_dhu=models.dhu;
         m_ds=models.ds;
-        ishutter=m_dh.parameters.ishutter;
-        ismacayealpattyn=m_dhu.parameters.ismacayealpattyn;
+        ishutter=m_dhu.parameters.ishutter;
+        ismacayealpattyn=m_dh.parameters.ismacayealpattyn;
         isstokes=m_ds.parameters.isstokes;
         if ishutter,
-                response=MassFlux(m_dhu.elements,m_dhu.nodes,m_dhu.vertices,m_dhu.loads,m_dhu.materials,m_dhu.parameters,results.u_g);
+                if isnumeric(m_dhu.parameters.qmu_mass_flux_segments)
+                        response=MassFlux(m_dhu.elements,m_dhu.nodes,m_dhu.vertices,m_dhu.loads,m_dhu.materials,m_dhu.parameters,results.u_g);
+                else
+                        save=m_dhu.parameters.qmu_mass_flux_segments;
+                        m_dhu.parameters.qmu_mass_flux_segments=m_dhu.parameters.qmu_mass_flux_segments{indx};
+                        response=MassFlux(m_dhu.elements,m_dhu.nodes,m_dhu.vertices,m_dhu.loads,m_dhu.materials,m_dhu.parameters,results.u_g);
+                        m_dhu.parameters.qmu_mass_flux_segments=save;
+                        clear save
+                end
         elseif ismacayealpattyn,
-                response=MassFlux(m_dh.elements,m_dh.nodes,m_dh.vertices,m_dh.loads,m_dh.materials,m_dh.parameters,results.u_g);
+
+% for now, separate all segments from double array for parallel to make cells
+        if (length(m_dh.parameters.qmu_mass_flux_num_segments) > 1)
+                segments=m_dh.parameters.qmu_mass_flux_segments;
+                m_dh.parameters.qmu_mass_flux_segments=cell(size(m_dh.parameters.qmu_mass_flux_num_segments));
+                ipt=1;
+                for i=1:length(m_dh.parameters.qmu_mass_flux_num_segments)
+                        if m_dh.parameters.qmu_mass_flux_num_segments(i)
+                                m_dh.parameters.qmu_mass_flux_segments{i}=segments(ipt:ipt+m_dh.parameters.qmu_mass_flux_num_segments(i)-1,:);
+                                ipt=ipt+m_dh.parameters.qmu_mass_flux_num_segments(i);
+                        end
+                end
+                clear segments
+        end
+
+                if isnumeric(m_dh.parameters.qmu_mass_flux_segments)
+                        response=MassFlux(m_dh.elements,m_dh.nodes,m_dh.vertices,m_dh.loads,m_dh.materials,m_dh.parameters,results.u_g);
+                else
+                        save=m_dh.parameters.qmu_mass_flux_segments;
+                        m_dh.parameters.qmu_mass_flux_segments=m_dh.parameters.qmu_mass_flux_segments{indx};
+                        response=MassFlux(m_dh.elements,m_dh.nodes,m_dh.vertices,m_dh.loads,m_dh.materials,m_dh.parameters,results.u_g);
+                        m_dh.parameters.qmu_mass_flux_segments=save;
+                        clear save
+                end
         elseif isstokes,
-                response=MassFlux(m_ds.elements,m_ds.nodes,m_ds.vertices,m_ds.loads,m_ds.materials,m_ds.parameters,results.u_g);
+                if isnumeric(m_ds.parameters.qmu_mass_flux_segments)
+                        response=MassFlux(m_ds.elements,m_ds.nodes,m_ds.vertices,m_ds.loads,m_ds.materials,m_ds.parameters,results.u_g);
+                else
+                        save=m_ds.parameters.qmu_mass_flux_segments;
+                        m_ds.parameters.qmu_mass_flux_segments=m_ds.parameters.qmu_mass_flux_segments{indx};
+                        response=MassFlux(m_ds.elements,m_ds.nodes,m_ds.vertices,m_ds.loads,m_ds.materials,m_ds.parameters,results.u_g);
+                        m_ds.parameters.qmu_mass_flux_segments=save;
+                        clear save
+                end
         else
                 error('qmuresponse error message: unsupported analysis type for mass_flux computation!');
         end
 else
-        error(['qmuresponse error message: unknow descriptor ' descriptor]);
+        error(['qmuresponse error message: unknown descriptor ' descriptor]);
 end
+