#pragma rtGlobals=1		// Use modern global access method.
#pragma IgorVersion = 4.0
#pragma version = 1.0
#include "CromerLieberman", version>=1.0

//#include  "Elements", version>=1.0

//Menu "Analysis"
//	Submenu "X-ray"
//		"X-ray atomic structure factor",GetfoXray("",NaN)
//	End
//End


Function GetfoXray(symb,K)	// compute fo(Q)
	String symb				// name of atom desired
	Variable K					// 2*pi/d	 (1/Angstroms)

	if (exists("root:Packages:Xray:focoefs")!=1)
		XrayDataInitPackage()
	endif

	symb = LowerStr(symb)
	symb[0,0] = UpperStr(symb[0])
	if (strlen(symb)<1 || numtype(K) || K<0)
		K = numtype(K) ? 0 : abs(K)
		SVAR symbs=root:Packages:Xray:symbs
		Prompt K, "Length of Q, 2pi/d (1/Angstrom)"
		Prompt symb, "symbol of atom", popup, symbs
		DoPrompt "pick atom and K", symb,K
		if (V_flag)
			return 1
		endif
	endif

	Variable fo = foXray(symb,K)
	if (numtype(fo))
		printf "  fo(%s,  K=%g) = Unknown symbol\r",symb,K
	else
		printf "  fo(%s,  K=%g) = %g  (electrons)\r",symb,K, fo
	endif
	return fo
End


Function foXray(symb,K)		// compute fo(Q)
	String symb				// name of atom desired
	Variable K					// 2*pi/d	 (1/Angstroms)

	SVAR symbs=root:Packages:Xray:symbs
	Wave focoefs=root:Packages:Xray:focoefs
	if (!WaveExists(focoefs) || !SVAR_Exists(symbs))
		print "you forgot to init this package, do a 'XrayDataInitPackage()'"
		Abort "you forgot to init this package, do a 'XrayDataInitPackage()'"
	endif
	symb = LowerStr(symb)
	symb[0,0] = UpperStr(symb[0])
	Variable m = WhichListItem(symb,symbs)	// index into symbs
	if (m<0)
		return NaN			// bad symbol
	endif

	Wave foLast=root:Packages:Xray:foLast
	Wave KLast=root:Packages:Xray:KLast
	Wave/T symbLast=root:Packages:Xray:symbLast
	NVAR nlast=root:Packages:Xray:nlast
	Variable Nbuf=numpnts(root:Packages:Xray:foLast), i,j
	for (i=0;i<Nbuf;i+=1)
		j = mod(i+nlast,Nbuf)
		if (stringmatch(symbLast[j],symb) && abs(KLast[j]-K)<.0001)		// user buffered value
			return foLast[j]
		endif
	endfor

	Variable s2=(abs(K)/4/PI)^2		// (sin(theta)/lambda)^2
	Variable fo=focoefs[8][m]			// the answer
	for(i=0;i<=3;i+=1)
		fo += (focoefs[i][m])*exp(-(focoefs[i+4][m])*s2)
	endfor
	//	focoefs[0]=a0,  focoefs[1]=a1,  focoefs[2]=a2,  focoefs[3]=a3
	//	focoefs[4]=b0,  focoefs[5]=b1,  focoefs[6]=b2,  focoefs[7]=b3,  
	//	focoefs[8]=c

	nLast = mod(nLast+1,Nbuf)			// save last value for likely reuse, for efficiency
	foLast[nLast] = fo
	KLast[nLast] = K
	symbLast[nLast] = symb
	return fo
End


Function XrayDataInitPackage()
	NewDataFolder /O root:Packages
	NewDataFolder /O root:Packages:Xray

	String/G root:Packages:Xray:symbs
	SVAR symbs=root:Packages:Xray:symbs
	symbs="H;H-1;He;Li;Li+1;Be;Be+2;B;C;N;O;O-1;F;F-1;Ne;Na;Na+1;Mg;Mg+2;Al;Al+3;"
	symbs += "Si;Si+4;P;S;Cl;Cl-1;Ar;K;K+1;Ca;Ca+2;Sc;Sc+3;Ti;Ti+3;Ti+4;V;V+2;V+3;V+5;"
	symbs += "Cr;Cr+2;Cr+3;Mn;Mn+2;Mn+3;Mn+4;Fe;Fe+2;Fe+3;Co;Co+2;Co+3;Ni;Ni+2;Ni+3;"
	symbs += "Cu;Cu+1;Cu+2;Zn;Zn+2;Ga;Ga+3;Ge;As;Se;Br;Br-1;Kr;Rb;Rb+1;Sr;Sr+2;Y;Y+3;"
	symbs += "Zr;Zr+4;Nb;Nb+3;Nb+5;Mo;Mo+3;Mo+5;Mo+6;Tc;Ru;Ru+3;Ru+4;Rh;Rh+3;Rh+4;"
	symbs += "Pd;Pd+2;Pd+4;Ag;Ag+1;Ag+2;Cd;Cd+2;In;In+3;Sn;Sn+2;Sn+4;Sb;Sb+3;Sb+5;Te;"
	symbs += "I;I-1;Xe;Cs;Cs+1;Ba;Ba+2;La;La+3;Ce;Ce+3;Ce+4;Pr;Pr+3;Pr+4;Nd;Nd+3;"
	symbs += "Pm;Pm+3;Sm;Sm+3;Eu;Eu+2;Eu+3;Gd;Gd+3;Tb;Tb+3;Dy;Dy+3;Ho;Ho+3;"
	symbs += "Er;Er+3;Tm;Tm+3;Yb;Yb+2;Yb+3;Lu;Lu+3;Hf;Ta;W;Re;Os;Ir;Pt;Au;Hg;"
	symbs += "Tl;Pb;Bi;Po;At;Rn;Fr;Ra;Ac;Th;Pa;U;Np;Pu;Am;Cm;Bk;Cf;Es;Fm;Md;No;Lw"
