source: issm/trunk/src/m/solutions/jpl/diagnostic.m@ 4002

function md=diagnostic(md);
%DIAGNOSTIC - compute the velocity field of a model
%
%   Usage:
%      md=diagnostic(md)
%
        %timing
        t1=clock;

        analyses=[DiagnosticHorizAnalysisEnum,DiagnosticVertAnalysisEnum,DiagnosticStokesAnalysisEnum,DiagnosticHutterAnalysisEnum,SlopecomputeAnalysisEnum];

        displaystring(md.verbose,'%s',['create fem model']);
        femmodel=CreateFemModel(md,analyses);

        %compute solution
        if ~femmodel.parameters.QmuAnalysis,
                if femodel.parameters.control_analysis,
                        
                        %launch core of control solution.
                        md.results.DiagnosticAnalysis=control_core(femmodel);

                else,
                        
                        %launch core of diagnostic solution.
                        md.results.DiagnosticAnalysis=diagnostic_core(femmodel);

                end
        else
                %launch dakota driver for diagnostic core solution
                Qmu(femmodel);
        end

        %stop timing
        t2=clock;
        displaystring(md.verbose,'\n%s\n',['solution converged in ' num2str(etime(t2,t1)) ' seconds']);