Printing zone abundances

[Bradley S. Meyer]

Question copied from this post:

Dear Prof. Meyer,
Can I ask you a question?
If we want to get the abundance as the function of nucleon number, we can type
"../analysis/print_abundances_vs_nucleon_number my_new_output.xml z "" > yz.txt"
or "../analysis/print_abundances_vs_nucleon_number my_new_output.xml a "" > ya.txt"
or "../analysis/print_abundances_vs_nucleon_number my_new_output.xml a "" > ya.txt"
Is it possible to get the two dimentional relationship, we can print the abundance for each nuclei. for example
atomic number neutron number abundance
Thank you very much for your help!!!

[Bradley S. Meyer]

Thanks for your post. I believe you are asking whether you can print out the Z, N, and abundance for a zone or zones (i.e., a timestep). Try using print_zone_abundances in examples/analysis. For example, to print the abundances in the last 5 zones, type

../analysis/print_zone_abundances my_new_output.xml "[position() > last() - 5]"

The output is Z, A, Y, X, where Z is the atomic number, A is the mass number, Y is the abundance, and X is the mass fraction. Of course the neutron number N is A - Z. For each zone (timestep) there is a header giving information and then some diagnostics at the end.

[Anonymous]

Dear Prof. Meyer,

Thank you very much!
It works very well!

Best wishes!