Added some basics to store/plot IASI ash mass (BTA-BTB method

bin/IASITools.py

@@ -1215,8 +1215,20 @@ class Scene(object):
             #             levels=[self.SEVIRI_dBTLimit]
             contour_colors = ['g', 'r']
 
-        elif  Type=='SEVIRI_AshML':
+        elif Type=='SEVIRI_AshML':
             data = self.MassLoading[:]
+            #            data = self.AshML[:]*1000  # Convert from g/m**3 to g/m**2 assuming a 1 km thick cloud
+            vmin=0
+            vmax= 10.0 #2.5 #0.007*1000 # 2 # 20 #data.max()*0.5#
+            step= 1.0 #0.8 # vmax/5. #0.01
+            cblabel='Ash Mass Loading (g/m$^2$)'
+            #cmap = plt.get_cmap('gist_ncar_r')
+            cmap = plt.get_cmap('jet')
+            #levels=[.5, 2.0]
+            levels=[.1]
+
+        elif Type=='IASI_AshML':
+            data = self.MassLoadingIASI[:]
 #            data = self.AshML[:]*1000  # Convert from g/m**3 to g/m**2 assuming a 1 km thick cloud
             vmin=0
             vmax= 10.0 #2.5 #0.007*1000 # 2 # 20 #data.max()*0.5#
@@ -1400,6 +1412,11 @@ class Scene(object):
                 except:
                     pass
 
+                try:
+                    self.MassLoadingIASI = ncfile.variables['MassLoadingIASI'][:]
+                except:
+                    pass
+
                 self.npoints = len(self.lats)
                 # These variables are possibly missing
                 try:

bin/plotCombined.py

@@ -39,6 +39,9 @@ def plotCombined(file_name,
     elif Type == "SEVIRI_AshML":
         plotmethod = 'SEVIRI_AshML'
         str1 = 'IASI_dBT'
+    elif Type == "IASI_AshML":
+        plotmethod = 'IASI_AshML'
+        str1 = 'IASI_dBT'        
     elif Type == "SEVIRI_dBT_AF":
         plotmethod = "SEVIRI_dBT_AF"
         str1 = 'IASI_dBT'
@@ -97,7 +100,8 @@ if __name__ == "__main__":
     parser.add_argument('-t', '--type', 
                         dest='Type',
                         default ='SEVIRI_dBT',
-                        help="plot type [IASI_dBT, SEVIRI_dBT, SEVIRI_AshML, SEVIRI_dBT_AF, IASI_dBT_AF]")
+                        help=("plot type [IASI_dBT, SEVIRI_dBT, SEVIRI_AshML, IASI_AshML, "
+                              "SEVIRI_dBT_AF, IASI_dBT_AF]"))
     parser.add_argument('-m', '--map_projection', 
                         dest='map_projection',
                         default ='PS_NE',