A macroscopic observation of a magnetic flux distribution on the surface of a superconducting bulk material was performed. First, a single-domained thin rectangular YBCO sample (10 mmxl0mmxlmm) was field-cooled in liquid nitrogen. After the applied field was turned off, a field distribution in the area of lOmmxlOmm was mapped by scanning the sample surface at intervals of 0. 5" using a Hall probe fixed to a XY-stage. The obtained field distribution was visualized into a three-dimensional distribution map and a twodimensional contour graph. The overall shape of the field distribution had a four-fold symmetry reflected by the sample geometry. Next, assuming that the critical current density J , is independent of the field strength, a field distribution generated by square current loops was calculated numerically. The intragranular J, was estimated to be 2.6x1@A/cm2 (77K, OT) by fitting the calculated profile to the measured profile. Finally, it is found that the measurement of field distributions on the surface of a superconductor using a Hall probe is applicable to an estimation of an ntragranular J , for a superconductor without destroying the samples.