  The nature of the zero temperature ordering transition in the 3D Gaussian
random field Ising magnet is studied numerically, aided by scaling analyses. In
the ferromagnetic phase the scaling of the roughness of the domain walls,
$w\sim L^\zeta$, is consistent with the theoretical prediction $\zeta = 2/3$.
As the randomness is increased through the transition, the probability
distribution of the interfacial tension of domain walls scales as for a single
second order transition. At the critical point, the fractal dimensions of
domain walls and the fractal dimension of the outer surface of spin clusters
are investigated: there are at least two distinct physically important fractal
dimensions. These dimensions are argued to be related to combinations of the
energy scaling exponent, $\theta$, which determines the violation of
hyperscaling, the correlation length exponent $\nu$, and the magnetization
exponent $\beta$. The value $\beta = 0.017\pm 0.005$ is derived from the
magnetization: this estimate is supported by the study of the spin cluster size
distribution at criticality. The variation of configurations in the interior of
a sample with boundary conditions is consistent with the hypothesis that there
is a single transition separating the disordered phase with one ground state
from the ordered phase with two ground states. The array of results are shown
to be consistent with a scaling picture and a geometric description of the
influence of boundary conditions on the spins. The details of the algorithm
used and its implementation are also described.
