Psykrometri och Värmeöverföring inom termodynamik och beräkning av temperatur vid jämvikt med matlab

Hej, jag har fastnat lite på en inlämningsuppgift inom termodynamik som jag arbetar med just nu. Specifikt försöker jag beräkna yttemperaturen på vattenfångaren under natten vid jämvikt. Jag är lite osäker på hur man ska tänka med funktionen för jämvikt, jag får fel svar någonstans men vet inte om problemet ligger i min jämviktsfunktion eller tidigare beräkningar i koden. Bifogar både kod och uppgiftsbeskrivningen

function termo6

tdag = 40; %[Celsius]
tnatt = 16; %[Celsius]

epsilon = 0.92; %Vattenfångarens emmissionstal
A = 1; %Vattenfångarens area
r = 2480e3; %Vattnets ångbildningsvärme
l = 335e3; %Vattnets smältvärme


Mv = 18;
Ma = 29;

sigma = 5.67e-8; %boltzman            
phidag = 0.2; %Luftens relativ fuktighet vid dagtid
alpha = 5; %luftens konvektiva värmeövergångstal
F12 = epsilon; %utstrålning nattetid
pbisv = 0.073750e5;
pv = pbisv*phidag
ptot = 1.01325e5;

Xw = (Mv/Ma)*(pv/(ptot-pv)) %vatteninnehåll
Xbisw = (Mv/Ma)*(pbisv/(ptot-pv)); %vatteninnehåll mättad ånga

Xwm = (Xbisw+Xw)/2; 

CpA = 1.006e3;
Cpv = 1.86e3;

Cpx1 = CpA + (Xwm*Cpv);
sigV = alpha/Cpx1;
mdot = sigV*A*(Xbisw-Xw);

thimmel = -80; %[Celsius]

f = @(Tp) (alpha*(Tp - tnatt+273.15)) - (F12*sigma*(((Tp^4) - (thimmel+273.15)^4))) + mdot*(r);



Tpstart = tnatt+273.15-10;
Tpsvar = fzero(f, Tpstart);
tp = Tpsvar-273.15

end