function B = radialrxrmodel(A,vo,Ts)% radialrxrmodel.m A nonlinear equation model object constructor, derived%           from the neqmodel class. Called as%%          B = radialrxrmodel(A)%% INPUT PARAMETERS%%         A : a radialrxrgeom objectunpack(A)if nargin < 2    vo = 20;   % inlet velocity (m/s)endif nargin < 3    Ts = 1000; % susceptor temperature (K)endxAo = 0.1;k1  = 1e19;E1  = 3e4;xBo = 0;k2  = 1e15;E2  = 2e4;xCo = 0;kd  = 100;Tamb = 300; % ambient temperature, Kh = 1; % heat transfer coefficientT  = scalarfield(R,(Tamb+Ts)/2); % Temperature solution initial guessxA = scalarfield(R,xAo);         % Species A profile solution initial guessxB = scalarfield(R,xAo/2);       % Species B profile solution initial guessxC = scalarfield(R,xAo/2);       % Species C profile solution initial guess%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% set up an instance of the rxnmodel objectvar   = assocarray({ 'T' T 'xA' xA 'xB' xB 'xC' xC });param = assocarray({ 'R' R 'RSF' RSF 'nrp' nrp ...    'Rs' Rs 'Rp' Rp 'Rw' Rw 'Ro' Ro ...    'vo' vo 'Tamb' Tamb 'Ts' Ts 'h' h ...    'xAo' xAo 'k1' k1 'E1' E1 ...    'xBo' xBo 'k2' k2 'E2' E2 ...    'xCo' xCo 'kd' kd  ...    });A = naemodel(var,param);B = struct([]);B = class(B,'radialrxrmodel',A);