TheiaSim - strana 14

Možnosti témy
displej

Ak niekto zanedbá aj gravitačný účinok asteroidov na planéty, mohol by jednoducho spočítať asteroidy medzi testovacími telesami.

theiasim

Aby boli čísla správne: to je 4,4 e17 štvorcových kilometrov.

Pravdepodobnosť potom opäť sedí .

Pre mňa samozrejme v tejto chvíli áno.

BTW: Dnes som sa trochu načítal do kódu SCATR, pretože by som chcel vypočítať porovnania s našim kódom BS. Musíte však byť schopní dobre prečítať Fortran, aby ste mohli upravovať vstupné súbory tak, aby sa pre všetky planéty počítali „naše“ počiatočné podmienky (= 01.04.2011 alebo niečo podobné). Počet skúšobných telies možno najskôr nastaviť na nulu. Zaujímalo by ma, ktoré výsledky vypočítava SCATR napríklad za rok a aké veľké sú odchýlky od výsledkov s BS kódom.

Máte pravdu a v skutočnosti vám ďakujem za tento „krok do A.“ Ukazuje to, že o túto tému je určitý záujem. Aktuálny projekt sa pre mňa stal opäť oveľa zaujímavejším prostredníctvom kódu SCATR. Ak chcete aj Fortranský kód, rád vám ho pošlem. Je to naozaj dobre urobené a veľmi poučné.
Zdravím vás

Vie niekto, ako pochopiť slnečnú hmotu 2,959e-4? Kvôli ďalším informáciám by táto hodnota mala byť úmerná hmotnosti SI. Geometrické jednotky nie sú.

Najskôr som na dovolenke.
Bernhard, možno v slnečných hmotách?

To je nepochopené. Hmota slnka sa v programe používa s hodnotou 2,959e-4.

a) Štartovacie podmienky skúšobných telies sa riadia prostredníctvom súboru tp.in. Tento súbor však obsahuje množstvo parametrov, ktoré mi zatiaľ nie sú úplne jasné.

b) Súbor param.in riadi tok programu a v súčasnosti pravdepodobne poskytuje simuláciu viac ako 1 milión rokov. Vypočítal som niekoľko príkladov s termínmi viac ako 1 rok a 100 rokov.

c) Pre počiatočné podmienky planét sú heliocentrické súradnice uvedené v AU a rýchlosti v AU/deň v súbore pl.in.

Nemôžem na to prísť.

Všeobecne by som navrhol dôsledne používať jednotky SI, aj keď sa stanú objemnými. Pretože určite neurobíte žiadne chyby pri prevode.

V niektorých prípadoch musíte udržiavať čísla v určitom rozmedzí, aby ste sa vyhli chybám pri výpočte a zaokrúhľovaní. Skôr ako ten kód pochopíme lepšie, alebo my to zatiaľ necháme tak. Jednotka hmotnosti by však mala byť objasnená. Aby ste zbytočne neobťažovali Nathana Kaiba, dalo by sa tiež kontaktovať autora SWIFTu.

BTW: Momentálne uvažujem aj o možnostiach preloženia celého kódu do C/C ++. Linuxový nástroj f2c na niekoľkých pasážach zlyhá, ale mohlo by to obmedziť písanie. Takýto preklad by bol určite užitočný.
Zdravím vás