Dear ARTS Community,

I am experiencing some issues with extracting optical properties of scattering 
particles such as 
"ext_mat_spt<https://atmtools.github.io/arts-docs-master/docserver/variables/ext_mat_spt.html>"
 and 
"abs_vec_spt<https://atmtools.github.io/arts-docs-master/docserver/variables/abs_vec_spt.html>"
 using the 
"opt_prop_sptFromData<https://atmtools.github.io/arts-docs-master/docserver/methods/opt_prop_sptFromData.html>"
 workspace method.

The Runtime Issue:

=> The python kernel crashes {"Aborted (core dumped)"} during runtime when the 
"opt_prop_sptFromData<https://atmtools.github.io/arts-docs-master/docserver/methods/opt_prop_sptFromData.html>"
 method is executed.

=> I am encountering two types of runtime errors based on za_grid and za_index.

=> During traceback, the runtime errors seem to be caused by checks on 
"za_index" and "aa_index".

=> When za_grid or aa_grid is empty===runtime error=>
python: /arts/src/matpackI.h:509: Numeric ConstVectorView::operator[](Index) 
const: Assertion `n < mrange.mextent' failed.

=> When za_grid or aa_grid is not empty===runtime error=>
python: /arts/src/array.h:215: base& Array<base>::operator[](Index) [with base 
= PropagationMatrix; Index = long int]: Assertion `n < nelem()' failed.

The Scattering Data:

=> The scattering data is extracted from the "totally_random" scattering 
particles "LiquidSphere_TotRand" from the ARTS single-scattering properties 
database.

=> ws.scat_data.za_grid = array([0.0,0.2,.....,90,....,179.8,180]).  Zenith 
angle grid from scat_data.

=> ws.scat_data.aa_grid = array([ ])  Empty azimuth angle grid from scat_data.

My Questions:

=> What could be the reason for the runtime errors?

=> Should I assign the same "za_grid" from "scat_data" when using the 
"opt_prop_sptFromData<https://atmtools.github.io/arts-docs-master/docserver/methods/opt_prop_sptFromData.html>"
 method? or is the "za_grid" selection dependent on something else?

=> Should I leave the aa_grid empty as mentioned in the ARTS user guide for 
totally random particles when using 
"opt_prop_sptFromData<https://atmtools.github.io/arts-docs-master/docserver/methods/opt_prop_sptFromData.html>"?.
 How would I do this without errors?

=> How do I select za_index intuitively? Should I loop through all elements or 
do all the index values give game results for "totally_random" scattering 
particles. An example or illustration would help. The lab and scattering 
coordinate transformation is making the selection of za_grid and za_index vague 
to me.

Misc Info:

=> ARTS Version: arts-2.5.0
=>  Pyarts

Thanks very much!
Renish Thomas


_______________________________________________
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