Changeset 3295 for trunk/GSASIIphsGUI.py


Ignore:
Timestamp:
Feb 22, 2018 11:01:07 AM (4 years ago)
Author:
vondreele
Message:

remove superfluous calls to SpcGroup? after SGData has been defined
add 3rd row magnetic form factors
add contour plots for projections of mustrain, size, % 3D pole distribution
improve quality of these surfaces

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIphsGUI.py

    r3291 r3295  
    24992499                ci = colLabels.index('I/A')
    25002500                ID = atomData[r][ci+8]
     2501                SGData = generalData['SGData']
    25012502                if Atoms.GetColLabelValue(c) in ['x','y','z']:
    25022503                    ci = colLabels.index('x')
     
    25062507                    SScol = colLabels.index('site sym')
    25072508                    Mulcol = colLabels.index('mult')
    2508                     E,SGData = G2spc.SpcGroup(generalData['SGData']['SpGrp'])
    25092509                    Sytsym,Mult = G2spc.SytSym(XYZ,SGData)[:2]
    25102510                    atomData[r][SScol] = Sytsym
     
    25562556                    ci = colLabels.index('I/A')
    25572557                    atomData[r][c] = atomData[r][c].replace(rbAtmDict.get(atomData[r][ci+8],''),'')
     2558                elif Atoms.GetColLabelValue(c) in ['Mx','My','Mz']:
     2559                    value = atomData[r][c]
     2560                    cx = colLabels.index('x')
     2561                    SpnFlp = generalData['SGData']['SpnFlp']
     2562                    SytSym,Mul,Nop,dupDir = G2spc.SytSym(atomData[r][cx:cx+3],SGData)
     2563                    CSI = G2spc.GetCSpqinel(SytSym,SpnFlp,dupDir)
     2564                    iM = CSI[0][c-colLabels.index('Mx')]
     2565                    for i in range(3):
     2566                        if iM == CSI[0][i]:
     2567                            atomData[r][i+colLabels.index('Mx')] = value*CSI[1][i]
    25582568                if 'Atoms' in data['Drawing']:
    25592569                    ci = colLabels.index('I/A')
     
    27202730                    if not SGData['SGGray']:
    27212731                        CSI = G2spc.GetCSpqinel(SytSym,SpnFlp,dupDir)
    2722 #                    print (SytSym,Nop,SpnFlp[Nop],CSI,dupDir)
     2732                    print (SytSym,Nop,SpnFlp[Nop],CSI,dupDir)
    27232733                    for i in range(3):
    27242734                        ci = i+colM
     
    27282738                            Atoms.SetCellStyle(row,ci,WHITE,False)
    27292739                            Atoms.SetCellTextColour(row,ci,BLACK)
    2730                            
    27312740                if 'X' in rbExcl:
    27322741                    for c in range(0,colX+3):
     
    28372846        generalData = data['General']
    28382847        atId = ran.randint(0,sys.maxsize)
    2839         E,SGData = G2spc.SpcGroup(generalData['SGData']['SpGrp'])
     2848        SGData = generalData['SGData']
    28402849        Sytsym,Mult = G2spc.SytSym([x,y,z],SGData)[:2]
    28412850        if generalData['Type'] == 'macromolecular':
     
    30863095        atomData = data['Atoms']
    30873096        generalData = data['General']
    3088         E,SGData = G2spc.SpcGroup(generalData['SGData']['SpGrp'])
     3097        SGData = generalData['SGData']
    30893098        Sytsym,Mult = G2spc.SytSym([x,y,z],SGData)[:2]
    30903099        atId = ran.randint(0,sys.maxsize)
     
    31783187        atomData = data['Atoms']
    31793188        generalData = data['General']
     3189        ifMag = False
    31803190        colLabels = [Atoms.GetColLabelValue(c) for c in range(Atoms.GetNumberCols())]
    31813191        ci = colLabels.index('I/A')
     
    31833193        if generalData['Type'] == 'magnetic':
    31843194            choices += ['Mx','My','Mz',]
     3195            ifMag = True
    31853196        dlg = wx.SingleChoiceDialog(G2frame,'Select','Atom parameter',choices)
    31863197        parm = ''
     
    31913202        dlg.Destroy()
    31923203        if parm in ['Type']:
    3193             dlg = G2elemGUI.PickElement(G2frame)
     3204            dlg = G2elemGUI.PickElement(G2frame,ifMag=ifMag)
    31943205            if dlg.ShowModal() == wx.ID_OK:
    31953206                if dlg.Elem not in ['None']:
Note: See TracChangeset for help on using the changeset viewer.