Added density dependent NSE, but it is still buggy.