Module nrtype
		integer, parameter :: dp = kind(1.d0)
		integer, parameter :: sp = kind(1.0)
end module nrtype