Changeset 989 for trunk/exports


Ignore:
Timestamp:
Jul 13, 2013 10:21:12 PM (8 years ago)
Author:
toby
Message:

add constraint derivs for single xtal; allow cellFill to work without esds; update docs; fix controls copy bug

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/exports/G2cif.py

    r981 r989  
    99import GSASIIgrid as G2gd
    1010import GSASIIstrIO as G2stIO
     11#reload(G2stIO)
    1112#import GSASIImapvars as G2mv
    1213import GSASIImath as G2mth
     14reload(G2mth)
    1315import GSASIIlattice as G2lat
    1416import GSASIIspc as G2spg
    1517#reload(G2spg)
    16 reload(G2mth)
    1718
    1819def getCallerDocString(): # for development
     
    585586            refprx = '_refln_' # normal
    586587
    587             print histblk.keys()
    588             #            for key in histblk:
    589             #                print key
    590             print inst
    591             print self.parmDict.keys()
    592             print self.sigDict.keys()
    593588            WriteCIFitem('\n# SCATTERING FACTOR INFO')
    594589            if 'Lam1' in inst:
     
    618613                WriteCIFitem('_diffrn_radiation_wavelength',G2mth.ValEsd(lam1,slam1))
    619614
    620             raise Exception, "testing"
    621615
    622616            if not oneblock:
     
    629623                                 '\n\t_pd_phase_block_id' +
    630624                                 '\n\t_pd_phase_mass_%')
     625                    wtFrSum = 0.
    631626                    for phasenam in phasebyhistDict.get(histlbl):
    632                         pass
    633 
    634             WriteCIFitem('_pd_proc_ls_prof_R_factor','?')
    635             WriteCIFitem('_pd_proc_ls_prof_wR_factor','?')
    636             WriteCIFitem('_pd_proc_ls_prof_wR_expected','?')
    637             WriteCIFitem('_refine_ls_R_Fsqd_factor','?')
     627                        hapData = self.Phases[phasenam]['Histograms'][histlbl]
     628                        General = self.Phases[phasenam]['General']
     629                        wtFrSum += hapData['Scale'][0]*General['Mass']
     630
     631                    for phasenam in phasebyhistDict.get(histlbl):
     632                        hapData = self.Phases[phasenam]['Histograms'][histlbl]
     633                        General = self.Phases[phasenam]['General']
     634                        wtFr = hapData['Scale'][0]*General['Mass']/wtFrSum
     635                        pfx = str(self.Phases[phasenam]['pId'])+':'+str(hId)+':'
     636                        if pfx+'Scale' in self.sigDict:
     637                            sig = self.sigDict[pfx+'Scale']*wtFr/hapData['Scale'][0]
     638                        else:
     639                            sig = -0.0001
     640                        WriteCIFitem(
     641                            '  '+
     642                            str(self.Phases[phasenam]['pId']) +
     643                            '  '+datablockidDict[phasenam]+
     644                            '  '+G2mth.ValEsd(wtFr,sig)
     645                            )
     646
     647            # this will need help from Bob
     648            # WriteCIFitem('_pd_proc_ls_prof_R_factor','?')
     649            # WriteCIFitem('_pd_proc_ls_prof_wR_factor','?')
     650            # WriteCIFitem('_pd_proc_ls_prof_wR_expected','?')
     651            # WriteCIFitem('_refine_ls_R_Fsqd_factor','?')
     652
     653            phasenam = self.Phases.keys()[0]
     654            for key in self.Phases[phasenam]['Histograms']:
     655                print key
     656                print '------------'
     657                print self.Phases[phasenam]['Histograms'][key]
     658            raise Exception, "testing"
     659            print histblk.keys()
     660            for key in histblk:
     661                print key,histblk[key]
     662            print inst
     663            #print self.parmDict.keys()
     664            #print self.sigDict.keys()
    638665           
    639666            #WriteCIFitem('_pd_meas_2theta_fixed',text)
     
    10501077                    WriteCIFitem('_pd_block_id',datablockidDict[hist])
    10511078                    WritePowderTemplate()
    1052                     WritePowderData(key1)
     1079                    WritePowderData(hist)
    10531080                elif hist.startswith("HKLF"):
    10541081                    WriteCIFitem('\ndata_'+self.CIFname+"_sx_"+str(i))
     
    10581085                    WriteCIFitem('_pd_block_id',datablockidDict[hist])
    10591086                    WriteSnglXtalTemplate()
    1060                     WriteSingleXtalData(key1)
     1087                    WriteSingleXtalData(hist)
    10611088
    10621089        WriteCIFitem('#--' + 15*'eof--' + '#')
Note: See TracChangeset for help on using the changeset viewer.