F.E. solution of a vertically averaged model for regional land subsidence