diff --git a/GalProf/detcnts.E01.ALL.resid.gal.fits b/GalProf/detcnts.E01.ALL.resid.gal.fits
index 503ac4e..0d688d8 100644
Binary files a/GalProf/detcnts.E01.ALL.resid.gal.fits and b/GalProf/detcnts.E01.ALL.resid.gal.fits differ
diff --git a/GalProf/detcnts.E01.ALL.resid.galprof.fits b/GalProf/detcnts.E01.ALL.resid.galprof.fits
index 4aae75a..102a3ed 100644
Binary files a/GalProf/detcnts.E01.ALL.resid.galprof.fits and b/GalProf/detcnts.E01.ALL.resid.galprof.fits differ
diff --git a/scripts/03_grxe_galprof.py b/scripts/03_grxe_galprof.py
index 888b921..ebabb76 100755
--- a/scripts/03_grxe_galprof.py
+++ b/scripts/03_grxe_galprof.py
@@ -58,14 +58,14 @@ maxiters=10
 
 with open(ignored_rev_file, 'rb') as fp:
     ignored_rev = pickle.load(fp)
-    print(ignored_rev)
+    #print(ignored_rev)
     print("{} orbits ignored".format(len(ignored_rev)))
 
 ign=ignored_rev.tolist()
 
 
 glon, step = np.linspace(-lon_max, lon_max, num=lon_nbin, endpoint=False,retstep=True)
-print(glon,step)
+#print(glon,step)
 
 
 #
@@ -99,6 +99,8 @@ lat_arr=[]
 for i in range(lon_nbin):
     dkey="{:04d}".format(i)
 
+    #if(i != 5):
+    #    continue
     
     df0 = df.query('LON > {} & LON <= {} & LAT > {} & LAT < {} & REV != @ign'.format(glon[i],glon[i]+step,-lat_max,lat_max))
 
@@ -124,7 +126,7 @@ for i in range(lon_nbin):
         lon_arr.append(row['LON'])
         lat_arr.append(row['LAT'])
 
-
+#sys.exit()
 
         
 if not os.path.exists(profdir):
diff --git a/scripts/03_grxe_galprof.sh b/scripts/03_grxe_galprof.sh
index ee04dd2..4afd58a 100644
--- a/scripts/03_grxe_galprof.sh
+++ b/scripts/03_grxe_galprof.sh
@@ -1,4 +1,4 @@
-./03_grxe_galprof.py E01 71
+./03_grxe_galprof.py E01 51
 ./03_grxe_galprof.py E13 17
 ./03_grxe_galprof.py E14 21
 
diff --git a/scripts/03_grxe_galprof_plot.py b/scripts/03_grxe_galprof_plot.py
index 2ffd023..747e601 100755
--- a/scripts/03_grxe_galprof_plot.py
+++ b/scripts/03_grxe_galprof_plot.py
@@ -4,7 +4,7 @@ __author__      = "Roman Krivonos"
 __copyright__   = "Space Research Institute (IKI)"
 
 from astropy.table import Table, Column
-
+from matplotlib import ticker
 import matplotlib.pyplot as plt
 import sys
 
@@ -51,30 +51,38 @@ ax1.ticklabel_format(style='sci', axis='y', scilimits=(0,0))
 ax2.ticklabel_format(style='sci', axis='y', scilimits=(0,0))
 ax3.ticklabel_format(style='sci', axis='y', scilimits=(0,0))
 
+ax1.xaxis.set_minor_locator(ticker.MultipleLocator(10))
+ax2.xaxis.set_minor_locator(ticker.MultipleLocator(10))
+ax3.xaxis.set_minor_locator(ticker.MultipleLocator(10))
+
+ax1.yaxis.set_minor_locator(ticker.MultipleLocator(0.2))
+ax2.yaxis.set_minor_locator(ticker.MultipleLocator(0.2))
+ax3.yaxis.set_minor_locator(ticker.MultipleLocator(0.2))
+
 lon=(df1['LON1']+df1['LON2'])/2
 ax1.errorbar(lon, df1['GRXE_SIM_FLUX']/scale,
              yerr=df1['GRXE_SIM_ERROR']/scale,
              xerr=(df1['LON2']-df1['LON1'])/2,
-             fmt='o' )
+             fmt='o', color='black')
 ax1.grid(visible=True)
 
 lon=(df2['LON1']+df2['LON2'])/2
 ax2.errorbar(lon, df2['GRXE_SIM_FLUX']/scale,
              yerr=df2['GRXE_SIM_ERROR']/scale,
-             xerr=(df2['LON2']-df2['LON1'])/2,fmt='o' )
+             xerr=(df2['LON2']-df2['LON1'])/2,fmt='o' , color='black')
 ax2.grid(visible=True)
 
 lon=(df3['LON1']+df3['LON2'])/2
 ax3.errorbar(lon, df3['GRXE_SIM_FLUX']/scale,
              yerr=df3['GRXE_SIM_ERROR']/scale,
-             xerr=(df3['LON2']-df3['LON1'])/2, fmt='o' )
+             xerr=(df3['LON2']-df3['LON1'])/2, fmt='o' , color='black')
 ax3.grid(visible=True)
 
 #x = np.arange(-10, 10, 0.001)
 #plot normal distribution with mean 0 and standard deviation 1
 #plt.plot(x, norm.pdf(x, 0, 1), color='red', linewidth=2)
 	    
-plt.xlabel('Galactic Longitude',fontsize=14, fontweight='normal')
+plt.xlabel('Galactic Longitude, deg.',fontsize=14, fontweight='normal')
 ax2.set_ylabel('GRXE flux, x100 mCrab',fontsize=14, fontweight='normal')
 
 #plt.xscale('linear')