various fixes to make --equilib-IC work. However, it appears that --equilib-IC 3mode doesn't provide the correct solution. Need to check on that