Added check that topo arrays cover full domain using new recursive routine topoarea in topo_module.f90