Changeset 3961 for trunk/GSASIIplot.py
- Timestamp:
- May 8, 2019 6:16:01 PM (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/GSASIIplot.py
r3959 r3961 7675 7675 7676 7676 # PlotStructure initialization here 7677 global mcsaXYZ,mcsaTypes,mcsaBonds,txID,contourSet 7677 global mcsaXYZ,mcsaTypes,mcsaBonds,txID,contourSet,Zslice 7678 7678 global cell, Vol, Amat, Bmat, A4mat, B4mat 7679 7679 txID = 0 … … 8904 8904 RenderPlane(plane,color) 8905 8905 if drawingData.get('showSlice',False): 8906 global contourSet 8906 global contourSet,Zslice 8907 8907 if len(D4mapData.get('rho',[])): #preferentially select 4D map if there 8908 8908 modQ = np.array(generalData['SuperVec'][0]) … … 8926 8926 if FourD: 8927 8927 SXYZT = np.vstack((SXYZ.T,np.inner(SXYZ,modQ)+G2frame.tau)).T 8928 Z = np.reshape(map_coordinates(rho,(SXYZT%1.*rho.shape).T,order=1,mode='wrap'),(npts,npts))8928 Zslice = np.reshape(map_coordinates(rho,(SXYZT%1.*rho.shape).T,order=1,mode='wrap'),(npts,npts)) 8929 8929 else: 8930 Z = np.reshape(map_coordinates(rho,(SXYZ%1.*rho.shape).T,order=1,mode='wrap'),(npts,npts))8931 Z = np.where(Z <=Rmax,Z,Rmax)8930 Zslice = np.reshape(map_coordinates(rho,(SXYZ%1.*rho.shape).T,order=1,mode='wrap'),(npts,npts)) 8931 Z = np.where(Zslice<=Rmax,Zslice,Rmax) 8932 8932 plt.rcParams['figure.facecolor'] = (1.,1.,1.,.5) 8933 8933 oldSize = plt.rcParams['figure.figsize']
Note: See TracChangeset
for help on using the changeset viewer.