Re: [arts-users] calculate absorption coefficients from abs_coefCalc

2016-02-16 Thread Pengwang Zhai
Thanks, Richard.  I tried your suggestion.  Unfortunately, it does not work. 
The error message is:

arts TestAbsCoeff.arts

Version: arts-2.2.57 (compiled Fri Feb 12 21:05:27 2016)
Executing Arts
{…
…
…
}
This run took 11.10s (11.07s CPU time)
Run-time error in controlfile: TestAbsCoeff.arts
Run-time error in method: AgendaExecute
Run-time error in method: abs_xsec_per_speciesInit
Method abs_xsec_per_speciesInit needs input variable: abs_species_active
Stopping ARTS execution.
Goodbye.

Here is the script:




#DEFINITIONS:  -*-sh-*-
#
# An example ARTS controlfile that calculates absorption
# coefficients. 
# SAB 16.06.2000

Arts2 {

INCLUDE "general/general.arts"
INCLUDE "general/continua.arts"
INCLUDE "general/agendas.arts"
INCLUDE "general/planet_earth.arts"

# Agenda for scalar gas absorption calculation
Copy(abs_xsec_agenda, abs_xsec_agenda__noCIA)
# on-the-fly absorption
Copy( propmat_clearsky_agenda, propmat_clearsky_agenda__OnTheFly )

abs_linesReadFromHitran( abs_lines,
   "/RT/HITRAN/HITRAN2012.par",
3.880627143214052e+14,
3.961610390054004e+14 )

abs_speciesSet( species=[ "H2O-PWR98",
  "O2-PWR93" ] )

# This separates the lines into the different tag groups and creates
# the workspace variable `abs_lines_per_species':
abs_lines_per_speciesCreateFromLines

# Dimensionality of the atmosphere
AtmosphereSet1D

VectorNLogSpace( p_grid, 200, 10, 10 )
 
# Atmospheric profiles
AtmRawRead( t_field_raw, z_field_raw, vmr_field_raw, abs_species, 
"testdata/tropical" )

# Now interpolate all the raw atmospheric input onto the pressure 
# grid and create the atmospheric variables `t_field', `z_field', `vmr_field'
AtmFieldsCalc

# Initialize the input variables of abs_coefCalc from the Atm fields:
AbsInputFromAtmFields

# Create the frequency grid `f_grid':
VectorNLinSpace( f_grid, 25000,3.880627143214052e+14, 3.961610390054004e+14 )

IndexSet(stokes_dim, 1)

# Calculate absorption coefficients, both total (`abs_coef') and
# separately for each tag group (`abs_coef_per_species'):

atmfields_checkedCalc
abs_xsec_agenda_checkedCalc
AgendaExecute(abs_xsec_agenda)
abs_coefCalcFromXsec

# Optionally write these to files:
WriteXML( output_file_format,abs_coef_per_species)
}



