Go to the U of M home page
School of Physics & Astronomy
School of Physics and Astronomy Wiki

User Tools


classes:2008:fall:astro4001.001:weblog

This is an old revision of the document!


Question and Answers - Web Log

How do I start?

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:

Change to directory 'src':

cd src

Open file dispatch.f:

emacs dispatch.f

You need to change this file so it looks as follows:

    PROGRAM dispatch
    USE demo1
    USE demo2
    USE ez_example
    CALL Do_Example
    END PROGRAM dispatch

Save and close it

Open ez_example.f:

emacs ez_example.f &

In there you find the lines to change metalicity and the mass:

For z = 0.02:

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

After you did the changes save and go to the directory 'make'. Type:

make

Now it should compile everything.

When it's done go to the directory 'run' and type:

./EZ

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.

Repeat for z = 0.001:

metals_dir = '../metals/z001' ! metallicity – use z0001, z0003, z001, z004, z01, z02, or z03

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 &

Open the plotting program gnuplot:

gnuplot

To plot, type:

    set terminal jpeg
    set output "stemperature.jpeg"
    set xrange [*:*]
    set yrange [*:*]
    set xlabel "Log age (yr)"
    set ylabel "Log surface temperature (K)"
    plot '0.02.txt' using 7:6 with lines, '0.001.txt' using 7:6 with lines
    reset

In this case, '0.02.txt' is the textfile for z = 0.02 and '0.001.txt' for z = 0.001. 'using 7:6' specifies which columns you want to plot. Column 7 corresponds to the age and column 6 to the surface temperature. For other plots change the 6 to whatever column you want to plot over time.

You will find the plot 'stemperature.jpeg' in the directory you are in.

classes/2008/fall/astro4001.001/weblog.1224002383.txt.gz · Last modified: 2008/10/14 11:39 by heger