Changeset 4322
- Timestamp:
- Feb 18, 2020 7:10:02 AM (3 years ago)
- Location:
- trunk
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/GSASIIphsGUI.py
r4321 r4322 2568 2568 Bstd = np.std(Bonds['Obonds']) 2569 2569 title = '%s-%s Octahedral bond lengths'%(Oatoms,Tatoms) 2570 G2plt.PlotBarGraph(G2frame,Bonds['Obonds'],Xname=r'$Bond, \AA$',Title=title,PlotName=' Bond')2570 G2plt.PlotBarGraph(G2frame,Bonds['Obonds'],Xname=r'$Bond, \AA$',Title=title,PlotName='Oct-Bond') 2571 2571 Tilts['Otilts'] = np.array(Tilts['Otilts']) 2572 2572 Tmean = np.mean(Tilts['Otilts']) 2573 2573 Tstd = np.std(Tilts['Otilts']) 2574 2574 G2plt.PlotBarGraph(G2frame,Tilts['Otilts'],Xname='Tilts, deg', 2575 Title='Octahedral %s tilts'%Oatoms,PlotName=' Tilts')2575 Title='Octahedral %s tilts'%Oatoms,PlotName='Oct-Tilts') 2576 2576 dVects['Ovec'] = np.reshape(np.array(dVects['Ovec']),(-1,3)) 2577 2577 for ix,aX in enumerate(['X','Y','Z']): 2578 2578 G2plt.PlotBarGraph(G2frame,dVects['Ovec'].T[ix],Xname=r'$%s%s, \AA$'%(GkDelta,aX), 2579 Title='%s Octahedral distortion'%Oatoms,PlotName=' %s-Delta'%aX)2579 Title='%s Octahedral distortion'%Oatoms,PlotName='Oct %s-Delta'%aX) 2580 2580 Vects['Ovec'] = np.array(Vects['Ovec']) #3D plot of tilt vectors 2581 2581 X = Vects['Ovec'].T[0] … … 2593 2593 Bstd = np.std(Bonds['Tbonds']) 2594 2594 title = '%s-%s Terahedral bond lengths'%(Oatoms,Tatoms) 2595 G2plt.PlotBarGraph(G2frame,Bonds['Tbonds'],Xname=r'$Bond, \AA$',Title=title,PlotName=' Bond')2595 G2plt.PlotBarGraph(G2frame,Bonds['Tbonds'],Xname=r'$Bond, \AA$',Title=title,PlotName='Tet-Bond') 2596 2596 Tilts['Ttilts'] = np.array(Tilts['Ttilts']) 2597 2597 Tmean = np.mean(Tilts['Ttilts']) 2598 2598 Tstd = np.std(Tilts['Ttilts']) 2599 2599 G2plt.PlotBarGraph(G2frame,Tilts['Ttilts'],Xname='Tilts, deg', 2600 Title='Tetrahedral %s tilts'%Oatoms,PlotName='T ilts')2600 Title='Tetrahedral %s tilts'%Oatoms,PlotName='Tet-Tilts') 2601 2601 dVects['Tvec'] = np.reshape(np.array(dVects['Tvec']),(-1,3)) 2602 2602 for ix,aX in enumerate(['X','Y','Z']): 2603 2603 G2plt.PlotBarGraph(G2frame,dVects['Tvec'].T[ix],Xname=r'$%s%s, \AA$'%(GkDelta,aX), 2604 Title='%s Tetrahedral distortion'%Oatoms,PlotName=' %s-Delta'%aX)2604 Title='%s Tetrahedral distortion'%Oatoms,PlotName='Tet %s-Delta'%aX) 2605 2605 Vects['Tvec'] = np.array(Vects['Tvec']) 2606 2606 X = Vects['Tvec'].T[0] -
trunk/GSASIIpwd.py
r4270 r4322 2261 2261 cell = eval(atcode[1].split('+')[1]) 2262 2262 fl.write('%6d%4s [%s]%19.15f%19.15f%19.15f%6d%4d%4d%4d\n'%( 2263 nat,atom[1] ,atcode[0],atom[3],atom[4],atom[5],(iat)%Natm+1,cell[0],cell[1],cell[2]))2263 nat,atom[1].strip(),atcode[0],atom[3],atom[4],atom[5],(iat)%Natm+1,cell[0],cell[1],cell[2])) 2264 2264 fl.close() 2265 2265 return fname,reset
Note: See TracChangeset
for help on using the changeset viewer.