get_hazard_2d function

Get estimated (log-)hazard surface with 2 time scales