Will vary gamma as a function of density, pressure and temperature
Will feed output of hydro code (time, temperature and density) into network code
Will use routine in network code to calculate initial abundance purely as a function of temperature, density and electron fraction.
iNSE lines commented out and fix makefile to include sed command
comment out not the NSE lines
Need MEUS for r process calculations