Fix a bug in construction of h_loc

The quadratic terms are now OK. The code still doesn't work
because deltat is empty.