Fixed a possible divide by zero and undefined density due to a bad call to find_rho_2d.  This caused optimized FPEs, sometimes.  Seems to run OK now.