Hi Jana,
thanks a lot. By changing the spectroscopic line data from the default
lines.xml to the cat folder in arts-xml-data/spectroscopy/, I correctly obtain
the cross section values for all the abs species now. Thank you again.
Mattia
Da: Jana Mendrok
Inviato: lunedì 14 febbraio 2022 12:10
A: Mattia Sabatini
Cc: Pengwang Zhai; Lemke, Oliver; ARTS Users List
Oggetto: Re: [arts-users] R: Cannot create abs_lookup with arts 2.5
Hi Mattia, Pengwang,
are you using an appropriate line catalog file (H2O-PWR98 is a so-called "full
absorption model" with its parameters implemented within ARTS, while "H2O"
triggers a line-by-line calculation that needs the line parameters as input). I
am not sure whether that specific test case comes with (or applies) a suitable
line file...
Best,
Jana
On Mon, Feb 14, 2022 at 11:00 AM Mattia Sabatini
<mattia.sabat...@artov.ismar.cnr.it> wrote:
Hello again Pengwang,
I tested controlfile TestAbs.arts with your edit:
AgendaSet( abs_xsec_agenda ){
abs_xsec_per_speciesInit
abs_xsec_per_speciesAddConts
lbl_checkedCalc
abs_xsec_per_speciesAddLines
}
And with the following abs species:
abs_speciesSet( species=[ "H2O", "H2O-PWR98", "O2-PWR93",
"N2-SelfContStandardType" ] ).
I thought that by adding abs_xsec_per_speciesAddLines in abs_xsec_agenda, just
like Richard suggested me to do in my post
(https://www.mail-archive.com/arts_users.mi@lists.uni-hamburg.de/msg00471.html),
would solve the problem you were experiencing of having a lookup table filled
with zeros. Unfortunately this is still an issue for me, even with the above
settings: the H2O cross section values are all zeros, while values for
"H2O-PWR98", "O2-PWR93", "N2-SelfContStandardType” are not. Also, during the
execution I had this “Deprecated function warning”:
abs_xsec_per_speciesAddLines is deprecated since 2021-07-13
This function is no longer up to date. It only exists to satisfy lookup table
calculations before these are updated.
Once the lookup table calculations are up-to-date, this function is fully
replaced with propmat_clearskyAddLines, with better functionality
Apparently, abs_xsec_per_speciesAddConts is doing its job but
abs_xsec_per_speciesAddLines does not. Did you had the chance to test it too?
Mattia
Thanks, Oliver and Mattia, for you help.
After studying Mattia and Richard's posts, I figured out the following:
Comment out line 15 of TestAbs.arts:
#Copy(abs_xsec_agenda, abs_xsec_agenda__noCIA)
And insert the following to line 16:
AgendaSet( abs_xsec_agenda ){
abs_xsec_per_speciesInit
abs_xsec_per_speciesAddConts
lbl_checkedCalc
abs_xsec_per_speciesAddLines
}
Now it seems working, though I have not yet got the chance to check the outputs
as it is still running.
Yours
Pengwang
> On Feb 11, 2022, at 2:06 AM, Lemke, Oliver <oliver.le...@uni-hamburg.de>
> wrote:
>
> Hi Mattia, hi Pengwang,
>
> Thanks Mattia for helping out. :-)
>
> Since you were referring to Richard's earlier post, I'll take this
> opportunity to point out the searchable archive of this list, which might
> come in handy at times:
>
> https://www.mail-archive.com/arts_users.mi@lists.uni-hamburg.de/
>
> Here is the post Mattia was referring to:
>
> https://www.mail-archive.com/arts_users.mi@lists.uni-hamburg.de/msg00471.html
>
> Cheers,
> Oliver
>
>
>> On 10 Feb 2022, at 17:38, mattia.sabat...@artov.ismar.cnr.it wrote:
>>
>> Pengwang,
>>
>> sorry for my partial answer. Recently I compiled ARTS 2.5 and I had a
>> similar issue concerning the absorption lookup table. Richard Larsson,
>> helped me replying to my mail sent in this mailing list, suggesting to add
>> abs_xsec_per_speciesAddLines into my cross-section agenda.
>>
>> He also added:
>>
>> "We are currently in a transition in the 2.5-branch of moving away from
>> abs_xsec_agenda entirely. There are some lingering problems when using
>> lookup table calculations at this time. Mainly, since the line calculations
>> should happen in propmat_clearsky_agenda now, they are not part of any of
>> the default cross-section agendas as in the past. The lookup generation
>> interface is currently under active development, so this interface might
>> change in the near future."
>>
>> Best regards,
>> Mattia
>>
>>
>>
>>
>>
>> Pengwang Zhai <pwz...@umbc.edu> ha scritto:
>>
>>> Thanks, Mattia. Would you advise how TestAbs.arts can be revised to include
>>> line-by-line calculation?
>>>
>>> Note that the behavior of TestAbs.arts is different in arts2.3, which does
>>> calculate line-by-line absorption coefficients.
>>>
>>> I do not want continua, as I am mainly interested in the visible spectra,
>>> which is out of most of those continuum models. Based on my experience with
>>> previous arts version, arts does not check the spectrum limits for the
>>> continuum models. I am not sure whether arts 2.5 has been improved on this.
>>>
>>> Pengwang
>>>
>>>
>>>
>>>
>>>
>>>
>>>> On Feb 10, 2022, at 9:46 AM, mattia.sabat...@artov.ismar.cnr.it wrote:
>>>>
>>>> Hi Pengwang,
>>>>
>>>> the controlfile TestAbs.arts, in line 15, has the following:
>>>> Copy(abs_xsec_agenda, abs_xsec_agenda__noCIA)
>>>>
>>>> I checked in controlfiles/general/agendas.arts for what
>>>> abs_xsec_agenda__noCIA does:
>>>> AgendaCreate( abs_xsec_agenda__noCIA )
>>>> AgendaSet( abs_xsec_agenda__noCIA ){
>>>> abs_xsec_per_speciesInit
>>>> abs_xsec_per_speciesAddConts
>>>> }
>>>>
>>>> Therefore the agenda calculates absorption for continua tag only (see
>>>> https://atmtools.github.io/arts-docs-master/docserver/methods/abs_xsec_per_speciesAddConts.html).
>>>> You selected H2O, and by doing this my guess is that you are not
>>>> considering continuum, as it is written in lines 36-37 of TestAbs.arts.
>>>>
>>>> I hope that this will help you,
>>>>
>>>> Mattia
>>>>
>>>>
>>>> Pengwang Zhai <pwz...@umbc.edu> ha scritto:
>>>>
>>>>> Hello, ARTS community,
>>>>>
>>>>> I downloaded and compiled the latest version of arts. Now I tested the
>>>>> the creation of abs_lookup with the example control file located in:
>>>>>
>>>>> arts/controlfiles/artscomponents/absorption/TestAbs.arts
>>>>>
>>>>> I only modified two occurrences of:
>>>>>
>>>>> abs_speciesSet( species=[ "H2O-PWR98",
>>>>> "O2-PWR93",
>>>>> "N2-SelfContStandardType" ] )
>>>>>
>>>>> to:
>>>>>
>>>>> abs_speciesSet( species=[ "H2O" ] )
>>>>>
>>>>> and run
>>>>>
>>>>> arts TestAbs.arts
>>>>>
>>>>> The resultant abs_lookup are all ZEROs.
>>>>>
>>>>> Any help?
>>>>>
>>>>> More background information: I used arts 2.3 to create abs_lookup for
>>>>> H2O, CO2, etc. in the visible by reading from HITRAN 2012. The baseline
>>>>> example was TestAbs.arts provided by the arts installation. Now I need to
>>>>> revisit the calculation based on HITRAN 2020, and the lookup table
>>>>> calculation seems not working with arts 2.5. I greatly appreciate it if
>>>>> you could provide a working example on how to create a abs_lookup with
>>>>> arts 2.5.
>>>>>
>>>>> Pengwang
>
_______________________________________________
arts_users.mi mailing list
arts_users.mi@lists.uni-hamburg.de
https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi
--
Jana Mendrok, Ph.D.
Deutscher Wetterdienst
Offenbach am Main, Germany
+49 (0)69 8062 3139
_______________________________________________
arts_users.mi mailing list
arts_users.mi@lists.uni-hamburg.de
https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi
