Changeset 3888 for trunk/GSASIIlattice.py

Timestamp:

Apr 11, 2019 3:59:48 PM (4 years ago)

Author:

vondreele

Message:

provide contour plot (not finished yet) on structure drawing

File:

• trunk/GSASIIlattice.py (modified) (3 diffs)

trunk/GSASIIlattice.py

@@ -891 +891 @@
     return rdsq
     
-def PlaneIntercepts(Amat,Bmat,H,phase,stack):
+def PlaneIntercepts(Amat,H,phase,stack):
     ''' find unit cell intercepts for a stack of hkl planes
     '''
@@ -906 +906 @@
                for j in [0,1,2,3]:
                     hx = [0,0,0]
-                    intcpt = (phase/360.+step-H[h]*Ux[j,0]-H[k]*Ux[j,1])/H[l]
+                    intcpt = ((phase)/360.+step-H[h]*Ux[j,0]-H[k]*Ux[j,1])/H[l]
                     if 0. <= intcpt <= 1.:                        
                         hx[h] = Ux[j,0]
@@ -924 +924 @@
             A[1:] = [np.dot(DX[1],dx) for dx in DX[1:]]
             HX = HX[np.argsort(A)]
-#            GSASIIpath.IPyBreak()
             Stack.append(HX)
     return Stack