Skip to content

Commit 9b57b71

Browse files
committed
New SC delt-F 3D plotting: Size of spot proportional to Fo^2 & color proportional to Fo^2-Fc^2.
1 parent 56ee7e5 commit 9b57b71

1 file changed

Lines changed: 9 additions & 4 deletions

File tree

GSASII/GSASIIplot.py

Lines changed: 9 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -1526,6 +1526,8 @@ def OnKey(event): #on key UP!!
15261526
[uBox[2],uBox[3]],[uBox[1],uBox[5]],[uBox[2],uBox[6]],[uBox[3],uBox[7]],
15271527
[uBox[4],uBox[5]],[uBox[5],uBox[6]],[uBox[6],uBox[7]],[uBox[7],uBox[4]]])
15281528
uColors = [Rd,Gr,Bl, Wt,Wt,Wt, Wt,Wt,Wt, Wt,Wt,Wt]
1529+
dFmin = np.min(hklRef.T[8+Super]-hklRef.T[9+Super])
1530+
dFmax = np.max(hklRef.T[8+Super]-hklRef.T[9+Super])
15291531

15301532
def FillHKLRC():
15311533
global HKL
@@ -1598,11 +1600,14 @@ def FillHKLRC():
15981600
elif Data['Type'] == 'dFsq':
15991601
dF = Fosq-Fcsq
16001602
if dF > 0:
1601-
R.append(dF)
1602-
C.append(Gr)
1603+
R.append(Fosq)
1604+
dw = int(255.*(1.0-(dF/dFmin)))
1605+
color = np.array([dw,255,0])
16031606
else:
1604-
R.append(-dF)
1605-
C.append(Rd)
1607+
R.append(Fosq)
1608+
dw = int(255.*(1.0+(dF/dFmax)))
1609+
color = np.array([255,dw,0])
1610+
C.append(color)
16061611
if len(R):
16071612
R = np.array(R)
16081613
if Data['Type'] == 'dFsq/sig':

0 commit comments

Comments
 (0)