source: issm/trunk-jpl/src/m/qmu/qmuresponse.m@ 13646

function response=qmuresponse(models,results,processedresults,descriptor)
%QMURESPONSE - compute response function from model results.

if strcmpi(descriptor,'max_vel'),
        response=max(processedresults.vel);
elseif strcmpi(descriptor,'min_vel'),
        response=min(processedresults.vel);
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 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_dhu.parameters.ishutter;
        ismacayealpattyn=m_dh.parameters.ismacayealpattyn;
        isstokes=m_ds.parameters.isstokes;
        if ishutter,

% for now, separate all segments from double array for parallel to make cells
                if (length(m_dhu.parameters.qmu_mass_flux_num_segments) > 1)
                        segments=m_dhu.parameters.qmu_mass_flux_segments;
                        m_dhu.parameters.qmu_mass_flux_segments=cell(size(m_dhu.parameters.qmu_mass_flux_num_segments));
                        ipt=1;
                        for i=1:length(m_dhu.parameters.qmu_mass_flux_num_segments)
                                if m_dhu.parameters.qmu_mass_flux_num_segments(i)
                                        m_dhu.parameters.qmu_mass_flux_segments{i}=segments(ipt:ipt+m_dhu.parameters.qmu_mass_flux_num_segments(i)-1,:);
                                        ipt=ipt+m_dhu.parameters.qmu_mass_flux_num_segments(i);
                                end
                        end
                        clear segments
                end

                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,

% 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,

% for now, separate all segments from double array for parallel to make cells
                if (length(m_ds.parameters.qmu_mass_flux_num_segments) > 1)
                        segments=m_ds.parameters.qmu_mass_flux_segments;
                        m_ds.parameters.qmu_mass_flux_segments=cell(size(m_ds.parameters.qmu_mass_flux_num_segments));
                        ipt=1;
                        for i=1:length(m_ds.parameters.qmu_mass_flux_num_segments)
                                if m_ds.parameters.qmu_mass_flux_num_segments(i)
                                        m_ds.parameters.qmu_mass_flux_segments{i}=segments(ipt:ipt+m_ds.parameters.qmu_mass_flux_num_segments(i)-1,:);
                                        ipt=ipt+m_ds.parameters.qmu_mass_flux_num_segments(i);
                                end
                        end
                        clear segments
                end

                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: unknown descriptor ' descriptor]);
end