Changeset 1952 for trunk/GSASIImath.py


Ignore:
Timestamp:
Aug 7, 2015 1:37:20 PM (7 years ago)
Author:
vondreele
Message:

Add modulation of thermal parameters to drawing

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIImath.py

    r1951 r1952  
    994994    for atom in atoms:   
    995995        atxyz = np.array(atom[cx:cx+3])
     996        atuij = np.array(atom[cia+2:cia+8])
    996997        waveType = atom[-1]['SS1']['waveType']
    997998        Spos = atom[-1]['SS1']['Spos']
     999        Sadp = atom[-1]['SS1']['Sadp']
    9981000        wave = np.zeros(3)
     1001        uwave = np.zeros(6)
    9991002        if len(Spos):
    10001003            scof = []
     
    10121015                    ccof.append(spos[0][3:])
    10131016            wave += np.sum(posFourier(tau,np.array(scof),np.array(ccof)),axis=1)
     1017        if len(Sadp):
     1018            scof = []
     1019            ccof = []
     1020            for i,sadp in enumerate(Sadp):
     1021                scof.append(sadp[0][:6])
     1022                ccof.append(sadp[0][6:])
     1023            uwave += np.sum(posFourier(tau,np.array(scof),np.array(ccof)),axis=1)
    10141024        indx = FindAtomIndexByIDs(drawAtoms,dci,[atom[cia+8],],True)
    10151025        for ind in indx:
     
    10171027            opr = drawatom[dcs-1]
    10181028            if atom[cia] == 'A':                   
    1019                 X,U = G2spc.ApplyStringOps(opr,SGData,atxyz+wave,atom[cia+2:cia+8])
     1029                X,U = G2spc.ApplyStringOps(opr,SGData,atxyz+wave,atuij+uwave)
    10201030                drawatom[dcx:dcx+3] = X
    1021 #                drawatom[dci-6:dci] = U
     1031                drawatom[dci-6:dci] = U
    10221032            else:
    10231033                X = G2spc.ApplyStringOps(opr,SGData,atxyz+wave)
Note: See TracChangeset for help on using the changeset viewer.