xr = randomu(s, 10000)*1024
yr = randomu(s, 10000)*1024

parLambda = readfits("/Users/colkin/Work04/Pluto_3um/Reduction/LambdaMap/LamMapParams4.fits")

.compile /Users/colkin/Work04/Pluto_3um/IDL_Pro/bivariateCubic.pro 

BIVARIATECubic, parLambda, xr, yr, lambdaR

;BIVARIATECubic, parLambda, FINDGEN(1024), FINDGEN(1024), lambdaOUT

parP = readfits('/Users/colkin/Work04/Pluto_3um/Reduction/PMap/pMapParamsCubic.fits')

BIVARIATECUBIC, parP, xr, yr, pR

ptsX=fltarr(10000,3)

ptsX[*, 0] = xr
ptsX[*, 1] = lambdaR
ptsX[*, 2] = pR

.compile /Users/colkin/Work04/Triton_NIRSPEC/IDL_Procedures/WaveCal.Pro

wavecal, ptsX, 3, 1024., 1024., parX, residX, chiSqX, xmap

ptsY=fltarr(10000,3)

ptsY[*, 0] = yr
ptsY[*, 1] = lambdaR
ptsY[*, 2] = pR

wavecal, ptsY, 4, 1024., 1024., parY, residY, chiSqY, ymap

dir = "/Users/colkin/Work04/Pluto_3um/Reduction/Rectify/"
writefits, dir+'xMapParams7.fits', parx
writefits, dir+'yMapParams7.fits', pary
writefits, dir+'xMap7.fits', xmap
writefits, dir+'yMap7.fits', ymap

np=300.
nlam = 1024.

lamg = replicate(1.0, np)##(2.7+(1./1024)*findgen(nlam))
pg = (findgen(np)+170.)##replicate(1.0, nlam)

BIVARIATECUBIC, parx, lamg, pg, xg
.compile '/Users/colkin/Work04/Pluto_3um/IDL_Pro/bivariateQuartic.pro'

BIVARIATEQUARTIC, pary, lamg, pg, yg

wavecal, ptsX, 4, 1024., 1024., parX4, residX4, chiSqX4, xmap4

writefits, dir+'xMapParams8.fits', parx4
writefits, dir+'xMap8.fits', xmap4