revised.modified.kostiakov: Parhi Revised modified kostiakov infiltration parameter...

In gowusu/vadose: The estimation of saturated hydraulic conductivity and soil water retention Curves in the vadose zone

Description

This function optimises Parhi et al (2007) revised modified kostiakov cumulative infiltration (I)

Usage

revised.modified.kostiakov(data = NULL, time, I, a1 = 0.1, a2 = 0.1,
  b1 = 0.1, b2 = 0.1, group = NULL)

## Default S3 method:
revised.modified.kostiakov(data = NULL, time, I, a1 = 0.1,
  a2 = 0.1, b1 = 0.1, b2 = 0.1, group = NULL)

## S3 method for class 'revised.modified.kostiakov'
predict(object, time = NULL, ...)

## S3 method for class 'revised.modified.kostiakov'
plot(x, xlab = "Time(Minutes)",
  ylab = "Cumulative (cm)", main = NULL, layout = NULL, ...)

## S3 method for class 'revised.modified.kostiakov'
summary(object, ...)

## S3 method for class 'revised.modified.kostiakov'
coef(object, ...)

Arguments

data	dataframe. It can contain data with column names of "time" and "I"
time	character or numeric. The name of time variable in the dataframe. If the "data" parameter contains "time", this will be ignored. The unit must be in seconds.
I	cumulative infiltration (I) or infiltration rate (i) [mm]
a1	constant
a2	constant
b1	constant
b2	constant
group	character. The name of the group variables if the data is from different areas.
object	Model output object
...	Any other graphical parameter
x	a return object of the function.
xlab	x label of the plot
ylab	y label of the plot
main	Title of the plot
layout	plot layout

Value

• a1: optimised constant

• a2: optimised constant

• b1optimised constant

• b2optimised constant

• output: output of the group simulation

Author(s)

George Owusu

References

Parhi, P. K., Mishra, S. K., & Singh, R. (2007). A modification to Kostiakov and modified Kostiakov infiltration modeles. Water Resour. Manage., 21, 1973-1989.

Examples

data=read.csv(system.file("ext","sys","exampleBEST.csv",package="vadose"))
revised.modified.kostiakov1<-revised.modified.kostiakov(data=data,time="time",I="I")
#print(gof(revised.modified.kostiakov1))
plot(revised.modified.kostiakov1)
predict(revised.modified.kostiakov1)
coef(revised.modified.kostiakov1)
#gof(revised.modified.kostiakov1)

#group simulation
data=read.csv(system.file("ext","sys","infiltration2.csv",package="vadose"))
assin_breko<- subset(data, ID=="41A20_1")
revised.modified.kostiakov1vg<-revised.modified.kostiakov(data=data,time="minutes",
I="CumInfil",group="ID")
coef(revised.modified.kostiakov1vg)

gowusu/vadose documentation built on May 17, 2019, 7:59 a.m.