[mesa-users] MESA code Questions!!
Farag, Ebraheem Khaled
farag.11 at lanl.gov
Tue Jul 7 11:22:53 EDT 2015
Hello,
I am trying to use the initial abundance from AGS05 and the code is returning this error:
------------------------------------------------------------------------------------------------------------
set_uniform_initial_composition
initial_h1 -1.0000000000000000D+00
initial_h2 -1.0000000000000000D+00
initial_he3 -1.0000000000000000D+00
initial_he4 -1.0000000000000000D+00
metals AGS05
when setting composition, need to provide values for h1, h2, he3, and he4
H1= -1.0000000000000000
H2= -1.0000000000000000
He3= -1.0000000000000000
He4= -1.0000000000000000
Z = 1.0000000000000000
set_uniform_initial_composition ierr -1
do_star_job_controls_after ierr -1
before_evolve_loop ierr -1
----------------------------------------------------------------------------------------------------
I tried to follow the star_job_defaults but it keeps returning the same error
This is all I have written:
---------------------------------------------------------------------------------------------------
! uniform initial composition of AGS05 abundances
set_uniform_initial_composition = .true.
initial_h1 = -1
!initial_h2 = -1
initial_he3 = -1
initial_he4 = -1
initial_zfracs = 5
------------------------------------------
thanks,
Ebraheem
________________________________
From: Robert FARMER [rjfarmer at asu.edu]
Sent: Tuesday, June 30, 2015 11:32 AM
To: Farag, Ebraheem Khaled; mesa-users at lists.sourceforge.net users
Subject: Re: [mesa-users] MESA code Questions!!
To keep the messages going to the users, when replying do a reply all and make sure mesa-users at lists.sourceforge.net<mailto:mesa-users at lists.sourceforge.net> is listed as a recipient.
Well your net will need to contain all the isotopes you think are important, if you add
show_net_species_info = .true.
To your star_job inlist then when mesa starts it will list the isotopes in the current net, then you can run models with different
nets and see which is most appropriate for your problem.
Rob
On Tue, Jun 30, 2015 at 8:23 AM, Farag, Ebraheem Khaled <farag.11 at lanl.gov<mailto:farag.11 at lanl.gov>> wrote:
Hello,
Thank you so much for your help! What net reccomendations would you have for me if I told you I was using new seperate opacity data files for elements z1-30. Will this affect the abundances at all in the evolution. I'm running this code specifically to see how these new opacity files affect the abundances I plot. Will these opacities affect my sounds speed plot?(I hope they do)- and do the abundances in the code depend at all on my input opacities?(I also hope they do). If these aren't the cases, how can I go about making them dependant on my new data.
p.s. with regards to keeping this on the mesa-users list, I'm not sure how to do that exactly except by responding to these emails. I tried emailing the other email on the website to get an account to no avail.
thanks for your help again!
-Ebraheem
________________________________
From: Robert FARMER [rjfarmer at asu.edu<mailto:rjfarmer at asu.edu>]
Sent: Monday, June 29, 2015 12:13 PM
To: Farag, Ebraheem Khaled; mesa-users at lists.sourceforge.net<mailto:mesa-users at lists.sourceforge.net> users
Subject: Re: [mesa-users] MESA code Questions!!
Please keep replies on the mesa-users list so if any else has this problem they can find the solution.
The nuclear reaction network ('net') tells mesa which isotopes to evolve, so by setting your job off to use the cno_extras_plus_fe56.net<http://cno_extras_plus_fe56.net> you will evolve the isotopes in this nets which will be fe56 plus the contents of the cno_extras.net<http://cno_extras.net> file (which might load another file..).
So by following this http://mesa.sourceforge.net/star_job_defaults.html#nuclear_reactions to set the net. Then in your star_job section
change_net = .true.
new_net_name = 'cno_extras_plus_fe56.net<http://cno_extras_plus_fe56.net>'
So now when you run mesa you should be able to add
Abundance_which_isos_to_show(1) = 'fe56'
to your abundance plot to show fe56
All of mesa options are in $MESA_DIR/defaults and on the web at http://mesa.sourceforge.net/ for the star_job, control and pgstar inlists.
Rob
On Mon, Jun 29, 2015 at 7:02 AM, Farag, Ebraheem Khaled <farag.11 at lanl.gov<mailto:farag.11 at lanl.gov>> wrote:
Hello,
I opened a file in the net directory called cno_extras_plus_fe56.net<http://cno_extras_plus_fe56.net>. I figure this gives info on how to add iron abundances into my plots but the information contained in the file just says :
include 'cno_extras.net<http://cno_extras.net>'
add_iso(fe56)
How do I manipulate this or put it into my inlist to allow me to plot iron abundances. I'm also wondering which inlist I should put it into?
thanks,
-Ebraheem
________________________________
From: Robert FARMER [rjfarmer at asu.edu<mailto:rjfarmer at asu.edu>]
Sent: Friday, June 26, 2015 3:34 PM
To: Farag, Ebraheem Khaled
Cc: mesa-users at lists.sourceforge.net<mailto:mesa-users at lists.sourceforge.net>
Subject: Re: [mesa-users] MESA code Questions!!
Hi
1) I suggest looking at http://mesa.sourceforge.net/pgstar_defaults.html#Abundance_window where you can see that you can set the isotope you want plotted as
Abundance_which_isos_to_show(1) = 'h1'
But that requires the isotope to be in your current nuclear network you are using, as MESA only evolves the isotopes you tell it too. Thus if there is no iron or silicon in your network then you wont be able to plot it. Without seeing your inlist i can't say what net your using but you can look in $MESA_DIR/data/net_data/net at the available ones. Then you could pick one with iron and silicon in which should then allow you to plot those abundances.
Though you'll have to work out which net is suitable for you problem, the bigger the net the slower mesa will run.
2) You can plot sound speed as a function of mass/radius etc look at it
http://mesa.sourceforge.net/pgstar_defaults.html#Profile_Panels
then by also plotting the isotope at the same time (using the _other_y plot settings) you can see the sound speed as a function of abundance.
3) Maybe there are, at this point your into doing science so you'll need to experiment to see if you think those parameters are suitable for the Sun. Don't just blindly trust the defaults.
Rob
On Fri, Jun 26, 2015 at 11:58 AM, Farag, Ebraheem Khaled <farag.11 at lanl.gov<mailto:farag.11 at lanl.gov>> wrote:
Hello,
I am a student work with the MESA code at Los Alamos National Lab, and i had a few questions needing answered.
1-Is there a method by which i can choose any element z1-30 to plot abundances for? It seems you can only use the 8 basic isotopes in the program, but I'm trying to plot abundances versus radaii for elemnts liek Iron and silicon too.. but if ISO those elements in my pgstar, nothing happens. where as if i remove H1 or He3 for example...then they either appear or disappear.
2- Is there a way to plot sound speed profiles for each element?
3- are there any specific parameters not included in the tutorial test_suite file that might be necessary to reference to ensure my models are closely replicating out sun?
Thanks,
Ebraheem Farag
------------------------------------------------------------------------------
Monitor 25 network devices or servers for free with OpManager!
OpManager is web-based network management software that monitors
network devices and physical & virtual servers, alerts via email & sms
for fault. Monitor 25 devices for free with no restriction. Download now
http://ad.doubleclick.net/ddm/clk/292181274;119417398;o
_______________________________________________
mesa-users mailing list
mesa-users at lists.sourceforge.net<mailto:mesa-users at lists.sourceforge.net>
https://lists.sourceforge.net/lists/listinfo/mesa-users
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20150707/dde7a6b2/attachment.html>
More information about the Mesa-users
mailing list