Campuses:
This shows you the differences between two versions of the page.
| Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
| classes:2008:fall:astro4001.001:weblog [2008/10/14 11:39] – heger | classes:2008:fall:astro4001.001:weblog [2008/10/20 13:20] (current) – x500_star0288 | ||
|---|---|---|---|
| Line 30: | Line 30: | ||
| **Open ez_example.f: | **Open ez_example.f: | ||
| - | emacs ez_example.f & | + | |
| **In there you find the lines to change metalicity and the mass:** | **In there you find the lines to change metalicity and the mass:** | ||
| Line 36: | Line 36: | ||
| **For z = 0.02:** | **For z = 0.02:** | ||
| - | metals_dir = ' | + | |
| **Change the mass:** | **Change the mass:** | ||
| - | mass = 25D0 ! initial mass -- anything from 0.3 to 100 | + | |
| + | |||
| + | **Change the output frequency to make nicer plots, and suppress repeated header** | ||
| + | |||
| + | Add the lines | ||
| + | |||
| + | head_cnt = 100000 | ||
| + | summary_cnt = 1 | ||
| + | |||
| + | after the line | ||
| + | |||
| + | CALL Init_Do_One_Utils(mass) | ||
| **After you did the changes save and go to the directory ' | **After you did the changes save and go to the directory ' | ||
| - | make | + | |
| **Now it should compile everything.** | **Now it should compile everything.** | ||
| Line 50: | Line 61: | ||
| **When it's done go to the directory ' | **When it's done go to the directory ' | ||
| - | ./EZ | + | |
| **You see an output on your screen. To redirect this output to a textfile:** | **You see an output on your screen. To redirect this output to a textfile:** | ||
| - | ./EZ > textfile.txt | + | |
| **This textfile contains a table of values which you will use to plot the evolution of the star.** | **This textfile contains a table of values which you will use to plot the evolution of the star.** | ||
| + | **if you'd like to see the output while copying it into he text file, you can use the `tee` command by " | ||
| + | ./EZ | tee textfile.txt | ||
| **Repeat for z = 0.001:** | **Repeat for z = 0.001:** | ||
| - | metals_dir = ' | + | |
| Line 72: | Line 85: | ||
| **Open textfile.txt to see which values are in which column.** | **Open textfile.txt to see which values are in which column.** | ||
| - | emacs textfile.txt & | + | |
| **Open the plotting program gnuplot:** | **Open the plotting program gnuplot:** | ||
| - | gnuplot | + | |
| **To plot, type:** | **To plot, type:** | ||
| Line 94: | Line 107: | ||
| **You will find the plot ' | **You will find the plot ' | ||
| + | |||
| + | |||
| + | ====Class project 1, Problem 2==== | ||
| + | I think I've found where to turn off the CNO cycle: line 246 of src/ | ||
| + | |||
| + | ===== Does anyone know how to turn off the CN, ON, or pp-chain cycles? ===== | ||
| + | |||
| + | --Mark Zastrow | ||
| + | |||
| + | In ez_nuclear_data.f you can find how the different reactions are named. You can then edit them in ez_nuclear.f | ||
| + | |||
| + | So for example to turn off the He3 + He3 reaction, change R33 = RHB*RRT(2)/ | ||
| + | |||
| + | --Andrea | ||
| + | |||