Changeset 2342


Ignore:
Timestamp:
Jun 24, 2016 3:36:47 PM (5 years ago)
Author:
toby
Message:

fix reading of CIF intensities with esds

Location:
trunk
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASII.py

    r2341 r2342  
    14941494                    var /= YI**2
    14951495                    rd.powderdata[2] = 1./var
    1496                 rd.powderdata[3] = np.zeros_like(rd.powderdata[0])                                       
    1497                 rd.powderdata[4] = np.zeros_like(rd.powderdata[0])                                       
    1498                 rd.powderdata[5] = np.zeros_like(rd.powderdata[0])                                       
     1496                rd.powderdata[3] = np.zeros_like(rd.powderdata[0])
     1497                rd.powderdata[4] = np.zeros_like(rd.powderdata[0])
     1498                rd.powderdata[5] = np.zeros_like(rd.powderdata[0])                 
    14991499            valuesdict = {
    15001500                'wtFactor':1.0,
  • trunk/imports/G2pwd_CIF.py

    r1168 r2342  
    273273                vl = np.zeros(len(x)) + 1.
    274274            else:
     275                vl = []
    275276                sucf = such[sui]
    276277                if sucf ==  '_pd_proc_ls_weight':
     
    293294                    for val in cf[blk].get(sucf,'?'):
    294295                        v,e = cif.get_number_with_esd(val)
    295                         if v is None: # not parsed
    296                             vl.append(0.0)
    297                         elif v <= 0:
    298                             vl.append(0.0)
    299                         else:
    300                             vl.append(1./(v*v))
     296                        if v is None or e is None: # not parsed or no ESD
     297                            vl.append(np.NaN)
     298                        elif e <= 0:
     299                            vl.append(1.)
     300                        else:
     301                            vl.append(1./(e*e))
    301302            w = np.array(vl)
    302303            # intensity modification factor
Note: See TracChangeset for help on using the changeset viewer.