Changeset 4350 for trunk/GSASIIimage.py

Timestamp:

Mar 8, 2020 11:04:21 PM (4 years ago)

Author:

toby

Message:

add warning for azimuths with no pixels

File:

• trunk/GSASIIimage.py (modified) (3 diffs)

trunk/GSASIIimage.py

@@ -1219 +1219 @@
     #make linear interpolators; outside limits give NaN
     H0 = []
-    for h0msk,h2msk in zip(H0msk,H2msk):
+    problemEntries = []
+    for i,(h0msk,h2msk) in enumerate(zip(H0msk,H2msk)):
         try:
             h0int = scint.interp1d(h2msk.compressed(),h0msk.compressed(),bounds_error=False)
@@ -1226 +1227 @@
             H0.append(h0)
         except ValueError:
+            problemEntries.append(i)
             H0.append(np.zeros(numChans,order='F',dtype=np.float32))
     H0 = np.nan_to_num(np.array(H0))        
@@ -1249 +1251 @@
     G2fil.G2Print ("Elapsed time:","%8.3fs"%(time.time()-tbeg))
     G2fil.G2Print ('Integration complete')
+    if problemEntries:
+        msg = ""
+        for i in problemEntries:
+            if msg: msg += ', '
+            msg += "{:.2f}".format((H1[i+1]+H1[i])/2.)
+        print('\nWarning, integrations have no pixels at these azimuth(s)\n\t',msg,'\n')
     if returnN:     #As requested by Steven Weigand
         return H0,H1,H2,NST,CancelPressed