fsd_build.f90 324 Bytes
Newer Older
Jérémy Baudry's avatar
Jérémy Baudry committed
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
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