test_interpolate2.f90
No OneTemporary
Actions

Size

1 KB

Subscribers

None

test_interpolate2.f90
View Options


	!-------------------------------------------------------------------
	program test_interpolate2
	!-------------------------------------------------------------------
	! This program was designed to test the subroutine interpolate_tabletype2
	! It is compiled and run using:
	! gfortran -L../HiggsBounds_KW/ -I../HiggsBounds_KW/ test_interpolate2.f90 -o test_interpolate2.exe -lHB
	! ./test_interpolate2.exe
	!-------------------------------------------------------------------
	use interpolate
	use S95tables
	use S95tables_type2
	use theory_BRfunctions, only : setup_BRSM

	implicit none
	integer :: elementnum,datcomp,id
	double precision :: x,y
	double precision :: interpol

	call setup_BRSM
	call setup_S95tables

	id=200
	!x=39.2D0
	!y=78.8D0
	!x=39.5D0
	!y=79.1D0
	x=39.5D0
	y=78.9D0
	datcomp=1

	elementnum=t2elementnumberfromid(S95_t2,id)

	call interpolate_tabletype2(x,y,S95_t2(elementnum),datcomp,interpol)

	write(,)'id=',id
	write(,)'x=',x
	write(,)'y=',y
	write(,)'interpol=',interpol

	end program test_interpolate2