Starting work on ash retrievals SEVIRI

bin/IASITools.py

@@ -547,14 +547,14 @@ class Scene(object):
             ylabel = 'Brightness temperature (K)'
         else:
             print "PlotSpectrum: unknown Type", Type
-            exit()
+            sys.exit()
         col = 'g'
         pl, = ax.plot(x, y, color=col)
         ax.set_xlabel(xlabel)
         ax.set_ylabel(ylabel)
         ax.set_title(title)
 #        plt.show()
-        exit()
+        sys.exit()
 
     def SelectSpectrum(self,SpectrumNr,verbose=False):
         latitude = self.lats[SpectrumNr]
@@ -715,7 +715,7 @@ class Scene(object):
         ax.set_ylabel(ylabel, fontsize = fsize)
 
 #        plt.show()
-        exit()
+        sys.exit()
         return
 
     def IfCrossDateLineCorrectLongitude(self, verbose=False):
@@ -796,7 +796,7 @@ class Scene(object):
         plt.scatter(x,y,marker='o',c='b',s=5)
 
 #        plt.show()
-        exit()
+        sys.exit()
         return
 
 
@@ -1203,7 +1203,7 @@ class Scene(object):
         if PlotFile==None:
             if PlotShow:
 #                plt.show()
-                exit()
+                sys.exit()
             else:
                 return ax
         else:

bin/ReadAndPlotIASI.py

@@ -33,6 +33,7 @@ IASI_PATH = "/sat_data_from_fred/eumetcast/EPS-11/"
 # for met.no
 NETCDF_OUT = "/viper/nadir/nvap/IASI/"
 PNG_OUT = "/viper/nadir/nvap/IASI/"
+args = None
 
 def at_exit(file_name):
     """ Delete lock file on exit.
@@ -59,23 +60,23 @@ def ReadAndPlotSingleIASIFile(fn, Projection='lcc', PlotFile=None, ncfFile=None,
         #    Scene.CalculateSEVIRIChannels()
 
 
-	if (args.write_png):
-        # Type is plot type
-        # dBT is brightness temperature difference between two IASI channels
-		Scene.PlotLatLonVals(Type='dBT', #'BT', # 'SEVIRIdBT',
-				     Projection=Projection,
-                             # llcrnrlon=181,
-                             # llcrnrlat=50,
-                             # urcrnrlon=120,
-                             # urcrnrlat=70,
-				     verbose=verbose,
-				     Title=True,
-				     test=True,
-				     PlotFile=PlotFile,
-				     PlotAshContour=False)
+	if getattr(args, 'write_png', True):
+            # Type is plot type
+            # dBT is brightness temperature difference between two IASI channels
+            Scene.PlotLatLonVals(Type='dBT', #'BT', # 'SEVIRIdBT',
+                                 Projection=Projection,
+                                 # llcrnrlon=181,
+                                 # llcrnrlat=50,
+                                 # urcrnrlon=120,
+                                 # urcrnrlat=70,
+                                 verbose=verbose,
+                                 Title=True,
+                                 test=True,
+                                 PlotFile=PlotFile,
+                                 PlotAshContour=False)
 
 
-	Scene.WriteNetcdf(fnncf)
+	Scene.WriteNetcdf(ncfFile)
 
         return Scene