function B = wafertempcontroller(N,Rw,Tset)

% wafertempcontroller.m A nonlinear equation model object constructor, derived
%                       from the neqmodel class. Called as
%
%            A = lampflux(N,Rw)
%
% INPUT PARAMETERS
%
%            N : number of quadrature points
%           Rw : wafer radius, (m)
%         Tset : mean wafer temperature setpoint


M = mwrmodel('cyln',N,'r',[0 Rw]);
S = get(M,'R');

power = 0.5; % initial guess for power input fraction
T     = scalarfield(S,Tset);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% set up an instance of the lampflux object

var   = assocarray({'power' power});
param = assocarray({ 'T' T,'Tset',Tset});

A = naemodel('wafertempcontroller',var,param);
B = struct([]);

B = class(B,'wafertempcontroller',A);