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
