Light shift per mk, extracted from MK mathematica

Ground state state for 850nm 

format as
sigma-
lin
sigma+

MF = ...-2 -1 0 1...

## S level ##
F=1
-1.99365*10^7	-2.0827*10^7	-2.17176*10^7
-2.0827*10^7	-2.0827*10^7	-2.0827*10^7
-2.17176*10^7	-2.0827*10^7	-1.99365*10^7

F=2
-2.26082*10^7	-2.17176*10^7	-2.0827*10^7	-1.99365*10^7	-1.90459*10^7
-2.0827*10^7	-2.0827*10^7	-2.0827*10^7	-2.0827*10^7	-2.0827*10^7
-1.90459*10^7	-1.99365*10^7	-2.0827*10^7	-2.17176*10^7	-2.26082*10^7

## P level ##
F'=0
9.52887*10^6
9.52887*10^6
9.52887*10^6

F'=1
4.70379*10^6	1.25405*10^7	1.13423*10^7
1.25405*10^7	3.50556*10^6	1.25405*10^7
1.13423*10^7	1.25405*10^7	4.70379*10^6

F'=2
2.89036*10^6	6.20961*10^6	9.52887*10^6	1.28481*10^7	1.61674*10^7
9.52887*10^6	9.52887*10^6	9.52887*10^6	9.52887*10^6	9.52887*10^6
1.61674*10^7	1.28481*10^7	9.52887*10^6	6.20961*10^6	2.89036*10^6

F'=3
3.33567*10^6	2.89036*10^6	3.95087*10^6	6.51721*10^6	1.05894*10^7	1.61674*10^7	2.32512*10^7
1.99973*10^6	9.52887*10^6	1.40464*10^7	1.55522*10^7	1.40464*10^7	9.52887*10^6	1.99973*10^6
2.32512*10^7	1.61674*10^7	1.05894*10^7	6.51721*10^6	3.95087*10^6	2.89036*10^6	3.33567*10^6

## Calculation of S level ##

P levels used
If[n=={5,1,1/2},v={2*Pi*2.998*10^8/(7949.8*10^-10),3.648*10^7},
If[n=={5,1,3/2},v={2*Pi*2.998*10^8/(7802.4*10^-10),3.843*10^7},
If[n=={6,1,1/2},v={2*Pi*2.998*10^8/(4215.524*10^-10),1.50*10^6},
If[n=={6,1,3/2},v={2*Pi*2.998*10^8/(4201.792*10^-10),1.77*10^6},
If[n=={7,1,1/2},v={2*Pi*2.998*10^8/(3591.572*10^-10),2.89*10^5},
If[n=={7,1,3/2},v={2*Pi*2.998*10^8/(3587.050*10^-10),3.96*10^5}]]]]]];
v]

v=shiftSJJ[5,1,1/2,J1,F1,mF1,q,I1]+
shiftSJJ[5,1,3/2,J1,F1,mF1,q,I1]+
shiftSJJ[6,1,1/2,J1,F1,mF1,q,I1]+
shiftSJJ[6,1,3/2,J1,F1,mF1,q,I1]

## Calcultion of P level ##
If[n=={5,0,1/2},v={2*Pi*2.998*10^8/(-7802.4*10^-10),3.843*10^7},
If[n=={6,0,1/2},v={2*Pi*2.998*10^8/(13668.7*10^-10),1.451*10^7},
If[n=={7,0,1/2},v={2*Pi*2.998*10^8/(7410.2*10^-10),4.536*10^6},
If[n=={8,0,1/2},v={2*Pi*2.998*10^8/(6161.3*10^-10),2.169*10^6},
If[n=={4,2,3/2},v={2*Pi*2.998*10^8/(15292.6*10^-10),1.864*10^6},
If[n=={4,2,5/2},v={2*Pi*2.998*10^8/(15293.7*10^-10),1.119*10^7},
If[n=={5,2,3/2},v={2*Pi*2.998*10^8/(7759.436*10^-10),4.161*10^5},
If[n=={5,2,5/2},v={2*Pi*2.998*10^8/(7757.651*10^-10),4.318*10^6},
If[n=={6,2,3/2},v={2*Pi*2.998*10^8/(6299.225*10^-10),3.876*10^5},
If[n=={6,2,5/2},v={2*Pi*2.998*10^8/(6298.325*10^-10),4.042*10^6}

v=shiftJJP3[5,0,1/2,J1,F1,mF1,q,I1]+
shiftJJP3[6,0,1/2,J1,F1,mF1,q,I1]+
shiftJJP3[7,0,1/2,J1,F1,mF1,q,I1]+
shiftJJP3[4,2,3/2,J1,F1,mF1,q,I1]+
shiftJJP3[4,2,5/2,J1,F1,mF1,q,I1]+
shiftJJP3[5,2,3/2,J1,F1,mF1,q,I1]+
shiftJJP3[5,2,5/2,J1,F1,mF1,q,I1]