Changeset 2474 for trunk/GSASIIstrIO.py


Ignore:
Timestamp:
Sep 20, 2016 3:35:05 PM (5 years ago)
Author:
vondreele
Message:

fix spin inversion plot & some refine list issues

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIstrIO.py

    r2473 r2474  
    938938        cx,ct,cs,cia = General['AtomPtrs']
    939939        cmx = cx+4
     940        AtInfo = dict(zip(General['AtomTypes'],General['Lande g']))
    940941        print >>pFile,'\n Magnetic moments:'
    941942        line = '   name    type  refine?  Mx        My        Mz    '
     
    943944        print >>pFile,135*'-'
    944945        for i,at in enumerate(Atoms):
    945             line = '%7s'%(at[ct-1])+'%7s'%(at[ct])+'%7s'%(at[ct+1])+'%10.5f'%(at[cmx])+'%10.5f'%(at[cmx+1])+ \
    946                 '%10.5f'%(at[cmx+2])
    947             print >>pFile,line
     946            if AtInfo[at[ct]]:
     947                line = '%7s'%(at[ct-1])+'%7s'%(at[ct])+'%7s'%(at[ct+1])+'%10.5f'%(at[cmx])+'%10.5f'%(at[cmx+1])+ \
     948                    '%10.5f'%(at[cmx+2])
     949                print >>pFile,line
    948950       
    949951               
     
    16281630        cx,ct,cs,cia = General['AtomPtrs']
    16291631        cmx = cx+4
     1632        AtInfo = dict(zip(General['AtomTypes'],General['Lande g']))
    16301633        print >>pFile,line
    16311634        print >>pFile,135*'-'
     
    16331636        noFXsig = {cmx:[10*' ','%10s'],cmx+1:[10*' ','%10s'],cmx+2:[10*' ','%10s']}
    16341637        for i,at in enumerate(Atoms):
    1635             name = fmt[0]%(at[ct-1])+fmt[1]%(at[ct])+':'
    1636             valstr = ' values:'
    1637             sigstr = ' sig   :'
    1638             for ind in range(cmx,cmx+3):
    1639                 sigind = str(i)+':'+str(ind)
    1640                 valstr += fmt[ind]%(at[ind])                   
    1641                 if sigind in atomsSig:
    1642                     sigstr += fmt[ind]%(atomsSig[sigind])
    1643                 else:
    1644                     sigstr += noFXsig[ind][1]%(noFXsig[ind][0])
    1645             print >>pFile,name
    1646             print >>pFile,valstr
    1647             print >>pFile,sigstr
     1638            if AtInfo[at[ct]]:
     1639                name = fmt[0]%(at[ct-1])+fmt[1]%(at[ct])+':'
     1640                valstr = ' values:'
     1641                sigstr = ' sig   :'
     1642                for ind in range(cmx,cmx+3):
     1643                    sigind = str(i)+':'+str(ind)
     1644                    valstr += fmt[ind]%(at[ind])                   
     1645                    if sigind in atomsSig:
     1646                        sigstr += fmt[ind]%(atomsSig[sigind])
     1647                    else:
     1648                        sigstr += noFXsig[ind][1]%(noFXsig[ind][0])
     1649                print >>pFile,name
     1650                print >>pFile,valstr
     1651                print >>pFile,sigstr
    16481652           
    16491653    def PrintWavesAndSig(General,Atoms,wavesSig):
Note: See TracChangeset for help on using the changeset viewer.