qpack uses some information from the Y structure to define the sensor, not
exactly sure which, but check that Y.F is equal to your backend frequencies in
Q.F_BACKEND. Furthermore I believe the zenith angle and position (altitude) of
the sensor is specified by Y so check that these are consistent.
If you don't find and discrepancies you can check the raw files generated by
qpack (and qarts) by running atmlab in debug mode (atmlab('DEBUG',1)). The
controlfiles used for the arts run is then stored on the harddisk, and you can
check that they are identical.
My 2 cents.
Ole Martin
________________________________
From: [email protected] [[email protected]] on behalf of
Pauline Martinet [[email protected]]
Sent: Friday, September 04, 2015 5:45 PM
To: [email protected]
Subject: [Qpack] differences between Qpack and arts simulation
Hi everybody,
I have noticed differences between the forward model calculations when I call
Qpack or the function arts_y.
I am simulating a microwave radiometer between 22 and 60 GHz.
I have a matlab routine where I define the structures Q and O.
Then I call the two functions:
TB=arts_y(Q);
and L2=qpack2(Q,O,Y);
with Y.Y=[] : from what I read in the Qpack documentation it should give the
same results.
But then I do not obtain the same values of brightness temperatures when
simulating the brightness temperature with Qpack or arts_y:
This is from Qpack:
37.6242
36.7163
32.3708
24.4462
22.0529
19.4580
18.6558
110.6579
152.7556
250.9260
279.8277
284.9573
285.2151
285.3064
This is from arts_y:
37.5430
36.6337
32.2927
24.3829
21.9958
19.4087
18.6113
110.5074
152.5850
250.8197
279.7959
284.9563
285.2181
285.3116
For the call to arts-y I had to define: Q.T_FIELD, Q.Z_FIELD, Q.VMR_FIELD
but I normally use the same values as the ones defined in Q.ABS_SPECIES and Q.T
for Qpack.
Do you have any idea on what I am doing wrong ?
Thank you very much.
Best regards,
Pauline
----- Météo-France -----
Dr. Pauline Martinet
Chercheur CNRM/GMEI/LISA
[email protected]
Fixe : +33 561079031
Site web: www.sites.google.com/site/martinetpauline31
<https://www.sites.google.com/site/martinetpauline31/>
_______________________________________________
qpack mailing list
[email protected]
https://www.sat.ltu.se/mailman/listinfo/qpack
