Dear all,
i updated our retrieval of the ozone line @142 GHz to Arts2/Qpack2 (Bern
radiometer)
and so far i used a constant measurement noise on our spectrum Y.Y.
The problem is, when i define a measurement noise with the length of my
frequency vector,
I get the following error message,
??? Undefined function or method 'tnoise_si' for input arguments of type
'double'.
Error in ==> qpack2>qp2_y2Q at 321
s(k) = prod( tnoise_si([i(k) j(k)]) ) * s(k);
Error in ==> qpack2 at 99
[Qm,Se] = qp2_y2Q( Q, Y, m );
Error in ==> qpack2_retrieval_45ch_C at 61
L2 = qpack2( Q, O, Y );
-> qpack2.m / qp2_y2Q.m (lines #315-324) :
if n == 1
Se = (Y(m).TNOISE^2) .* Q.TNOISE_C;
else
% Avoid to generate full matrix that is n x n:
[i,j,s] = find( Q.TNOISE_C );
for k = 1 : length(i)
s(k) = prod( tnoise_si([i(k) j(k)]) ) * s(k);
end
Se = sparse(i,j,s);
end
So far, i didn't encounter the problem since my "n was == 1", and Se set
correctly.
Am I doing something wrong?
Or how can I fix the problem, if I don't wanna change the thermal noise
correlation T_NOISE_C,
since I'm assuming that their is no correlation of thermal noise between
channels?
Thank you and
best regards,
Simone_______________________________________________
qpack mailing list
[email protected]
https://www.sat.ltu.se/mailman/listinfo/qpack
