function B = rooftemp(reactortype)% rooftemp.m A nonlinear algebraic equation model object constructor,%          derived from the naemodel class. Called as%%          A = rooftemp(reactortype)%% INPUT PARAMETERS%%         reactortype % 'aixtron','thomasswann'switch reactortype        case 'aixtron'        disp('aixtron geom')        Tsus = 1200+273; % K        kRoof = 0.06; % W/(m K) - value for carbon foam                      % which is approximated by the kAir        Lrxr = 2.54/100; % susceptor/roof spacing m        Lins = 2.54/100; % insulation thickness m        Esus = 0.91; % SiC emissivity        Eroof = Esus;    case 'thomasswann'        disp('thomas swann geom')        Tsus = 1040+273; % K        kRoof = 25; % W/(m K) - value for stainless steel        Lrxr = 1/100; % wafer/showerhead spacing m        Lins = 5/100; % showerhead effective thickness m        Esus = 0.91; % SiC emissivity        Eroof = 0.2; % stainless steel emissivityendsig = 5.67e-8;      % J/(sec*m^2*K^4)kH2 = 0.185;        % W/(m K)% Solution initial guessT = (Tsus+300)/2;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% set up an instance of the twodww objectvar   = assocarray({'T' T});param = assocarray({ 'Tsus' Tsus 'kRoof' kRoof 'Lrxr' Lrxr ...    'Lins' Lins 'Esus' Esus 'Eroof' Eroof 'sig' sig 'kH2' kH2 });A = naemodel('rooftemp',var,param);B = struct([]);B = class(B,'rooftemp',A);