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

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,numanalyses]=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']);
        end

        %Check consistency of the object
        if verLessThan('matlab', '7.6')
                md=checkconsistency(md.(field),md,solution,analyses);
        else
                md=md.(field).checkconsistency(md,solution,analyses);
        end
end

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

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



switch solutiontype,

        case DiagnosticSolutionEnum(),
                numanalyses=5;
                analyses=[DiagnosticHorizAnalysisEnum();DiagnosticVertAnalysisEnum();DiagnosticHutterAnalysisEnum();SurfaceSlopeAnalysisEnum();BedSlopeAnalysisEnum()];

        case SteadystateSolutionEnum(),
                numanalyses=7; 
                analyses=[DiagnosticHorizAnalysisEnum();DiagnosticVertAnalysisEnum();DiagnosticHutterAnalysisEnum();SurfaceSlopeAnalysisEnum();BedSlopeAnalysisEnum();ThermalAnalysisEnum();MeltingAnalysisEnum()];

        case ThermalSolutionEnum(),
                numanalyses=2; 
                analyses=[ThermalAnalysisEnum();MeltingAnalysisEnum()];

        case EnthalpySolutionEnum(),
                numanalyses=1; 
                analyses=[EnthalpyAnalysisEnum()];

        case PrognosticSolutionEnum(),
                numanalyses=1; 
                analyses=[PrognosticAnalysisEnum()];

        case BalancethicknessSolutionEnum(),
                numanalyses=1; 
                analyses=[BalancethicknessAnalysisEnum()];

        case SurfaceSlopeSolutionEnum(),
                numanalyses=1; 
                analyses=[SurfaceSlopeAnalysisEnum()];

        case BedSlopeSolutionEnum(),
                numanalyses=1; 
                analyses=[BedSlopeAnalysisEnum()];

        case TransientSolutionEnum(),
                numanalyses=9; 
                analyses=[DiagnosticHorizAnalysisEnum();DiagnosticVertAnalysisEnum();DiagnosticHutterAnalysisEnum();SurfaceSlopeAnalysisEnum();BedSlopeAnalysisEnum();ThermalAnalysisEnum();MeltingAnalysisEnum();EnthalpyAnalysisEnum();PrognosticAnalysisEnum()];

        case FlaimSolutionEnum(),
                numanalyses=1; 
                analyses=[FlaimAnalysisEnum()];

        case HydrologySolutionEnum(),
                numanalyses=3; 
                analyses=[BedSlopeAnalysisEnum();SurfaceSlopeAnalysisEnum();HydrologyAnalysisEnum()];

        otherwise
                error('%s%s%s',' solution type: ',EnumToString(solutiontype),' not supported yet!');

        end
end % }}}