Renish and Freddy,

It is likely very long time since someone else used opt_prop_sptFromData!

Freddy: This method is most closely related to DOIT. I noticed that it is opt_prop_sptFromMonoData that is actually used by DOIT (for the version in main at least). Any comments about these methods? Should they still work in a stand-alone fashion?

Renish; That said, you should not get assert, but a proper error message. If you could send a complete example that would help to sort this out.

Bye,

Patrick



On 2022-02-18 06:03, Thomas,Renish wrote:
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
_______________________________________________
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