Changeset 526 for trunk/GSASIIIO.py


Ignore:
Timestamp:
Apr 2, 2012 2:30:28 PM (10 years ago)
Author:
vondreele
Message:

make SetNewPhase? routine in GSASIIIO.py to initialize phase info
begin Fourier map search routine
fix atom plotting if no map

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIIO.py

    r502 r526  
    10351035    print 'index peak list saved'
    10361036   
     1037def SetNewPhase(Name='New Phase',SGData=G2spc.SpcGroup('P 1')[1],cell=[1.0,1.0,1.0,90.,90,90.,1.]):
     1038    phaseData = {
     1039        'General':{
     1040            'Name':Name,
     1041            'Type':'nuclear',
     1042            'SGData':SGData,
     1043            'Cell':[False,]+cell,
     1044            'Pawley dmin':1.0,
     1045            'Data plot type':'Mustrain',
     1046            'SH Texture':{
     1047                'Order':0,
     1048                'Model':'cylindrical',
     1049                'Sample omega':[False,0.0],
     1050                'Sample chi':[False,0.0],
     1051                'Sample phi':[False,0.0],
     1052                'SH Coeff':[False,{}],
     1053                'SHShow':False,
     1054                'PFhkl':[0,0,1],
     1055                'PFxyz':[0,0,1],
     1056                'PlotType':'Pole figure'}},
     1057        'Atoms':[],
     1058        'Drawing':{},
     1059        'Histograms':{},
     1060        'Pawley ref':[],
     1061        'Models':{},
     1062        }
     1063    return phaseData
     1064   
    10371065def ReadEXPPhase(G2frame,filename):
    10381066    shModels = ['cylindrical','none','shear - 2/m','rolling - mmm']
     
    11501178        textureData['SH Coeff'] = [False,shCoef]
    11511179       
    1152     Phase = {
    1153             'General':{
    1154                 'Name':PhaseName,
    1155                 'Type':Ptype,
    1156                 'SGData':SGData,
    1157                 'Cell':[False,]+abc+angles+[Volume,],
    1158                 'Pawley dmin':1.0,
    1159                 'SH Texture':textureData},
    1160             'Atoms':Atoms,
    1161             'Drawing':{},
    1162             'Histograms':{},
    1163             'Pawley ref':[],
    1164             'Models':{},
    1165             }
     1180    Phase = SetNewPhase(Name=PhaseName,SGData=SGData,cell=abc+angles+[Volume,])
     1181    general = Phase['General']
     1182    general['Type'] = Ptype
     1183    general['SH Texture'] = textureData
     1184    Phase['Atoms'] = Atoms
    11661185    return Phase
    11671186       
     
    12411260    else:
    12421261        PhaseName = 'None'
    1243     Phase['General'] = {'Name':PhaseName,'Type':'macromolecular','SGData':SGData,
    1244         'Cell':[False,]+cell+[Volume,]}
     1262    Phase = SetNewPhase(Name=PhaseName,SGData=SGData,cell=cell+[Volume,])
     1263    Phase['General']['Type'] = 'macromolecular'
    12451264    Phase['Atoms'] = Atoms
    1246     Phase['Drawing'] = {}
    1247     Phase['Histograms'] = {}
    12481265   
    12491266    return Phase
     
    12771294        self.strictExtension = strictExtension
    12781295        # define a default Phase structure
    1279         self.Phase = {}
    1280         for i in 'General', 'Atoms', 'Drawing', 'Histograms':
    1281             self.Phase[i] = {}
    1282         self.Phase['General']['Name'] = 'default'
    1283         self.Phase['General']['Type'] = 'nuclear'
    1284         self.Phase['General']['SGData'] = SGData
    1285         self.Phase['General']['Cell'] = [
    1286             False, # refinement flag
    1287             1.,1.,1.,    # a,b,c
    1288             90.,90.,90., # alpha, beta, gamma
    1289             1.           # volume
    1290             ]
     1296        self.Phase = SetNewPhase(Name='new phase',SGData=SGData)
    12911297        self.warnings = ''
    12921298        self.errors = ''
Note: See TracChangeset for help on using the changeset viewer.