Re: [arts-users] Cannot create abs_lookup with arts 2.5
Dear all,
The deprecation warning is there because we will remove
abs_xsec_per_speciesAddLines before Arts 3.0 is released. It is just
warning you that some code will have to change in the future. Ignore it
for the purpose of what you want to do here.
About the speed. The new code should be comparable and often faster than
the old code (depending on your line-data inputs).
There are two caveats to that though.
1) Hitran2020 is much larger than previous copies.
2) You will not see parallelism unless you have enough temperature
perturbations. There is a pending pull-request to fix this.
As Oliver points out, the recommended data to use is the
arts-cat-data/lines/ data. This is based on Hitran2020 with a few
modifications we found necessary to improve the quality of the catalog.
//Richard
Den mån 14 feb. 2022 kl 13:49 skrev Pengwang Zhai :
> Hello, Mattia,
>
> Thanks for the update. I have to admit that I did not use the line.xml
> file. I used the H2O file downloaded from HITRAN 2020. Originally I said
> “it seems working”, today I still cannot definitely confirm whether it
> works from my side, as ARTS 2.5 did not finish the calculation after
> several days. Maybe I am trying to calculate too many frequency lines,
> though I believe ARTS 2.3 only took several hours to build the same lookup
> table.
>
> I did see the same warning that you saw regarding
> abs_xsec_per_speciesAddLines.
>
> Yours,
>
> Pengwang
>
>
> Pengwang Zhai (he/his/him)
> Graduate Program Director, Atmospheric Physics,
> Associate Professor, Physics Department, UMBC
> pwz...@umbc.edu
>
>
>
>
>
>
> > On Feb 14, 2022, at 5: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
Re: [arts-users] Cannot create abs_lookup with arts 2.5
Hello, Mattia,
Thanks for the update. I have to admit that I did not use the line.xml file. I
used the H2O file downloaded from HITRAN 2020. Originally I said “it seems
working”, today I still cannot definitely confirm whether it works from my
side, as ARTS 2.5 did not finish the calculation after several days. Maybe I am
trying to calculate too many frequency lines, though I believe ARTS 2.3 only
took several hours to build the same lookup table.
I did see the same warning that you saw regarding abs_xsec_per_speciesAddLines.
Yours,
Pengwang
Pengwang Zhai (he/his/him)
Graduate Program Director, Atmospheric Physics,
Associate Professor, Physics Department, UMBC
pwz...@umbc.edu
> On Feb 14, 2022, at 5:00 AM, Mattia Sabatini
> 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
> > 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 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
Re: [arts-users] Cannot create abs_lookup with arts 2.5
Hi Mattia,
For ARTS 2.5, the line catalog has been updated and moved from arts-xml-data
into its own package arts-cat-data. It's best to use this updated version with
the development version of ARTS. You can obtain it from SVN:
svn co https://arts.mi.uni-hamburg.de/svn/rt/arts-cat-data/trunk/ arts-cat-data
You find the catalog in the "lines" subdirectory.
Cheers,
Oliver
> On 14 Feb 2022, at 12:23, Mattia Sabatini
> wrote:
>
> 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
> 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
> > 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
> >>
> >>
[arts-users] R: R: Cannot create abs_lookup with arts 2.5
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
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
> 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 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
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
> 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 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
[arts-users] R: Cannot create abs_lookup with arts 2.5
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
> 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 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 ha scritto:
> Hello, ARTS community,
[arts-users] PhD position in sub-mm observations of clouds
Hi all, We have now a new PhD student position at Chalmers University of Technology in sub-millimetre observations of clouds This position is related to two completely new satellite instruments ICI (Ice Cloud Imager)and AWS (Arctic Weather Satellite). You will take part in the final pre-launch preparations and analyse the first data from the satellites. See further: https://www.chalmers.se/en/about-chalmers/Working-at-Chalmers/Vacancies/Pages/default.aspx?rmpage=job=10254=UK Contact: Patrick Eriksson, patrick.eriks...@chalmers.se There is also a new PhD position in radar remote sensing of forests https://www.chalmers.se/en/about-chalmers/Working-at-Chalmers/Vacancies/Pages/default.aspx?rmpage=job=10255=UK ___ arts_users.mi mailing list arts_users.mi@lists.uni-hamburg.de https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi
