show_cv.py: allow plotting multiple VIIRS granules

cris-viirs/show_cv.py

@@ -10,22 +10,24 @@ def main():
 
     f = h5py.File(sys.argv[1], 'r')
     g = f['All_Data']['CrIS-SDR-GEO_All']
-    lat = g['Latitude'][:,:,4]
-    lon = g['Longitude'][:,:,4]
+    lat = g['Latitude'][2:-2,:,4]
+    lon = g['Longitude'][2:-2,:,4]
     b = Basemap(projection='ortho', lon_0=mid(lon), lat_0=mid(lat))
-    b.plot(*b(mask(outline(lon)), mask(outline(lat))))
+    b.plot(*b(mask(outline(lon)), mask(outline(lat))), color='b')
 
-    f = h5py.File(sys.argv[2], 'r')
-    g = f['All_Data']['VIIRS-MOD-GEO-TC_All']
-    lat = g['Latitude'][:]
-    lon = g['Longitude'][:]
-    b.plot(*b(mask(outline(lon)), mask(outline(lat))))
+    for gmtco_file in sys.argv[2:]:
+        f = h5py.File(gmtco_file, 'r')
+        g = f['All_Data']['VIIRS-MOD-GEO-TC_All']
+        lat = g['Latitude'][:]
+        lon = g['Longitude'][:]
+        b.plot(*b(mask(outline(lon)), mask(outline(lat))), color='g')
 
     b.drawcoastlines()
     b.drawmapboundary()
     plt.show()
 
 def mid(a):
+
     idx = tuple(np.array(a.shape) / 2)
     return a[idx]
 
@@ -35,7 +37,8 @@ def outline(a):
 
 def mask(a):
 
-    return np.ma.masked_array(a, a == -9999.0)
+    return np.ma.masked_array(a, a < -200.0)
 
 if __name__ == '__main__':
     main()
+