My humble codes for the simulation of the impact of satellite constellations on astronomical optical observations is available here for download. Feel free to download, use and modify it at your will. Simply, please, acknowledge the source in case of preparing any publication based on my work.
The programs are stand-alone and they are written in fortran. They are not an example of good programming practice, but they work and they may be useful to somebody.
For the moment I can only offer the code and some indications on how to use each piece. In the near future I may add more detailed documentation, getting into the details of what and how the programs do. For now, this is what you have:
A) Montecarlo simulations: satellite counts above a specified elevation, and generation of simulated orbital elements in several formats, including TLE
This is your software if you want to count how many satellites are there above a specific elevation (for instance, above the horizon), and which of them are illuminated by the Sun. Also, this program allows you to generate simulated orbital elements in TLE format, or in the format required by ESA’s program PROOF.
| handbook_A.pdf | The handbook for section A |
| constellation4.f | The source code in fortran, for Linux systems |
| starlink.dat | A sample of constellation profile file |
| averator | The script for iterative execution |
| doaverage.f | One of the two mini-programs required by averator |
| doaversim.f | The second of the two mini-programs required by averator |
B) Montecarlo simulations: satellite trails in specific shots
You may be interested in knowing how many trails you can expect if you point your telescope to a specific part of the sky and you perform an integration, given a constellation profile.
| handbook_B.pdf | The handbook for section B |
| constellation9b.f | The source code in fortran, for Linux systems |
| starlink.dat | A sample of constellation profile file |
C) Analytic simulations: satellite trails in specific shots
A sophisticated turn to the same problem solved in section B. Use the analytic probability density function to save computer time to get accurate results in a fraction of a second.
| handbook_C.pdf | The handbook for section C |
| constellation10.f | The source code in fortran, should work for all operating systems |
| starlink.dat | A sample of constellation profile file |