> On Feb 16, 2016, at 3:15 AM, Richard Larsson <ric.lars...@gmail.com> wrote:
> 
> Dear Pengwang Zhai,
> 
> I think we still have what you are looking for in ARTS.
> 
> The variable you are actually interested in is not abs_coef, but 
> abs_coef_per_species.
> 
> The way to create this is to first run abs_xsec_agenda and then 
> abs_coefCalcFromXsec.  
> 
> I think your code should only change from:
> 
> abs_coefCalc
> 
> to
> 
> abs_xsec_agendaExecute
> abs_coefCalcFromXsec
> 
> Now abs_coef_per_species is an array of matrices.  The array-size is the 
> species and the matrices are the size of f_grid and p_grid.  Save this and do 
> what you want with it.
> 
> This said, there is potentially a few other changes depending on what version 
> you are using and what species you are interested in.  You will want to run 
> jacobianOff, for instance, if you are using the dev-branch.
> 
> Cheers,
> //Richard
> 
> 2016-02-15 22:11 GMT+01:00 Pengwang Zhai <pwz...@umbc.edu>:
> Thanks, Jana.
> 
> > Note that for deriving abs_coef for this, you need to sum up over all 
> > species (i.e. over the 0th dimension of propmat_clearsky_field.
> 
> This is exactly what I wanted to do within arts to get “abs_coef”.  If it is 
> not possible with arts, I will use matlab to do this.
> 
> The three dimensions are: [species, f_grid, p_grid]
> 
> Cheers,
> 
> PZ
> 
> 
> > On Feb 15, 2016, at 3:47 PM, Jana Mendrok <jana.mend...@gmail.com> wrote:
> >
> > Hi,
> >
> > what do you intend to do with abs_coeff_user?
> > Do you really process this further within ARTS itself? only then it makes 
> > sense to process propmat_clearsky_field into a Tensor3*, i think.
> > Else, I strongly recommend to do any re-shaping / reduction outside of 
> > ARTS; other programming languages are much better suited for this kind of 
> > task (e.g. you can use the matlab-interface atmlab to have easy access to 
> > ARTS output. or the python interface typhon).
> >
> > abs_coef was a Matrix of dimension [f_grid, abs_p] (and 
> > abs_coef_per_species an Array holding one abs_coef matrix per defined 
> > abs_species).
> > propmat_clearsky_field is a Tensor7 of dimension [species, f_grid, 
> > stokes_dim, stokes_dim, p_grid, lat_grid, lon_grid].
> >
> > That is, for reducing propmat_clearsky_field to what was 
> > abs_coef_pe

Re: [arts-users] Compatibility with HITRAN2016

2018-01-10 Thread Pengwang Zhai
Thanks, Oliver. I tried your recipe and it works fine to me.

Another small thing (or bug) is that arts LUT function sometimes generates 
small negative values for absorption coefficient, when the spectral regions is 
slightly outside an absorption band. If the developers are interested in fixing 
this problem, I would be happy to provide more info that I have.

Yours,

Pengwang
 

> On Jan 9, 2018, at 3:01 AM, Oliver Lemke <oliver.le...@uni-hamburg.de> wrote:
> 
> Hi Pengwang,
> 
> ARTS hasn't been tested extensively with HITRAN 2016 and thus we don't 
> officially support it yet. However, in general it works. You can work around 
> the error you're seeing by renaming your catalog file and prepend HITRAN2012_ 
> to the filename. In your case for example rename it to 
> HITRAN2012_5a4fe828.par and ARTS will read it.
> 
> When using HITRAN 2016 make sure you're using the latest version from the 
> Subversion repository.
> 
> Cheers,
> Oliver
> 
> 
>> On 6 Jan 2018, at 22:21, Pengwang Zhai <pwz...@umbc.edu> wrote:
>> 
>> Dear ARTS developers:
>> 
>> ARTS is an invaluable tool for the RT community and thanks very much for 
>> your great contribution.
>> 
>> I have recently run into a problem with HITRAN 2016, the newest version of 
>> HITRAN. I have downloaded the O2 par file from the HITRAN website. It is 
>> named as: 5a4fe828.par. I used the example script to calculate the 
>> absorption Look-Up-Table. One function is:
>> 
>> abs_linesReadFromHitran( abs_lines,
>>  “~/Oxygen_HITRAN2016/5a4fe828.par",
>>   1.3e+14,
>>   1.0e+15 )
>> 
>> Then I got the following error message:
>> 
>> Run-time error in method: abs_linesReadFromHitran
>> ‘~/Oxygen_HITRAN2016/5a4fe828.par'
>> does not appear to be a HITRAN 2012 catalogue. The catalog filename
>> must start with HITRAN2012. If you intend to use a HITRAN 2008 catalog
>> run 'cmake -DWITH_HITRAN2008=1 ..' and recompile arts
>> Stopping ARTS execution.
>> 
>> Any suggestion to make this work? Have you attempted to confirm the 
>> compatibility with HITRAN 2016?
>> 
>> My arts version is: arts-2.3.657.
>> 
>> Thanks very much and have a great 2018!
>> 
>> Yours,
>> 
>> Pengwang
> 

___
arts_users.mi mailing list
arts_users.mi@lists.uni-hamburg.de
https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi


[arts-users] Compatibility with HITRAN2016

2018-01-06 Thread Pengwang Zhai
Dear ARTS developers:

ARTS is an invaluable tool for the RT community and thanks very much for your 
great contribution.

I have recently run into a problem with HITRAN 2016, the newest version of 
HITRAN. I have downloaded the O2 par file from the HITRAN website. It is named 
as: 5a4fe828.par. I used the example script to calculate the absorption 
Look-Up-Table. One function is:

abs_linesReadFromHitran( abs_lines,
   “~/Oxygen_HITRAN2016/5a4fe828.par",
1.3e+14,
1.0e+15 )

Then I got the following error message:

Run-time error in method: abs_linesReadFromHitran
‘~/Oxygen_HITRAN2016/5a4fe828.par'
does not appear to be a HITRAN 2012 catalogue. The catalog filename
must start with HITRAN2012. If you intend to use a HITRAN 2008 catalog
run 'cmake -DWITH_HITRAN2008=1 ..' and recompile arts
Stopping ARTS execution.

Any suggestion to make this work? Have you attempted to confirm the 
compatibility with HITRAN 2016?

My arts version is: arts-2.3.657.

Thanks very much and have a great 2018!

Yours,

Pengwang

___
arts_users.mi mailing list
arts_users.mi@lists.uni-hamburg.de
https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi


[arts-users] Cannot create abs_lookup with arts 2.5

2022-02-10 Thread Pengwang Zhai
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


Re: [arts-users] Cannot create abs_lookup with arts 2.5

2022-02-14 Thread 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 
>  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 regar

[arts-users] read in HITRAN par file

2022-02-03 Thread Pengwang Zhai
Hello, ARTS Community,

I used arts 2.3 to generate gas absorption coefficient lookup table for o2, 
co2, etc. One key function was: "abs_linesReadFromHitran”. Now I need to 
revisit the calculation with HITRAN2020. I thought I would use arts 2.4, a 
newer version. After successfully compiling arts2.4, and run the same .arts 
file I used before (a modification of TestAbs.arts), arts tells me that 
“abs_linesReadFromHitran” is not available any more. I went over the 
documentation, and found “READHITRAN”. I tried the following:

ReadHITRAN( abs_lines,
   “./h2o_61fbd3ce.par",
1.3e+14,
1.0e+15 )

arts reports:

"Run-time error in method: ReadHITRAN
Error parsing quantum number Sym”.

Is there any way we can read in HITRAN par file and create a lookup table? Do 
you suggest me going back to arts 2.3?

Thanks for your advice,

Pengwang

____
Pengwang Zhai (he/his/him)
Graduate Program Director, Atmospheric Physics,
Associate Professor, Physics Department, UMBC
pwz...@umbc.edu






___
arts_users.mi mailing list
arts_users.mi@lists.uni-hamburg.de
https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi


[arts-users] cannot compile the latest version of ARTS

2022-02-08 Thread Pengwang Zhai
Hello,

I had some issue with arts 2.4 in reading HITRAN par file. I was advised to use 
the latest version of arts. I did the following:

git clone https://github.com/atmtools/arts.git
cd arts
mkdir build
cd build
cmake -DENABLE_NETCDF=1 ..
make

The compilation stoped at 25%, with the following error message:

/arts/src/nc_io.h:37:10: fatal error: 'abs_species_tags.h' file not found
#include "abs_species_tags.h"
 ^~~~
1 error generated.
make[2]: *** [src/CMakeFiles/artscore.dir/auto_md.cc.o] Error 1
make[1]: *** [src/CMakeFiles/artscore.dir/all] Error 2
make: *** [all] Error 2

Any suggestion here on how to fix this? 

Note that my purpose is to create an absorption lookup table from HITRAN 2020. 
I need the NetCDF support to make the resultant absorption lookup table 
manageable.

Appreciate it very much,

Yours,

Pengwang






___
arts_users.mi mailing list
arts_users.mi@lists.uni-hamburg.de
https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi


Re: [arts-users] Cannot create abs_lookup with arts 2.5

2022-02-11 Thread Pengwang Zhai
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,
>>>>> 
>>>>> 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


Re: [arts-users] Cannot create abs_lookup with arts 2.5

2022-02-10 Thread Pengwang Zhai
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,
>> 
>> 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
> 
> 
> 

___
arts_users.mi mailing list
arts_users.mi@lists.uni-hamburg.de
https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi