Ignore:
Timestamp:
Apr 11, 2019 3:59:48 PM (3 years ago)
Author:
vondreele
Message:

provide contour plot (not finished yet) on structure drawing

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIlattice.py

    r3801 r3888  
    891891    return rdsq
    892892   
    893 def PlaneIntercepts(Amat,Bmat,H,phase,stack):
     893def PlaneIntercepts(Amat,H,phase,stack):
    894894    ''' find unit cell intercepts for a stack of hkl planes
    895895    '''
     
    906906               for j in [0,1,2,3]:
    907907                    hx = [0,0,0]
    908                     intcpt = (phase/360.+step-H[h]*Ux[j,0]-H[k]*Ux[j,1])/H[l]
     908                    intcpt = ((phase)/360.+step-H[h]*Ux[j,0]-H[k]*Ux[j,1])/H[l]
    909909                    if 0. <= intcpt <= 1.:                       
    910910                        hx[h] = Ux[j,0]
     
    924924            A[1:] = [np.dot(DX[1],dx) for dx in DX[1:]]
    925925            HX = HX[np.argsort(A)]
    926 #            GSASIIpath.IPyBreak()
    927926            Stack.append(HX)
    928927    return Stack
Note: See TracChangeset for help on using the changeset viewer.