Added commented code in microphysical_ec that allows for switching between 100 lmp nucs and all abundant nucs with ni56 weak rate.
Now loops over all abundant nuclei
Full opacity table run with ni56 rates for nuclei with unknown rates.
Added extrapolation in weak rates out to lrhoYe=15.0
Changed table bounds to allow for extrapolation
Inprogress..
Added deltaye binary that outputs A,Z,dye,normalized_dye,mass_fraction,normalized_dye/mass_fraction to file for each step in time
Currently producing dye_i/sum(dye_i) and binning based on A from 0:5, 5:25, 45:65, ... , 385:405.
Added patch to qec_solver that returns after extrapolation
Merge branch 'qec_extrap_patch' into analysis_deltaye
Added time integrated analysis loop
Reproduction of Ye(t) and total contribution of each nucleus is completed.
Added analytic rates, performing test over range of electron density and temperature
Test program was modified to output an analysis of the approximate electron capture rate equation detailed in arxiv:astro-ph/0302459v1 for pf-shell nuclei.
Updated rate interpolation to return zero when outside of range
LMP rate tables list a reaction Q value with the reaction products and not the parent nuclei, adjusted read routine to match.
WIP: analytic integration
Added weakrates_density_extrapolation flag to nulib for easy toggling of extrapolation in electron density to log10(rhoYe)=15.0. Rate tables extend only to 12.5. This branch contains most recent patches and will be merged with master.
Updated Makefile, ready for merge