Changeset 3800


Ignore:
Timestamp:
Jan 26, 2019 3:46:41 PM (3 years ago)
Author:
toby
Message:

remove duplicate _pd_proc_ls_profile_function entry

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/exports/G2export_CIF.py

    r3796 r3800  
    7171            if name.strip():
    7272                fp.write(name+'\n')
    73             fp.write('; '+value+'\n')
    74             fp.write('; '+'\n')
     73            fp.write(';\n'+value+'\n')
     74            fp.write(';'+'\n')
    7575        elif " " in value:
    7676            if len(name)+len(value) > 65:
     
    753753                mustrain = hapData['Mustrain']
    754754                hstrain = hapData['HStrain']
    755                 if len(self.powderDict) > 1:
    756                     if s:
    757                         s += '\n'
    758                     else:
    759                         s += '  Crystallite size model "%s" for %s (microns)\n  '%(size[0],phasenam)
    760                     s += '  Parameters for histogram #'+str(hId)+' '+str(histogram)+'\n'
    761                 else:
    762                     s += '  Crystallite size model "%s" for %s (microns)\n  '%(size[0],phasenam)
    763 
     755                if s: s += '\n'
     756                if len(self.powderDict) > 1: # if one histogram, no ambiguity
     757                    s += '  Parameters for histogram #{:} {:} & phase {:}\n'.format(
     758                        str(hId),str(histogram),phasenam)
     759                s += '  Crystallite size in microns with "%s" model:\n  '%(size[0])
    764760                names = ['Size;i','Size;mx']
    765761                if 'uniax' in size[0]:
     
    787783                        s += G2mth.ValEsd(size[1][i],sig)+', '
    788784                        i = 2    #skip the aniso value
    789                 s += '\n  Mustrain model "%s" for %s (10^6^)\n  '%(mustrain[0],phasenam)
     785                s += '\n  Microstrain, "%s" model (10^6^ * delta Q/Q)\n  '%(mustrain[0])
    790786                names = ['Mustrain;i','Mustrain;mx']
    791787                if 'uniax' in mustrain[0]:
     
    14671463            #CALL WRVAL(IUCIF,'_gsas_exptl_extinct_corr_T_max',TEXT(11:20))
    14681464
    1469             if not oneblock:                 # instrumental profile terms go here
    1470                 WriteCIFitem(self.fp, '_pd_proc_ls_profile_function',
    1471                     FormatInstProfile(histblk["Instrument Parameters"],histblk['hId']))
     1465            # code removed because it is causing duplication in histogram block 1/26/19 BHT
     1466            #if not oneblock:                 # instrumental profile terms go here
     1467            #    WriteCIFitem(self.fp, '_pd_proc_ls_profile_function',
     1468            #        FormatInstProfile(histblk["Instrument Parameters"],histblk['hId']))
    14721469
    14731470            #refprx = '_refln.' # mm
     
    22982295            # report the phase info
    22992296            WritePhaseInfo(phasenam)
    2300             if hist.startswith("PWDR"):
     2297            if hist.startswith("PWDR"):  # this is invoked for single-block CIFs
    23012298                # preferred orientation
    23022299                SH = FormatSH(phasenam)
     
    23082305                else:
    23092306                    WriteCIFitem(self.fp, '_pd_proc_ls_pref_orient_corr', 'none')
    2310                     # report profile, since one-block: include both histogram and phase info
     2307                # report profile, since one-block: include both histogram and phase info (N.B. there is only 1 of each)
    23112308                WriteCIFitem(self.fp, '_pd_proc_ls_profile_function',
    23122309                    FormatInstProfile(histblk["Instrument Parameters"],histblk['hId'])
    23132310                    +'\n'+FormatPhaseProfile(phasenam))
     2311
    23142312                histblk = self.Histograms[hist]["Sample Parameters"]
    23152313                writeCIFtemplate(histblk,'powder',histblk['InstrName']) # write powder template
     
    24052403                    else:
    24062404                        WriteCIFitem(self.fp, '_pd_proc_ls_pref_orient_corr', 'none')
    2407                 # report sample profile terms
     2405                # report sample profile terms for all histograms with current phase
    24082406                PP = FormatPhaseProfile(phasenam)
    24092407                if PP:
Note: See TracChangeset for help on using the changeset viewer.