* 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