Activity for Bradley S. Meyer
-
Bradley S. Meyer posted a comment on discussion Getting Help with Nucnet Tools
Sorry for this very late reply. The make data and make weak_data should now work. I will be working on the other errors. Best wishes.
-
Thanks for your email. Sourceforge has upgraded its security requirements. I have updated the links, so you should be able to cd nucnet-tools-code svn up cd examples/network make data Let me know if this doesn't work. Also, if you find other broken links, please let me know. Thanks again, and best wishes.
-
Thanks for your email. Sourceforge has upgraded it security requirements. I have updated the links, so you should be able to cd nucnet-tools-code svn up cd examples/network make data Let me know if this doesn't work. Also, if you find other broken links, please let me know. Thanks again, and best wishes.
-
Use https.
-
Compile
-
Home
-
Home
-
Home
-
Create release 0.8 tag.
-
Update for SPARSKIT2_F95 v1.1.0.
-
Restore missing header lines.
-
Fix download links. Update release date.
-
Thanks for your post. munuekT is the chemical potential of the electron-type neutrinos divided by kT. The value -infinity means that they escape the system. When the neutrinos are trapped, their chemical potential can be different from -infinity, which affects weak rates (either neutrinos are captured by nuclei or the neutrino degeneracy hinders decays). In short, we keep track of this quantity to compute weak interaction rates. I hope this helps. Best wishes.
-
-
Add make.
-
Add flags.
-
Add.
-
Bump version.
-
Remove files. Update for f95.
-
Fix typo.
-
Update for f595.
-
I should have added that we have a Python code that checks rates (forward and reverse) over a temperature range: https://osf.io/6ygax/wiki/home/ This is also linked from webnucleo.org: https://webnucleo.readthedocs.io/en/latest/scripts.html
-
Thanks for your post. Historically the fit parameters included Tlowfit and Thighfit, which gave the range of temperature over which the fit holds. Typically, these were T9=0.01 to 10. libnucnet XML by default uses this range and the rate outside this range as either at Tlowfit for T9 < Tlowfit and at Thighfit for T9 > Thighfit (that is, no extrapolation). I'm not sure what the current status of the valid range of fit parameters for JINA rates is, but fits can certainly misbehave outside the valid...
-
I'm glad to hear it worked for you. I don't have any i-process posts yet. I'd like to do that soon, but for now, you could try the specify abundances recipe but use neutron densities in the range 10^12 to 10^15 per cc. Timescales on the order of hours to months would be interesting. Recent papers should give some guidance. I hope that helps. Best wishes.
-
Thanks for your post. In (n,g)-(g,n) equilibrium, you are assuming that neutron-capture and disintegration reaction rate are fast (especially faster than beta decay rates). The equilibrium establishes itself, then the nuclei in an isotope chain beta decay. Those beta decays move nuclei up in Z. You can compute the abundance of an isotope chain Y_Z = \sum_A Y(Z, A), where Y(Z, A) is the (n,g)-(g,n) equilibrium abundance of (Z, A). Then define p(Z, A) = Y(Z, A) / Y_Z. Then define lambda_beta_Z = \sum_A...
-
I'm very, very late to respond to this--sorry for that. The problem is that sourceforge is no longer allowing codes to access the schemas or stylesheets. I hope to finally get around to fixing this via XML Catalogs. You can turn off validation inside the code for now.
-
Curious. Could you tell me what compiler you are using? If gcc, you could type gcc -v and let me know the result.
-
Sorry for the trouble. With all the updates to Mac Os, I haven't had time to fix everything. There are blog posts that try to address it. You might also try setting the NNT_NO_OPENMP flag: export NNT_NO_OPENMP=1 and then compiling. If you get totally frustrated, your best bet is to run with Docker: https://webnucleo-cookbooks.readthedocs.io/en/latest/cookbooks/installation/Docker.html Let me know if you still have trouble. Best wishes.
-
To run the nucnet tools codes from Docker, do docker run -it webnucleo/single_zone Then, once the container starts, you are in the directory /my-projects/single_zone. To run the nucnet-tools-code examples, go to that directory cd ../nucnet-tools-code/examples and proceed.
-
Sorry for the trouble. With all the updates to Mac Os, I haven't had time to fix everything. There are blog posts that try to address it. You might also try setting the NNT_NO_OPENMP flag: $ export NNT_NO_OPENMP=1 and then compiling. If you get totally frustrated, your best bet is to run with Docker: https://webnucleo-cookbooks.readthedocs.io/en/latest/cookbooks/installation/Docker.html Let me know if you still have trouble. Best wishes.
-
Thanks for your post. You'll need to compile with make. The make file downloads the required libraries, etc. You might also be interested in https://webnucleo-cookbooks.readthedocs.io/en/latest/cookbooks/cookbooks.html Let me know if you still have trouble. Best wishes.
-
Thanks for your posts. Obviously there are several issues here. Let me first point you to https://bitbucket.org/mbradle/simple_sn2/src/master/ That code suite is more up-to-date and better maintained. Second, for making figures, I suggest looking into using wnutils: https://wnutils.readthedocs.io/en/latest/ Regarding the problem with large XML and decaying abundances, I have another code for that which I willl be working on. It will be updated on https://webnucleo-cookbooks.readthedocs.io/en/latest/...
-
Use phi instead of exp.
-
Return if no decay reactions.
-
Thanks for the last two posts. For (n,g)-(g,n) equilibrium, you can look at the following Jupyter notebooks, which are also linked from webnucleo.org. You might also be interested in this paper. The codes use webnucleo XML. We have notebooks to create and manipulate such data. It would be useful to have a code that computes the r process evolution from (n,g)-(g,n) equilibrium. That's an old-fashioned approach but instructive. We'll put it on our todo list.
-
Thanks for your email. I'm not sure exactly what you'd like to do. Perhaps you can communicate with me directly at mbradle at clemson dot edu. Best wishes.
-
Thanks for your note. You might also be interested in our new version of things at: https://webnucleo-cookbooks.readthedocs.io/en/latest/cookbooks/cookbooks.html Best wishes. Brad Meyer
-
Thanks for your post. Your compiler does not like the string comparison. You can try changing line 579 from: if( S_SOLVER == nnt::s_ARROW ) to if( strcmp(S_SOLVER, nnt::s_ARROW) == 0 ) and recompiling. Best wishes.
-
Thanks for your post. I'll look at this soon. Just so I know, what compiler are you using? In particular, if you type gcc -v what is the response? Thanks, and best wishes.
-
Thanks for your post. I think one of the tarballs got corrupted. You can try make clean_all which will clean everything up and then make run_single_zone You could instead also just type ./examples_make I hope that helps. Best wishes.
-
Thanks for your post. It is not a problem. The conditions for this calculation are for those like at the center of the Sun. The temperature and density are not high enough to create nuclei above z about equal to 4. Try, for example, changing t9_0 to 1.0 and rho_0 to 1.e3 in zone.xml and run the calculation again with "[z <= 10]" and "[a <= 280]". You'll see a difference. Best wishes.
-
Thanks for your post. I just addressed this in a tardy reply from the post above, but I repeat here for completeness. The issue is that SourceForge no longer allows direct parsing of stylesheets over the web. I've moved the examples.xsl to the local directory, so compilation should now work. Try cd nucnet-tools-code svn up cd build make -f Makefile.libnuceq all_libnuceq That should now work. Let me know if not. Best wishes.
-
Thanks for your post, and I'm very sorry for the tardy reply. The issue is that SourceForge no longer allows direct parsing of stylesheets over the web. I've moved the examples.xsl to the local directory, so compilation should now work. Try cd nucnet-tools-code svn up cd build make -f Makefile.libnuceq all_libnuceq That should now work. Let me know if not. Best wishes.
-
Move examples stylesheet to local directory.
-
Thanks for your posts, and I'm glad the problem is resolved. You might also check out the webnucleo cookbooks. Best wishes.
-
Thanks for your post. Does this help? I also see this post. If these don't work, I'd try further searching on the error message. I hope one of those suggestions works.. Best wishes.
-
-
Thanks for your post. The code will adjust the timestep according to the change in abundances, so you will typically not have exactly the same timesteps in the two different calculations. Also it may not be good to force the two calculations to have the same timesteps, unless they are small, because other abundances may be changing too fast. You can compare most easily by plotting the two calculations against each other. I suggest you use wnutils. You can use the Multi_Xml class to do this most easily....
-
Thanks for your post. The code will adjust the timestep according to the change in abundances, so you will typically not have exactly the same timesteps in the two different calculations. Also it may not be good to force the two calculations to have the same timesteps, unless they are small, because other abundances may be changing too fast. You can compare most easily by plotting the two calculations against each other. I suggest you use wnutils. You can use the Multi_Xml class to do this most easily....
-
Change / to in to avoid hdf5 conflict.
-
Thanks for your post. The default procedure is for the code to use detailed balance in which case the forward flow is balanced by the reverse flow in the case of nuclear statistical equilibrium, which is a standard approach. This technical report gives some details, including how to include corrections to the reverse factors. It is linked from the libnucnet home page. I hope that helps. Best wishes.
-
Update ignores.
-
Update for new release. Update to align with libxml2 API.
-
Sorry for the delayed response. As I say, I'm not an expert on plasma effects. The point is to compute the correction to the chemical potential of the nuclear species. Those can be accommodated by the code. If it is different from their Eq. (12), then the correction term would need to be different.
-
Add release 0.20.
-
Fix typo.
-
Update regression test. Fix typos and add output.
-
Add file and ignore.
-
Add matrix by matrix multiplication. Update documentation and release notes.
-
Sorry for the delayed response. I'm somehow just seeing it. I'm glad things are working. Best wises.
-
Thanks for your post. I'm not sure exactly what you are asking since the Coulomb corrections make the plasma non ideal. In any event, the corrections I used are those in the Bravo and Garcia-Senz paper, which has the details in case you are interested. Best wishes.
-
Fix typos.
-
Thanks for your post. I am now recommending that people use wnutils to update data. You can read the data easily and update. You can check that the data are getting updated. Let me know if you still have trouble. Best wishes.
-
Update.
-
Thanks for your post. The value -1.1209e-7 for dY(15N)/dt comes from -1.1204e-7 - 4.8204e-11 = -1.1209e-7. If there is a non-zero value for the gamma reverse reaction , you would add the rate * Y(16O) to dY(15N)/dt and dY(1H)/dt and subtract from the dY(16O)/dt. I hope that helps. Best wishes.
-
Thanks for your post. Yes, you may certainly do this. You will need to convert your data to XML. This post shows how. We now have Jupyter notebooks that make this even easier. See webnucleo.org, and, in particular, the create_reaction_xml notebook at webnucleo_xml. You can merge your data into the network file, or you can use XInclude. Once you've created the XML with your data, try plotting it with the plot_reaction_xml project. Let me know if you have any trouble. Best wishes.
-
I believe I have communicated with you previously, but I will answer this post for completeness. You might want to try using wnutils. See the Creating and Writing XML Data tutorial. We also have a Jupyter notebook that can be of use. Best wishes.
-
Thanks for your post, and apologies for the late reply. It seems you are not able to link to the sourceforge servers, perhaps because you are behind a firewall. You can download the wn_matrix, libnucnet, libnuceq, and libstatmech tarbars directly to nucnet-tools-code/vendor and then compile. for example. svn checkout svn://svn.code.sf.net/p/nucnet-tools/code/trunk nucnet-tools-code mkdir vendor Now download the various tarballs to the vendor directory. Do this from their websites: wn_matrix, libnucnet,...
-
Sorry for the very late reply. It looks like gsl is not installed. Could you try to install it? You should get information returned when you type: gsl-config Best wishes.
-
My apologies for the very late reply. I recommend using wnutils. Please see the tutorials. Best wishes.
-
Sorry for the very late reply. It looks like you have a very old version of Boost. You might want to update that. Best wishes.
-
Sorry for the very late reply. If the abundance per nucleon of species (Z, A) is Y(Z, A), then the mass fraction is A * Y(Z, A). The mass fraction is the number of grams of species (Z, A) in a sample per total gram of material in the sample. Thus, if the mass fraction in the Sun is X(Z, A), the number of Solar masses of (Z, A) in the Sun is X(Z, A) * M_sun, so X(Z, A) solar masses. I hope that helps. Best wishes.
-
Thanks for your post, and sorry for the late reply. The value 2.431E+10 is the number of H per 106 Si atoms. The 0.7110 is the H mass fraction. We converted the abundances per 106 Si atoms to abundances per nucleon and the computed mass fractions to get the Lodders.xml and Lodders_0.xml files. I hope that answers the question. Best wishes.
-
Sorry for the very late reply. It looks like you haven't installed xsltproc. Could you install xsltproc and trying compiling again? Best wishes.
-
Sorry for the late reply. I used IDL, but I now recommend wnutils.
-
Sorry for the very late response. The scaling option on ./flow_graph.sh allows the flow scaling to be logarithmic. That would help show small flows. Type ./flow_graph.sh to see the options. I hope that helps. Best wishes.
-
Thanks for your post. "[a = 1]" XPath selects all nuclei with A = 1 (neutrons and protons). "[a = 1 or a = 30]" selects nuclei with A = 1 or A = 30, thus, neutrons, protons, and A = 30 nuclei. You need the neutrons and protons to calculate binding energies., which are the differences between the mass of the species and the constituent nucleons. Imagine pulling the nucleus apart into all its neutrons and protons--how much energy would that require? Thats the binding energy. I obviously need to know...
-
Sorry for the very late reply. Please see the post about two below here. Best wishes.
-
Thanks for your post, and sorry for the late reply. Have you looked at this blog entry? It might help. if not, let me know. Best wishes.
-
Thanks for your post, and apologies for the late reply. It seems you are not able to link to the sourceforge servers, perhaps because you are behind a firewall. You can download the wn_matrix, libnucnet, libnuceq, and libstatmech tarbars directly to nucnet-tools-code/vendor and then compile. for example. svn checkout svn://svn.code.sf.net/p/nucnet-tools/code/trunk nucnet-tools-code mkdir vendor Now download the various tarballs to the vendor directory. Do this from their websites: wn_matrix, libnucnet,...
-
Thanks for your post, and apologies for the reply. It seems you are not able to link to the sourceforge servers, perhaps because you are behind a firewall. You can download the wn_matrix, libnucnet, libnuceq, and libstatmech tarbars directly to nucnet-tools-code/vendor and then compile. for example. svn checkout svn://svn.code.sf.net/p/nucnet-tools/code/trunk nucnet-tools-code mkdir vendor Now download the various tarballs to the vendor directory. Do this from their websites: wn_matrix, libnucnet,...
-
Update libncnet.
-
The calculations for the blog post used the original webnucleo data file. More recent data files have faster beta decay from 273Ac. This competes strongly with the fission, which leads to little fission during the calculation. With the more recent files, you can use the XPath expression "[z < 90]", or even "[z < 90 and a <= 273]" to force 273Ac to fission, like in this post.
-
This seems to be a firewall issue. Are you able to download libnucnet.tar.gz directly: https://sourceforge.net/projects/libnucnet/ Click on the Download link. Let me know what happens. Best wishes.
-
It looks like a firewall issue. You might be able to use wget -dv ... to debug. Type something like: wget -dv -O /dev/null --timeout=300 --tries=1 http://sourceforge.net/projects//libnucnet/files/libnucnet/libnucnet_0.40.tar.gz as described here Let me know how it goes. Best wishes.
-
Thanks for your post. Make sure that the name run_single_zone_new is among the target executables in the Makefile in your new directory. Thus, change NETWORK_EXEC = run_constant_entropy \ run_entropy \ run_energy_generation \ run_multiple_zone_omp \ run_single_zone to NETWORK_EXEC = run_constant_entropy \ run_entropy \ run_energy_generation \ run_multiple_zone_omp \ run_single_zone. \ run_single_zone_new Let me know if that doesn't work. Best wishes.
-
Thanks for your post. Are you behind a firewall? It looks like the compilation failed to download libstatmech correctly. You should be able to try ./project_make again if you are not behind a firewall. Best wishes.
-
Thanks for your post. Could you look at this discussion? Best wishes.
-
Thanks for your post. 0.27 and 0.73 are the approximate Solar System mass fractions from Lodders. I tweaked them slightly to ensure the sum was unity. Best wishes.
-
Thanks for your post. The code was not originally written to allow the gamma as a reactant. That should be possible, so I've updated the code. Do svn up in your talys_to_webnucleo directory, compile again, and try running. It should work. If not, let me know. Best wishes.
-
Extract substrings more generally. Allow for reactant gamma.
-
Thanks for your post. The download of libstatmech is timing out. Are you perhaps behind a firewall? What happens if you simply type wget http://sourceforge.net/projects//libstatmech/files/libstatmech/libstatmech_0.10.tar.gz (all on one line)? Best wishes.
-
Interested users may find this Cygwin cheat sheet useful.
-
Sorry for the delayed response. Are you still getiing these errors? Best wishes.
-
Thanks for your post. To be honest, I've lost track of where the nova profile came from. Sorry about that. On the second question, there are discussions going on to create a trajectory library. That will take some time, but I think it would be a great resource. I'll put a note to myself to announce it if it ever does get constructed. Best wishes.
-
Thanks for your post. It is a logical question, and we do have some new codes coming on line that do that. Unfortunately, those are not yet released and documented. I will try to announce those here when we do release them. For now, I guess you can evolve to time 1, use the output from that calculation in a new calculation to time 2, and make the integrated current diagram for that. That's very crude, but at least it gives you something to go on. Best wishes.
-
Thanks for your post. I recommend you use wnutils on the output from your calculations. There are tutorials that demonstrate how to do that as well as some Jupyter notebooks illustrating how to use wnutils. Best wishes.
-
Those are the input mass fractions.
-
Simplify.
1 >