TheiaSim - страница 14
Опции за темата
дисплей
Ако някой пренебрегне и гравитационния ефект на астероидите върху планетите, може просто да преброи астероидите сред тестовите тела.
За да са правилни числата: Това е 4.4e17 квадратни километра.
След това вероятността отново се вписва .
За мен, разбира се, да на това в момента.
BTW: Прочетох малко SCATR кода днес, защото бих искал да изчисля сравнения с нашия BS код. Но трябва да можете да четете Fortran доста добре, за да можете да модифицирате входните файлове по такъв начин, че "нашите" начални условия (= 01.04.2011 или нещо подобно) да бъдат изчислени за всички планети. Броят на тестовите тела може първо да бъде зададен на нула. Би ме интересувало какви резултати SCATR изчислява например за една година и колко големи са отклоненията от резултатите с BS кода.
Прав си и всъщност ти благодаря за тази „стъпка към А.“ Това показва, че има определен интерес към темата. Настоящият проект отново стана много по-интересен за мен чрез кода на SCATR. Ако искате и кода на Fortran, ще се радвам да ви го изпратя. Това е наистина добре направено и много образователно.
Поздравления
Някой знае ли как да разберем слънчевата маса от 2.959e-4? Поради другата информация тази стойност трябва да бъде пропорционална на масата на SI. Геометричните единици не са.
Първо съм на почивка.
Бернхард, може би в слънчеви маси?
Това е неразбрано. Масата на слънцето се използва в програмата със стойност 2.959e-4.
а) Условията за стартиране на тестовите тела се контролират чрез файла tp.in. Този файл обаче съдържа редица параметри, които все още не са ми напълно ясни.
б) Файлът param.in контролира програмната последователност и понастоящем вероятно осигурява симулация за 1 милион години. Изчислил съм няколко примера със срокове над 1 година и 100 години.
в) За началните условия на планетите хелиоцентричните координати са дадени в AU, а скоростите в AU/ден във файла pl.in.
Не мога да разбера това.
Като цяло бих препоръчал последователно да се използват единици SI, дори и да станат обемисти. Защото определено няма да допускате грешки при преобразуването.
В някои случаи трябва да поддържате числа в определен диапазон, за да избегнете грешки при изчисление и закръгляване. Преди ние или аз да не разбираме кода по-добре, засега можете да го оставите така. Единицата за масата обаче трябва да бъде изяснена. За да не притеснявате Нейтън Кайб ненужно, евентуално бихте могли да се свържете и с автор на SWIFT.
BTW: В момента също обмислям възможностите да преведа целия код в C/C ++. Инструментът за Linux f2c се проваля в няколко пасажа, но може да намали въвеждането. Такъв превод определено би бил полезен.
Поздравления