Changeset 743


Ignore:
Timestamp:
Aug 29, 2012 12:48:04 PM (10 years ago)
Author:
vondreele
Message:

fix PO refinement output & calcs for SH

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIstruct.py

    r742 r743  
    12891289            line += ' equatorial:%12.1f'%(hapData[1][0])
    12901290            if mustrainSig[0][0]:
    1291                 line += ', sig: %8.1f'%(mustrainSig[0][0])
     1291                line += ', sig:%8.1f'%(mustrainSig[0][0])
    12921292                refine = True
    12931293            if hapData[0] == 'uniaxial':
     
    13441344            print >>pFile,sigstr
    13451345       
    1346     def PrintSHPOAndSig(hapData,POsig):
     1346    def PrintSHPOAndSig(pdx,hapData,POsig):
    13471347        print >>pFile,'\n Spherical harmonics preferred orientation: Order:'+str(hapData[4])
    13481348        ptlbls = ' names :'
     
    13521352            ptlbls += '%12s'%(item)
    13531353            ptstr += '%12.3f'%(hapData[5][item])
    1354             if item in POsig:
    1355                 sigstr += '%12.3f'%(POsig[item])
     1354            if pfx+item in POsig:
     1355                sigstr += '%12.3f'%(POsig[pfx+item])
    13561356            else:
    13571357                sigstr += 12*' '
     
    14711471                            print >>pFile,' March-Dollase PO: %10.4f, sig %10.4f'%(hapData['Pref.Ori.'][1],PhFrExtPOSig[pfx+'MD'])
    14721472                    else:
    1473                         PrintSHPOAndSig(hapData['Pref.Ori.'],PhFrExtPOSig)
     1473                        PrintSHPOAndSig(pfx,hapData['Pref.Ori.'],PhFrExtPOSig)
    14741474                    PrintSizeAndSig(hapData['Size'],SizeMuStrSig[pfx+'Size'])
    14751475                    PrintMuStrainAndSig(hapData['Mustrain'],SizeMuStrSig[pfx+'Mustrain'],SGData)
     
    22332233            POcorr = sumMD/len(refl[11])
    22342234    else:   #spherical harmonics
    2235         if calcControls[pfx+'SHord']:
     2235        if calcControls[phfx+'SHord']:
    22362236            POcorr = SHPOcal(refl,g,phfx,hfx,SGData,calcControls,parmDict)
    22372237    return POcorr
     
    22532253        POderv[phfx+'MD'] = sumdMD/len(refl[11])
    22542254    else:   #spherical harmonics
    2255         if calcControls[pfx+'SHord']:
     2255        if calcControls[phfx+'SHord']:
    22562256            POcorr,POderv = SHPOcalDerv(refl,g,phfx,hfx,SGData,calcControls,parmDict)
    22572257    return POcorr,POderv
Note: See TracChangeset for help on using the changeset viewer.