dwdataclass.m

% Class to combine double wedge geometry data
classdef dwdataclass
    properties
        theta1 double {mustBeReal, mustBeNonnegative, mustBeFinite}
        theta2 double {mustBeReal, mustBeNonnegative, mustBeFinite}
        L1 double {mustBeReal, mustBeNonnegative, mustBeFinite}
        L2 double {mustBeReal, mustBeNonnegative, mustBeFinite}
        xu double {mustBeReal, mustBeNonnegative, mustBeFinite}
    end
    methods
        function obj = dwdataclass(varargin)
            if nargin == 0 
                % Choose default values
                obj.theta1 = 9*pi/180;          %(rad) First wedge angle
                obj.theta2 = 45*pi/180;         %(rad) Second wedge angle
                obj.L1 = 0.2;                   %(m) First wedge length
                obj.L2 = obj.L1;                %(m) Second wedge length
                obj.xu = obj.L1/4;              %(m) Upstream x-length
            elseif nargin == 1
                % Choose default values for wall and given values for
                % domain boundary
                obj.theta1 = 9*pi/180;          %(rad) First wedge angle
                obj.theta2 = 45*pi/180;         %(rad) Second wedge angle
                obj.L1 = 0.2;                   %(m) First wedge length
                obj.L2 = obj.L1;                %(m) Second wedge length
                obj.xu = varargin;              %(m) Upstream length
            elseif nargin == 5
                obj.theta1 = varargin{1};       %(rad) First wedge angle
                obj.theta2 = varargin{2};       %(rad) Second wedge angle
                obj.L1 = varargin{3};           %(m) First wedge length
                obj.L2 = varargin{4};           %(m) Second wedge length
                obj.xu = varargin{5};           %(m) Upstream length
            end
        end
    end
end