Added get_meta.py

get_meta.py

+import sys
+
+from pyhdf.SD import SD, SDC
+
+# verify command line inputs
+if len(sys.argv) != 3:
+    print 'usage: python %s <airs_file> <output_file>' % sys.argv[0]
+    sys.exit(1)
+airs_filename = sys.argv[1]
+output_filename = sys.argv[2]
+
+# get "metadata" that wer'e interested in from the AIRS file
+sd = SD(airs_filename)
+lat = sd.select('Latitude').get()
+lon = sd.select('Longitude').get()
+scan_ang = sd.select('scanang').get()
+sol_zen = sd.select('solzen').get()
+
+# output the selected arrays to a new HDF file
+output_sd = SD(output_filename, SDC.WRITE | SDC.CREATE | SDC.TRUNC)
+lat_sds = output_sd.create('Latitude', SDC.FLOAT64, lat.shape)
+lon_sds = output_sd.create('Longitude', SDC.FLOAT64, lon.shape)
+scan_ang_sds = output_sd.create('Scan Angle', SDC.FLOAT32, scan_ang.shape)
+sol_zen_sds = output_sd.create('Solar Zenith Angle',
+                               SDC.FLOAT32,
+                               sol_zen.shape)
+lat_sds[:] = lat
+lon_sds[:] = lon
+scan_ang_sds[:] = scan_ang
+sol_zen_sds[:] = sol_zen