working simudark comparison plot_simudark_corrections.py.