-- crystal field operator in Oh 
-- symmetry acting on a d shell
NF = 10
NB = 0
IndexDn_3d = {0,2,4,6,8}
IndexUp_3d = {1,3,5,7,9}
tenDq = 1.1
Akm = {{4, 0,(21/10)*tenDq},
       {4,-4,(3/2)*math.sqrt(7/10)*tenDq},
       {4, 4,(3/2)*math.sqrt(7/10)*tenDq}}
OppCF = NewOperator("CF", NF, IndexUp_3d, IndexDn_3d, Akm)