Modelar la formación de planetas es un problema extremadamente difícil y se realiza por etapas. Para las primeras etapas necesitas hidrocódigos. Para las últimas etapas, necesitas códigos de N cuerpos. Las cosas entre las personas todavía están tratando de averiguar qué hacer. Estoy más familiarizado con los códigos N-body utilizados para la acumulación de etapa tardía. Los principales actores son Swift y Symba de Hal Levison y Mercury de John Chambers. Esos están escritos en Fortran. La mayoría de los códigos hyrdo que conozco también están escritos en Fortran o C. Python es un lenguaje horrible para el trabajo de simulación numérica, especialmente para este tipo de problemas que realmente requieren muchos cálculos rápidos de números.
En cuanto a la gravedad dependiente del tiempo, puede modificar cualquiera de los anteriores para incluirlo siempre que se sienta cómodo editando el código Fortran.
- ¿Cuál fue el estado inicial de la tierra?
- ¿Cuánto tiempo duraría un viaje a Wolf 1061c?
- ¿Están los vientos solares quitando la atmósfera de la Tierra y, de ser así, a qué velocidad?
- ¿A qué temperatura debe enfriarse el Sol para que muera la vida humana en la Tierra?
- ¿Cuán sensible es la trayectoria de la órbita de la Tierra a su masa, es decir, cuánto aumento de la masa de la Tierra puede causar un cambio en su trayectoria orbital?