Changeset 793


Ignore:
Timestamp:
Oct 24, 2012 9:38:52 AM (10 years ago)
Author:
vondreele
Message:

more iparm fixes

Location:
trunk
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASII.py

    r792 r793  
    557557            S = File.readline()               
    558558        File.close()
    559         return dict(zip(newItems,zip(newVals,newVals,len(newVals)*[False,])))
     559        inst = dict(zip(newItems,zip(newVals,newVals,len(newVals)*[False,])))
     560        for item in inst:
     561            inst[item] = list(inst[item])
     562        return inst
    560563       
    561564    def ReadPowderIparm(self,instfile,bank,databanks,rd):
     
    655658                    data.extend([0.0,0.0,0.002,azm])                                      #OK defaults if fxn #3 not 1st in iprm file
    656659                codes.extend([0,0,0,0,0,0,0])
    657                 return dict(zip(names,zip(data,data,codes)))
     660                inst = dict(zip(names,zip(data,data,codes)))
     661                for item in inst:
     662                    inst[item] = list(inst[item])
     663                return inst
    658664            elif 'T' in DataType:
    659665                names = ['Type','2-theta','difC','difA','Zero','alpha','beta-0','beta-1','var-inst','X','Y','Azimuth']
     
    682688                        s = Iparm['INS  1PRCF12'].split()
    683689                        data.extend([G2IO.sfloat(s[0]),0.0,0.0,azm])
    684                 return dict(zip(names,zip(data,data,codes)))
     690                inst = dict(zip(names,zip(data,data,codes)))
     691                for item in inst:
     692                    inst[item] = list(inst[item])
     693                return inst
    685694
    686695        # stuff we might need from the reader
     
    11191128                data = ['PKS',Cuka,0.0]
    11201129                names = ['Type','Lam','Zero']
    1121                 codes = [0,0]
    1122                 self.PatternTree.SetItemPyData(self.PatternTree.AppendItem(Id,text='Instrument Parameters'),dict(zip(names,zip(data,data,codes))))
     1130                codes = [0,0,0]
     1131                inst = dict(zip(names,zip(data,data,codes)))
     1132                for item in inst:
     1133                    inst[item] = list(inst[item])
     1134            self.PatternTree.SetItemPyData(self.PatternTree.AppendItem(Id,text='Instrument Parameters'),inst)
    11231135                self.PatternTree.SetItemPyData(self.PatternTree.AppendItem(Id,text='Comments'),comments)
    11241136                self.PatternTree.SetItemPyData(self.PatternTree.AppendItem(Id,text='Index Peak List'),peaks)
     
    14301442        SumList = []
    14311443        Names = []
    1432         Inst = []
     1444        Inst = {}
    14331445        SumItemList = []
    14341446        Comments = ['Sum equals: \n']
  • trunk/GSASIIIO.py

    r792 r793  
    496496    Comments = G2frame.PatternTree.GetItemPyData(G2gd.GetPatternTreeItemId(G2frame,Id, 'Comments'))
    497497    names = ['Type','Lam','Zero','Polariz.','U','V','W','X','Y','SH/L','Azimuth']
    498     codes = [0 for i in range(11)]
     498    codes = [0 for i in range(12)]
    499499    LRazm = data['LRazimuth']
    500500    Azms = []
     
    522522            G2frame.PatternTree.SetItemPyData(G2gd.GetPatternTreeItemId(G2frame,Id,'Background'),[['chebyschev',1,3,1.0,0.0,0.0],
    523523                            {'nDebye':0,'debyeTerms':[],'nPeaks':0,'peaksList':[]}])
    524             G2frame.PatternTree.SetItemPyData(G2gd.GetPatternTreeItemId(G2frame,Id,'Instrument Parameters'),dict(zip(names,zip(parms,parms,codes))))
     524            inst = dict(zip(names,zip(parms,parms,codes)))
     525            for item in inst:
     526                inst[item] = list(inst[item])
     527            G2frame.PatternTree.SetItemPyData(G2gd.GetPatternTreeItemId(G2frame,Id,'Instrument Parameters'),inst)
    525528            G2frame.PatternTree.SetItemPyData(G2gd.GetPatternTreeItemId(G2frame,Id,'Peak List'),[])
    526529            G2frame.PatternTree.SetItemPyData(G2gd.GetPatternTreeItemId(G2frame,Id,'Index Peak List'),[])
     
    554557        Inst = G2frame.PatternTree.GetItemPyData(G2gd.GetPatternTreeItemId(G2frame, \
    555558            PickId, 'Instrument Parameters'))
    556         print Inst['Type']
    557559        prm.write( '            123456789012345678901234567890123456789012345678901234567890        '+'\n')
    558560        prm.write( 'INS   BANK      1                                                               '+'\n')
    559561        prm.write(('INS   HTYPE   %sR                                                              '+'\n')%(Inst['Type'][0]))
    560562        if 'Lam1' in Inst:              #Ka1 & Ka2
    561             prm.write(('INS  1 ICONS%10.7f%10.7f    0.0000               0.990    0     0.500   '+'\n')%(Inst['Lam1'],Inst['Lam2'][0]))
     563            prm.write(('INS  1 ICONS%10.7f%10.7f    0.0000               0.990    0     0.500   '+'\n')%(Inst['Lam1'][0],Inst['Lam2'][0]))
    562564        elif 'Lam' in Inst:             #single wavelength
    563             prm.write(('INS  1 ICONS%10.7f%10.7f    0.0000               0.990    0     0.500   '+'\n')%(Inst['Lam'][0],0.0))
     565            prm.write(('INS  1 ICONS%10.7f%10.7f    0.0000               0.990    0     0.500   '+'\n')%(Inst['Lam'][1],0.0))
    564566        prm.write( 'INS  1 IRAD     0                                                               '+'\n')
    565567        prm.write( 'INS  1I HEAD                                                                    '+'\n')
     
    567569        prm.write(('INS  1DETAZM%10.3f                                                          '+'\n')%(Inst['Azimuth'][0]))
    568570        prm.write( 'INS  1PRCF1     3    8   0.00100                                                '+'\n')
    569         prm.write(('INS  1PRCF11     %15.6g%15.6g%15.6g%15.6g   '+'\n')%(Inst['U'][0],Inst['V'][0],Inst['W'][0],0.0))
    570         prm.write(('INS  1PRCF12     %15.6g%15.6g%15.6g%15.6g   '+'\n')%(Inst['X'][0],Inst['Y'][0],Inst['SH/L'][0]/2.,Inst['SH/L'][0]/2.))
     571        prm.write(('INS  1PRCF11     %15.6g%15.6g%15.6g%15.6g   '+'\n')%(Inst['U'][1],Inst['V'][1],Inst['W'][1],0.0))
     572        prm.write(('INS  1PRCF12     %15.6g%15.6g%15.6g%15.6g   '+'\n')%(Inst['X'][1],Inst['Y'][1],Inst['SH/L'][1]/2.,Inst['SH/L'][1]/2.))
    571573        prm.close()
    572574        file = open(filename,'w')
Note: See TracChangeset for help on using the changeset viewer.