Dear Mattia,
I leave it to Simon to comment on the covariance matrix part.
The error you report now is different. The iterations give you a water
vapour VMR that is out of the range of your lookup table. If you look at
the (first) error message, you can see that you get a H2O VMR of 2.9% at
1 Pa. This is not realistic. The reference value is 3.1 ppm, which looks
OK.
In the text you indicate that you retrieve H2O only at 101540, 101420,
101170, 100890 Pa. I don't know what measurements you use, but this
seems odd. Anyhow, please note that the value you retrieve at the end
point of the pressure retrieval grid, is applied all the way to the top
of the atmosphere.
Bye,
Patrick
On 2023-06-05 17:04, Mattia Sabatini wrote:
Hi Simon,
as you suggested, I tried adding the blocks separately. This time I've
also changed the variables to retrieve: they are the surface temperature
and the water vapour concetration at 4 pressure levels (101540, 101420,
101170, 100890 Pa).
-------------------------------------------
MatrixCreate(m)
# Surface temperature block
retrievalAddSurfaceQuantity(g1=lat_grid, g2=lon_grid, quantity="Skin
temperature")
ReadXML(m, "block1.xml")
covmat_sxAddBlock(block=m, i=0, j=0)
# Water vapour block
retrievalAddAbsSpecies(species="H2O, H2O-SelfContCKDMT400,
H2O-ForeignContCKDMT400", unit="vmr", g1=p_ret_grid, g2=lat_grid,
g3=lon_grid)
ReadXML(m, "block2.xml")
covmat_sxAddBlock(block=m, i=1, j=1)
# Upper off-diagonal block
ReadXML(m, "block3.xml")
covmat_sxAddBlock(block=m, i=0, j=1)
-------------------------------------------
Now I get this error:
-------------------------------------------
Run-time error in method: OEM
Run-time error in agenda: inversion_iterate_agenda
Run-time error in method: yCalc
User Error: failed
Error is found at:
/home/Mattia.Sabatini/arts-2.5.8/src/m_rte.cc:2774
Please follow these instructions to correct your error:
User Error: failed
Error is found at:
/home/Mattia.Sabatini/arts-2.5.8/src/rte.cc:1749
Please follow these instructions to correct your error:
Run-time error in function: iyb_calc
Run-time error in agenda: iy_main_agenda
Run-time error in method: iyEmissionStandard
User Error: do_abort
Error is found at:
/home/Mattia.Sabatini/arts-2.5.8/src/m_rte.cc:1696
Please follow these instructions to correct your error:
Error messages from failed cases:
Runtime-error in source calculation at index 0:
Run-time error in agenda: propmat_clearsky_agenda
Run-time error in method: propmat_clearskyAddFromLookup
Problem with gas absorption lookup table.
VMR for H2O (species 0) is outside the range covered by the lookup table.
Your VMR was 0.02937 at a pressure of 1.0002 Pa.
The reference VMR value there is 3.05987e-06
The fractional VMR relative to the reference value is 9598.46.
The allowed range is -4.46834 to 49.1517.
The fractional VMR perturbation grid range in the table is 0 to 44.6834.
We allow a bit of extrapolation, but NOT SO MUCH! Runtime-error in
source calculation at index 54:
Run-time error in agenda: propmat_clearsky_agenda
Run-time error in method: propmat_clearskyAddFromLookup
Problem with gas absorption lookup table.
VMR for H2O (species 0) is outside the range covered by the lookup table.
Your VMR was 0.02937 at a pressure of 7111.9 Pa.
The reference VMR value there is 4.42132e-06
The fractional VMR relative to the reference value is 6642.81.
The allowed range is -4.46834 to 49.1517.
The fractional VMR perturbation grid range in the table is 0 to 44.6834.
We allow a bit of extrapolation, but NOT SO MUCH! Runtime-error in
source calculation at index 18:
Run-time error in agenda: propmat_clearsky_agenda
Run-time error in method: propmat_clearskyAddFromLookup
Problem with gas absorption lookup table.
VMR for H2O (species 0) is outside the range covered by the lookup table.
Your VMR was 0.02937 at a pressure of 209.65 Pa.
The reference VMR value there is 5.3185e-06
The fractional VMR relative to the reference value is 5522.23.
The allowed range is -4.46834 to 49.1517.
The fractional VMR perturbation grid range in the table is 0 to 44.6834.
We allow a bit of extrapolation, but NOT SO MUCH! Runtime-error in
source calculation at index 36:
Run-time error in agenda: propmat_clearsky_agenda
Run-time error in method: propmat_clearskyAddFromLookup
Problem with gas absorption lookup table.
VMR for H2O (species 0) is outside the range covered by the lookup table.
Your VMR was 0.02937 at a pressure of 1986.8 Pa.
The reference VMR value there is 4.17148e-06
The fractional VMR relative to the reference value is 7040.67.
The allowed range is -4.46834 to 49.1517.
The fractional VMR perturbation grid range in the table is 0 to 44.6834.
We allow a bit of extrapolation, but NOT SO MUCH! Runtime-error in
source calculation at index 71:
Run-time error in agenda: propmat_clearsky_agenda
Run-time error in method: propmat_clearskyAddFromLookup
Problem with gas absorption lookup table.
VMR for H2O (species 0) is outside the range covered by the lookup table.
Your VMR was 0.02937 at a pressure of 17119 Pa.
The reference VMR value there is 7.45595e-05
The fractional VMR relative to the reference value is 393.914.
The allowed range is -4.46834 to 49.1517.
The fractional VMR perturbation grid range in the table is 0 to 44.6834.
We allow a bit of extrapolation, but NOT SO MUCH! Runtime-error in
source calculation at index 88:
Run-time error in agenda: propmat_clearsky_agenda
Run-time error in method: propmat_clearskyAddFromLookup
Problem with gas absorption lookup table.
VMR for H2O (species 0) is outside the range covered by the lookup table.
Your VMR was 0.02937 at a pressure of 36743 Pa.
The reference VMR value there is 0.000457873
The fractional VMR relative to the reference value is 64.1444.
The allowed range is -4.46834 to 49.1517.
The fractional VMR perturbation grid range in the table is 0 to 44.6834.
We allow a bit of extrapolation, but NOT SO MUCH!
Stopping ARTS execution.
Goodbye.
-------------------------------------------
I can still open an issue with other details about the controlfile, if
it is better. Thank you!
Mattia
Il 02/06/2023 16:43, Pfreundschuh,Simon ha scritto:
Hi Mattia,
Sorry for the delayed reply. I forgot to subscribe to the mailing list
after changing
positions.
This does look like a bug. Would you mind creating an issue at
https://github.com/atmtools/arts <https://github.com/atmtools/arts?>?
It should be an easy fix, so I could have this ready early next week.
If you don't want to update ARTS, the appropriate work around would be
to add the blocks,
i.e., the two on the diagonal and the upper off-diagonal block,
separately.
Let me know if you have any further questions.
Kind regards,
Simon Pfreundschuh
Postdoctoral researcher
Department of Atmospheric Science
Colorado State University
