This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision Next revision Both sides next revision | ||
nu_rnet [2014/06/05 13:54] schutrumpg |
nu_rnet [2014/06/05 14:57] schutrumpg |
||
---|---|---|---|
Line 1: | Line 1: | ||
- | **[Group 7]** | + | ===$\nu$Rnet=== |
- | + | ||
- | *** Group name:** $\nu$Rnet | + | |
*** Members:** | *** Members:** | ||
Line 33: | Line 31: | ||
{{ :out.jpg?500 |}} | {{ :out.jpg?500 |}} | ||
+ | To set the initial electron fraction $Y_e$ which is equivalent to the proto fraction due to charge neutrality, we ran the nse version of the xnet code. This calculates the equilibrium initial abundances for a given electron fraction, density and temperature. Therefore the lines in the code tagged by !NSE have to be uncommented and the lines tagged by !NOTNSE have to be commented. This can be done quickly with the following commands in Unix: | ||
+ | <code> | ||
+ | sed -e 's/^\!\(.*\!NSE\)/ \1/' conditions.f90 >conditions.f90- | ||
+ | sed -e 's/^\!\(.*\!NSE\)/ \1/' net_mpi.f90 >net_mpi.f90- | ||
+ | sed -e 's/^\!\(.*\!NSE\)/ \1/' net.f90 >net.f90- | ||
+ | sed -e 's/^\!\(.*\!NSE\)/ \1/' data_distribute_mpi.f90 >data_distribute_mpi.f90- | ||
+ | mv conditions.f90- conditions.f90 | ||
+ | mv net_mpi.f90- net_mpi.f90 | ||
+ | mv net.f90- net.f90 | ||
+ | mv data_distribute_mpi.f90- data_distribute_mpi.f90 | ||
+ | sed -e 's/^\!\(.*\!NOTNSE\)/!\1/' net.f90 >net.f90- | ||
+ | mv net.f90- net.f90 | ||
+ | </code> | ||
+ | Additionally the line | ||
+ | <code>NSE_OBJ = nse.o </code> | ||
+ | in the Makefile has to be uncommented. | ||
+ | |||
+ | It came out that at high Temperatures around 10 GK mainly free neutrons and free protons exist. The other abundances were found to be very low. Thus, for our final calculations we did not use the nse, but changed the proton fraction directly by varying the proton to neutron rate in the initial abundances and neglecting the all other abundances. |