solve_T_profile

Author:

Keith Cherkauer

Date:

June 11, 1997

Description:

This subroutine iteratively solves for the soil temperature profile using a numerical difference equation. The difference equation is second order in space, and first order in time. If the maximum number of iterations is exceeded the program will dump the last temperature profile to stderr, and exit.

The program makes use of soil thermal properties (conductivity, heat capacity), which are assinged to each of the Tlayer solution depths, based on soil moisture conditions.

References

Input Variables:

Variable	Type	Units	Description
T0	double *	C	Temperature profile at previous time step
dz	double *	m	Soil layer thickness
kappa	double *	W/m/K	Soil layer thermal conductivity
Cs	double *	J/m^3/K	Soil layer volumetric heat capacity
deltat	double	s	Model time step
Tlayer	int	N/A	Total number of layers in soil profile

Output Variables:

Variable	Type	Units	Description
T	double *	C	Temperature profile at current time step

Subroutines Called:

Modifications and Limitations

VIC Home Page

VIC Code Page

Page Last Modified on 6/12/97