Comments (3)
The underlying fortran code can become unstable when presented with very large particles. It's a known issue of the code. There has been some discussion between various people about fixing it, but not sure anyone has had the time to take it on.
from pytmatrix.
You could increase value of the attribute ndgs in the Scatterer class. By default ndgs=2, try increasing it to 100 (play around with this value). This attribute controls convergence of the T-matrix computations.
from pytmatrix.
I'll comment now on this long-standing issue since there hasn't been any progress: Unfortunately this is a problem with the core Fortran code that PyTMatrix uses (a code originally from NASA GISS).
The best solution to this problem would be to rewrite the T-Matrix core in Python/NumPy. This would also eliminate the compiling problems that many people have had with PyTMatrix. However, this is a serious effort and I don't really expect to have the time to do it in the foreseeable future, especially as my own career is gravitating away from scattering. I'll leave this issue open just in case some brave soul wants to tackle it.
from pytmatrix.
Related Issues (20)
- vertical incidence HOT 1
- Conda packages seem to be out of date HOT 2
- Calculate Nw HOT 1
- License mismatch HOT 2
- How to get VSF with an incident radiant of none polarization
- Error in drop shape relationship function HOT 1
- Is there something wrong with this code? HOT 4
- The unit of the scattering amplitude the code defines is 'km' ?
- How to calculate snow particles'reflectivity?
- Can the psd funtion be added in the psd.py?
- How to get RCS of single particle?
- Roundoff error during test_adaptive_orient HOT 2
- How large size parameters can this code calculate? HOT 2
- What is the limitation of this code? HOT 1
- Something wrong with the computation of Av? HOT 13
- When calculating the reflectance factor in the WRF scene by PSD
- numpy issue with python >= 3.12 HOT 14
- scatterer.psd_integrator.geometries = (geom_back, geom_forw) HOT 1
- I cant get the gfortran installed in my Jypyter Lab
- The example is working but my code with different parameters isn't working
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from pytmatrix.