my_qromo

Author:

modified by Xu Liang, and Keith Cherkauer

Date:

Description:

This routine was modified from the Numerical Recipes routine qromo.c so that the function used is that declared in funct_design.c. qromo.c applies Romberg integration on an open interval. It returns ss as the integral of the function from a to b, using the midpoint rule integration technique of midpnt.c. Configured in this fashion the routine will solve an improper integral.

References:

See qromo.c or qromo.f in Numerical Recipes (pg 137 of Numerical Recipes in Fortran, 2^nd Edition)

Input Variables:

Variable	Type	Units	Description
a	double	N/A	Lower bound of integral
b	double	N/A	Upper bound of integral
ppt	double	mm	Average precipitation
mu	double	fraction	Fraction of grid cell that receives precipitation
im	double	N/A	Maximum infiltration capacity
i0	double	N/A
bi	double	N/A	Infiltration parameter
bp_index	double	N/A	Not used
func_index	int	N/A	Index for function type, not currently used

Output Variables:

Variable	Type	Units	Description
k	int *	N/A	Number of iterations for solution

Subroutines Called:

Modifications and Limitations:

VIC Home Page

VIC Code Page

Page Last Modified on 2/8/98