fixed setaux.f90 to extrapolate topo properly to ghost cells in the case where interior grid does not start at (1,1).  Seems to work even on case where small patch has mx=-1.