Changeset 793
- Timestamp:
- Oct 24, 2012 9:38:52 AM (10 years ago)
- Location:
- trunk
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/GSASII.py
r792 r793 557 557 S = File.readline() 558 558 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 560 563 561 564 def ReadPowderIparm(self,instfile,bank,databanks,rd): … … 655 658 data.extend([0.0,0.0,0.002,azm]) #OK defaults if fxn #3 not 1st in iprm file 656 659 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 658 664 elif 'T' in DataType: 659 665 names = ['Type','2-theta','difC','difA','Zero','alpha','beta-0','beta-1','var-inst','X','Y','Azimuth'] … … 682 688 s = Iparm['INS 1PRCF12'].split() 683 689 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 685 694 686 695 # stuff we might need from the reader … … 1119 1128 data = ['PKS',Cuka,0.0] 1120 1129 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) 1123 1135 self.PatternTree.SetItemPyData(self.PatternTree.AppendItem(Id,text='Comments'),comments) 1124 1136 self.PatternTree.SetItemPyData(self.PatternTree.AppendItem(Id,text='Index Peak List'),peaks) … … 1430 1442 SumList = [] 1431 1443 Names = [] 1432 Inst = []1444 Inst = {} 1433 1445 SumItemList = [] 1434 1446 Comments = ['Sum equals: \n'] -
trunk/GSASIIIO.py
r792 r793 496 496 Comments = G2frame.PatternTree.GetItemPyData(G2gd.GetPatternTreeItemId(G2frame,Id, 'Comments')) 497 497 names = ['Type','Lam','Zero','Polariz.','U','V','W','X','Y','SH/L','Azimuth'] 498 codes = [0 for i in range(1 1)]498 codes = [0 for i in range(12)] 499 499 LRazm = data['LRazimuth'] 500 500 Azms = [] … … 522 522 G2frame.PatternTree.SetItemPyData(G2gd.GetPatternTreeItemId(G2frame,Id,'Background'),[['chebyschev',1,3,1.0,0.0,0.0], 523 523 {'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) 525 528 G2frame.PatternTree.SetItemPyData(G2gd.GetPatternTreeItemId(G2frame,Id,'Peak List'),[]) 526 529 G2frame.PatternTree.SetItemPyData(G2gd.GetPatternTreeItemId(G2frame,Id,'Index Peak List'),[]) … … 554 557 Inst = G2frame.PatternTree.GetItemPyData(G2gd.GetPatternTreeItemId(G2frame, \ 555 558 PickId, 'Instrument Parameters')) 556 print Inst['Type']557 559 prm.write( ' 123456789012345678901234567890123456789012345678901234567890 '+'\n') 558 560 prm.write( 'INS BANK 1 '+'\n') 559 561 prm.write(('INS HTYPE %sR '+'\n')%(Inst['Type'][0])) 560 562 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])) 562 564 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)) 564 566 prm.write( 'INS 1 IRAD 0 '+'\n') 565 567 prm.write( 'INS 1I HEAD '+'\n') … … 567 569 prm.write(('INS 1DETAZM%10.3f '+'\n')%(Inst['Azimuth'][0])) 568 570 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.)) 571 573 prm.close() 572 574 file = open(filename,'w')
Note: See TracChangeset
for help on using the changeset viewer.