An Eulerian‐Lagrangian approach with an adaptively corrected method of characteristics to simulate variably saturated water flow