rmf_create_huf creates an RMODFLOW huf object.
rmf_create_huf( dis, ihufcb = 0, hdry = -888, nhuf, iohufheads = 0, iohufflows = 0, lthuf = rmfi_ifelse0(dis$nlay == 1, 1, c(1, rep(0, dis$nlay - 1))), laywt = rep(0, dis$nlay), wetfct = 0.1, iwetit = 1, ihdwet = 0, wetdry = NULL, hgunam, top, thck, hguhani = rep(1, nhuf), hguvani = rep(0, nhuf), parameters = NULL, hk = rep(1e-04, nhuf), hani = rep(1, nhuf), vk = hk, ss = rep(1e-05, nhuf), sy = rep(0.15, nhuf), sytp = NULL )
| dis | RMODFLOW dis object |
|---|---|
| ihufcb | flag and unit number for writing cell-by-cell flow terms; defaults to 0 |
| hdry | head assigned to cells that are converted to dry cells; defaults to -888 |
| nhuf | number of hydrogeological units (hgu's) |
| iohufheads | flag and unit number indicating whether interpolated heads should be written. Defaults to 0 (interpolated heads are not written). See details. |
| iohufflows | flag and unit number indicating whether interpolated cell-by-cell flows should be written. Defaults to 0 (interpolated flows are not written). See details. |
| lthuf | vector length |
| laywt | vector length |
| wetfct | is a factor that is included in the calculation of the head that is initially established at a cell when it is converted from dry to wet; defaults to 0.1 |
| iwetit | is the iteration interval for attempting to wet cells; defaults to 1 |
| ihdwet | is a flag that determines which equation is used to define the initial head at cells that become wet; defaults to 0 |
| wetdry | 3d array with a wetting threshold and flag indicating which neighboring cells can cause a cell to become wet; defaults to NULL. If not read for a specific layer, set all values in that layer to NA. |
| hgunam | character vector of length nhuf specifying the names of the hgu's |
| top | 3d array with the top elevations of each hgu (hgu is the third dimension) |
| thck | 3d array with the thicknesses of each hgu (hgu is the third dimension) |
| hguhani | vector of flags indicating whether HANI parameters (0) or hguhani (> 0) is used to specify horizontal anisotropies for each hgu; defaults to 1 for each hgu |
| hguvani | vector of flags indicating whether vk parameters will be specified (0) or HGUVANI is used for vertical anisotropy (> 0). VANI parameters may overwrite the the latter value; defaults to 0 for each hgu |
| parameters | list of |
| hk | (named) vector with hydraulic conductivity along rows; one value per hgu specified in the order of |
| hani | (named) vector with the ratio of hydraulic conductivity along columns to that along rows; one value per hgu specified in the order of |
| ss | (named) vector with specific storage; only required when there are transient stress periods; one value per hgu specified in the order of |
| sy | (named) vector with specific yield; only required when there are transient stress periods; one value per hgu specified in the order of |
| sytp | numeric value defining sytp; only required when there are transient stress periods and all layers are confined; defaults to NULL. See details. |
| vka | (named) vector with vertical hydraulic conductivity or the ratio of horizontal to vertical conductivity depending on the value of hguvani; one value per hgu specified in the order of |
Object of class huf
If iohufheads (iohufflows) > 0, heads (cell-by-cell flows) will be interpolated to each hgu and saved, resulting in a single value per hgu at every row-column combination.
Flow variables are any of HK, HANI, VK, VANI, SS, SY and SYTP. Variables can be specified either through the use of parameters or direct vector input through hk, hani, vk, ss, sy and sytp arguments.
Specifying directly with parameters is useful if the parametrization is complex. If only a single value per hgu per partyp is needed, the direct vector input can be used.
If any parameters are specified, all flow variables have to be specified by parameters. If hk, hani, vk, ss, sy or sytp are not named, the hgu number will be added as a suffix to the type (e.g. 'HK_1').
Flow variable SYTP can not be specified per hgu, only for all hgu's. Using direct vector input, VANI can not be specified, only VK. VANI can be set through hguvani.
\code{\link{rmf_convert_huf_to_grid}} can be used to convert parameters defined on the HUF grid to the numerical grid