Hi Michie,

It looks like there are just no lines available from HITRAN (which our catalog 
is based on) in that region (I attached a zoomed in plot in case others want to 
chime in). Can you link to a plot that shows the feature you're expecting?

Cheers,
Oliver

> On 16. Jun 2025, at 19:30, De Vera, Michie Vianca Dichaves via arts_users.mi 
> <[email protected]> wrote:
> 
> Good day!
>  I am new to ARTS and have been using the newest version on Python. I’m 
> trying to plot the absorption cross section for ozone as a function of 
> wavelength. However, I noticed that I’m only getting some of the features 
> right.  As seen in the image attached, for example, there should be a dip in 
> the values around 6.3 microns based on textbook images, but I don’t see that 
> here. Other wavelengths should also have dips that are not seen in the image. 
> I have tried to increase the spectral resolution, but I still do not get 
> these features.
>  Below is part of the Python code I’m using to get the absorption coefficient 
> values. I tried to follow one of the example files, but I would appreciate 
> any help to see if I’m missing something.
>  ws = pyarts.workspace.Workspace()
> ws.abs_speciesSet(species=[‘O3’])
> ws.abs_lines_per_speciesReadSpeciesSplitCatalog(basename="/dir/arts/arts-cat-data-2.6.14/lines/")
> ws.lbl_checkedCalc()
> ws.propmat_clearsky_agendaAuto()
> ws.jacobian_quantities = []
> ws.select_abs_species = [] 
> ws.f_grid = frq_grid
> ws.rtp_mag = [] 
> ws.rtp_los = [] 
> ws.rtp_pressure = 100681
> ws.rtp_temperature = 285.77 
> ws.rtp_nlte = pyarts.arts.EnergyLevelMap() 
> ws.rtp_vmr = [0.2676793e-07] 
> ws.stokes_dim = 1 
> ws.AgendaExecute(a=ws.propmat_clearsky_agenda)
>   Thank you so much!
>  Sincerely,
> Michie<arts_O3.png>




Reply via email to