source: issm/trunk-jpl/src/m/consistency/ismodelselfconsistent.m@ 21530

function ismodelselfconsistent(md),
%ISMODELSELFCONSISTENT - check that model forms a closed form solvable problem.
%
%   Usage:
%      ismodelselfconsistent(md),

%initialize consistency as true
md.private.isconsistent=true;

%Get solution and associated analyses
solution=md.private.solution;
[analyses]=AnalysisConfiguration(solution);

%Go through a model field, check that it is a class, and call checkconsistency
fields=properties('model');
for i=1:length(fields),
        field=fields{i};

        %Some properties do not need to be checked
        if ismember(field,{'results' 'debug' 'radaroverlay'}),
                continue;
        end

        %Check that current field is an object
        if ~isobject(md.(field))
                md=checkmessage(md,['field ''' char(field) ''' is not an object']);
                continue;
        end

        %Check consistency of the object
        md=checkconsistency(md.(field),md,solution,analyses);
end

%error message if mode is not consistent
if md.private.isconsistent==false,
        error('Model not consistent, see messages above');
end
end

function [analyses]=AnalysisConfiguration(solutiontype), % {{{
%ANALYSISCONFIGURATION - return type of analyses, number of analyses 
%
%   Usage:
%      [analyses]=AnalysisConfiguration(solutiontype);

        if strcmp(solutiontype,'StressbalanceSolution')
                analyses={'StressbalanceAnalysis','StressbalanceVerticalAnalysis','StressbalanceSIAAnalysis','L2ProjectionBaseAnalysis'};
        elseif strcmp(solutiontype,'SteadystateSolution')
                analyses={'StressbalanceAnalysis','StressbalanceVerticalAnalysis','StressbalanceSIAAnalysis','L2ProjectionBaseAnalysis','ThermalAnalysis','MeltingAnalysis','EnthalpyAnalysis'};
        elseif strcmp(solutiontype,'ThermalSolution')
                analyses={'EnthalpyAnalysis','ThermalAnalysis','MeltingAnalysis'};
        elseif strcmp(solutiontype,'MasstransportSolution')
                analyses={'MasstransportAnalysis'};
        elseif strcmp(solutiontype,'BalancethicknessSolution')
                analyses={'BalancethicknessAnalysis'};
        elseif strcmp(solutiontype,'Balancethickness2Solution')
                analyses={'Balancethickness2Analysis'};
        elseif strcmp(solutiontype,'BalancethicknessSoftSolution')
                analyses={'BalancethicknessAnalysis'};
        elseif strcmp(solutiontype,'BalancevelocitySolution')
                analyses={'BalancevelocityAnalysis'};
        elseif strcmp(solutiontype,'SurfaceSlopeSolution')
                analyses={'L2ProjectionBaseAnalysis'};
        elseif strcmp(solutiontype,'BedSlopeSolution')
                analyses={'L2ProjectionBaseAnalysis'};
        elseif strcmp(solutiontype,'GiaIvinsSolution')
                analyses={'GiaIvinsAnalysis'};
        elseif strcmp(solutiontype,'EsaSolution')
                analyses={'EsaAnalysis'};
        elseif strcmp(solutiontype,'TransientSolution')
                analyses={'StressbalanceAnalysis','StressbalanceVerticalAnalysis','StressbalanceSIAAnalysis','L2ProjectionBaseAnalysis','ThermalAnalysis','MeltingAnalysis','EnthalpyAnalysis','MasstransportAnalysis','HydrologySommersAnalysis'};
        elseif strcmp(solutiontype,'SealevelriseSolution')
                analyses={'SealevelriseAnalysis'};
        elseif strcmp(solutiontype,'HydrologySolution')
                analyses={'L2ProjectionBaseAnalysis','HydrologyShreveAnalysis','HydrologyDCInefficientAnalysis','HydrologyDCEfficientAnalysis'};
        elseif strcmp(solutiontype,'DamageEvolutionSolution')
                analyses={'DamageEvolutionAnalysis'};
        else
                error(' solution type: %s' , solutiontype, ' not supported yet!');
        end
end % }}}