Differences

This shows you the differences between two versions of the page.

--- classes:2008:fall:astro4001.001:weblog [2008/10/13 22:33] – x500_mehne002
+++ classes:2008:fall:astro4001.001:weblog [2008/10/20 13:20] (current) – x500_star0288
@@ Line 1: / Line 1: @@
-====== Question and Answers - Web Log ======
+======= Question and Answers - Web Log =======
-===== Title of First Question =====
-====== How do I start? ======
-----
+===== How do I start? =====
-I assume everyone has the code installed. I also assume that everyone knows how to get to the directory EZ_3_6_3.
-===== 1 problem: =====
+Let's assume everyone has the code installed. I also assume that everyone knows how to get to the directory EZ_3_6_3.
+===== 1st Problem: =====
@@ Line 30: / Line 30: @@
 **Open ez_example.f:**
-emacs ez_example.f &
+      emacs ez_example.f &
 **In there you find the lines to change metalicity and the mass:**
@@ Line 36: / Line 36: @@
 **For z = 0.02:**
-metals_dir = '../metals/z02' ! metallicity -- use z0001, z0003, z001, z004, z01, z02, or z03
+      metals_dir = '../metals/z02' ! metallicity -- use z0001, z0003, z001, z004, z01, z02, or z03
 **Change the mass:**
-mass = 25D0  ! initial mass -- anything from 0.3 to 100
+      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 'make'. Type:**
-make
+      make
 **Now it should compile everything.**
@@ Line 50: / Line 61: @@
 **When it's done go to the directory 'run' and type:**
-./EZ
+      ./EZ
 **You see an output on your screen. To redirect this output to a textfile:**
-./EZ > textfile.txt
+      ./EZ > textfile.txt
 **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 "pipelining" the output as input to the `tee` command using `|`**
+      ./EZ | tee textfile.txt
 **Repeat for z = 0.001:**
-metals_dir = '../metals/z001' ! metallicity -- use z0001, z0003, z001, z004, z01, z02, or z03
+      metals_dir = '../metals/z001' ! metallicity -- use z0001, z0003, z001, z004, z01, z02, or z03
-===== Title of Second Question =====
-====== How do I plot? ======
+===== How do I plot? =====
-----
 **Stay in the directory where your textfile.txt is.**
 **Open textfile.txt to see which values are in which column.**
-emacs textfile.txt &
+      emacs textfile.txt &
 **Open the plotting program gnuplot:**
-gnuplot
+      gnuplot
 **To plot, type:**
-set terminal jpeg
+      set terminal jpeg
-set output "stemperature.jpeg"
+      set output "stemperature.jpeg"
-set xrange [*:*]
+      set xrange [*:*]
-set yrange [*:*]
+      set yrange [*:*]
-set xlabel "Log age (yr)"
+      set xlabel "Log age (yr)"
-set ylabel "Log surface temperature (K)"
+      set ylabel "Log surface temperature (K)"
-plot '0.02.txt' using 7:6 with lines, '0.001.txt' using 7:6 with lines
+      plot '0.02.txt' using 7:6 with lines, '0.001.txt' using 7:6 with lines
-reset
+      reset
 **In this case, '0.02.txt' is the textfile for z = 0.02 and '0.001.txt' for z = 0.001.
@@ Line 92: / Line 107: @@
 **You will find the plot 'stemperature.jpeg' in the directory you are in.**
+====Class project 1, Problem 2====
+I think I've found where to turn off the CNO cycle: line 246 of src/star_controls.f.
+===== 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)/2D0 into R33 = 0 !RHB*RRT(2)/2D0
+--Andrea

User Tools

Differences

Page Tools