//	symbs="H;H-1;HE;LI;LI+1;BE;BE+2;B;C;N;O;O-1;F;F-1;NE;NA;NA+1;MG;MG+2;AL;AL+3;"
//	symbs += "SI;SI+4;P;S;CL;CL-1;AR;K;K+1;CA;CA+2;SC;SC+3;TI;TI+3;TI+4;V;V+2;V+3;V+5;"
//	symbs += "CR;CR+2;CR+3;MN;MN+2;MN+3;MN+4;FE;FE+2;FE+3;CO;CO+2;CO+3;NI;NI+2;NI+3;"
//	symbs += "CU;CU+1;CU+2;ZN;ZN+2;GA;GA+3;GE;AS;SE;BR;BR-1;KR;RB;RB+1;SR;SR+2;Y;Y+3;"
//	symbs += "ZR;ZR+4;NB;NB+3;NB+5;MO;MO+3;MO+5;MO+6;TC;RU;RU+3;RU+4;RH;RH+3;RH+4;"
//	symbs += "PD;PD+2;PD+4;AG;AG+1;AG+2;CD;CD+2;IN;IN+3;SN;SN+2;SN+4;SB;SB+3;SB+5;TE;"
//	symbs += "I;I-1;XE;CS;CS+1;BA;BA+2;LA;LA+3;CE;CE+3;CE+4;PR;PR+3;PR+4;ND;ND+3;"
//	symbs += "PM;PM+3;SM;SM+3;EU;EU+2;EU+3;GD;GD+3;TB;TB+3;DY;DY+3;HO;HO+3;"
//	symbs += "ER;ER+3;TM;TM+3;YB;YB+2;YB+3;LU;LU+3;HF;TA;W;RE;OS;IR;PT;AU;HG;"
//	symbs += "TL;PB;BI;PO;AT;RN;FR;RA;AC;TH;PA;U;NP;PU;AM;CM;BK;CF;ES;FM;MD;NO;LW"

	Make/N=(9,182)/O root:Packages:Xray:focoefs
	Wave focoefs=root:Packages:Xray:focoefs
	//	focoefs[0]=a0,  focoefs[1]=a1,  focoefs[2]=a2,  focoefs[3]=a3
	//	focoefs[4]=b0,  focoefs[5]=b1,  focoefs[6]=b2,  focoefs[7]=b3,  
	//	focoefs[8]=c
  	focoefs[][000] = {0.39875, 0.31285, 0.2144,  0.07135, 58.3331, 14.7175, 236.7147, 3.4848, 0.00125}
  	focoefs[][001] = {0.7975,  0.6257,  0.4288,  0.1427,  58.3331, 14.7175, 236.7147, 3.4848, 0.0025}
  	focoefs[][002] = {0.76844, 0.72694, 0.27631, 0.21572, 10.9071, 4.30779, 1.33127, 25.6848, 0.01249}
  	focoefs[][003] = {0.99279, 0.87402, 0.84240, 0.23131, 4.33979, 1.26006, 98.7088, 212.088, 0.05988}
  	focoefs[][004] = {6.08475, 0.86773, 0.80588, 0.17720, 0.00498, 1.53730, 4.28524, 9.81413,-5.93560}
  	focoefs[][005] = {2.22744, 1.55249, 1.40060, 0.58290, 0.04965, 42.9165, 1.66379, 100.361,-1.76339}
  	focoefs[][006] = {5.69034, 1.19706, 1.03057, 0.20150,-0.01336, 0.39000, 1.97441, 4.90642,-6.11950}
  	focoefs[][007] = {2.03876, 1.41491, 1.11609, 0.73273, 23.0888, 0.97848, 59.8985, 0.08538,-0.30409}
  	focoefs[][008] = {1.93019, 1.87812, 1.57415, 0.37108, 12.7188, 28.6498, 0.59645, 65.0337, 0.24637}
  	focoefs[][009] = {12.7913, 3.28546, 1.76483, 0.54709, 0.02064, 10.7018, 30.7773, 1.48044,-11.3926}
  	focoefs[][010] = {2.95648, 2.45240, 1.50510, 0.78135, 13.8964, 5.91765, 0.34537, 34.0811, 0.30413}
  	focoefs[][011] = {3.22563, 3.01717, 1.42553, 0.90525, 18.4991, 6.65680, 0.40589, 61.1889, 0.42362}
  	focoefs[][012] = {3.30393, 3.01753, 1.35754, 0.83645, 11.2651, 4.66504, 0.33760, 27.9898, 0.48398}
  	focoefs[][013] = {3.63220, 3.51057, 1.26064, 0.94071, 5.27756, 14.7353, 0.44226, 47.3437, 0.65340}
  	focoefs[][014] = {3.71272, 3.52631, 1.19237, 0.8308, 3.91091, 9.63126, 0.40483, 23.9546, 0.73728}
  	focoefs[][015] = {5.26400, 2.17549, 1.36690, 1.08859, 4.02579, 10.4796, 0.84222, 133.617, 1.09912}
  	focoefs[][016] = {3.99479, 3.37245, 1.13877, 0.65118, 3.11047, 7.14318, 0.40692, 15.7319, 0.84267}
  	focoefs[][017] = {5.59229, 2.68206, 1.72235, 0.73055, 4.41142, 1.36549, 93.4885, 32.5281, 1.26883}
  	focoefs[][018] = {4.30491, 3.14719, 1.12859, 0.49034, 2.55961, 5.60660, 0.41574, 11.4840, 0.92893}
  	focoefs[][019] = {5.35047, 2.92451, 2.27309, 1.16531, 3.48665, 1.20535, 42.6051, 107.170, 1.28489}
  	focoefs[][020] = {4.17448, 3.38760, 1.20296, 0.52814, 1.93816, 4.14553, 0.22875, 8.28524, 0.70679}
  	focoefs[][021] = {5.79411, 3.22390, 2.42795, 1.32149, 2.57104, 34.1775, 0.86937, 85.3410, 1.23139}
  	focoefs[][022] = {4.43918, 3.20345, 1.19453, 0.41653, 1.64167, 3.43757, 0.21490, 6.65365, 0.74630}
  	focoefs[][023] = {6.92073, 4.14396, 2.01697, 1.53860, 1.83778, 27.0198, 0.21318, 67.1086, 0.37870}
  	focoefs[][024] = {7.18742, 5.88671, 5.15858, 1.64403, 1.43280, 0.02865, 22.1101, 55.4651,-3.87732}
  	focoefs[][025] = {9.83957, 7.53181, 6.07100, 1.87128,-0.00053, 1.11119, 18.0846, 45.3666,-8.31430}
  	focoefs[][026] = {18.0842, 7.47202, 6.46337, 2.43918, 0.00129, 1.12976, 19.3079, 59.0633,-16.4654}
  	focoefs[][027] = {16.8752, 8.32256, 6.91326, 2.18515,-0.01456, 0.83310, 14.9177, 37.2256,-16.2972}
  	focoefs[][028] = {8.11756, 7.48062, 1.07795, 0.97218, 12.6684, 0.76409, 211.222, 37.2727, 1.35009}
  	focoefs[][029] = {9.70659, 7.37245, 5.67228, 1.90668, 0.59947, 11.8765,-0.08359, 26.7668,-6.65819}
  	focoefs[][030] = {8.60272, 7.50769, 1.75117, 0.96216, 10.2636, 0.62794, 149.301, 60.2274, 1.17430}
  	focoefs[][031] = {13.2063, 11.0586, 7.73221, 1.72057, 0.39466,-0.08204, 9.62976, 20.3341,-15.7176}
  	focoefs[][032] = {9.06482, 7.55526, 2.05017, 1.28745, 8.77431, 0.53306, 123.880, 36.8890, 1.03849}
  	focoefs[][033] = {13.4008, 8.02730, 1.65943, 1.57936, 0.29854, 7.96290,-0.28604, 16.0662,-6.66668}
  	focoefs[][034] = {9.54969, 7.60067, 2.17223, 1.75438, 7.60579, 0.45899, 109.099, 27.5715, 0.91762}
  	focoefs[][035] = {17.7344, 8.73816, 5.25691, 1.92134, 0.22061, 7.04716,-0.15762, 15.9768,-14.6519}
  	focoefs[][036] = {19.5114, 8.23473, 2.01341, 1.52080, 0.17885, 6.67018,-0.29263, 12.9464,-13.2803}
  	focoefs[][037] = {10.0661, 7.61420, 2.23551, 2.23170, 6.67721, 0.40322, 98.5954, 22.5720, 0.84574}
  	focoefs[][038] = {9.34513, 7.68833, 2.94531, 0.26998, 6.49985, 0.39491, 15.9868, 41.0832, 0.75143}
  	focoefs[][039] = {9.43141, 7.74190, 2.15343, 0.01686, 6.39535, 0.38335, 15.1908, 63.9690, 0.65657}
  	focoefs[][040] = {15.6887, 8.14208, 2.03081,-9.57602, 0.67900, 5.40135, 9.97278, 0.94046, 1.71430}
  	focoefs[][041] = {10.4757, 7.51402, 3.50115, 1.54902, 6.01658, 0.37426, 19.0654, 97.4599, 0.95226}
  	focoefs[][042] = {9.54034, 7.75090, 3.58274, 0.50911, 5.66078, 0.34426, 13.3075, 32.4224, 0.61690}
  	focoefs[][043] = {9.68090, 7.81136, 2.87603, 0.11357, 5.59463, 0.33439, 12.8288, 32.8761, 0.51827}
  	focoefs[][044] = {11.2519, 7.36935, 3.04107, 2.27703, 5.34818, 0.34373, 17.4089, 84.2139, 1.05195}
  	focoefs[][045] = {9.78094, 7.79153, 4.18544, 0.72736, 4.98303, 0.30421, 11.4399, 27.7750, 0.51454}
  	focoefs[][046] = {9.84521, 7.87194, 3.56531, 0.32361, 4.91797, 0.29439, 10.8171, 24.1281, 0.39397}
  	focoefs[][047] = {9.96253, 7.97057, 2.76067, 0.05445, 4.84850, 0.28330, 10.4852, 27.5730, 0.25188}
  	focoefs[][048] = {11.9185, 7.04848, 3.34326, 2.27228, 4.87394, 0.34023, 15.9330, 79.0339, 1.40818}
  	focoefs[][049] = {10.1270, 7.78007, 4.71825, 0.89547, 4.44133, 0.27418, 10.1451, 24.8302, 0.47888}
  	focoefs[][050] = {10.0333, 7.90625, 4.20562, 0.55048, 4.36007, 0.26250, 9.35847, 20.4105, 0.30429}
  	focoefs[][051] = {12.6158, 6.62642, 3.57722, 2.25644, 4.48994, 0.35459, 14.8402, 74.7352, 1.91452}
  	focoefs[][052] = {10.5942, 7.67791, 5.15947, 1.01440, 4.00858, 0.25410, 9.21931, 22.7516, 0.55358}
  	focoefs[][053] = {10.3380, 7.88173, 4.76795, 0.72559, 3.90969, 0.23867, 8.35583, 18.3491, 0.28667}
  	focoefs[][054] = {13.3239, 6.18746, 3.74792, 2.23195, 4.17742, 0.38682, 14.0123, 71.1195, 2.49899}
  	focoefs[][055] = {11.1650, 7.45636, 5.51106, 1.09496, 3.65944, 0.24397, 8.52556, 21.1647, 0.77218}
  	focoefs[][056] = {10.7806, 7.75868, 5.22746, 0.84711, 3.54770, 0.22314, 7.64468, 16.9673, 0.38604}
  	focoefs[][057] = {13.9352, 5.84833, 4.64221, 1.44753, 3.97779, 0.44555, 13.3971, 74.1605, 3.11686}
  	focoefs[][058] = {12.4655, 6.63111, 5.76679, 1.34230, 3.54270, 0.28920, 9.31140, 26.9799, 1.79285}
  	focoefs[][059] = {11.8168, 7.11181, 5.78135, 1.14523, 3.37484, 0.24408, 7.98760, 19.8970, 1.14431}
  	focoefs[][060] = {14.6744, 5.62816, 3.92540, 2.16398, 3.71486, 0.50033, 12.8862, 65.4071, 3.59838}
  	focoefs[][061] = {12.5225, 6.68507, 5.98382, 1.17317, 3.13961, 0.25431, 7.55544, 18.8543, 1.63497}
  	focoefs[][062] = {15.3412, 5.74150, 3.10733, 2.52764, 3.63868, 0.65640, 16.0719, 70.7609, 4.26842}
  	focoefs[][063] = {12.6920, 6.69883, 6.06692, 1.00660, 2.80262, 0.22789, 6.36441, 14.4122, 1.53545}
  	focoefs[][064] = {15.4378, 6.00432, 3.05158, 2.93572, 3.39715, 0.73097, 18.9533, 63.7969, 4.56068}
  	focoefs[][065] = {15.4043, 6.13723, 3.74679, 3.01390, 3.07517, 0.74113, 21.0014, 57.7446, 4.69149}
  	focoefs[][066] = {15.5372, 5.98288, 4.83996, 2.93549, 2.71530, 0.68962, 21.0079, 52.4308, 4.70026}
  	focoefs[][067] = {15.9934, 6.02439, 5.51599, 2.88716, 2.35651, 19.7393, 0.58143, 47.3323, 4.57602}
  	focoefs[][068] = {15.4080, 6.78083, 6.00715, 2.99332, 2.43532, 22.0832, 0.68621, 64.9193, 4.80234}
  	focoefs[][069] = {16.8494, 7.19790, 4.92564, 2.91606, 2.01856, 18.0409, 0.39741, 42.5054, 4.10864}
  	focoefs[][070] = {11.4809, 9.46904, 9.16981, 1.42608, 1.08140, 18.2800, 2.38825, 185.293, 5.43921}
  	focoefs[][071] = {17.8943, 8.59341, 7.91428, 2.47499, 1.71750, 0.09258, 15.4484, 32.5110,-.087756}
  	focoefs[][072] = {11.6164, 9.73009, 8.68081, 2.60986, 1.85574, 14.6109, 0.89852, 139.830, 5.34841}
  	focoefs[][073] = {18.2430, 8.90811, 1.69192,-32.1118, 1.51215, 13.6536, 27.8238,-0.01488, 39.2691}
  	focoefs[][074] = {19.0567, 6.50783, 4.81524, 2.84786, 1.24615, 9.68019, 18.8903, 121.353, 5.76121}
  	focoefs[][075] = {18.4207, 9.75213, 1.05270,-33.4755, 1.34457, 12.0631, 25.1684,-0.01023, 40.2513}
  	focoefs[][076] = {19.2273, 10.1378, 2.48177, 2.42892, 1.15488, 10.7877, 120.126, 33.3722, 5.71886}
  	focoefs[][077] = {19.1301, 10.1098, 0.98896,-0.00004, 1.16051, 10.4084, 20.7214,-3.20442, 5.77164}
  	focoefs[][078] = {19.3496, 10.8737, 3.47687, 1.64516, 1.06626, 10.5977, 32.6174, 120.397, 5.65073}
  	focoefs[][079] = {19.1248, 18.2989, 11.0121, 2.04325, 1.07235, 0.00315, 10.3385, 25.9292,-12.4799}
  	focoefs[][080] = {19.0175, 10.7591, 1.09900, 0.48469, 1.06028, 9.36239, 0.03765, 20.9764, 4.64045}
  	focoefs[][081] = {19.3885, 11.8308, 3.75919, 1.46772, 0.97877, 10.0885, 31.9738, 117.932, 5.55047}
  	focoefs[][082] = {19.6761, 18.0893, 11.7086, 2.50624, 0.95118,-0.00669, 9.61097, 24.0356,-12.9813}
  	focoefs[][083] = {19.6054, 17.9292, 11.3451, 1.04247, 0.94029,-0.00795, 8.76715, 19.3690,-12.9217}
  	focoefs[][084] = {19.4800, 17.6328, 11.0940, 0.37154, 0.94043,-0.00723, 8.29745, 18.9700,-12.5778}
  	focoefs[][085] = {19.3597, 12.8087, 3.41372, 1.99926, 0.89356, 9.27497, 32.3513, 107.406, 5.41556}
  	focoefs[][086] = {19.4316, 13.7309, 4.26537, 1.28720, 0.82092, 8.97737, 28.2621, 111.501, 5.28192}
  	focoefs[][087] = {20.8024, 13.2995, 3.27542, 2.21026, 0.74711, 8.36626, 20.6179,-0.14664, 1.41087}
  	focoefs[][088] = {41.5821, 12.9936, 2.71276,-24.2593, 0.61466, 7.99801, 18.1564, 0.43857, 6.97025}
  	focoefs[][089] = {19.4524, 14.6845, 4.50240, 1.24740, 0.75019, 8.42622, 26.1564, 107.780, 5.11007}
  	focoefs[][090] = {25.0958, 14.1510, 3.64428,-12.5768, 0.61346, 7.80244, 19.0932, 0.13532, 11.6838}
  	focoefs[][091] = {41.5236, 13.8272, 3.07969,-25.9694, 0.52905, 7.49419, 16.9498, 0.32686, 8.53824}
  	focoefs[][092] = {19.5123, 15.3800, 5.38330, 0.81015, 0.68583, 7.95714, 23.1808, 65.9295, 4.91427}
  	focoefs[][093] = {19.4652, 15.5805, 4.04748, 0.02216, 0.68159, 7.80880, 20.9573, 110.020, 4.88510}
  	focoefs[][094] = {51.1288, 14.6979, 3.41607,-38.2678, 0.43734, 7.03139, 15.8623, 0.26589, 11.0241}
  	focoefs[][095] = {19.5284, 16.5811, 4.99150, 1.21404, 0.62387, 7.39504, 22.2282, 100.226, 4.68114}
  	focoefs[][096] = {19.5416, 16.4239, 5.12995, 0.24053, 0.62273, 7.39663, 20.5530, 59.0604, 4.66470}
  	focoefs[][097] = {19.5152, 16.4852, 4.32525, 0.02777, 0.62050, 73.0347, 19.3673, 92.9184, 4.64695}
  	focoefs[][098] = {19.5528, 17.5717, 4.47374, 1.98562, 0.56604, 6.79630, 21.2907, 85.2777, 4.41158}
  	focoefs[][099] = {19.5901, 17.3740, 4.62594, 0.03770, 0.56389, 6.83082, 17.8856, 76.2909, 4.37269}
  	focoefs[][100] = {19.5872, 18.7169, 4.02722, 2.51452, 0.51510, 6.29430, 22.7308, 88.5675, 4.14542}
  	focoefs[][101] = {19.6698, 18.1942, 4.09851, 0.00365, 0.50926, 6.28098, 15.4189, 160.227, 4.03396}
  	focoefs[][102] = {19.6527, 19.5108, 3.86895, 3.14764, 0.46604, 5.76321, 24.0627, 78.1533, 3.81227}
  	focoefs[][103] = {19.7166, 18.9265, 3.79775, 1.86248, 0.46027, 5.66448, 17.7248, 42.8086, 3.69648}
  	focoefs[][104] = {19.7914, 18.9162, 3.64761, 0.00000, 0.45879, 5.76682, 13.3733, 0.00000, 3.64494}
  	focoefs[][105] = {20.0755, 19.7766, 4.30389, 3.44952, 5.24328, 0.41858, 26.0178, 70.1646, 3.38881}
  	focoefs[][106] = {19.8617, 19.5199, 3.73465, 1.61027, 0.41409, 5.18292, 16.8529, 35.1406, 3.27356}
  	focoefs[][107] = {19.9613, 19.5889, 3.24333, 0.00000, 0.41262, 5.30028, 11.7603, 0.00000, 3.20701}
  	focoefs[][108] = {20.4608, 20.0336, 5.38664, 3.33079, 4.74225, 0.37041, 27.3458, 65.0573, 2.78462}
  	focoefs[][109] = {20.7492, 20.5640, 6.86158, 2.97589, 4.27091, 0.31960, 27.3186, 61.5375, 1.84739}
  	focoefs[][110] = {20.8307, 20.4454, 7.52618, 3.18616, 4.29514, 0.32402, 29.8990, 81.4344, 2.00513}
  	focoefs[][111] = {21.6679, 21.0085, 8.43382, 2.62265, 0.26422, 3.83526, 26.2297, 58.4830, 0.26635}
  	focoefs[][112] = {22.3163, 21.1792, 10.7382, 1.46163, 0.23092, 3.49464, 25.1864, 232.829,-0.70709}
  	focoefs[][113] = {23.9649, 21.2204, 9.76727, 1.61550, 0.20448, 3.43876, 23.4941, 49.7057,-2.56728}
  	focoefs[][114] = {27.7489, 21.3777, 11.0400, 2.68186, 0.15152, 3.09817, 20.6774, 178.819,-6.85854}
  	focoefs[][115] = {29.2996, 21.4669, 10.9209, 0.80126, 0.14047, 3.08785, 20.8818, 46.8842,-8.48753}
  	focoefs[][116] = {33.2109, 21.7181, 11.6222, 3.17239, 0.11040, 2.83641, 19.3886, 144.438,-12.7404}
  	focoefs[][117] = {43.6346, 21.7192, 11.7264, 0.32945, 0.07854, 2.78360, 18.4930, 49.2222,-23.4085}
  	focoefs[][118] = {29.4100, 22.2428, 11.9818, 3.19259, 0.12335, 2.74837, 18.3794, 139.603,-8.84560}
  	focoefs[][119] = {49.1105, 22.3499, 11.8399, 0.67455, 0.06535, 2.67229, 17.2040, 38.1904,-28.9739}
  	focoefs[][120] = {66.7693, 21.8563, 12.2486, 0.09617, 0.04464, 2.53711, 16.4477, 64.4675,-46.9691}
  	focoefs[][121] = {22.9220, 22.2518, 12.2269, 2.72431, 2.78604, 0.18015, 17.6663, 160.915,-1.13930}
  	focoefs[][122] = {49.4655, 22.9705, 11.8015, 1.12179, 0.06197, 2.57634, 16.0371, 32.3673,-29.3586}
  	focoefs[][123] = {62.6752, 22.4952, 12.4946, 0.20294, 0.04586, 2.45900, 15.5713, 46.5889,-42.8667}
  	focoefs[][124] = {23.4069, 19.7073, 12.5016, 2.72850, 2.71587, 0.20950, 16.9122, 156.556, 1.64038}
  	focoefs[][125] = {49.4292, 23.6116, 11.6190, 1.68986, 0.05936, 2.48611, 14.9366, 28.4515,-29.3493}
  	focoefs[][126] = {23.8480, 17.5535, 12.7324, 2.72975, 2.65746, 0.24780, 16.2463, 152.682, 4.12018}
  	focoefs[][127] = {49.2699, 24.2700, 11.3481, 2.32869, 0.05709, 2.40059, 13.9124, 25.6906,-29.2165}
  	focoefs[][128] = {24.2242, 15.9132, 12.9238, 2.72836, 2.60993, 0.29475, 15.6554, 149.221, 6.19355}
  	focoefs[][129] = {36.3271, 24.8202, 11.3426, 2.62300, 0.07823, 2.33602, 13.1872, 24.3996,-16.1429}
  	focoefs[][130] = {24.5148, 14.8058, 13.0799, 2.72477, 2.57225, 0.34930, 15.1280, 146.103, 7.85731}
  	focoefs[][131] = {25.6516, 23.9387, 10.5738, 4.05853, 2.36073, 0.13260, 12.6495, 25.0026,-3.22358}
  	focoefs[][132] = {33.2862, 25.5041, 11.1494, 3.13496, 0.08350, 2.26275, 12.3883, 22.8351,-13.0748}
  	focoefs[][133] = {24.4004, 14.0308, 13.1754, 3.24472, 2.47491, 0.40238, 14.4670, 119.738, 9.12488}
  	focoefs[][134] = {29.0290, 26.1387, 11.0510, 3.52244, 0.09521, 2.19696, 11.7141, 21.6929,-8.74150}
  	focoefs[][135] = {24.3736, 13.8649, 13.2510, 3.24435, 2.46637, 0.47517, 14.0424, 117.446, 10.2420}
  	focoefs[][136] = {26.7821, 25.9463, 10.9724, 3.88172, 2.13333, 0.10597, 11.0974, 20.7042,-5.58307}
  	focoefs[][137] = {24.6193, 14.2735, 13.3567, 2.70316, 2.52208, 0.54556, 13.8487, 138.385, 11.0290}
  	focoefs[][138] = {27.3805, 22.2062, 10.9975, 4.10030, 2.07832, 0.12643, 10.5960, 19.9671,-1.68516}
  	focoefs[][139] = {24.3162, 14.9012, 13.3895, 2.69309, 2.52724, 0.61572, 13.5041, 136.246, 11.6817}
  	focoefs[][140] = {27.9956, 19.9560, 11.0106, 4.33205, 2.02324, 0.14275, 10.1165, 19.2589, 0.70499}
  	focoefs[][141] = {23.8201, 15.8796, 13.3938, 2.68190, 2.54419, 0.68445, 13.1932, 134.282, 12.2062}
  	focoefs[][142] = {28.5315, 17.4316, 11.1113, 4.43156, 1.97796, 0.17182, 9.73821, 18.7294, 3.49325}
  	focoefs[][143] = {23.1386, 17.1707, 13.3703, 2.66981, 2.57320, 0.74948, 12.9126, 132.468, 12.6322}
  	focoefs[][144] = {29.0215, 15.6168, 11.2288, 4.49403, 1.93707, 0.20467, 9.40342, 18.2607, 5.63812}
  	focoefs[][145] = {22.3028, 18.7202, 13.3200, 2.65701, 2.61393, 0.80868, 12.6590, 130.783, 12.9818}
  	focoefs[][146] = {29.1313, 13.5855, 11.4132, 4.69659, 1.99979, 0.32335, 9.59277, 20.3507, 9.17182}
  	focoefs[][147] = {29.4761, 14.4357, 11.3446, 4.54681, 1.89879, 0.23793, 9.09408, 17.8206, 7.19600}
  	focoefs[][148] = {21.1866, 20.1760, 13.0532, 3.21190, 0.88654, 2.68610, 12.2746, 107.128, 13.3489}
  	focoefs[][149] = {29.8480, 13.6268, 11.4750, 4.56009, 1.86596, 0.27623, 8.82479, 17.4364, 8.48923}
  	focoefs[][150] = {24.6725, 17.2295, 12.8069, 3.55970, 0.97400, 2.89038, 12.2897, 93.4381, 13.7049}
  	focoefs[][151] = {28.1757, 14.4288, 12.6412, 3.74436, 1.04034, 3.20784, 12.5054, 85.0183, 13.9824}
  	focoefs[][152] = {31.0935, 12.5273, 12.3769, 3.79138, 1.07885, 12.8331, 3.63298, 79.7647, 14.1842}
  	focoefs[][153] = {33.2961, 12.3497, 11.2819, 3.72367, 1.09315, 13.2559, 4.16736, 76.6562, 14.3239}
  	focoefs[][154] = {34.8667, 11.9524, 11.1851, 3.56436, 1.08840, 13.8042, 4.79179, 75.1399, 14.4097}
  	focoefs[][155] = {35.9454, 11.9980, 11.2501, 3.34312, 1.06924, 5.43443, 14.4983, 74.7918, 14.4449}
  	focoefs[][156] = {36.8102, 13.0747, 11.3323, 2.31421, 1.04422, 6.07340, 15.7018, 73.8375, 14.4526}
  	focoefs[][157] = {37.3027, 14.9306, 10.3425, 2.01229, 1.00810, 6.52550, 16.5100, 76.9117, 14.3992}
  	focoefs[][158] = {37.5186, 17.0353, 8.51121, 2.63340, 0.96455, 6.65786, 16.8438, 76.7228, 14.2911}
  	focoefs[][159] = {37.6947, 19.7195, 6.38290, 3.00960, 0.92263, 6.78248, 19.2435, 85.9267, 14.1800}
  	focoefs[][160] = {37.7383, 21.3394, 5.17527, 3.71604, 0.87755, 6.58964, 21.2437, 78.8094, 14.0203}
  	focoefs[][161] = {37.7143, 22.4542, 4.84549, 4.14816, 0.83222, 6.27051, 24.4693, 72.1558, 13.8301}
  	focoefs[][162] = {37.6297, 23.1323, 5.59203, 4.04218, 0.78640, 5.86644, 27.8678, 68.1617, 13.5991}
  	focoefs[][163] = {37.4971, 23.5635, 7.15953, 3.45924, 0.74012, 5.42694, 29.8350, 66.3564, 13.3183}
  	focoefs[][164] = {37.3308, 23.8933, 9.02222, 2.77349, 0.69354, 4.98696, 30.0338, 65.5799, 12.9796}
  	focoefs[][165] = {37.1902, 24.1306, 11.5026, 1.47980, 0.65303, 4.61305, 29.2597, 257.965, 12.6868}
  	focoefs[][166] = {36.9820, 24.2495, 11.8719, 2.72428, 0.60394, 4.17857, 24.3782, 200.024, 12.1642}
  	focoefs[][167] = {36.8705, 24.7131, 12.3889, 3.26501, 0.56458, 3.88776, 23.1506, 161.726, 11.7484}
  	focoefs[][168] = {36.7754, 25.2506, 13.0681, 3.63791, 0.52510, 3.61658, 22.3410, 139.164, 11.2497}
  	focoefs[][169] = {37.1457, 25.2998, 13.7846, 3.29611, 0.52020, 3.66300, 20.6539, 150.973, 11.4561}
  	focoefs[][170] = {37.2808, 25.6563, 14.3501, 3.30732, 0.50239, 3.58562, 19.6342, 146.633, 11.3864}
  	focoefs[][171] = {37.3968, 26.0671, 14.8366, 3.31586, 0.48676, 3.52325, 18.7419, 142.798, 11.3632}
  	focoefs[][172] = {37.6407, 26.5603, 15.4492, 2.79814, 0.47976, 3.57178, 17.9814, 165.232, 11.5358}
  	focoefs[][173] = {37.6909, 27.1436, 15.7842, 2.79600, 0.46617, 3.52195, 17.3069, 161.931, 11.5685}
  	focoefs[][174] = {37.5543, 27.6657, 15.8858, 3.32758, 0.44932, 3.38713, 16.6498, 133.547, 11.5431}
  	focoefs[][175] = {37.5273, 28.3202, 16.1181, 3.32793, 0.43930, 3.35014, 16.1000, 131.027, 11.6823}
  	focoefs[][176] = {37.6111, 29.2465, 16.4566, 2.78216, 0.43255, 3.39285, 15.6791, 153.766, 11.8853}
  	focoefs[][177] = {37.4979, 30.0495, 16.5881, 2.77596, 0.42353, 3.35234, 15.2381, 151.474, 12.0698}
  	focoefs[][178] = {37.3380, 30.8936, 16.6818, 2.76929, 0.41562, 3.31193, 14.8362, 149.344, 12.2983}
  	focoefs[][179] = {37.1301, 31.7721, 16.7422, 2.76232, 0.40883, 3.27132, 14.4683, 147.353, 12.5741}
  	focoefs[][180] = {36.8731, 32.6784, 16.7732, 2.75513, 0.40324, 3.23045, 14.1302, 145.481, 12.9008}
  	focoefs[][181] = {36.3813, 33.1999, 16.6469, 3.31406, 0.40165, 3.13608, 13.7255, 119.377, 13.4313}

	Variable Nbuf=8
	Make/N=(Nbuf)/O root:Packages:Xray:foLast, root:Packages:Xray:KLast	// used to buffer calculation
	Make/N=(Nbuf)/T/O root:Packages:Xray:symbLast
	if (exists("root:Packages:Xray:nLast")!=2)
		Variable/G root:Packages:Xray:nLast=Nbuf-1	// number to fill next
	endif
End