save_gpf seems to be working well. Note sure the test is ready though. It is there but uncertain if filecmp.cmp with shallow=True will allow it to pass since they are not exactly the same (e.g. number of decimal places). Also cat file be written too?