wellanalyse: Analyse receptors in a DRW model output

Description Usage Arguments Value

Description

Extract the flux to individual receptors and calculate the abstracted concentration to receptors in the output of a DRW model.

Usage

1
2
3
wellJ(dr, welref)

wellC(dr, WEL, DIS, MFt0, welref = WEL)

Arguments

dr

a DRWmodel S4 object,

welref

data.table (or object coercible to one) with named columns C, R and L giving the CRL references of receptors to be analysed, integer [3], NULL, WEL.MFpackage (see read.WEL), or (for wellJ only) various other MFpackage objects; NULL implies all cell references to which a flux has been recorded by DRW

WEL

WEL.MFpackage or NetCDF object of MODFLOW outputs; if a NetCDF object, there must be a "Wells" array in the data set

DIS

DIS.MFpackage (see read.DIS); only required if WEL is a WEL.MFpackage object

MFt0

numeric [1]; only required if WEL is a WEL.MFpackage object start time of the MODFLOW model, which should be the same as for the NetCDF MODFLOW data set used for the DRW model

Value

wellJ: list of fluxes to each receptor specified in welref, with a value for each time step in dr (dr@time)
wellC: as for wellJ but concentrations, determined by diluting the flux with the transient abstraction rate

In each case, the elements of the list are named giving the cell reference as a character string (e.g. "C30R20L1")


CJBarry/DRW documentation built on May 6, 2019, 9:25 a.m.