fsd_build.f90 358 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)
	elseif (FSD_scheme.eq.1) then
	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