estimate_ksat.f

	subroutine estimate_ksat(perc_clay,esti_ksat)

!!	This subroutine calculates ksat value for a soil layer
!!	given the % of clay in the soil layer

!!	Background: Published work of Walter Rawls- calculated 
!!	ksat values  based on soil texture (sand,silt and clay)
!!	idea: there exists a relationship between % clay and Ksat
!!	Equations used in this subroutine are based on the above
!!	idea (Jimmy willimas). 

!!	NK June 28,2006	

	implicit none

	integer::i,eof
	real::esti_ksat,perc_clay,xc,exksat
	exksat=5.0

!	print *,"Enter the % clay in the soil layer"
!	read *,perc_clay

	xc=100.0-perc_clay
	esti_ksat=12.7*xc/(xc+exp(11.45-0.097*xc))+1.0

!	print *,"The estimated ksat value is ",min(esti_ksat,exksat)

	return

	end subroutine estimate_ksat