def get_collocation_points(x, t, z, eta_meas_max, N_pde, N_bc_sur, N_bc_bot, N_pb): d = np.max(np.abs(z)) # water depth # domain bounds total lb = np.array([np.min(x ...