docs: def median(ma):
      """ do it row by row, to save memory...."""
      _median = 0*ma[0].filled(fill_value=0)
      for i in range(ma.shape[-1]):
  	t=xmedian(ma[:,:,i])
      	_median[:,i]=t
  	t=None
      return _median