compute Jacobian of observation operator at ub
Dh = JObsOp(ub)