I hope that PISMHydrology::max_timestep() is closer to the right thing.  And it now reads option -input_to_bed_reference_year.  But attempt to run causes hang; see issue soon.