Convert real world coordinates to modflow coordinates

rmf_convert_xyz_to_grid(
  dis,
  x,
  y,
  z = NULL,
  prj = rmf_get_prj(dis),
  output = "xyz"
)

Arguments

dis

RMODFLOW dis object

x

real world x coordinate

y

real world y coordinate

z

real world z coordinate; optional

prj

optional RMODFLOW prj object

output

character; containing 'xyz','ijk' and/or 'off' for the return of x, y, z, i, j, k, roff, coff and loff modflow coordinates

Value

data frame with modflow coordinates

Details

If prj is not provided, x, y and/or z input coordinates already represent modflow coordinates. If z is not provided, no third dimension coordinates are returned. If the xyz coordinate falls on a boundary of two cells, the minimum ijk indices are returned.

If the z coordinate falls within a Quasi-3D confining bed, the layer index of the overlying model layer is returned. The loff value then represents the fractional distance from the center of the overlying model layer.

The coordinate system has its origin in the lower left corner in the top layer. Row indices (i) increase from front to back (with decreasing Y values), column indices (j) increase from left to right (with increasing X values), layer indices (k) increase from top to bottom (with decreasing Z values).

See also