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