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 
> <mattia.sabat...@artov.ismar.cnr.it> 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 
> <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

Reply via email to