A Mathematica based Version of the CKMfitter Package

D.3. LoadLUT 69
else
dTableauO2 = 0D0
k = int( (dble(n-1)*ObsPred +xmax -dble(n)*xmin)/(xmax -xmin) )
SplinePred = ObsPred -tabsave(fn,0,k)
dTableauO2 = tabsave(fn,3,k) +
> 2D0 *tabsave(fn,4,k) *SplinePred +
> 3D0 *tabsave(fn,5,k) *SplinePred**2
end if
End
D.3 LoadLUT
dTableauO2 = dTableauO2 /2D0
!-----------------------------------------------------------------------
C
C subroutine to load LUT’s for Tableau function
C
subroutine LoadLUT( nbOfLUT, LUTfname)
Implicit None
include ’dimarray.f’
integer nbOfLUT
character*200 LUTfname(maxLUT)
integer i, fn, lengthTable(maxLUT)
double precision tabsave(maxLUT,0:5,maxGranularity)
common /LUT/ lengthTable, tabsave
do fn = 1, nbOfLUT
open(22, file = LUTfname(fn), status = ’old’)
read(22,*) lengthTable(fn)