Do wyboru:

Modelowanie numeryczne krzywych blasku

Cel ćwiczenia: Modelowania numeryczne krzywej blasku ciasnego układu podwójnego z misji Kepler za pomocą kodu bazującego na kodzie Wilsona-Devinney'a.

Dane: Można pracować na obiekcie wybranym do ćwiczenia z efektem O'Connella. Do modelowania należy sobie wybrać jedną krzywą blaksu z zestawów z poprzednich ćwiczeń. Najlepiej będzie wybrać krzywą blasku z silnym efektem O'Connella.

Pliki z zajęć znajdują się w katalogu kopiuj_mnie oraz są dostępne jako paczka tar.

UWAGA! Aby uproscic proces modelowania zmodyfikowalem plik bat2new. Jego nowa wersja nazywa sie failsafe.sh i znajduje sie w paczkach pod powyzszymi linkami oraz tutaj. Nalezy skopiowac ten plik do swojego katalogu z plikami do modelowania i uzywac go zamiast bat2new.

Uzywajac pliku failsafe.sh mozna uruchomic modelowanie na wlasna reke. Nalezy mnie koniecznie poinformowac o fakcie rozpoczecia modelowania. Prosze wyslac mi mail z informacja oraz numerem PID (wyswietla sie po uruchomieniu programu). Kazdy niezgloszony proces moze byc przeze mnie wstrzymany ze wzgledow bezpieczenstwa.

Plan ramowy zadania:

1. Wybrać krzywą blaksu z efektem O'Connella
2. Przygotować pliki wsadow do programu sp.prog
3. Rozpocząć proces modelowania
4. Przygotować plik wsadowy do programu lc-generator
5. Narysować krzywą blasku
6. Wykonać obraz za pomocą rysownika
7. Całość opisać w raporcie i dołączyć obrazy

Termin oddawania raportów: 24 grudnia 2015 (sześć tygodni po wprowadzeniu ćwiczenia na zajęciach (+1 tydzień przez awarię!))
Termin oddawania spóźnionych raportów: do sesji, 28 stycznia 2016. Na zaliczenie trzeba oddać raport powiększony o dodatkowy rozdział - szczegóły niżej.

---

Raport musi się składać z następujących części:

1. Wstęp teoretyczny
• Ogólne przedstawienie ćwiczenia
• Krótki opis geometrii Roche'a, fillout factor, potencjałów i innych waznych terminów używanych w ćwiczeniu
• Krótki opis ogólny obiektu
• Krótki opis działania programu (m.in. metoda Monte Carlo, uzywane parametry itd.)
• W przypadku raportu wysłanego po 24 grudnia: rozdział opisujący, czym jest pociemnienie brzegowe w gwiazdach; jakie są (co najmniej) trzy podstawowe "prawa" analityczne opisujące rozkład jasności na widomej tarczy gwiazdy.
• W przypadku raportu wysłanego po 24 grudnia: rozdział opisujący, czym jest pociemnienie (pojaśnienie?) grawitacyjne; koniecznie napisac o "prawie von Zeipela" i wspomnieć, jaki wykładnik teoretycznie powinien opisywać rozkład temperatury efektywnej w przypadku powierzchni konwektywnych.
2. Przygotowanie danych
• Napisać, jakie sie wybrało dane: obiekt + epoka
• Krótko opisać proces przygotowania plików wsadowych (do sp.prog i do lc-generator)
3. Modelowanie numeryczne
• Zwięźle przedstawić sposób działania programu sp.prog
• Zamieścić tabelę z zakresami modelowania oraz z wynikami i niepewnościami
4. Podsumowanie
• Przedstawić rysunek z danymi obserwacyjnymi i naniesiona krzywą symulowaną
• Przedstawić obraz z rysownika
• Skomentować wyniki: czy dopasowanie jest dobre, co mogłoby być dopadowane lepiej, o ile cokolwiek być może, co można wysnuć z modelowania itd.

---

Podpowiedzi:

• Modelowanie odbywa sie programem sp.prog.
• Program sp.prog jest wolany przez skrypt failsafe.sh.
• To znaczy, ze moze trzeba bedzie podac sciezke do tego programu!
• ...program dolaczam w paczce z plikami.
• Wyniki modelowania sprawdza sie programem disp2err.
• Uruchamianie procesu w tle: nohup ./failsafe.sh > log&
• Podgladanie pliku "na zywo": tail -f log&
• Wynik rysowania krzywej bedzie w pliku plt.dat.
• Rysowanie ukladu odbywa sie skryptem rysownik.sh.

Struktura pliku wsadowego do programu lc-generator