Remove unused code.
Remove PismIntMask.

If we assume that a B.C. location mask contains zeros and
ones *only*, "(int)floor(bc_mask(i,j) + 0.5) == 1" is equivalent to
"bc_mask(i,j) > 0.5".