Numerical study of the F model with domain-wall boundaries