Hi, Fran

(I'm adding the qpack mailing list here as well)

could you attach your Q file and perhaps the script you use to execute the 
retrievals? I assume this problem could be due to some mismatch between the 
number of channels specified in you Y.F (qpack Y struckture) and 
Q.SENSOR_RESPONSE.

Ole Martin
________________________________________
From: [email protected] [[email protected]] on 
behalf of Fran Navas [[email protected]]
Sent: Wednesday, July 15, 2015 10:15 AM
To: ARTS Users List
Subject: [arts-users]  arts 2.2

Dear Patrick, Jana and all,
I am having some problems running our tropospheric temperature
retrievals (from filter bank measurements) in ARTS 2.2. These retrievals
were running in arts 2.0 until now. I have adapted the scripts in order
to run it in the last version and now everything look fine until the
different iterations start in the retrieval.

Below you can find the error message. I do not why but the dimensions of
the 'yf' that is calculated in the function: [R,yf,J] = calc_y_and_j(
comfun, Q, R, O, xnorm, x, iter ); are different to the ones from 'y'.
In my retrievals (y: [108*1]; yf: [126*1]). This difference in the
dimensions is the cause of the error. I guess that something is wrong in
the Q-structure but I don't find what can be. I would appreciate a lot
if someone could help me.

Thank you very much in advance!!!
Fran Navas


Error in ==> oem>calc_cost at 945
     dd       = y - yf;

Error in ==> oem>calc_cost_for_X at 956
     [cost,cost_y,cost_x] = calc_cost( y, yf, Seinv, xa, x, Sxinv, O,
xnorm );

Error in ==> oem at 558
           X = calc_cost_for_X( y, yf, Seinv, xa, x, Sxinv, O, xnorm, X,
iter );

Error in ==> qpack2 at 116
     [X,R] = oem( O, Qm, R, @arts_oem, Sx, Se, Sxinv, [], xa, Y(m).Y );





--
*****************************************************************************
Francisco Navas Guzmán, PhD

Microwave Physics Group
Division of Atmospheric radiometry
Institute of Applied Physics
University of Bern

Address:                                    | http://www.iap.unibe.ch
Sidlerstrasse 5                             | Ph:  +41-31 631 50 46
CH-3012 Bern (Switzerland)                  | Fax: +41-31 631 37 65
*****************************************************************************

_______________________________________________
arts-users mailing list
[email protected]
https://www.sat.ltu.se/mailman/listinfo/arts-users
_______________________________________________
qpack mailing list
[email protected]
https://www.sat.ltu.se/mailman/listinfo/qpack

Reply via email to