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 model
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 wro
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 controlfil
[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&rmjob=10254&rmlang=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&rmjob=10255&rmlang=UK ___ arts_users.mi mailing list arts_users.mi@lists.uni-hamburg.de https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi