Changeset 3266 for trunk/GSASIIstrIO.py


Ignore:
Timestamp:
Feb 8, 2018 1:31:55 PM (4 years ago)
Author:
vondreele
Message:

change modulated wave type assignments so one for each kind of wave (pos, frac, etc).
implement optional display of PNT data set positions on pole figure plots. Picker displays histo. name.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIstrIO.py

    r3255 r3266  
    987987                if len(Waves):
    988988                    pFile.write(' atom: %s, site sym: %s, type: %s wave type: %s:\n'%
    989                         (at[ct-1],at[cs],Stype,AtomSS['waveType']))
    990                 for iw,wave in enumerate(Waves):                   
     989                        (at[ct-1],at[cs],Stype,Waves[0]))
     990                else:
     991                    continue
     992                for iw,wave in enumerate(Waves[1:]):                   
    991993                    line = ''
    992                     if AtomSS['waveType'] in ['Block','ZigZag'] and Stype == 'Spos' and not iw:
     994                    if Waves[0] in ['Block','ZigZag'] and Stype == 'Spos' and not iw:
    993995                        for item in names[Stype][6:]:
    994996                            line += '%8s '%(item)                       
     
    12441246                if General.get('Modulated',False):
    12451247                    AtomSS = at[-1]['SS1']
    1246                     waveType = AtomSS['waveType']
    1247                     phaseDict[pfx+'waveType:'+str(i)] = waveType
    12481248                    for Stype in ['Sfrac','Spos','Sadp','Smag']:
    12491249                        Waves = AtomSS[Stype]
     1250                        if len(Waves):
     1251                            waveType = Waves[0]
     1252                        else:
     1253                            continue
     1254                        phaseDict[pfx+Stype[1].upper()+'waveType:'+str(i)] = waveType
    12501255                        nx = 0
    1251                         for iw,wave in enumerate(Waves):
     1256                        for iw,wave in enumerate(Waves[1:]):
    12521257                            if not iw:
    12531258                                if waveType in ['ZigZag','Block']:
    12541259                                    nx = 1
    1255                                 CSI = G2spc.GetSSfxuinel(waveType,1,at[cx:cx+3],SGData,SSGData)
     1260                                CSI = G2spc.GetSSfxuinel(waveType,Stype,1,at[cx:cx+3],SGData,SSGData)
    12561261                            else:
    1257                                 CSI = G2spc.GetSSfxuinel('Fourier',iw+1-nx,at[cx:cx+3],SGData,SSGData)
     1262                                CSI = G2spc.GetSSfxuinel('Fourier',Stype,iw+1-nx,at[cx:cx+3],SGData,SSGData)
    12581263                            uId,uCoef = CSI[Stype]
    12591264                            stiw = str(i)+':'+str(iw)
     
    16851690        for i,at in enumerate(Atoms):
    16861691            AtomSS = at[-1]['SS1']
    1687             waveType = AtomSS['waveType']
    16881692            for Stype in ['Sfrac','Spos','Sadp','Smag']:
    16891693                Waves = AtomSS[Stype]
    16901694                if len(Waves):
     1695                    waveType = Waves[0]
     1696                else:
     1697                    continue
     1698                if len(Waves):
    16911699                    pFile.write(' atom: %s, site sym: %s, type: %s wave type: %s:\n'%
    16921700                        (at[ct-1],at[cs],Stype,waveType))
    1693                     for iw,wave in enumerate(Waves):
     1701                    for iw,wave in enumerate(Waves[1:]):
    16941702                        stiw = ':'+str(i)+':'+str(iw)
    16951703                        namstr = '  names :'
     
    20222030                if General.get('Modulated',False):
    20232031                    AtomSS = at[-1]['SS1']
    2024                     waveType = AtomSS['waveType']
    20252032                    for Stype in ['Sfrac','Spos','Sadp','Smag']:
    20262033                        Waves = AtomSS[Stype]
    2027                         for iw,wave in enumerate(Waves):
     2034                        if len(Waves):
     2035                            waveType = Waves[0]
     2036                        else:
     2037                            continue
     2038                        for iw,wave in enumerate(Waves[1:]):
    20282039                            stiw = str(i)+':'+str(iw)
    20292040                            if Stype == 'Spos':
     
    20472058                                    'MXcos:'+stiw,'MYcos:'+stiw,'MZcos:'+stiw]
    20482059                            for iname,name in enumerate(names):
    2049                                 AtomSS[Stype][iw][0][iname] = parmDict[pfx+name]
     2060                                AtomSS[Stype][iw+1][0][iname] = parmDict[pfx+name]
    20502061                                if pfx+name in sigDict:
    20512062                                    wavesSig[name] = sigDict[pfx+name]
Note: See TracChangeset for help on using the changeset viewer.