added Trevor MacDougall's neutral density code. Helped sophie compile it today and the steps were not easy. Had to fix typos in makefile and update it to use gfortran. All fortran code compiles successfully now. Haven't tested it in MATLAB yet