TheiaSim - pagina 14
Opțiuni teme
afişa
Dacă se neglijează și efectul gravitațional al asteroizilor asupra planetelor, s-ar putea număra pur și simplu asteroizii printre corpurile testate.
Pentru a avea numerele corecte: asta înseamnă 4,4e17 kilometri pătrați.
Probabilitatea se potrivește apoi din nou .
De la mine, desigur, un da pentru asta în acest moment.
BTW: Astăzi am citit puțin în codul SCATR, pentru că aș dori să calculez comparațiile cu codul nostru BS. Dar trebuie să puteți citi Fortran destul de bine pentru a putea modifica fișierele de intrare în așa fel încât condițiile de pornire „noastre” (= 01.04.2011 sau ceva similar) să fie calculate pentru toate planetele. Numărul corpurilor de testare poate fi setat mai întâi la zero. Aș fi interesat să aflu ce rezultate calculează SCATR pentru un an, de exemplu, și cât de mari sunt abaterile de la rezultatele cu codul BS.
Ai dreptate și, de fapt, îți mulțumesc pentru acest „pas în A.” Arată că există un anumit interes pentru subiect. Proiectul actual a redevenit mult mai interesant pentru mine prin codul SCATR. Dacă doriți și codul Fortran, vă voi trimite cu plăcere. Este foarte bine făcut și foarte educativ.
Salutari
Știe cineva cum să înțeleagă o masă solară de 2.959e-4? Datorită celorlalte informații, această valoare ar trebui să fie proporțională cu masa SI. Unitățile geometrice nu sunt.
Mai întâi sunt în vacanță.
Bernhard, poate în mase solare?
Acest lucru este înțeles greșit. Masa soarelui este utilizată în program cu valoarea 2.959e-4.
a) Condițiile de pornire ale corpurilor de testare sunt controlate prin intermediul fișierului tp.in. Cu toate acestea, acest fișier conține o serie de parametri care nu sunt încă complet clari pentru mine.
b) Fișierul param.in controlează secvența programului și probabil oferă în prezent o simulare de peste 1 milion de ani. Am calculat câteva exemple cu termeni mai mari de 1 an și 100 de ani.
c) Pentru condițiile de pornire ale planetelor, coordonatele heliocentrice sunt date în UA și viteza în UA/zi în fișierul pl.in.
Nu pot să-mi dau seama.
În general, aș sugera utilizarea unităților SI în mod consecvent, chiar dacă devin voluminoase. Pentru că cu siguranță nu veți face erori de conversie.
În unele cazuri, trebuie să păstrați numerele într-un anumit interval pentru a evita erorile de calcul și de rotunjire. Înainte ca noi sau să nu înțelegem mai bine codul, îl puteți lăsa așa pentru moment. Unitatea pentru masă ar trebui însă clarificată. Pentru a nu-l deranja în mod inutil pe Nathan Kaib, s-ar putea contacta și un autor SWIFT.
BTW: În prezent, mă gândesc și la posibilitățile de a traduce întregul cod în C/C ++. Instrumentul Linux f2c eșuează pe câteva pasaje, dar ar putea reduce tastarea. O astfel de traducere ar fi cu siguranță utilă.
Salutari