# this allows you to put comments here #You have to indent commands, with the name of the macro at the left. #all is the name of this macro all data starmodl_05a.dat #change the file name to what you used for 0.5 solar mass model lines 24 449 #don't want to plot the top lines read x 1 #the independent variable is in column 1 read y 5 #the dependent variable for pressure is in column 5 set z = x/6.24E+10 #this changes the radius to the fraction of the total radius set w = lg(y) #puts y into logarithmic form limits 0 1.05 0 18 #hard set the limits box #draw the box with the set limits ltype 0 #set the line type connect z w #connect the dots # data starmodl_10a.dat #new file -- model for 1 solar mass lines 24 387 #there are a different number of lines in this model read x 1 #the independent variable; same limits as above read y 5 #the dependent variable set a = x/2.36E+11 #this star has a different radius set b = lg(y) #pressure in log form ltype 1 #change the line type connect a b #connect xlabel r/R_{star} #label the x axis ylabel log10(pressure) #label the y axis # data starmodl_3a.dat #new file with different mass lines 24 425 #different number of total lines read x read y 5 set c = x/1.11E+11 #different radius -- but still get r/R set d = lg(y) ltype 3 connect c d relocate 0.2 10 #this is like xyouts in IDL, but uses set limits for coordinates ltype 1 #set line type draw 0.3 10 #draws a line from 0.2, 10 to 0.3, 10 label 10 M_{sun}i #and puts in a label right where the "cursor" left off relocate 0.2 8 ltype 0 draw 0.3 8i #ditto label 0.5 M_{sun} relocate 0.2 9 ltype 3 draw 0.3 9 label 3 M_{sun} relocate 0.2 5 label X = 0.7, Y = 0.292, Z = 0.008 # #Add additional data from files with models for different masses, just adjusting the radius and the line type. #FOR HELP IN SUPERMONGO, use help topic; example, help ltype