fsd_build.f90 324 Bytes
subroutine fsd_build
use parameters

	implicit none
	double precision :: coeff

	if (Dmax(i).eq.Dmin) then
	Dave(i)=Dmin
	elseif (Dmax(i).eq.D0) then
	Dave(i)=Dmax(i)
	else
	coeff=1/((1/(1-gam))*(Dmax(i)**(1-gam)-Dmin**(1-gam)))
	Dave(i)=coeff*(1/(2-gam))*(Dmax(i)**(2-gam)-Dmin**(2-gam))
	end if


end subroutine fsd_build