function B = wafer2d% wafer2d.m A nonlinear algebraic equation model object constructor,%          derived from the naemodel class. Called as%%          A = wafer2d%%         thickness = 350; % microns%         diameter = 2;    % inRTobj = rooftemp('thomasswann');RTobj = newton(RTobj);Troof = gettemp(RTobj)unpack(RTobj,'param') % to get kH2, Lrxr, TsusLrxr = Lrxr*1000; % convert to mmkH2 = kH2/1000;   % J/(sec*mm*K)ESiC = 0.91;         % unitlesssig = 5.67e-8;       % J/(sec*m^2*K^4)kSiC = 300;          % J/(sec*m*K) (approx)kSapp = 10;          % J/(sec*m*K) (8 at 600C; 9 at 700C; about 11 at 800C)betaSapp = 5.6e-6;   % expansion coeff  1/Kthickness = 350;     % wafer thickness, micronsLw = thickness/10^6; % convert to malph = 0.75;         % 1-transmittivityRad = 2*2.54/100;    % wafer radiuskwaf = kSapp;Esat  = ESiC;Eroof = 0.4;% Number of quadrature pointsnr = 20;nz = 10;R = quadgrid('cyln',nr,'r',[0 Rad]);S = R*quadgrid('slab',nz,'z',[0 Lw]);Dr  = linearoperator(S,'d', 'r');DDr = linearoperator(S,'dd','r');Dz  = linearoperator(S,'d', 'z');DDz = linearoperator(S,'dd','z');% Solution initial guessT = scalarfield(S,rand(nr,nz)+Tsus);y = scalarfield(R,0);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% set up an instance of the twodww objectvar   = assocarray({'T' T 'y' y});param = assocarray({ 'Dr' Dr 'DDr' DDr 'Dz' Dz 'DDz' DDz ...                     'Troof' Troof 'Tsus' Tsus 'sig' sig ...                     'alph' alph 'Ewaf' ESiC ...                     'Esat' Esat 'Eroof' Eroof ...                     'kwaf' kwaf 'kH2' kH2 'Lrxr' Lrxr ...                     'Lw' Lw 'betaSapp' betaSapp ...                 });A = naemodel('wafer2d',var,param);B = struct([]);B = class(B,'wafer2d',A);