get_hazard_2d_LMM function

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